Ceci est une ancienne révision du document !

FreeFem++

• installé dans /local/apps/freefem++-3.13-2
• un exemple

[henry@octopus ~]$ /local/apps/freefem++-3.13-2/bin/FreeFem++ /local/apps/src/freefem++-3.13-2/examples++-tutorial/LaplaceP1.edp 
 EXEC of the plot  : ffglut
-- FreeFem++ v  3.130002 (date Mar 28 jui 2011 09:22:37 CEST)
 Load: lg_fem lg_mesh lg_mesh3 eigenvalue 
    1 :  mesh Th=square(10,10);
    2 :  fespace Vh(Th,P1);     // P1 FE space
    3 :  Vh uh,vh;              // unkown and test function. 
    4 :  func f=1;                 //  right hand side function 
    5 :  func g=0;                 //  boundary condition function
    6 :  
    7 :  problem laplace(uh,vh,solver=GMRES,tgv=1e5) =                    //  definion of  the problem 
    8 :     int2d(Th)( dx(uh)*dx(vh) + dy(uh)*dy(vh) ) //  bilinear form
    9 :   + int1d(Th,1)( uh*vh) 
   10 :   - int1d(Th,1)( vh)
   11 :   - int2d(Th)( f*vh )                          //  linear form
   12 :   + on(2,3,4,uh=g) ;                      //  boundary condition form
   13 : 
   14 :   laplace; // solve the problem plot(uh); // to see the result
   15 :   plot(uh,ps="LaplaceP1.eps",value=true);
   16 :  sizestack + 1024 =1896  ( 872 )

  -- Square mesh : nb vertices  =121 ,  nb triangles = 200 ,  nb boundary edges 40
GMRES converges: 27 1.116e-07 0.154071 7.24338e-07 < 1e-06
  -- Solve : 
          min 1.62405e-12  max 0.361312
times: compile 0s, execution 0s,  mpirank:0
Err ReadOnePlot 1
 CodeAlloc : nb ptr  2299,  size :305528 mpirank: 0
Bien: On a fini Normalement