-- FreeFem++ v4.14 (mer. 06 mars 2024 16:59:04 CET - git v4.14-1-g2b2052ae) file : mshmet.edp Load: lg_fem lg_mesh lg_mesh3 eigenvalue 1 : // mshmet.edp 2 : // 3 : // Discussion: 4 : // 5 : // The domain is a 3D version of the L shaped region, a cube with 6 : // a corner removed. 7 : // 8 : // Morice's mshmet plugin computes an isotropic adapted mesh. 9 : // 10 : // Licensing: 11 : // 12 : // This code is distributed under the MIT license. 13 : // 14 : // Modified: 15 : // 16 : // 26 May 2020 17 : // 18 : // Reference: 19 : // 20 : // Frederic Hecht, 21 : // New development in FreeFem++, 22 : // Journal of Numerical Mathematics, 23 : // Volume 20, Number 3-4, 2012, pages 251-265. 24 : // 25 : cout << "\n"; 26 : cout << "mshmet():\n"; 27 : cout << " FreeFem++ version\n"; 28 : 29 : load "msh3" 30 : load "tetgen" load: tetgen 31 : load "mshmet" load: mshmet 32 : load "medit" 33 : 34 : int n = 6; 35 : int[int] l1 = [1,1,1,1],l01=[0,1],l11=[1,1]; 36 : // 37 : // Label numbering for boundary condition. 38 : // 39 : mesh3 Th3 = buildlayers ( square ( n, n, region = 0, label = l1 ), 40 : n, zbound = [0,1], labelmid = l11, labelup = l01, 41 : labeldown = l01 ); 42 : // 43 : // Remove the [0.5,1.0]^3 cube. 44 : // 45 : Th3 = trunc ( Th3, ( x < 0.5 ) | ( y < 0.5 ) | ( z < 0.5 ), 46 : label = 1 ); 47 : 48 : fespace Vh ( Th3, P1 ); 49 : Vh u, v, usol; 50 : 51 : macro Grad(u) [ dx(u), dy(u), dz(u) ] ) // EOM 52 : 53 : solve Poisson ( u, v, solver = CG ) = 54 : int3d ( Th3 ) ( 1.0e-6*u*v + Grad(u) [ dx(u), dy(u), dz(u) ] ' * Grad(v) [ dx(v), dy(v), dz(v) ] ) 55 : - int3d ( Th3 ) ( ( z - y / 2 - x / 2 ) * v ); 56 : 57 : real errm = 1.0e-02; 58 : 59 : for ( int ii = 0; ii < 5; ii++ ) 60 : { 61 : Vh h; 62 : h[] = mshmet ( Th3, u, normalization = true, aniso = false, nbregul = 1, 63 : hmin = 1.0E-03, hmax = 0.3, err = errm ); 64 : 65 : errm = errm * 0.6; 66 : 67 : Th3 = tetgreconstruction ( Th3, switch = "raAQ", 68 : sizeofvolume = h*h*h / 6.0 ); 69 : 70 : Poisson; 71 : // 72 : // A 3D plot can't be saved to a PostScript file. 73 : // 74 : plot ( u, wait = true, nbiso = 15 ); 75 : } 76 : // 77 : // Invoke medit to display the solution. 78 : // 79 : medit ( "U3", Th3, u ); 80 : // 81 : // Terminate. 82 : // 83 : cout << "\n"; 84 : cout << "mshmet():\n"; 85 : cout << " Normal end of execution.\n"; 86 : 87 : exit ( 0 ); 88 : 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 after 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. 06 mars 2024 16:53:00 CET NUMBER OF GIVEN VERTICES 316 NUMBER OF GIVEN TETRAHEDRA 1134 NUMBER OF GIVEN DATA 316 -- DATA READING COMPLETED. 0.20 sec. &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& MODULE MSHMET-LJLL : 3.0a (Feb. 15, 2010) &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& -- PHASE 1 : ANALYSIS -- PHASE 1 COMPLETED. 0.20 sec. -- PHASE 2 : METRIC -- PHASE 2 COMPLETED. 0.21 sec. &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& END OF MODULE MSHMET &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& -- CPU REQUIREMENTS in/out 32.92 % 1. calls, 0.20 sec/call analysis 33.06 % 1. calls, 0.20 sec/call metric 34.03 % 1. calls, 0.21 sec/call total 100.00 % 1. calls, 0.60 sec/call ELAPSED TIME 0.60 SEC. (0.60) &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& 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 :: 370 1446 498 Finish Mesh3 tetgen :: Vertex, Element, Border370 1446 498 FreeFEM: End check mesh given by tetgen finish reconstruction FreeFEM: End check mesh given by tetgen -- FESpace: Nb of Nodes 370 Nb of DoF 370 GC: converge after 57 g=3.82162e-16 rho= 1.37891 gamma= 0.670151 -- Solve : min -0.101861 max 0.0898152 Plot:: Sorry no ps version for this type of plot 6 -- MSHMET, Release 3.0a (Feb. 15, 2010) Copyright (c) LJLL, 2007-10 mer. 06 mars 2024 16:53:00 CET NUMBER OF GIVEN VERTICES 370 NUMBER OF GIVEN TETRAHEDRA 1446 NUMBER OF GIVEN DATA 370 -- DATA READING COMPLETED. 0.58 sec. &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& MODULE MSHMET-LJLL : 3.0a (Feb. 15, 2010) &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& -- PHASE 1 : ANALYSIS -- PHASE 1 COMPLETED. 0.58 sec. -- PHASE 2 : METRIC -- PHASE 2 COMPLETED. 0.60 sec. &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& END OF MODULE MSHMET &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& -- CPU REQUIREMENTS in/out 33.03 % 1. calls, 0.58 sec/call analysis 33.10 % 1. calls, 0.58 sec/call metric 33.87 % 1. calls, 0.60 sec/call total 100.00 % 1. calls, 1.77 sec/call ELAPSED TIME 1.77 SEC. (1.77) &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& END OF MODULE MSHMET &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& refine tetgen: Vertex Triangle Border 370 1446 498 Before reconstruction: nbhole=0 nbregion=0 3D RemplissageSurf3D:: Vertex triangle2 border 370 1446 498 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 :: 944 4140 1004 Finish Mesh3 tetgen :: Vertex, Element, Border944 4140 1004 FreeFEM: End check mesh given by tetgen finish reconstruction FreeFEM: End check mesh given by tetgen -- FESpace: Nb of Nodes 944 Nb of DoF 944 GC: converge after 74 g=1.82716e-16 rho= 1.44226 gamma= 0.543498 -- Solve : min -0.100542 max 0.0921181 Plot:: Sorry no ps version for this type of plot 6 -- MSHMET, Release 3.0a (Feb. 15, 2010) Copyright (c) LJLL, 2007-10 mer. 06 mars 2024 16:53:00 CET NUMBER OF GIVEN VERTICES 944 NUMBER OF GIVEN TETRAHEDRA 4140 NUMBER OF GIVEN DATA 944 -- DATA READING COMPLETED. 1.12 sec. &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& MODULE MSHMET-LJLL : 3.0a (Feb. 15, 2010) &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& -- PHASE 1 : ANALYSIS -- PHASE 1 COMPLETED. 1.13 sec. -- PHASE 2 : METRIC -- PHASE 2 COMPLETED. 1.13 sec. &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& END OF MODULE MSHMET &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& -- CPU REQUIREMENTS in/out 33.28 % 1. calls, 1.12 sec/call analysis 33.29 % 1. calls, 1.13 sec/call metric 33.43 % 1. calls, 1.13 sec/call total 100.00 % 1. calls, 3.38 sec/call ELAPSED TIME 3.38 SEC. (3.38) &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& END OF MODULE MSHMET &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& refine tetgen: Vertex Triangle Border 944 4140 1004 Before reconstruction: nbhole=0 nbregion=0 3D RemplissageSurf3D:: Vertex triangle2 border 944 4140 1004 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 :: 2320 11248 1922 Finish Mesh3 tetgen :: Vertex, Element, Border2320 11248 1922 FreeFEM: End check mesh given by tetgen finish reconstruction FreeFEM: End check mesh given by tetgen -- FESpace: Nb of Nodes 2320 Nb of DoF 2320 GC: converge after 105 g=1.10947e-16 rho= 1.48412 gamma= 0.70557 -- Solve : min -0.0997535 max 0.0936481 Plot:: Sorry no ps version for this type of plot 6 -- MSHMET, Release 3.0a (Feb. 15, 2010) Copyright (c) LJLL, 2007-10 mer. 06 mars 2024 16:53:00 CET NUMBER OF GIVEN VERTICES 2320 NUMBER OF GIVEN TETRAHEDRA 11248 NUMBER OF GIVEN DATA 2320 -- DATA READING COMPLETED. 1.38 sec. &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& MODULE MSHMET-LJLL : 3.0a (Feb. 15, 2010) &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& -- PHASE 1 : ANALYSIS -- PHASE 1 COMPLETED. 1.38 sec. -- PHASE 2 : METRIC -- PHASE 2 COMPLETED. 1.39 sec. &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& END OF MODULE MSHMET &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& -- CPU REQUIREMENTS in/out 33.18 % 1. calls, 1.38 sec/call analysis 33.21 % 1. calls, 1.38 sec/call metric 33.62 % 1. calls, 1.39 sec/call total 100.00 % 1. calls, 4.15 sec/call ELAPSED TIME 4.15 SEC. (4.15) &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& END OF MODULE MSHMET &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& refine tetgen: Vertex Triangle Border 2320 11248 1922 Before reconstruction: nbhole=0 nbregion=0 3D RemplissageSurf3D:: Vertex triangle2 border 2320 11248 1922 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 :: 5565 28656 3672 Finish Mesh3 tetgen :: Vertex, Element, Border5565 28656 3672 FreeFEM: End check mesh given by tetgen finish reconstruction FreeFEM: End check mesh given by tetgen -- FESpace: Nb of Nodes 5565 Nb of DoF 5565 GC: converge after 136 g=6.37043e-17 rho= 1.60095 gamma= 0.693046 -- Solve : min -0.0999705 max 0.0938937 Plot:: Sorry no ps version for this type of plot 6 -- MSHMET, Release 3.0a (Feb. 15, 2010) Copyright (c) LJLL, 2007-10 mer. 06 mars 2024 16:53:00 CET NUMBER OF GIVEN VERTICES 5565 NUMBER OF GIVEN TETRAHEDRA 28656 NUMBER OF GIVEN DATA 5565 -- DATA READING COMPLETED. 2.10 sec. &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& MODULE MSHMET-LJLL : 3.0a (Feb. 15, 2010) &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& -- PHASE 1 : ANALYSIS -- PHASE 1 COMPLETED. 2.10 sec. -- PHASE 2 : METRIC -- PHASE 2 COMPLETED. 2.13 sec. &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& END OF MODULE MSHMET &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& -- CPU REQUIREMENTS in/out 33.15 % 1. calls, 2.10 sec/call analysis 33.19 % 1. calls, 2.10 sec/call metric 33.67 % 1. calls, 2.13 sec/call total 100.00 % 1. calls, 6.33 sec/call ELAPSED TIME 6.33 SEC. (6.33) &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& END OF MODULE MSHMET &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& refine tetgen: Vertex Triangle Border 5565 28656 3672 Before reconstruction: nbhole=0 nbregion=0 3D RemplissageSurf3D:: Vertex triangle2 border 5565 28656 3672 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 :: 12922 69494 6794 Finish Mesh3 tetgen :: Vertex, Element, Border12922 69494 6794 FreeFEM: End check mesh given by tetgen finish reconstruction FreeFEM: End check mesh given by tetgen -- FESpace: Nb of Nodes 12922 Nb of DoF 12922 GC: converge after 182 g=3.64607e-17 rho= 1.08649 gamma= 0.728681 -- Solve : min -0.100407 max 0.0939255 Plot:: Sorry no ps version for this type of plot 6 -- Medit, Release 3.0a (Nov. 30, 2007) Copyright (c) LJLL, 1999-2007. compiled: mer. 06 mars 2024 17:04:49 CET (with ff++ 4.14). 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.08 Rendering scene(s) Total running seconds: 0.10 Thank you for using Medit. version de medit ffmedit -popen -filebin -addsol 1 U3 mshmet(): Normal end of execution. current line = 87 exit(0) err code 0 , mpirank 0 CodeAlloc : nb ptr 4100, size :542104 mpirank: 0 Ok: Normal End