Chapter 1 parallel sparse solver - freeFEM.org

18 CHAPTER 1. PARALLEL SPARSE SOLVER
9: iparm[1]=1; // PCG as Krylov method
10:iparm[4]=0; // Matrix are symmetric
11:iparm[5]=1; // Pattern are also symmetric
12: iparm[9]=1; // Scale matrix
13:dparm[0]=1e-13; // Tolerance to convergence
14: dparm[1]=5e-4; // Threshold in ILUT
15: dparm[2]=5e-4; // Threshold for Schur preconditionner
16: mesh Th2=square(nn,nn);
17: fespace Vh2(Th2,P2);
18: Vh2 ux,uz,p2;
19: int[int] rup=[0,2], rdown=[0,1], rmid=[1,1,2,1,3,1,4,1];
20:real zmin=0,zmax=1;
21: mesh3 Th=buildlayers(Th2,nn,
zbound=[zmin,zmax],
reffacemid=rmid,
reffaceup = rup,
reffacelow = rdown);
22: savemesh(Th,"copie.mesh");
23: mesh3 Th3("copie.mesh");
24: fespace Vh(Th,P2);
25: func ue = 2*x*x + 3*y*y + 4*z*z + 5*x*y+6*x*z+1;
26: func uex= 4*x+ 5*y+6*z;
27: func uey= 6*y + 5*x;
28: func uez= 8*z +6*x;
29: func f= -18. ;
30: Vh uhe = ue; //
31: cout