unimplemented code fix

Author: Angelyr

src/pmpo_c.cpp

@@ -53,9 +53,10 @@ void polympo_deleteMPMesh_f(MPMesh_ptr p_mpmesh) {
 }
 
 void polympo_setMPICommunicator_f(MPMesh_ptr p_mpmesh, MPI_Fint fcomm){
-    MPI_Comm comm = MPI_Comm_f2c(fcomm);
-    auto p_MPs = ((polyMPO::MPMesh*)p_mpmesh)->p_MPs;
-    p_MPs->setMPIComm(comm);
+  checkMPMeshValid(p_mpmesh);
+  MPI_Comm comm = MPI_Comm_f2c(fcomm);
+  auto p_MPs = ((polyMPO::MPMesh*)p_mpmesh)->p_MPs;
+  p_MPs->setMPIComm(comm);
 }
 
 void polympo_createMPs_f(MPMesh_ptr p_mpmesh,
@@ -453,18 +454,27 @@ void polympo_setMPStrainRate_f(MPMesh_ptr p_mpmesh, const int nComps, const int
   p_MPs->parallel_for(setMPStrainRate, "setMPStrainRate");
 }
 void polympo_getMPStrainRate_f(MPMesh_ptr p_mpmesh, const int nComps, const int numMPs, double* mpStrainRateHost){
+  checkMPMeshValid(p_mpmesh);
+  std::cerr << "Error: This routine is not implemented yet\n";
+  exit(1);
   (void)p_mpmesh;
   (void)nComps;
   (void)numMPs;
   (void)mpStrainRateHost;
 }
 void polympo_setMPStress_f(MPMesh_ptr p_mpmesh, const int nComps, const int numMPs, const double* mpStressIn){
+  checkMPMeshValid(p_mpmesh);
+  std::cerr << "Error: This routine is not implemented yet\n";
+  exit(1);
   (void)p_mpmesh;
   (void)nComps;
   (void)numMPs;
   (void)mpStressIn;
 }
 void polympo_getMPStress_f(MPMesh_ptr p_mpmesh, const int nComps, const int numMPs, double* mpStressHost){
+  checkMPMeshValid(p_mpmesh);
+  std::cerr << "Error: This routine is not implemented yet\n";
+  exit(1);
   (void)p_mpmesh;
   (void)nComps;
   (void)numMPs;
@@ -949,12 +959,16 @@ void polympo_setReconstructionOfVel_f(MPMesh_ptr p_mpmesh, const int order, cons
 
 void polympo_setReconstructionOfStrainRate_f(MPMesh_ptr p_mpmesh, const int order, const int meshEntType){
   checkMPMeshValid(p_mpmesh);
+  std::cerr << "Error: This routine is not implemented yet\n";
+  exit(1);
   (void)order;
   (void)meshEntType;
 }
 
 void polympo_setReconstructionOfStress_f(MPMesh_ptr p_mpmesh, const int order, const int meshEntType){
   checkMPMeshValid(p_mpmesh);
+  std::cerr << "Error: This routine is not implemented yet\n";
+  exit(1);
   (void)order;
   (void)meshEntType;
 }

test/testFortranMPAdvection.f90

@@ -117,17 +117,17 @@ subroutine runReconstructionTest(mpMesh, numMPs, numPush, nCells, nVertices, mp2
       call polympo_setReconstructionOfMass(mpMesh,0,polympo_getMeshFVtxType())
       call polympo_setReconstructionOfVel(mpMesh,0,polympo_getMeshFVtxType())
       call polympo_push(mpMesh)
-      ! call polympo_getMeshElmMass(mpMesh,nCells,c_loc(meshElmMass))
-      ! call polympo_getMeshVtxMass(mpMesh,nVertices,c_loc(meshVtxMass))
-      ! call polympo_getMPCurElmID(mpMesh,numMPs,c_loc(mp2Elm))
+      call polympo_getMeshElmMass(mpMesh,nCells,c_loc(meshElmMass))
+      call polympo_getMeshVtxMass(mpMesh,nVertices,c_loc(meshVtxMass))
+      call polympo_getMPCurElmID(mpMesh,numMPs,c_loc(mp2Elm))
 
-      ! do i = 1, numMPs
-      !   vID = verticesOnCell(1,mp2Elm(i))
-      !   call assert(meshVtxMass(vID) < TEST_VAL+TOLERANCE .and. meshVtxMass(vID) > TEST_VAL-TOLERANCE, "Error: wrong vtx mass")
+      do i = 1, numMPs
+        vID = verticesOnCell(1,mp2Elm(i))
+        call assert(meshVtxMass(vID) < TEST_VAL+TOLERANCE .and. meshVtxMass(vID) > TEST_VAL-TOLERANCE, "Error: wrong vtx mass")
 
-      !   call assert(meshElmMass(mp2Elm(i)) < TEST_VAL+TOLERANCE &
-      !         .and. meshElmMass(mp2Elm(i)) > TEST_VAL-TOLERANCE, "Error: wrong elm mass")
-      ! end do
+        call assert(meshElmMass(mp2Elm(i)) < TEST_VAL+TOLERANCE &
+              .and. meshElmMass(mp2Elm(i)) > TEST_VAL-TOLERANCE, "Error: wrong elm mass")
+      end do
     end do
 
     deallocate(mpMass)