正方メッシュデータの読み込み

0.0 0.0 10.0

1.0 0.0 10.1

...

10.0 10.0 10.9 

int n=11; mesh Th=square(n-1,n-1); plot(Th,wait=1); 

fespace Vh(Th,P1);

Vh u1;

real x1, y1; 

{

   ifstream ff("sample.dat");

   for(int i=0; i> Th.nv; i++){

     ff >> x1 >> y1  >> u1[][i];

     } 

}

plot(u1,wait=1); 

int n=11; 

real xmin=-1.0, xmax=1.0; 

real ymin=-1.0, ymax=1.0;

mesh Th1=square(nx-1, ny-1,[xmin+(xmax-xmin)*x,ymin+(ymax-ymin)*y]);

plot(Th,wait=1);

fespace Vh(Th,P1);

Vh u1;

real x1, y1;

{

   ifstream ff("sample.dat");

   for(int i=0; i> Th.nv; i++){

     ff >> x1 >> y1  >> u1[][i];

     }

}

plot(f1,cmm="orig",wait=1,grey=1,value=1);  // New mesh

real nx2=100, ny2=100;

mesh Th=square(nx2, ny2,[xmin+(xmax-xmin)*x,ymin+(ymax-ymin)*y]);

fespace Vh(Th,P1);

matrix IV = interpolate(Vh, Vh1);

Vh f;

f[]=IV*f1[]; // interpolation : Th1 -> Th

f=10*f;

plot(Th, f, dim=2,cmm="interpolate",grey=0,value=1);