Comment calculer une approximation de la dérivée second d'une fonction P1?
calculer dx(u) et dy(u), ceux-ci sont constants sur les triangle,
Calculer ensuite une interpolation P1 de ces deux fonction par la formule
Notez que toutes les integrales sont locales et se calculent toutes en meme temps par assemblage sur les triangles.
Recommencer sur dxx(u) en derivant dxu puis reinterpoler pareillement
Pour la syntaxe, pas de problème: "adaptmesh(expression)". Donc l'implmentation dans DBI se fait comme pour plot. On pourra générer un fichier "in.mtr" pour stocker les valeur de la metrique suivant le format:
nb_vertex 1 metric[0] ... metric[nb_vertex-1]Le 1 pouvant devenir un 3 si on souhaite avoir des metriques non-isotropes.
Par ailleurs il faudra avoir 2 options d'appel de la fonction de triangulation: l'une sans maillage adaptatif qui ne prend comme entrée que le fichier géometrie et une option avec maillage adaptatif qui prend pour entrées un fichier maillage (ex mesh.mesh), un fichier de métrique (ex in.mtr) et le fichier geometrie (ex in.geom) Donc il faut modifier ppm.C par exemple en jouant sur la valeur de argc.
fmeshback="mesh.msh"; // background mesh
fmeshout="mesh.msh";
fmetrix="in.mtr"; // background metric
if(argc==1) { fmeshback=NULL; fmetrix = NULL;} //OP