Hello everyone,
I’m having a problem with make check in the installation of FreeFEM on Ubuntu 23.10. When I do the command make -j4 check, I get the following error:
Making check in 3dSurf
make[2]: Entering directory '/home/sgau/lib/FREEFEM/examples/3dSurf'
make check-TESTS
make[3]: Entering directory '/home/sgau/lib/FREEFEM/examples/3dSurf'
make[4]: Entering directory '/home/sgau/lib/FREEFEM/examples/3dSurf'
PASS: buildMeshSfromMesh3.edp
PASS: glu3DSurf.edp
PASS: Laplace3dSP2.edp
PASS: LaplacianSurf.edp
PASS: operatorsOnMeshS.edp
PASS: extractSurfFromMesh3.edp
../../bin/test-driver-ff: line 126: 54197 Segmentation fault (core dumped) ${TEST_FFPP} ${FLAGS_FFPP_B} "$@" ${FLAGS_FFPP_A} > $log_file 2>&1
FAIL: testvtk.edp
PASS: LapP1withflux.edp
PASS: cauchyproblem.edp
PASS: periodicS.edp
PASS: Connectivite-S.edp
PASS: LaplaceRT1Surf.edp
PASS: buildmeshS.edp
PASS: HeatTorus.edp
PASS: Pinocchio.edp
============================================================================
Testsuite summary for FreeFEM 4.14
============================================================================
# TOTAL: 15
# PASS: 14
# SKIP: 0
# XFAIL: 0
# FAIL: 1
# XPASS: 0
# ERROR: 0
============================================================================
See examples/3dSurf/test-suite.log
Please report to frederic.hecht@sorbonne-universite.fr
============================================================================
In the testsuite.log file there is written the following:
==================================================
FreeFEM 4.14: examples/3dSurf/test-suite.log
==================================================
# TOTAL: 15
# PASS: 14
# SKIP: 0
# XFAIL: 0
# FAIL: 1
# XPASS: 0
# ERROR: 0
.. contents:: :depth: 2
PASS: buildMeshSfromMesh3.edp
=============================
-- FreeFem++ v4.14 (Tue Dec 26 09:55:30 AM CET 2023 - git no git)
file : ./buildMeshSfromMesh3.edp
Load: lg_fem lg_mesh lg_mesh3 eigenvalue
1 : load "msh3"
2 : //load "medit"
3 :
4 : //load "gmsh"
5 : //load "iovtk"
6 :
7 : int n= 10;
8 : int nvb = (n+1)^3 - (n-1)^3;// Nb boundary vertices
9 : int ntb = n*n*12; // Nb of Boundary triangle
10 : mesh3 Th=cube(n,n,n);
11 : Th = buildBdMesh(Th); // build the surface mesh
12 :
13 :
14 : int[int] ref2=[1,0,2,0];
15 : mesh3 Thm=change(Th,label=ref2);
16 :
17 :
18 : meshS Th1 = Th.Gamma;
19 :
20 : assert(Th1.nv==602);
21 : assert(Th1.nt==1200);
22 : assert(Th1.nbe==120);
23 :
24 : plot(Th1);
25 : savemesh(Th1, "Th1.mesh");
26 :
27 : meshS Th2 = readmeshS("Th1.mesh");
28 :
29 : assert(Th1.nv==Th2.nv);
30 : assert(Th1.nt==Th2.nt);
31 : assert(Th1.nbe==Th2.nbe);
32 : sizestack + 1024 =1896 ( 872 )
-- Cube nv=1331 nt=6000 nbe=1200 kind= 6
Build meshS from mesh3....
Plot:: Sorry no ps version for this type of plot 50
read meshS ok 0surface Mesh, num Triangles:= 1200, num Vertice:= 602 num boundary Edges:= 120
-- MeshS : Th1.mesh, space dimension 3, num Triangle elts 1200, num Vertice 602 num Bordary elts 120
times: compile 0.117124s, execution 0.210865s, mpirank:0
CodeAlloc : nb ptr 3945, size :527640 mpirank: 0
Ok: Normal End
PASS: Laplace3dSP2.edp
======================
-- FreeFem++ v4.14 (Tue Dec 26 09:55:30 AM CET 2023 - git no git)
file : ./Laplace3dSP2.edp
Load: lg_fem lg_mesh lg_mesh3 eigenvalue
1 : load "gmsh"
2 :
3 : load "msh3"
4 : int nn=20;
5 : mesh Th2=square(nn,nn,region=0);
6 : fespace Vh2(Th2,P2);
7 : Vh2 ux,uz,p2;
8 : int[int] rup=[0,2], rdown=[0,1], rmid=[1,1,2,1,3,1,4,1];
9 : real zmin=0,zmax=1;
10 : mesh3 Th=buildlayers(Th2,nn,
11 : zbound=[zmin,zmax],
12 : labelmid=rmid,
13 : reffaceup = rup,
14 : reffacelow = rdown);
15 :
16 : Th = buildBdMesh(Th); // build the surface mesh
17 : meshS ThS= Th.Gamma;
18 :
19 : fespace Vh(ThS,P2);
20 : Vh u,v;
21 :
22 : func ue = 2*x*x + 3*y*y + 4*z*z + 5*x*y+6*x*z+1;
23 : func uex= 4*x+ 5*y+6*z;
24 : func uey= 6*y + 5*x;
25 : func uez= 8*z +6*x;
26 : func f= -18. ;
27 :
28 : macro Grad3(u) [dx(u),dy(u),dz(u)] ) // EOM
29 : //
30 : problem Lap3dS(u,v,solver=CG) =
31 : int2d(ThS)(Grad3(v) [dx(v),dy(v),dz(v)] ' *Grad3(u) [dx(u),dy(u),dz(u)] ) //') for emacs
32 : + int1d(ThS,2)(u*v)
33 : - int2d(ThS)(f*v)
34 : - int1d(ThS,2) ( ue*v + (uex*N.x +uey*N.y +uez*N.z)*v )
35 : + on(1,u=ue);
36 : Lap3dS;
37 :
38 : plot(u,wait=1);
39 :
40 :
41 : savegmsh(ThS,"ThS");
42 :
43 :
44 :
45 : sizestack + 1024 =7088 ( 6064 )
-- Square mesh : nb vertices =441 , nb triangles = 800 , nb boundary edges 80
Build meshS from mesh3....
-- Build Nodes/DF on mesh : n.v. 2402, n. elmt. 4800, n b. elmt. 240
nb of Nodes 9602 nb of DoF 9602 DFon=1100
-- FESpace: Nb of Nodes 9602 Nb of DoF 9602
GC: converge after 118 g=2.9998e-08 rho= 1.56891 gamma= 0.657623
-- Solve :
min 0.989202 max 21
Plot:: Sorry no ps version for this type of plot 8
times: compile 0.119846s, execution 1.63546s, mpirank:0
CodeAlloc : nb ptr 4158, size :540088 mpirank: 0
Ok: Normal End
FAIL: testvtk.edp
=================
-- FreeFem++ v4.14 (Tue Dec 26 09:55:30 AM CET 2023 - git no git)
file : ./testvtk.edp
Load: lg_fem lg_mesh lg_mesh3 eigenvalue
1 : // basic example io vtk for meshS/meshL
2 : load "iovtk" load: iovtk
3 : load "msh3"
4 :
5 : mesh Th = square(10,10);
6 : savevtk("Th.vtk",Th,bin=1);
7 : mesh Th2=vtkload("Th.vtk")Load mesh given by VTK
;
8 : meshS ThS = square3(10,10);
9 :
10 : savevtk("ThS.vtk",ThS,bin=0);
11 : savevtk("ThS.vtu",ThS,bin=1);
12 : meshS ThS2=vtkloadS("ThS.vtk")Load mesh given by VTK
I do not know how to deal with this and I’m looking forward to your help.
Thank you so much for reading and for your help.
Best wishes,
Gabriel