Test memcpy by changing SD field to contiguous RightLayout

Author: dhyan1272

src/pmpo_MPMesh.hpp

@@ -56,12 +56,14 @@ class MPMesh{
 
       timer.reset();
       std::vector<double> fieldData1(nEntities*numEntries);
-      pumipic::RecordTime("SD: STL alloc" + std::to_string(self), timer.seconds());
+      pumipic::RecordTime("SD: STL alloc-" + std::to_string(self), timer.seconds());
       timer.reset();
+      /*
       for (int i=0;i<nEntities;i++)
         for (int j=0;j<numEntries;j++)
           fieldData1[i*numEntries+j]=reconVals_host(i,j);
-      //std::memcpy(fieldData1.data(), reconVals_host.data(), nEntities*numEntries*sizeof(double));
+      */
+      std::memcpy(fieldData1.data(), reconVals_host.data(), nEntities*numEntries*sizeof(double));
       pumipic::RecordTime("SD: STL alloc_copy-" + std::to_string(self), timer.seconds());
 
 /*

src/pmpo_mesh.cpp

@@ -54,7 +54,7 @@ namespace polyMPO{
 
     auto stressDivergenceEntry = meshFields2TypeAndString.at(MeshF_StressDivergence);
     PMT_ALWAYS_ASSERT(stressDivergenceEntry.first == MeshFType_VtxBased);
-    stressDivergence_ = MeshFView<MeshF_StressDivergence>(stressDivergenceEntry.second, numVtxs_);
+    stressDivergence_ = MeshFView<MeshF_StressDivergence>(stressDivergenceEntry.second, numVtxs_, 2);
 
     auto solveVelocityEntry = meshFields2TypeAndString.at(MeshF_SolveVelocity);
     PMT_ALWAYS_ASSERT(solveVelocityEntry.first == MeshFType_VtxBased);

src/pmpo_mesh.hpp

@@ -66,7 +66,7 @@ template <> struct meshFieldToType < MeshF_TanLatVertexRotatedOverRadius > { usi
 template <> struct meshFieldToType < MeshF_SolveStress       > { using type = IntView; };
 template <> struct meshFieldToType < MeshF_SolveVelocity     > { using type = IntView; };
 template <> struct meshFieldToType < MeshF_InteriorVertex    > { using type = IntView; };
-template <> struct meshFieldToType < MeshF_StressDivergence  > { using type = Kokkos::View<vec2d_t*>; };
+template <> struct meshFieldToType < MeshF_StressDivergence  > { using type = Kokkos::View<double**, Kokkos::LayoutRight>; };
 template <> struct meshFieldToType < MeshF_TotalMassVtx      > { using type = Kokkos::View<doubleSclr_t*>; };
 template <> struct meshFieldToType < MeshF_AirStress         > { using type = Kokkos::View<vec2d_t*>; };
 template <> struct meshFieldToType < MeshF_SurfaceTilt       > { using type = Kokkos::View<vec2d_t*>; };