Discussion :

[REDK131]

Bonjour,

Je viens de commencer la programmation en Fortran et j'ai comme projet d'élaborer une méthode de résolution Newton afin de l'intégrer à une simulation de réacteur.

Néanmoins, je bloque toujours dans la première partie (le Newton) car lors de l'éxecution le programme fait un Break Point et je ne comprend pas pourquoi ni comment y remédier , je suis toujours en train d'y travailler mais si quelqu'un pourrait m'éclaircie entre temps, je serai vraiment reconnaissante.

Merci d'avance,

Wissale

[REDK131]

J'ai réussi à connaitre l'origine du problème décrit au-dessus. Cependant un autre se présente: j'ai des valeurs aberrantes en exécutant la boucle de mon Newton - à la fin j'ai meme des valeurs de l'ordre de 10^30-. Encore, si quelqu'un pourrait m'éclaircir j'en serai reconnaissante.

[Ehouarn]

Bonjour,

Je n'ai pas regardé ton code.

Comme tu précise que tu débutes:
- Penses-tu bien à compiler ton code avec les options "debug" (recherche de variables non-initialisées, de valeurs aberrantes, de dépassements de tableau, etc.) pour t'assurer qu'il fait bien ce que tu veux qu'il fasse?
- La méthode de Newton ne converge pas nécessairement, ça dépend du problème à résoudre et de la condition initiale qu'on se donne. As-tu testé sur un cas simple (que tu peux suivre à la main, analytiquement) qu'elle est correctement appliquée?

Bon courage.

[REDK121]

Bonjour,

Désolée pour le délai de réponse (j'avais un problème avec ma license de Visual Basic et j'ai perdu mon mot de passe c'est pour ça que je vous répons d'un autre compte):

-En ce qui concerne ce que vous m'avez indiqué, j'avais fait un suivi des valeurs en affichant les différentes solutions de mon systeme linéaire (le deltaX de la méthode Newton-Raphson) avec un Print, en ayant fait le calcul à la main je me suis rendu compte que pour la première itératon c était bon mais dès la seconde les valeurs deviennent incroyablement grandes -elles atteignent des 10^30 très loin des solution recherchées qu avoisinnent 1 et racine de 3-, je me suis dit que c est peut-etre un problème d'initialisation du coup j ai pris l initiative de tout mettre à 0 avant le début des calculs sans succès.

-J'ai aussi vérifié mes autres modules et le calcul des valeurs e la fonction et du jacobien anaytique vont comme il faut, la résolution du système linéaire aussi d'ailleurs dans la premiere itération (pour cette dernière)

Le programme sort tjrs avec un Break point: l array solution (c est un de mes vecteurs) à une valeur de 16 qui est plus grande que 15 l upper-bound

Je compend ce que veut dire cet erreur mais je n'arrive pas à percevoir son origine.:?

Merci Monsieur pour votre réponse et merci d'avance pour tous ceux qui liront ou répondront éventuellement.

[REDK121]

J'ai résolu mon problème de Break Point, c'était à cause d'un compteur mal placé et dont l'expression n'était pas juste et que j'associais au tableau solution qui causait l'erreur, cependant j'ai tjrs des valeurs aberrantes pour mes calculs.:)

Je remets le code mis à jour pour qui voudrais bien y jeter un coup d'oeil.

[REDK121]

Quelqu'un pourrait-il m'aider? :? Ou jeter un coup d'oeil à mon code