Add NodeByCellData::operator=

performs a shallow copy

src/main.cpp

@@ -197,7 +197,6 @@ int main(int argc, char *argv[])
           typedef TinyVector<MeshType::dimension> Rd;
 
           const Kokkos::View<const double*> Vj = mesh_data.Vj();
-          const NodeByCellData<Rd> Cjr = mesh_data.Cjr();
 
           const double tmax=0.2;
           double t=0;

src/mesh/MeshData.hpp

@@ -214,12 +214,8 @@ public:
           m_Cjr(j,0)=-1;
           m_Cjr(j,1)= 1;
         });
-      // in 1d njr=Cjr
-#warning set equality between njr and Cjr
-      Kokkos::parallel_for(m_mesh.numberOfCells(), KOKKOS_LAMBDA(const int& j) {
-          m_njr(j,0)=-1;
-          m_njr(j,1)= 1;
-        });
+      // in 1d njr=Cjr (here performs a shallow copy)
+      m_njr=m_Cjr;
       Kokkos::parallel_for(m_mesh.numberOfCells(), KOKKOS_LAMBDA(const int& j) {
           m_ljr(j,0)= 1;
           m_ljr(j,1)= 1;

src/scheme/NodeByCellData.hpp

@@ -96,6 +96,8 @@ class NodeByCellData
     }
   };
 
+  NodeByCellData& operator=(const NodeByCellData&) = default;
+
   KOKKOS_INLINE_FUNCTION
   DataType& operator()(const size_t& j, const size_t& r)
   {