-- FreeFem++ v4.6 (Thu Apr 2 15:47:38 CEST 2020 - git v4.6) Load: lg_fem lg_mesh lg_mesh3 eigenvalue 1 : // Location: 2 : // 3 : // http://people.sc.fsu.edu/~jburkardt/freefem_src/mshmet/mshmet.edp 4 : // 5 : // Discussion: 6 : // 7 : // The domain is a 3D version of the L shaped region, a cube with 8 : // a corner removed. 9 : // 10 : // Morice's mshmet plugin computes an isotropic adapted mesh. 11 : // 12 : // Licensing: 13 : // 14 : // This code is distributed under the GNU LGPL license. 15 : // 16 : // Modified: 17 : // 18 : // 26 May 2020 19 : // 20 : // Reference: 21 : // 22 : // Frederic Hecht, 23 : // New development in FreeFem++, 24 : // Journal of Numerical Mathematics, 25 : // Volume 20, Number 3-4, 2012, pages 251-265. 26 : // 27 : cout << "\n"; 28 : cout << "mshmet:\n"; 29 : cout << " FreeFem++ version\n"; 30 : 31 : load "msh3" 32 : load "tetgen" load: tetgen 33 : load "mshmet" load: mshmet 34 : load "medit" 35 : 36 : int n = 6; 37 : int[int] l1 = [1,1,1,1],l01=[0,1],l11=[1,1]; 38 : // 39 : // Label numbering for boundary condition. 40 : // 41 : mesh3 Th3 = buildlayers ( square ( n, n, region = 0, label = l1 ), 42 : n, zbound = [0,1], labelmid = l11, labelup = l01, 43 : labeldown = l01 ); 44 : // 45 : // Remove the [0.5,1.0]^3 cube. 46 : // 47 : Th3 = trunc ( Th3, ( x < 0.5 ) | ( y < 0.5 ) | ( z < 0.5 ), 48 : label = 1 ); 49 : 50 : fespace Vh ( Th3, P1 ); 51 : Vh u, v, usol; 52 : 53 : macro Grad(u) [ dx(u), dy(u), dz(u) ] ) // EOM 54 : 55 : solve Poisson ( u, v, solver = CG ) = 56 : int3d ( Th3 ) ( 1.0e-6*u*v + Grad(u) [ dx(u), dy(u), dz(u) ] ' * Grad(v) [ dx(v), dy(v), dz(v) ] ) 57 : - int3d ( Th3 ) ( ( z - y / 2 - x / 2 ) * v ); 58 : 59 : real errm = 1.0e-02; 60 : 61 : for ( int ii = 0; ii < 5; ii++ ) 62 : { 63 : Vh h; 64 : h[] = mshmet ( Th3, u, normalization = 1, aniso = 0, nbregul = 1, 65 : hmin = 1.0E-03, hmax = 0.3, err = errm ); 66 : 67 : errm = errm * 0.6; 68 : 69 : Th3 = tetgreconstruction ( Th3, switch = "raAQ", 70 : sizeofvolume = h*h*h / 6.0 ); 71 : Poisson; 72 : 73 : plot ( u, wait = 1, nbiso = 15, ps = "fuck.ps" ); 74 : } 75 : // 76 : // Invoke medit to display the solution. 77 : // 78 : medit ( "U3", Th3, u ); 79 : // 80 : // Terminate. 81 : // 82 : cout << "\n"; 83 : cout << "mshmet\n"; 84 : cout << " Normal end of execution.\n"; 85 : 86 : sizestack + 1024 =2624 ( 1600 ) mshmet: FreeFem++ version -- Square mesh : nb vertices =49 , nb triangles = 72 , nb boundary edges 24 -- FESpace: Nb of Nodes 316 Nb of DoF 316 GC: converge in 27 g=2.27536e-16 rho= 1.19229 gamma= 0.248966 -- Solve : min -0.101113 max 0.0873718 -- MSHMET, Release 3.0a (Feb. 15, 2010) Copyright (c) LJLL, 2007-10 Mer 25 avr 2012 09:30:16 CEST NUMBER OF GIVEN VERTICES 316 NUMBER OF GIVEN TETRAHEDRA 1134 NUMBER OF GIVEN DATA 316 -- DATA READING COMPLETED. 0.04 sec. &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& MODULE MSHMET-LJLL : 3.0a (Feb. 15, 2010) &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& -- PHASE 1 : ANALYSIS -- PHASE 1 COMPLETED. 0.04 sec. -- PHASE 2 : METRIC -- PHASE 2 COMPLETED. 0.04 sec. &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& END OF MODULE MSHMET &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& -- CPU REQUIREMENTS in/out 33.17 % 1. calls, 0.04 sec/call analysis 33.24 % 1. calls, 0.04 sec/call metric 33.60 % 1. calls, 0.04 sec/call total 100.00 % 1. calls, 0.12 sec/call ELAPSED TIME 0.12 SEC. (0.12) &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& END OF MODULE MSHMET &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& refine tetgen: Vertex Triangle Border 316 1134 432 Before reconstruction: nbhole=0 nbregion=0 3D RemplissageSurf3D:: Vertex triangle2 border 316 1134 432 tetgenio: vertex tetrahedrons lecture des facettes tetgen: before tetrahedralize( , &in, &out); numberof regions 0 numberof hole 0 tetgen: after tetrahedralize( , &in, &out); Th3 :: Vertex Element Border :: 500 1751 756 Finish Mesh3 tetgen :: Vertex, Element, Border500 1751 756 FreeFEM: End check mesh given by tetgen finish reconstruction FreeFEM: End check mesh given by tetgen -- FESpace: Nb of Nodes 500 Nb of DoF 500 GC: converge in 56 g=3.02957e-16 rho= 1.46253 gamma= 0.49334 -- Solve : min -0.101373 max 0.0908572 Plot:: Sorry no ps version for this type of plot 6 -- MSHMET, Release 3.0a (Feb. 15, 2010) Copyright (c) LJLL, 2007-10 Mer 25 avr 2012 09:30:16 CEST NUMBER OF GIVEN VERTICES 500 NUMBER OF GIVEN TETRAHEDRA 1751 NUMBER OF GIVEN DATA 500 -- DATA READING COMPLETED. 0.07 sec. &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& MODULE MSHMET-LJLL : 3.0a (Feb. 15, 2010) &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& -- PHASE 1 : ANALYSIS -- PHASE 1 COMPLETED. 0.07 sec. -- PHASE 2 : METRIC -- PHASE 2 COMPLETED. 0.07 sec. &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& END OF MODULE MSHMET &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& -- CPU REQUIREMENTS in/out 32.99 % 1. calls, 0.07 sec/call analysis 33.07 % 1. calls, 0.07 sec/call metric 33.94 % 1. calls, 0.07 sec/call total 100.00 % 1. calls, 0.21 sec/call ELAPSED TIME 0.21 SEC. (0.21) &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& END OF MODULE MSHMET &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& refine tetgen: Vertex Triangle Border 500 1751 756 Before reconstruction: nbhole=0 nbregion=0 3D RemplissageSurf3D:: Vertex triangle2 border 500 1751 756 tetgenio: vertex tetrahedrons lecture des facettes tetgen: before tetrahedralize( , &in, &out); numberof regions 0 numberof hole 0 tetgen: after tetrahedralize( , &in, &out); Th3 :: Vertex Element Border :: 1000 3918 1394 Finish Mesh3 tetgen :: Vertex, Element, Border1000 3918 1394 FreeFEM: End check mesh given by tetgen finish reconstruction FreeFEM: End check mesh given by tetgen -- FESpace: Nb of Nodes 1000 Nb of DoF 1000 GC: converge in 97 g=2.36864e-16 rho= 1.25369 gamma= 0.514847 -- Solve : min -0.0747503 max 0.115848 Plot:: Sorry no ps version for this type of plot 6 -- MSHMET, Release 3.0a (Feb. 15, 2010) Copyright (c) LJLL, 2007-10 Mer 25 avr 2012 09:30:16 CEST NUMBER OF GIVEN VERTICES 1000 NUMBER OF GIVEN TETRAHEDRA 3918 NUMBER OF GIVEN DATA 1000 -- DATA READING COMPLETED. 0.14 sec. &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& MODULE MSHMET-LJLL : 3.0a (Feb. 15, 2010) &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& -- PHASE 1 : ANALYSIS -- PHASE 1 COMPLETED. 0.14 sec. -- PHASE 2 : METRIC -- PHASE 2 COMPLETED. 0.15 sec. &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& END OF MODULE MSHMET &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& -- CPU REQUIREMENTS in/out 33.00 % 1. calls, 0.14 sec/call analysis 33.09 % 1. calls, 0.14 sec/call metric 33.91 % 1. calls, 0.15 sec/call total 100.00 % 1. calls, 0.43 sec/call ELAPSED TIME 0.43 SEC. (0.43) &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& END OF MODULE MSHMET &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& refine tetgen: Vertex Triangle Border 1000 3918 1394 Before reconstruction: nbhole=0 nbregion=0 3D RemplissageSurf3D:: Vertex triangle2 border 1000 3918 1394 tetgenio: vertex tetrahedrons lecture des facettes tetgen: before tetrahedralize( , &in, &out); numberof regions 0 numberof hole 0 tetgen: after tetrahedralize( , &in, &out); Th3 :: Vertex Element Border :: 2519 11119 2988 Finish Mesh3 tetgen :: Vertex, Element, Border2519 11119 2988 FreeFEM: End check mesh given by tetgen finish reconstruction FreeFEM: End check mesh given by tetgen -- FESpace: Nb of Nodes 2519 Nb of DoF 2519 GC: converge in 277 g=9.67044e-17 rho= 1.32481 gamma= 0.698231 -- Solve : min -0.0933285 max 0.0963402 Plot:: Sorry no ps version for this type of plot 6 -- MSHMET, Release 3.0a (Feb. 15, 2010) Copyright (c) LJLL, 2007-10 Mer 25 avr 2012 09:30:16 CEST NUMBER OF GIVEN VERTICES 2519 NUMBER OF GIVEN TETRAHEDRA 11119 NUMBER OF GIVEN DATA 2519 -- DATA READING COMPLETED. 0.33 sec. &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& MODULE MSHMET-LJLL : 3.0a (Feb. 15, 2010) &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& -- PHASE 1 : ANALYSIS -- PHASE 1 COMPLETED. 0.33 sec. -- PHASE 2 : METRIC -- PHASE 2 COMPLETED. 0.34 sec. &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& END OF MODULE MSHMET &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& -- CPU REQUIREMENTS in/out 32.88 % 1. calls, 0.33 sec/call analysis 33.02 % 1. calls, 0.33 sec/call metric 34.10 % 1. calls, 0.34 sec/call total 100.00 % 1. calls, 1.00 sec/call ELAPSED TIME 1.00 SEC. (1.00) &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& END OF MODULE MSHMET &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& refine tetgen: Vertex Triangle Border 2519 11119 2988 Before reconstruction: nbhole=0 nbregion=0 3D RemplissageSurf3D:: Vertex triangle2 border 2519 11119 2988 tetgenio: vertex tetrahedrons lecture des facettes tetgen: before tetrahedralize( , &in, &out); numberof regions 0 numberof hole 0 tetgen: after tetrahedralize( , &in, &out); Th3 :: Vertex Element Border :: 6938 33865 6526 Finish Mesh3 tetgen :: Vertex, Element, Border6938 33865 6526 FreeFEM: End check mesh given by tetgen finish reconstruction FreeFEM: End check mesh given by tetgen -- FESpace: Nb of Nodes 6938 Nb of DoF 6938 GC: converge in 238 g=5.41763e-17 rho= 1.08687 gamma= 0.758779 -- Solve : min -0.189084 max 0.00408403 Plot:: Sorry no ps version for this type of plot 6 -- MSHMET, Release 3.0a (Feb. 15, 2010) Copyright (c) LJLL, 2007-10 Mer 25 avr 2012 09:30:16 CEST NUMBER OF GIVEN VERTICES 6938 NUMBER OF GIVEN TETRAHEDRA 33865 NUMBER OF GIVEN DATA 6938 -- DATA READING COMPLETED. 0.79 sec. &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& MODULE MSHMET-LJLL : 3.0a (Feb. 15, 2010) &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& -- PHASE 1 : ANALYSIS -- PHASE 1 COMPLETED. 0.79 sec. -- PHASE 2 : METRIC -- PHASE 2 COMPLETED. 0.82 sec. &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& END OF MODULE MSHMET &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& -- CPU REQUIREMENTS in/out 32.86 % 1. calls, 0.79 sec/call analysis 33.01 % 1. calls, 0.79 sec/call metric 34.14 % 1. calls, 0.82 sec/call total 100.00 % 1. calls, 2.39 sec/call ELAPSED TIME 2.39 SEC. (2.39) &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& END OF MODULE MSHMET &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& refine tetgen: Vertex Triangle Border 6938 33865 6526 Before reconstruction: nbhole=0 nbregion=0 3D RemplissageSurf3D:: Vertex triangle2 border 6938 33865 6526 tetgenio: vertex tetrahedrons lecture des facettes tetgen: before tetrahedralize( , &in, &out); numberof regions 0 numberof hole 0 tetgen: after tetrahedralize( , &in, &out); Th3 :: Vertex Element Border :: 7607 37027 7220 Finish Mesh3 tetgen :: Vertex, Element, Border7607 37027 7220 FreeFEM: End check mesh given by tetgen finish reconstruction FreeFEM: End check mesh given by tetgen -- FESpace: Nb of Nodes 7607 Nb of DoF 7607 GC: converge in 239 g=4.9243e-17 rho= 1.46535 gamma= 0.760915 -- Solve : min -0.190787 max 0.00229083 Plot:: Sorry no ps version for this type of plot 6 -- Medit, Release 3.0a (Nov. 30, 2007) Copyright (c) LJLL, 1999-2007. compiled: Thu Apr 2 15:52:17 CEST 2020 (with ff++ 4.6). medit with binary version of popen : mesh(es) and solution(s) mesh_name= U3 Loading data file(s) End of mesh .sol: Dimension 3 (mesh)3 (lecture)1 Input seconds: 0.01 medit1() Building scene(s) Creating scene 1 Loading default options Scene seconds: 0.05 Rendering scene(s) Total running seconds: 0.07 Thank you for using Medit. version de medit ffmedit -popen -filebin -addsol 1 U3 mshmet Normal end of execution. times: compile 0.007603s, execution 1.33116s, mpirank:0 ######## We forget of deleting 6 Nb pointer, 0Bytes , mpirank 0, memory leak =40064 CodeAlloc : nb ptr 3733, size :489792 mpirank: 0 Ok: Normal End