massp0.edp

FreeFem++可以生成高质量的有限元网格。可以用于流体力学

mesh Th=square(2,2);
fespace Mh(Th,P0);
varf a(u,v)=int2d(Th)(u*v);
matrix A=a(Mh,Mh);
cout << A << endl;