Bonjour,
Je dois programmer une fonction renvoyant la matrice du laplacien discret en 2D (dans la but d'avoir ses valeurs et vecteurs propres). évidemment, mon maillage est constitué de triangles plus ou moins quelconque, je vous ai mis en image l'énoncé où se trouve la formule à implémenter.
J'ai testé sur un rectangle pour comparer avec la formule analytique. Bien évidement mon programme ne marche pas... les valeurs propres sont délirantes, mais bizarrement, les vecteurs propres semblent corrects
Je ne vois vraiment pas ce qui ne marche pas, j'ai mis mon code en pièce jointe.
Le laplacien discret est la fonction "M" dans le fichier "fonction_in_curve"
-----