-- FreeFem++ v4.6 (Thu Apr 2 15:47:38 CEST 2020 - git v4.6) Load: lg_fem lg_mesh lg_mesh3 eigenvalue 1 : // migration.edp 2 : // 3 : // Discussion: 4 : // 5 : // Migration and proliferation of biological cells. 6 : // 7 : // The region is a 3x2 rectangle, with periodic conditions 8 : // at the left and right ends. 9 : // 10 : // Location: 11 : // 12 : // http://people.sc.fsu.edu/~jburkardt/freefem++/migration/migration.edp 13 : // 14 : // Modified: 15 : // 16 : // 19 June 2015 17 : // 18 : // Author: 19 : // 20 : // Florian De Vuyst 21 : // 22 : // Reference: 23 : // 24 : // Florian De Vuyst, 25 : // Numerical modeling of transport problems using freefem++ software - 26 : // with examples in biology, CFD, traffic flow and energy transfer, 27 : // HAL id: cel-00842234 28 : // https://cel.archives-ouvertes.fr/cel-00842234 29 : // 30 : real dt = 0.05; 31 : real uf = 1.0; 32 : real rhoc = 100.0; 33 : // 34 : // Define the mesh. 35 : // 36 : real Lx = 3.0; 37 : real Ly = 2.0; 38 : mesh Th = square ( 60, 40, [x*Lx,y*Ly] ); 39 : // 40 : // Define the finite element spaces using periodicity. 41 : // 42 : fespace Uh ( Th, P2, periodic=[[3,x],[1,x]] ); 43 : fespace Vh ( Th, P1, periodic=[[3,x],[1,x]] ); 44 : // 45 : // Declare variables. 46 : // 47 : Uh rho, lrho, rhoold, rhoh, c, ch, cold; 48 : Vh u1, u2, u, n1, n2, v1, v2, rhop1, cp1; 49 : // 50 : // Set the function rho(x,y). 51 : // 52 : rho = rhoc * exp ( -40.0*(x-Lx/4.0)^2 - 40.0*(y-Ly/2)^2 ) 53 : + rhoc * exp ( -40.0*(x-3.0*Lx/4.0)^2 - 40.0*(y-Ly/2)^2 ) 54 : + rhoc * exp ( -40.0*(x-0.55*Lx)^2 - 40.0*(y-Ly/2)^2 ); 55 : // 56 : // Adapt the mesh according to rho and periodicity. 57 : // 58 : Th = adaptmesh ( Th, rho, periodic=[[3,x],[1,x]] ); 59 : // 60 : // Plot the adapted mesh. 61 : // 62 : plot ( Th, ps = "migration_mesh_initial.ps" ); 63 : // 64 : // Transfer rho to the new mesh. 65 : // 66 : rho = rho; 67 : c = rho / rhoc; 68 : rhoold = rho; 69 : cold = c; 70 : // 71 : // Plot the initial density. 72 : // 73 : rhop1 = rho; 74 : plot ( Th, rhop1, nbiso=50, fill=0, value=1, ps="migration_rho_initial.ps" ); 75 : // 76 : // Define the weak form of the migration equations. 77 : // 78 : problem migration ( [rho, c], [rhoh, ch]) = 79 : int2d ( Th ) ( rho*rhoh/dt ) 80 : - int2d ( Th ) ( convect([v1,v2], -dt,rhoold)*rhoh/dt ) 81 : + int2d ( Th ) ( dx(v1)*rho*rhoh+dy(v2)*rho*rhoh ) 82 : + int2d ( Th ) ( 0.01*dx(rho)*dx(rhoh)+0.01*dy(rho)*dy(rhoh) ) 83 : - int2d ( Th ) ( 0.01*rho*(rhoc-rhoold)*rhoh ) 84 : + int2d ( Th ) ( c*ch/dt ) 85 : - int2d ( Th ) ( cold*ch/dt ) 86 : + int2d ( Th ) ( dx(c)*dx(ch)+dy(c)*dy(ch) ) 87 : - int2d ( Th ) ( 10.0*(rho/rhoc-c)*ch ); 88 : // 89 : // Time loop. 90 : // 91 : for ( int it = 0; it < 20; it++ ) 92 : { 93 : for ( int substep = 0; substep < 2; substep++ ) 94 : { 95 : u1 = - dx ( cold ); 96 : u2 = - dy ( cold ); 97 : v1 = 0.5 * u1 * rhoold / rhoc; 98 : v2 = 0.5 * u2 * rhoold / rhoc; 99 : 100 : migration; 101 : 102 : Th = adaptmesh ( Th, rho, periodic=[[3,x],[1,x]] ); 103 : 104 : rho=rho; 105 : c=c; 106 : rhoold = rho; 107 : cold = c; 108 : } 109 : // 110 : // Display the current state. 111 : // 112 : rhop1 = rho; 113 : cp1 = c; 114 : plot ( Th, rhop1, nbiso=50, fill=0, value=1 ); 115 : } 116 : // 117 : // Plot the final mesh. 118 : // 119 : plot ( Th, ps="migration_mesh_final.ps" ); 120 : // 121 : // 122 : // Plot the final state. 123 : // 124 : plot ( Th, rhop1, nbiso=50, fill=0, value=1, ps="migration_rho_final.ps" ); 125 : // 126 : // Terminate. 127 : // 128 : cout << "\n"; 129 : cout << "MIGRATION:\n"; 130 : cout << " Normal end of execution.\n"; 131 : sizestack + 1024 =12184 ( 11160 ) -- Square mesh : nb vertices =2501 , nb triangles = 4800 , nb boundary edges 200 -- adaptmesh Regulary: Nb triangles 4184 , h min 0.00821544 , h max 0.673875 area = 6 , M area = 1888.98 , M area/( |Khat| nt) 1.04264 infiny-regularity: min 0.0854602 max 2.44384 anisomax 11.6093, beta max = 2.06694 min 0.72004 -- mesh: Nb of Triangles = 4184, Nb of Vertices 2104 kk 2 6 : -- Solve : min -6.7928e-05 max 117.274 min 0.00122453 max 0.326769 -- adaptmesh Regulary: Nb triangles 14530 , h min 0.00383049 , h max 0.851889 area = 6 , M area = 6340.68 , M area/( |Khat| nt) 1.00779 infiny-regularity: min 0.317327 max 2.42576 anisomax 9.46768, beta max = 1.3855 min 0.711912 -- mesh: Nb of Triangles = 14530, Nb of Vertices 7276 kk 2 6 : -- Solve : min -0.000442129 max 102.64 min 0.00334607 max 0.251032 -- adaptmesh Regulary: Nb triangles 14079 , h min 0.00326008 , h max 0.716436 area = 6 , M area = 6526.87 , M area/( |Khat| nt) 1.07061 infiny-regularity: min 0.30911 max 2.49005 anisomax 11.185, beta max = 1.40992 min 0.676495 -- mesh: Nb of Triangles = 14079, Nb of Vertices 7052 kk 2 6 : -- Solve : min -0.000441292 max 84.2349 min 0.00580619 max 0.243774 -- adaptmesh Regulary: Nb triangles 12552 , h min 0.00337164 , h max 0.621118 area = 6 , M area = 5689.84 , M area/( |Khat| nt) 1.04686 infiny-regularity: min 0.252767 max 2.4781 anisomax 17.1219, beta max = 1.55861 min 0.659329 -- mesh: Nb of Triangles = 12552, Nb of Vertices 6289 kk 2 6 : -- Solve : min -0.000614512 max 71.5845 min 0.00814333 max 0.244429 -- adaptmesh Regulary: Nb triangles 11244 , h min 0.00386548 , h max 0.582638 area = 6 , M area = 5171.9 , M area/( |Khat| nt) 1.06225 infiny-regularity: min 0.348698 max 2.43199 anisomax 12.7184, beta max = 1.36904 min 0.685688 -- mesh: Nb of Triangles = 11244, Nb of Vertices 5635 kk 2 6 : -- Solve : min -0.000858661 max 63.5281 min 0.0100205 max 0.246831 -- adaptmesh Regulary: Nb triangles 10628 , h min 0.00501437 , h max 0.631605 area = 6 , M area = 4835.06 , M area/( |Khat| nt) 1.05063 infiny-regularity: min 0.303166 max 2.38865 anisomax 12.1152, beta max = 1.50112 min 0.686403 -- mesh: Nb of Triangles = 10628, Nb of Vertices 5328 kk 2 6 : -- Solve : min -0.00113775 max 58.5974 min 0.0114927 max 0.249729 -- adaptmesh Regulary: Nb triangles 10066 , h min 0.00573217 , h max 0.636616 area = 6 , M area = 4527.92 , M area/( |Khat| nt) 1.03882 infiny-regularity: min 0.289171 max 2.49688 anisomax 10.573, beta max = 1.53399 min 0.679723 -- mesh: Nb of Triangles = 10066, Nb of Vertices 5047 kk 2 6 : -- Solve : min -0.000651053 max 55.4328 min 0.0127097 max 0.252853 -- adaptmesh Regulary: Nb triangles 9508 , h min 0.00616478 , h max 0.582947 area = 6 , M area = 4242.8 , M area/( |Khat| nt) 1.03053 infiny-regularity: min 0.252424 max 2.43666 anisomax 10.9308, beta max = 1.66618 min 0.697881 -- mesh: Nb of Triangles = 9508, Nb of Vertices 4769 kk 2 6 : -- Solve : min -0.00078041 max 53.3045 min 0.0138391 max 0.256109 -- adaptmesh Regulary: Nb triangles 9029 , h min 0.00783382 , h max 0.501573 area = 6 , M area = 4000.58 , M area/( |Khat| nt) 1.02325 infiny-regularity: min 0.226527 max 2.48353 anisomax 11.9663, beta max = 1.64231 min 0.69446 -- mesh: Nb of Triangles = 9029, Nb of Vertices 4531 kk 2 6 : -- Solve : min -0.000699449 max 51.8073 min 0.0147802 max 0.259483 -- adaptmesh Regulary: Nb triangles 8640 , h min 0.00789619 , h max 0.522256 area = 6 , M area = 3785.7 , M area/( |Khat| nt) 1.01189 infiny-regularity: min 0.258532 max 2.64644 anisomax 9.43904, beta max = 1.55655 min 0.656429 -- mesh: Nb of Triangles = 8640, Nb of Vertices 4337 kk 2 6 : -- Solve : min -0.00058835 max 50.7569 min 0.0156655 max 0.263091 -- adaptmesh Regulary: Nb triangles 8259 , h min 0.00917386 , h max 0.534203 area = 6 , M area = 3592.84 , M area/( |Khat| nt) 1.00464 infiny-regularity: min 0.256077 max 2.59576 anisomax 9.64279, beta max = 1.54026 min 0.655593 -- mesh: Nb of Triangles = 8259, Nb of Vertices 4146 kk 2 6 : -- Solve : min -0.000906277 max 50.0866 min 0.016532 max 0.26687 -- adaptmesh Regulary: Nb triangles 7797 , h min 0.00902247 , h max 0.503552 area = 6 , M area = 3432.51 , M area/( |Khat| nt) 1.01668 infiny-regularity: min 0.338344 max 2.75844 anisomax 9.9766, beta max = 1.40182 min 0.653845 -- mesh: Nb of Triangles = 7797, Nb of Vertices 3915 kk 2 6 : -- Solve : min -0.00268222 max 49.6613 min 0.0173757 max 0.270798 -- adaptmesh Regulary: Nb triangles 7525 , h min 0.00997442 , h max 0.492841 area = 6 , M area = 3280.06 , M area/( |Khat| nt) 1.00664 infiny-regularity: min 0.305369 max 2.51608 anisomax 11.2621, beta max = 1.52086 min 0.672255 -- mesh: Nb of Triangles = 7525, Nb of Vertices 3778 kk 2 6 : -- Solve : min -0.000671962 max 49.3625 min 0.0182275 max 0.274906 -- adaptmesh Regulary: Nb triangles 7289 , h min 0.0106883 , h max 0.542028 area = 6 , M area = 3150.48 , M area/( |Khat| nt) 0.998177 infiny-regularity: min 0.321296 max 2.63427 anisomax 10.2891, beta max = 1.53264 min 0.70814 -- mesh: Nb of Triangles = 7289, Nb of Vertices 3659 kk 2 6 : -- Solve : min -0.00109957 max 49.1651 min 0.0191011 max 0.279145 -- adaptmesh Regulary: Nb triangles 7031 , h min 0.0108904 , h max 0.568813 area = 6 , M area = 3049.43 , M area/( |Khat| nt) 1.00162 infiny-regularity: min 0.273908 max 2.63506 anisomax 9.2399, beta max = 1.35646 min 0.676008 -- mesh: Nb of Triangles = 7031, Nb of Vertices 3531 kk 2 6 : -- Solve : min -0.0010269 max 49.0415 min 0.0199022 max 0.283537 -- adaptmesh Regulary: Nb triangles 6864 , h min 0.012124 , h max 0.536081 area = 6 , M area = 2957.42 , M area/( |Khat| nt) 0.995026 infiny-regularity: min 0.3212 max 2.60069 anisomax 9.26384, beta max = 1.39856 min 0.668783 -- mesh: Nb of Triangles = 6864, Nb of Vertices 3447 kk 2 6 : -- Solve : min -0.00100484 max 48.9844 min 0.0208192 max 0.288132 -- adaptmesh Regulary: Nb triangles 6694 , h min 0.0110674 , h max 0.548742 area = 6 , M area = 2873.84 , M area/( |Khat| nt) 0.991461 infiny-regularity: min 0.341672 max 2.44219 anisomax 9.23596, beta max = 1.42071 min 0.677707 -- mesh: Nb of Triangles = 6694, Nb of Vertices 3363 kk 2 6 : -- Solve : min -0.000862431 max 48.9884 min 0.0218221 max 0.292804 -- adaptmesh Regulary: Nb triangles 6546 , h min 0.0115735 , h max 0.545979 area = 6 , M area = 2805.13 , M area/( |Khat| nt) 0.989639 infiny-regularity: min 0.345755 max 2.78291 anisomax 8.35223, beta max = 1.34724 min 0.664419 -- mesh: Nb of Triangles = 6546, Nb of Vertices 3288 kk 2 6 : -- Solve : min -0.000866927 max 49.0304 min 0.022709 max 0.297579 -- adaptmesh Regulary: Nb triangles 6442 , h min 0.0114682 , h max 0.558185 area = 6 , M area = 2734.69 , M area/( |Khat| nt) 0.980361 infiny-regularity: min 0.299742 max 2.5971 anisomax 10.2891, beta max = 1.48026 min 0.670852 -- mesh: Nb of Triangles = 6442, Nb of Vertices 3238 kk 2 6 : -- Solve : min -0.000936798 max 49.0953 min 0.0239104 max 0.30245 -- adaptmesh Regulary: Nb triangles 6328 , h min 0.011613 , h max 0.498617 area = 6 , M area = 2675.82 , M area/( |Khat| nt) 0.976538 infiny-regularity: min 0.270968 max 2.53049 anisomax 9.58936, beta max = 1.42207 min 0.675516 -- mesh: Nb of Triangles = 6328, Nb of Vertices 3182 kk 2 6 : -- Solve : min -0.00105059 max 49.1927 min 0.0250232 max 0.307429 -- adaptmesh Regulary: Nb triangles 6177 , h min 0.0126719 , h max 0.573954 area = 6 , M area = 2622.64 , M area/( |Khat| nt) 0.980529 infiny-regularity: min 0.278809 max 2.65005 anisomax 9.36569, beta max = 1.4318 min 0.676467 -- mesh: Nb of Triangles = 6177, Nb of Vertices 3104 kk 2 6 : -- Solve : min -0.000805203 max 49.3122 min 0.0261406 max 0.312502 -- adaptmesh Regulary: Nb triangles 6065 , h min 0.0118891 , h max 0.475803 area = 6 , M area = 2577.21 , M area/( |Khat| nt) 0.981338 infiny-regularity: min 0.292247 max 2.42956 anisomax 10.4357, beta max = 1.37217 min 0.671574 -- mesh: Nb of Triangles = 6065, Nb of Vertices 3050 kk 2 6 : -- Solve : min -0.000956053 max 49.4604 min 0.0273474 max 0.317671 -- adaptmesh Regulary: Nb triangles 5954 , h min 0.013923 , h max 0.505375 area = 6 , M area = 2532.4 , M area/( |Khat| nt) 0.982252 infiny-regularity: min 0.336949 max 2.46223 anisomax 9.11561, beta max = 1.4267 min 0.68512 -- mesh: Nb of Triangles = 5954, Nb of Vertices 2995 kk 2 6 : -- Solve : min -0.000935552 max 49.6255 min 0.0286266 max 0.32284 -- adaptmesh Regulary: Nb triangles 5872 , h min 0.0140478 , h max 0.499362 area = 6 , M area = 2482.4 , M area/( |Khat| nt) 0.976303 infiny-regularity: min 0.292592 max 2.77382 anisomax 9.74376, beta max = 1.40735 min 0.673741 -- mesh: Nb of Triangles = 5872, Nb of Vertices 2955 kk 2 6 : -- Solve : min -0.00104094 max 49.8147 min 0.0299676 max 0.328101 -- adaptmesh Regulary: Nb triangles 5747 , h min 0.0136705 , h max 0.51476 area = 6 , M area = 2443 , M area/( |Khat| nt) 0.981707 infiny-regularity: min 0.3194 max 2.46472 anisomax 7.61368, beta max = 1.38649 min 0.693327 -- mesh: Nb of Triangles = 5747, Nb of Vertices 2892 kk 2 6 : -- Solve : min -0.00170098 max 50.0163 min 0.0313854 max 0.333439 -- adaptmesh Regulary: Nb triangles 5598 , h min 0.0135016 , h max 0.523149 area = 6 , M area = 2403.85 , M area/( |Khat| nt) 0.991686 infiny-regularity: min 0.326148 max 2.71716 anisomax 9.81146, beta max = 1.4276 min 0.663932 -- mesh: Nb of Triangles = 5598, Nb of Vertices 2817 kk 2 6 : -- Solve : min -0.00138893 max 50.2339 min 0.0328171 max 0.33888 -- adaptmesh Regulary: Nb triangles 5510 , h min 0.0153115 , h max 0.566644 area = 6 , M area = 2364.62 , M area/( |Khat| nt) 0.991082 infiny-regularity: min 0.278727 max 2.58506 anisomax 11.3316, beta max = 1.34182 min 0.627724 -- mesh: Nb of Triangles = 5510, Nb of Vertices 2774 kk 2 6 : -- Solve : min -0.00152461 max 50.4731 min 0.0343457 max 0.344459 -- adaptmesh Regulary: Nb triangles 5461 , h min 0.0161817 , h max 0.506467 area = 6 , M area = 2332.16 , M area/( |Khat| nt) 0.986248 infiny-regularity: min 0.296355 max 2.46155 anisomax 8.33587, beta max = 1.39127 min 0.709703 -- mesh: Nb of Triangles = 5461, Nb of Vertices 2751 kk 2 6 : -- Solve : min -0.0013594 max 50.7214 min 0.0359909 max 0.350008 -- adaptmesh Regulary: Nb triangles 5377 , h min 0.0156632 , h max 0.468657 area = 6 , M area = 2299.88 , M area/( |Khat| nt) 0.987791 infiny-regularity: min 0.280471 max 2.6538 anisomax 8.12542, beta max = 1.42378 min 0.68131 -- mesh: Nb of Triangles = 5377, Nb of Vertices 2710 kk 2 6 : -- Solve : min -0.00160907 max 50.9784 min 0.0376891 max 0.355614 -- adaptmesh Regulary: Nb triangles 5295 , h min 0.0155529 , h max 0.49335 area = 6 , M area = 2267.08 , M area/( |Khat| nt) 0.988782 infiny-regularity: min 0.272521 max 2.43163 anisomax 8.93372, beta max = 1.48749 min 0.68182 -- mesh: Nb of Triangles = 5295, Nb of Vertices 2669 kk 2 6 : -- Solve : min -0.00269408 max 51.244 min 0.0394513 max 0.361263 -- adaptmesh Regulary: Nb triangles 5207 , h min 0.0169928 , h max 0.463872 area = 6 , M area = 2233 , M area/( |Khat| nt) 0.990375 infiny-regularity: min 0.284829 max 2.70305 anisomax 9.4966, beta max = 1.43816 min 0.662361 -- mesh: Nb of Triangles = 5207, Nb of Vertices 2626 kk 2 6 : -- Solve : min -0.00279914 max 51.5322 min 0.041298 max 0.36694 -- adaptmesh Regulary: Nb triangles 5124 , h min 0.0170475 , h max 0.503824 area = 6 , M area = 2204.43 , M area/( |Khat| nt) 0.993544 infiny-regularity: min 0.284297 max 2.51025 anisomax 8.58228, beta max = 1.46548 min 0.692827 -- mesh: Nb of Triangles = 5124, Nb of Vertices 2585 kk 2 6 : -- Solve : min -0.00272129 max 51.8385 min 0.0432543 max 0.372523 -- adaptmesh Regulary: Nb triangles 5042 , h min 0.0164 , h max 0.52532 area = 6 , M area = 2177.18 , M area/( |Khat| nt) 0.997218 infiny-regularity: min 0.286485 max 2.47438 anisomax 10.1184, beta max = 1.55189 min 0.657737 -- mesh: Nb of Triangles = 5042, Nb of Vertices 2544 kk 2 6 : -- Solve : min -0.00240368 max 52.1535 min 0.0452895 max 0.378296 -- adaptmesh Regulary: Nb triangles 4980 , h min 0.0173561 , h max 0.486167 area = 6 , M area = 2147.23 , M area/( |Khat| nt) 0.995744 infiny-regularity: min 0.19986 max 2.74638 anisomax 9.96642, beta max = 1.67323 min 0.66529 -- mesh: Nb of Triangles = 4980, Nb of Vertices 2513 kk 2 6 : -- Solve : min -0.00230127 max 52.4763 min 0.0474163 max 0.384201 -- adaptmesh Regulary: Nb triangles 4898 , h min 0.0163958 , h max 0.445988 area = 6 , M area = 2124.75 , M area/( |Khat| nt) 1.00181 infiny-regularity: min 0.170366 max 2.41561 anisomax 11.328, beta max = 1.57503 min 0.697036 -- mesh: Nb of Triangles = 4898, Nb of Vertices 2473 kk 2 6 : -- Solve : min -0.00242449 max 52.8051 min 0.049641 max 0.390104 -- adaptmesh Regulary: Nb triangles 4812 , h min 0.0168617 , h max 0.460787 area = 6 , M area = 2098.54 , M area/( |Khat| nt) 1.00714 infiny-regularity: min 0.205297 max 2.54822 anisomax 10.4994, beta max = 1.69202 min 0.686806 -- mesh: Nb of Triangles = 4812, Nb of Vertices 2431 kk 2 6 : -- Solve : min -0.00306772 max 53.1405 min 0.051975 max 0.396016 -- adaptmesh Regulary: Nb triangles 4748 , h min 0.0180374 , h max 0.443084 area = 6 , M area = 2073.91 , M area/( |Khat| nt) 1.00874 infiny-regularity: min 0.246832 max 2.58175 anisomax 8.92382, beta max = 1.48188 min 0.694117 -- mesh: Nb of Triangles = 4748, Nb of Vertices 2399 kk 2 6 : -- Solve : min -0.0032278 max 53.4825 min 0.0544529 max 0.401961 -- adaptmesh Regulary: Nb triangles 4689 , h min 0.016027 , h max 0.430196 area = 6 , M area = 2050.54 , M area/( |Khat| nt) 1.00992 infiny-regularity: min 0.280729 max 2.37901 anisomax 9.55729, beta max = 1.36492 min 0.657474 -- mesh: Nb of Triangles = 4689, Nb of Vertices 2370 kk 2 6 : -- Solve : min -0.0025422 max 53.835 min 0.0570139 max 0.408025 -- adaptmesh Regulary: Nb triangles 4647 , h min 0.0165108 , h max 0.415208 area = 6 , M area = 2031.82 , M area/( |Khat| nt) 1.00974 infiny-regularity: min 0.241301 max 2.32127 anisomax 9.75751, beta max = 1.64776 min 0.680943 -- mesh: Nb of Triangles = 4647, Nb of Vertices 2351 kk 2 6 : -- Solve : min -0.00292385 max 54.1941 min 0.059817 max 0.414144 -- adaptmesh Regulary: Nb triangles 4568 , h min 0.0177287 , h max 0.382049 area = 6 , M area = 2010.69 , M area/( |Khat| nt) 1.01653 infiny-regularity: min 0.229423 max 2.5146 anisomax 10.9525, beta max = 1.58238 min 0.646068 -- mesh: Nb of Triangles = 4568, Nb of Vertices 2312 kk 2 6 : -- Solve : min -0.00242452 max 54.5602 min 0.062685 max 0.420297 -- adaptmesh Regulary: Nb triangles 4512 , h min 0.017769 , h max 0.406604 area = 6 , M area = 1983.33 , M area/( |Khat| nt) 1.01514 infiny-regularity: min 0.271838 max 2.48144 anisomax 13.7132, beta max = 1.46959 min 0.692471 -- mesh: Nb of Triangles = 4512, Nb of Vertices 2285 MIGRATION: Normal end of execution. times: compile 0.005875s, execution 161.614s, mpirank:0 ######## We forget of deleting 0 Nb pointer, 0Bytes , mpirank 0, memory leak =106128 CodeAlloc : nb ptr 4010, size :498088 mpirank: 0 Ok: Normal End