Discussion :

[Butterfly83]

Je ne m y connais pas encore trop, j espère juste que je pose ma question au bon endroit....
Bref, actuellement je travaille sur un code assez complexe (en fortran, mais mes problème ne sont pas, en tout cas pour l instant de l ordre de la programmation), avec un gros problème de raideur.
Ma question c est qu à plusieurs reprises, dans le code sous certaines conditions (en générale c est une condition sur le nombre d itération) le résultat rendu est multiplié pas une certaine constante ou encore une fois par une fonction de mon nombre d itérations, ett je pense que c est une astuce numérique, mais je ne comprend pas trop le but de la manoeuvre...
J espère que j ai réussi à me faire comprendre, je vous remercie...

[vinzzzz]

Faudrait peut être poster un bout de code ou d'algorithme histoire d'avoir une idée plus claire.

Et de plus tu aura plus de chance de trouver ta réponse dans le forum algorithmique ou maths

[genteur slayer]

trop flou ta question... on résout quoi? c'est quoio cette raideur? un ressort? un gradient raide?

[Butterfly83]

En fait ce serrai plus une question de méthode numérique, alors encore une fois je ne sais pas trop si je pose cette question au bon endroit!!!

[souviron34]

la constante est-ele supérieure ou inférieure à 1 ?

[Butterfly83]

Je t ai envoyé un message privé l 'as tu ressu?

[souviron34]

non...

Et de toutes façons je ne réponds pas aux questions techniques par MP. Le forum est fait pour ça

[Butterfly83]

D accord, excuse moi je ne connais pas encore bien le fonctionnement...
Bref..
J utilise une méthode de Newton Raphson pour résoudre un système
mon vecteur solution est deltax=x(i+1)-x(i) a priori
donc sachant que l inconnue que je recherche est x(i+1) je devrais faire
x(i+1)=x(i)+deltax
cependant quand ma fameuse constante (disons C) est supérieure à 10:
P=exp(-fonction(C))
x(i+1)=x(i)+P*deltax

avec fonction(C)>0

Voilà, je te remercie.

[souviron34]

désolé je n'ai rien compris...

Je pense qu'il faudrait expliciter ceci :

Ma question c est qu à plusieurs reprises, dans le code sous certaines conditions (en générale c est une condition sur le nombre d itération) le résultat rendu est multiplié pas une certaine constante ou encore une fois par une fonction de mon nombre d itérations, ett je pense que c est une astuce numérique, mais je ne comprend pas trop le but de la manoeuvre...

qui était ta question initiale...

Et ce que je te demandais par rapport à ça c'est qu'elle était la valeur de ta constante... ou qu'elle était la formule..

[Butterfly83]

En fait c est le P de mon message précédent que je ne comprend pas.
Je simplifie assez, mais je voudrais juste savoir si cette manipulation est courante et comment l interpreter.
La méthode de NR dit
F'(x(n))(x(n+1)-x(n))=-F(x(n))
l inconnue est x(n+1) mais pour éviter d inverser la jacobienne à chaque itération (elle est déjà calculée à chaque itération ce qui engendre déjà un fort coût) notre vecteur inconnu est x(n+1)-x(n) = deltax
or ce que je cherche c est toujours x(n+1)
donc on dit que x(n+1)=deltax+x(n)
cependant il introduit cette constante P dont je ne comprend pas l utilité.

Voilà j espère t avoir éclairé....

[pseudocode]

cependant il introduit cette constante P dont je ne comprend pas l utilité.

P comme... "PAS" d'integration, ou "integration step" en anglais.

Car je rapelle que le principe de l'utilisation de la derivé c'est:

F(x+epsilon) = F(x) + epsilon*F'(x)

avec un "petit" epsilon.


De la à dire que P c'est epsilon, il n'y a qu'un pas (d'integration)...