assert spherical

Author: Angelyr

src/pmpo_MPMesh.cpp

@@ -8,6 +8,7 @@ namespace polyMPO{
 void printVTP_mesh(MPMesh& mpMesh, int printVTPIndex=-1);
 
 void MPMesh::calcBasis() {
+    assert(p_mesh->getGeomType() == geom_spherical_surf);
     auto MPsPosition = p_MPs->getPositions();
     auto mp_basis_field = p_MPs->getData<MPF_Basis_Vals>(); // we can implement MPs->getBasisVals() like MPs->getPositions()
     auto elm2VtxConn = p_mesh->getElm2VtxConn();

src/pmpo_createTestMPMesh.cpp

@@ -149,7 +149,7 @@ MaterialPoints* initTestMPs(Mesh* mesh, int testMPOption){
         default:
             fprintf(stderr,"TestMPOption not avaiable! return an empty one!");
     }
-    auto vtxCoords = mesh->getMeshField<polyMPO::MeshF_VtxCoords>();   
+    auto vtxCoords = mesh->getMeshField<MeshF_VtxCoords>();   
     auto elm2VtxConn = mesh->getElm2VtxConn();
     auto geomType = mesh->getGeomType();
 
@@ -245,8 +245,8 @@ Mesh* replicateMesh(Mesh* mesh, int replicateFactor){
     const int nCells = nCells_size*replicateFactor;
     const int nVertices = nVertices_size*replicateFactor;
 
-    auto v_array = mesh->getMeshField<polyMPO::MeshF_VtxCoords>(); 
-    MeshFView<polyMPO::MeshF_VtxCoords> vtxCoords("verticesCoordinates", nVertices);
+    auto v_array = mesh->getMeshField<MeshF_VtxCoords>(); 
+    MeshFView<MeshF_VtxCoords> vtxCoords("verticesCoordinates", nVertices);
     Kokkos::parallel_for("set vtxCoords", nVertices_size, KOKKOS_LAMBDA(const int i){
         for(int f=0; f<replicateFactor; f++){
             vtxCoords(i+f*nVertices_size,0) = v_array(i,0);