load "iovtk" string filename1="image1"; // Noms des fichiers textes contenant les images string filename2="image2"; mesh Sh; { Sh=triangulate(filename1); //On introduit un maillage uniforme sur l'image 1 } fespace Vh(Sh,P1); // On souhaite travailler avec les éléments finis P1 sur Sh Vh u1=0; Vh u2=0; // On définit les images { { ifstream file(filename1);real xx,yy;for(int i=0;i> xx >>yy >> u1[][i]; } // import de l'image 1 (la syntaxe est celle du C++) { ifstream file(filename2);real xx,yy;for(int i=0;i> xx >>yy >> u2[][i]; } // import de l'image 2 } plot(u1,wait=1,fill=1); // on dessine u1 mesh Th=adaptmesh(Sh,u1,u2,err=1.e-3,nbvx=10000000); // on adapte le maillage aux images 1 et 2 plot(Th,wait=1); // On dessine le maillage obtenu.