Bonjour,
j'ai un doute sur la discrétisation de l'équation d'Alvarez, Lions Morel, à savoir :
du/dt = g( |grad(G)*u| ) uEE
où G est une Gaussienne, g est par exemple la fonction inverse 1/(1+(x/K)²) et * le produit de convolution, et uEE la dérivée par rapport à la direction orthogonale au gradient.
Si c'est pas clair, vous pouvez voir ce lien là :
http://people.csail.mit.edu/sparis/old_project/PDE/
Donc en gros j'essaie de discrétiser cette équation.
Pour la Gaussienne, j'ai pas de soucis. Par contre je ne sais pas comment traiter la convolution grad(G)*u ?
Il faut que je calcule le gradient de G dans chacune des directions puis que je somme en valeur absolue ?
Merci.
Partager