Version imprimable
bonjour tout le monde,
le problème que je rencontre en écrivant un code en fortran qui intégrent des équations, est le suivant:: le nombre de pas (temps) que je dois faire est 8270000, ce qui est un nombre trop grand de pas, et en plus ma machine ne peut allouer en memoire que des tableaux de taille (Nx=200,Nt=210000)...je ne sais pas si vous avez des solutions envisageables à mon problème et merci...bonne journée!!
Bonjour,
première idée à froid, tu peux le découper en plusieurs fois non? genre tous les 1000 pas de temps et tu écris les résultats dans un fichier. Tu utilises juste en mémoire ce dont tu as besoin pour calculer...
Salut !
Sous cette forme, c'est sans espoir: les erreurs d'arrondis vont s'accumuler et à la fin, les résultats n'auront plus aucun sens. Je pense deviner pourquoi tu as un si grand nombre de pas: ton système est "stiff" et l'algorithme d'intégration se met à diverger si tu prends un pas plus grand. Fais-moi savoir si c'est bien ça.Citation:
le nombre de pas (temps) que je dois faire est 8270000
Mon expérience m'a montré que les problèmes intrinsèquement "stiff" sont plutôt rares; dans 9 cas sur 10, c'est par maladresse qu'on les a rendu "stiff". Alors explique-nous comment tu es arrivé à tes équations différentielles et on essaiera de voir ce que tu as fait sinon faux, du moins de manière maladroite.
Jean-Marc Blanc
si tu as une équa diff "stiff", tu peux utiliser le solveur CVode (de Sundials). Il est super bien !