-- FreeFem++ v4.6 (Thu Apr  2 15:47:38 CEST 2020 - git v4.6)
 Load: lg_fem lg_mesh lg_mesh3 eigenvalue 
    1 : //  Discussion:
    2 : //
    3 : //    Solve a simple boundary value problem on the unit square.
    4 : //    Repeat the process after halving h.
    5 : //    Report the errors and error convergence rates.
    6 : //
    7 : //  Location:
    8 : //
    9 : //    https://people.sc.fsu.edu/~jburkardt/freefem_src/convergence/convergence.edp  
   10 : //
   11 : //  Licensing:
   12 : //
   13 : //    This code is distributed under the GNU LGPL license.
   14 : //
   15 : //  Modified:
   16 : //
   17 : //    11 March 2021
   18 : //
   19 : //  Author:
   20 : //
   21 : //    John Burkardt
   22 : //
   23 : cout << endl;
   24 : cout << "convergence:" << endl;
   25 : cout << "  FreeFem++ version." << endl;
   26 : cout << "  Solve a boundary value problem in a square." << endl;
   27 : cout << "  Repeatedly refine the mesh, and report the L2 e
  ... : rror norm." << endl;
   28 : //
   29 : //  Declare some convergence variables here.
   30 : //
   31 : real a;
   32 : real aold;
   33 : real arate;
   34 : real b;
   35 : real bold;
   36 : real brate;
   37 : real c;
   38 : real cold;
   39 : real crate;
   40 : real h;
   41 : real hold;
   42 : //
   43 : //  Loop for h = 1/8, 1/16, ...
   44 : //
   45 : for ( int log2 = 3; log2 <= 7; log2++ )
   46 : {
   47 : //
   48 : //  n: number of elements in each spatial direction.
   49 : //
   50 :   int n = pow ( 2.0, log2 );
   51 : //
   52 : //  h: characteristic length of each element.
   53 : //
   54 :   if ( 3 < log2 )
   55 :   {
   56 :     hold = h;
   57 :   }
   58 :   h = 1.0 / n;
   59 : //
   60 : //  Th: the triangulation of the region.
   61 : //
   62 :   mesh Th = square ( n, n );
   63 : //
   64 : //  Plot the mesh the first time.
   65 : //
   66 :   if ( log2 == 3 )
   67 :   {
   68 :     plot ( Th, wait = False The Identifier False does not exist 

 Error line number 68, in file convergence.edp, before  token False

  current line = 68
Compile error : 
	line number :68, False
error Compile error : 
	line number :68, False
 code = 1 mpirank: 0