add P1ncS and add too in CHANGELOG.md, and update plugin/seq/Makefile.am

Author: frederichecht

AutoGeneratedFile.tar.gz

@@ -26,7 +26,7 @@ All notable changes to this project will be documented in this file.
 
 ### Added
 
-- Element P1nc3d  (crouzeix Raviart)  in plugin Element_P1nc_3d 
+- Element P1nc3d  (crouzeix Raviart)  in plugin Element_P1nc_3d in 3d and in Surface (meshS)
   and the test is testFE-P1nc3d.edp
 
 ### Changed

CHANGELOG.md

@@ -1,19 +1,13 @@
 load "Element_P1nc_3d"
+{ // 3d 
 mesh3 Th=cube(10,10,10);
 macro grad(u) [dx(u),dy(u),dz(u)]//
-
+func f = x+2*y+3*z;
 fespace Vh(Th,P1nc3d);
 fespace Wh(Th,P0face3d);
 cout << " Vh ndof "<< Vh.ndof <<  " == " << Wh.ndof << endl;
 Vh u,v;
-Vh X,Y,Z;
-/*X = x;
-Y = y;
-Z = z;
-cout << X << endl; 
-cout << Y << endl; 
-cout << Z << endl; */
-func f = x+2*y+3*z;
+
 
 u = f;
 cout << " err L1 : " <<  int3d(Th)( abs(f-u)) << endl;
@@ -23,4 +17,26 @@ plot(u,wait=1);
 int[int] labs = labels(Th);
 solve Pb(u,v)= int3d(Th)(grad(u)'*grad(v))- int3d(Th)(v)+ on(labs,u=0);
 plot(u,wait=1);
-cout << " u max "<< u[] .max << endl;
+cout << " u max "<< u[] .max << endl;
+}
+
+{ // Surf. 
+meshS Th=square3(10,10);
+macro grad(u) [dx(u),dy(u),dz(u)]//
+func f = x+ 2*y;
+fespace Vh(Th,P1ncS);
+fespace Wh(Th,P0edgeS);
+cout << " Vh ndof "<< Vh.ndof <<  " == " << Wh.ndof << endl;
+Vh u,v;
+
+
+u = f;
+cout << "2Surf err L1 : " <<  int2d(Th)( abs(f-u)) << endl;
+plot(u,wait=1);
+
+
+int[int] labs = labels(Th);
+solve Pb(u,v)= int2d(Th)(grad(u)'*grad(v))- int2d(Th)(v)+ on(labs,u=0);
+plot(u,wait=1);
+cout << " u max "<< u[].max << endl;
+}

examples/plugin/testFE-P1nc3d.edp

@@ -52,7 +52,9 @@ LIST_COMPILE=myfunction.$(DYLIB_SUFFIX) BernardiRaugel.$(DYLIB_SUFFIX) \
 	geophysics.$(DYLIB_SUFFIX) CircumCenter.$(DYLIB_SUFFIX) MatD-VFP0.$(DYLIB_SUFFIX) \
 	meshtools.$(DYLIB_SUFFIX) Helmholtz_FD.$(DYLIB_SUFFIX)  Element_P3pnc.$(DYLIB_SUFFIX) \
 	vortextools.$(DYLIB_SUFFIX)  MatrixMarket.$(DYLIB_SUFFIX) Element_P3nc.$(DYLIB_SUFFIX) \
-	plotPDF.$(DYLIB_SUFFIX) Element_P2pnc_3d.$(DYLIB_SUFFIX) Element_P3pnc_3d.$(DYLIB_SUFFIX)
+	plotPDF.$(DYLIB_SUFFIX) Element_P2pnc_3d.$(DYLIB_SUFFIX) Element_P3pnc_3d.$(DYLIB_SUFFIX) \
+	Element_P1nc_3d.$(DYLIB_SUFFIX)
+
 
 # FFCS - some libraries are skipped because the corresponding tool is deactivated.