How to avoid multiple Jacobian calls
1 vue (au cours des 30 derniers jours)
Afficher commentaires plus anciens
Hi all,
I have a problem: when solving a stiff differential equation system (with analytical jacobian provided by me), I noticed (using the profile viewer) that, over the total 33 hours needed to solve the system, 30 are dedicated to calculate the Jacobian analytical matrix.
Now, the Jacobian matrix is a 11000x11000 matrix and, to be computed it uses a lot of "if" and "for", for this reason it takes so long to be calculated.
The ode23s calls the function 32 times (so every call need approx. 1 hour to compute the Jacobian matrix).
Is it possible to avoid the calculation of the Jacobian matrix at every call? I was wondering if it wasn't possible to calculate it just one time and then exploit other functions or the global variables...
Thanks for your help
1 commentaire
Réponses (0)
Voir également
Catégories
En savoir plus sur Ordinary Differential Equations dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!