diff --git a/packages/stokhos/src/sacado/kokkos/vector/belos/Belos_TpetraAdapter_MP_Vector.hpp b/packages/stokhos/src/sacado/kokkos/vector/belos/Belos_TpetraAdapter_MP_Vector.hpp
index 9987be272f92..549971d79fb2 100644
--- a/packages/stokhos/src/sacado/kokkos/vector/belos/Belos_TpetraAdapter_MP_Vector.hpp
+++ b/packages/stokhos/src/sacado/kokkos/vector/belos/Belos_TpetraAdapter_MP_Vector.hpp
@@ -324,6 +324,38 @@ namespace Belos {
       return mv.isConstantStride ();
     }
 
+    template <class DstType, class SrcType>
+    static void deep_copy_2d_view_with_intercessory_space(
+      const DstType& dst,
+      const SrcType& src
+    )
+    {
+      if (std::is_same<typename SrcType::memory_space, Kokkos::HostSpace>::value)
+      {
+        auto dst_i = Kokkos::create_mirror(dst);
+        for (size_t i = 0; i < src.extent(0); i++)
+        {
+          for (size_t j = 0; j < src.extent(1); j++)
+          {
+            dst_i(i, j) = src(i, j);
+          }
+        }
+        Kokkos::deep_copy(dst, dst_i);
+      }
+      else
+      {
+        auto src_i = Kokkos::create_mirror(src);
+        Kokkos::deep_copy(src_i, src);
+        for (size_t i = 0; i < src.extent(0); i++)
+        {
+          for (size_t j = 0; j < src.extent(1); j++)
+          {
+            dst(i, j) = src_i(i, j);
+          }
+        }
+      }
+    }
+
     static void
     MvTimesMatAddMv (const dot_type& alpha,
                      const Tpetra::MultiVector<Scalar,LO,GO,Node>& A,
@@ -372,7 +404,15 @@ namespace Belos {
       typedef Kokkos::View<dot_type*, Kokkos::LayoutLeft, execution_space> b_1d_view_type;
       b_1d_view_type B_1d_view_dev(Kokkos::ViewAllocateWithoutInitializing("B"), numRowsB*numColsB);
       b_view_type B_view_dev( B_1d_view_dev.data(), numRowsB, numColsB);
-      Kokkos::deep_copy(B_view_dev, B_view_host);
+
+      if (Kokkos::SpaceAccessibility<Kokkos::HostSpace, typename execution_space::memory_space>::accessible)
+      {
+        Kokkos::deep_copy(B_view_dev, B_view_host);
+      }
+      else
+      {
+        deep_copy_2d_view_with_intercessory_space(B_view_dev, B_view_host);
+      }
 
       // Do local multiply
       {
@@ -499,16 +539,34 @@ namespace Belos {
       // reduce across processors -- could check for RDMA
       RCP<const Comm<int> > pcomm = A.getMap()->getComm ();
       if (pcomm->getSize () == 1)
-        Kokkos::deep_copy(C_view_host, C_view_dev);
-      else {
+      {
+        if (Kokkos::SpaceAccessibility<execution_space, Kokkos::HostSpace>::accessible)
+        {
+          Kokkos::deep_copy(C_view_host, C_view_dev);
+        }
+        else
+        {
+          deep_copy_2d_view_with_intercessory_space(C_view_host, C_view_dev);
+        }
+      }
+      else
+      {
         typedef Kokkos::View<dot_type*, Kokkos::LayoutLeft, Kokkos::HostSpace> c_1d_host_view_type;
         c_1d_host_view_type C_1d_view_tmp(Kokkos::ViewAllocateWithoutInitializing("C_tmp"), strideC*numColsC);
-        c_host_view_type C_view_tmp( C_1d_view_tmp.data(),
-                                     strideC, numColsC);
-        Kokkos::deep_copy(Kokkos::subview(C_view_tmp, 
-                                          Kokkos::pair<int,int>(0,numRowsC), 
-                                          Kokkos::pair<int,int>(0,numColsC)) , 
-                          C_view_dev);
+        c_host_view_type C_view_tmp_input( C_1d_view_tmp.data(), strideC, numColsC);
+        auto C_view_tmp = Kokkos::subview(C_view_tmp_input,
+                                          Kokkos::pair<int,int>(0,numRowsC),
+                                          Kokkos::pair<int,int>(0,numColsC));
+        
+        if (Kokkos::SpaceAccessibility<execution_space, Kokkos::HostSpace>::accessible)
+        {
+          Kokkos::deep_copy(C_view_tmp, C_view_dev);
+        }
+        else
+        {
+          deep_copy_2d_view_with_intercessory_space(C_view_tmp, C_view_dev);
+        }
+
         reduceAll<int> (*pcomm, REDUCE_SUM, strideC*numColsC,
                         C_view_tmp.data(),
                         C_view_host.data());
diff --git a/packages/stokhos/test/UnitTest/Stokhos_TpetraCrsMatrixMPVectorUnitTest.hpp b/packages/stokhos/test/UnitTest/Stokhos_TpetraCrsMatrixMPVectorUnitTest.hpp
index deb7e1e7c528..fd920498f37e 100644
--- a/packages/stokhos/test/UnitTest/Stokhos_TpetraCrsMatrixMPVectorUnitTest.hpp
+++ b/packages/stokhos/test/UnitTest/Stokhos_TpetraCrsMatrixMPVectorUnitTest.hpp
@@ -156,20 +156,20 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
   // Fill vectors
   RCP<Tpetra_Vector> x1 = Tpetra::createVector<Scalar>(map);
   RCP<Tpetra_Vector> x2 = Tpetra::createVector<Scalar>(map);
-  ArrayRCP<Scalar> x1_view = x1->get1dViewNonConst();
-  ArrayRCP<Scalar> x2_view = x2->get1dViewNonConst();
-  Scalar val1(VectorSize, BaseScalar(0.0)), val2(VectorSize, BaseScalar(0.0));
-  for (size_t i=0; i<num_my_row; ++i) {
-    const GlobalOrdinal row = myGIDs[i];
-    for (LocalOrdinal j=0; j<VectorSize; ++j) {
-      val1.fastAccessCoeff(j) = generate_vector_coefficient<BaseScalar,size_t>(nrow, VectorSize, row, j);
-      val2.fastAccessCoeff(j) = 0.12345 * generate_vector_coefficient<BaseScalar,size_t>(nrow, VectorSize, row, j);
+  {
+    auto x1_view = x1->getLocalViewHost(Tpetra::Access::OverwriteAll);
+    auto x2_view = x2->getLocalViewHost(Tpetra::Access::OverwriteAll);
+    Scalar val1(VectorSize, BaseScalar(0.0)), val2(VectorSize, BaseScalar(0.0));
+    for (size_t i=0; i<num_my_row; ++i) {
+      const GlobalOrdinal row = myGIDs[i];
+      for (LocalOrdinal j=0; j<VectorSize; ++j) {
+        val1.fastAccessCoeff(j) = generate_vector_coefficient<BaseScalar,size_t>(nrow, VectorSize, row, j);
+        val2.fastAccessCoeff(j) = 0.12345 * generate_vector_coefficient<BaseScalar,size_t>(nrow, VectorSize, row, j);
+      }
+      x1_view(i,0) = val1;
+      x2_view(i,0) = val2;
     }
-    x1_view[i] = val1;
-    x2_view[i] = val2;
   }
-  x1_view = Teuchos::null;
-  x2_view = Teuchos::null;
 
   // Add
   Scalar alpha = 2.1;
@@ -181,7 +181,7 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
   //             Teuchos::VERB_EXTREME);
 
   // Check
-  ArrayRCP<Scalar> y_view = y->get1dViewNonConst();
+  auto y_view = y->getLocalViewHost(Tpetra::Access::ReadOnly);
   Scalar val(VectorSize, BaseScalar(0.0));
   BaseScalar tol = 1.0e-14;
   for (size_t i=0; i<num_my_row; ++i) {
@@ -191,9 +191,9 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
         nrow, VectorSize, row, j);
       val.fastAccessCoeff(j) = alpha.coeff(j)*v + 0.12345*beta.coeff(j)*v;
     }
-    TEST_EQUALITY( y_view[i].size(), VectorSize );
+    TEST_EQUALITY( y_view(i,0).size(), VectorSize );
     for (LocalOrdinal j=0; j<VectorSize; ++j)
-      TEST_FLOATING_EQUALITY( y_view[i].fastAccessCoeff(j), val.fastAccessCoeff(j), tol );
+      TEST_FLOATING_EQUALITY( y_view(i,0).fastAccessCoeff(j), val.fastAccessCoeff(j), tol );
   }
 }
 
@@ -235,20 +235,20 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
   // Fill vectors
   RCP<Tpetra_Vector> x1 = Tpetra::createVector<Scalar>(map);
   RCP<Tpetra_Vector> x2 = Tpetra::createVector<Scalar>(map);
-  ArrayRCP<Scalar> x1_view = x1->get1dViewNonConst();
-  ArrayRCP<Scalar> x2_view = x2->get1dViewNonConst();
-  Scalar val1(VectorSize, BaseScalar(0.0)), val2(VectorSize, BaseScalar(0.0));
-  for (size_t i=0; i<num_my_row; ++i) {
-    const GlobalOrdinal row = myGIDs[i];
-    for (LocalOrdinal j=0; j<VectorSize; ++j) {
-      val1.fastAccessCoeff(j) = generate_vector_coefficient<BaseScalar,size_t>(nrow, VectorSize, row, j);
-      val2.fastAccessCoeff(j) = 0.12345 * generate_vector_coefficient<BaseScalar,size_t>(nrow, VectorSize, row, j);
+  {
+    auto x1_view = x1->getLocalViewHost(Tpetra::Access::OverwriteAll);
+    auto x2_view = x2->getLocalViewHost(Tpetra::Access::OverwriteAll);
+    Scalar val1(VectorSize, BaseScalar(0.0)), val2(VectorSize, BaseScalar(0.0));
+    for (size_t i=0; i<num_my_row; ++i) {
+      const GlobalOrdinal row = myGIDs[i];
+      for (LocalOrdinal j=0; j<VectorSize; ++j) {
+        val1.fastAccessCoeff(j) = generate_vector_coefficient<BaseScalar,size_t>(nrow, VectorSize, row, j);
+        val2.fastAccessCoeff(j) = 0.12345 * generate_vector_coefficient<BaseScalar,size_t>(nrow, VectorSize, row, j);
+      }
+      x1_view(i,0) = val1;
+      x2_view(i,0) = val2;
     }
-    x1_view[i] = val1;
-    x2_view[i] = val2;
   }
-  x1_view = Teuchos::null;
-  x2_view = Teuchos::null;
 
   // Dot product
   dot_type dot = x1->dot(*x2);
@@ -337,26 +337,26 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
   size_t ncol = 5;
   RCP<Tpetra_MultiVector> x1 = Tpetra::createMultiVector<Scalar>(map, ncol);
   RCP<Tpetra_MultiVector> x2 = Tpetra::createMultiVector<Scalar>(map, ncol);
-  ArrayRCP< ArrayRCP<Scalar> > x1_view = x1->get2dViewNonConst();
-  ArrayRCP< ArrayRCP<Scalar> > x2_view = x2->get2dViewNonConst();
-  Scalar val1(VectorSize, BaseScalar(0.0)), val2(VectorSize, BaseScalar(0.0));
-  for (size_t i=0; i<num_my_row; ++i) {
-    const GlobalOrdinal row = myGIDs[i];
-    for (size_t j=0; j<ncol; ++j) {
-      for (LocalOrdinal k=0; k<VectorSize; ++k) {
-        BaseScalar v =
-          generate_multi_vector_coefficient<BaseScalar,size_t>(
-            nrow, ncol, VectorSize, row, j, k);
-        val1.fastAccessCoeff(k) = v;
-        val2.fastAccessCoeff(k) = 0.12345 * v;
+  {
+    auto x1_view = x1->getLocalViewHost(Tpetra::Access::OverwriteAll);
+    auto x2_view = x2->getLocalViewHost(Tpetra::Access::OverwriteAll);
+    Scalar val1(VectorSize, BaseScalar(0.0)), val2(VectorSize, BaseScalar(0.0));
+    for (size_t i=0; i<num_my_row; ++i) {
+      const GlobalOrdinal row = myGIDs[i];
+      for (size_t j=0; j<ncol; ++j) {
+        for (LocalOrdinal k=0; k<VectorSize; ++k) {
+          BaseScalar v =
+            generate_multi_vector_coefficient<BaseScalar,size_t>(
+              nrow, ncol, VectorSize, row, j, k);
+          val1.fastAccessCoeff(k) = v;
+          val2.fastAccessCoeff(k) = 0.12345 * v;
+        }
+        x1_view(i,j) = val1;
+        x2_view(i,j) = val2;
       }
-      x1_view[j][i] = val1;
-      x2_view[j][i] = val2;
     }
   }
-  x1_view = Teuchos::null;
-  x2_view = Teuchos::null;
-
+    
   // Add
   Scalar alpha = 2.1;
   Scalar beta = 3.7;
@@ -367,7 +367,7 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
   //             Teuchos::VERB_EXTREME);
 
   // Check
-  ArrayRCP< ArrayRCP<Scalar> > y_view = y->get2dViewNonConst();
+  auto y_view = y->getLocalViewHost(Tpetra::Access::ReadOnly);
   Scalar val(VectorSize, BaseScalar(0.0));
   BaseScalar tol = 1.0e-14;
   for (size_t i=0; i<num_my_row; ++i) {
@@ -378,9 +378,9 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
           nrow, ncol, VectorSize, row, j, k);
         val.fastAccessCoeff(k) = alpha.coeff(k)*v + 0.12345*beta.coeff(k)*v;
       }
-      TEST_EQUALITY( y_view[j][i].size(), VectorSize );
+      TEST_EQUALITY( y_view(i,j).size(), VectorSize );
       for (LocalOrdinal k=0; k<VectorSize; ++k)
-        TEST_FLOATING_EQUALITY( y_view[j][i].fastAccessCoeff(k),
+        TEST_FLOATING_EQUALITY( y_view(i,j).fastAccessCoeff(k),
                                 val.fastAccessCoeff(k), tol );
     }
   }
@@ -425,25 +425,25 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
   size_t ncol = 5;
   RCP<Tpetra_MultiVector> x1 = Tpetra::createMultiVector<Scalar>(map, ncol);
   RCP<Tpetra_MultiVector> x2 = Tpetra::createMultiVector<Scalar>(map, ncol);
-  ArrayRCP< ArrayRCP<Scalar> > x1_view = x1->get2dViewNonConst();
-  ArrayRCP< ArrayRCP<Scalar> > x2_view = x2->get2dViewNonConst();
-  Scalar val1(VectorSize, BaseScalar(0.0)), val2(VectorSize, BaseScalar(0.0));
-  for (size_t i=0; i<num_my_row; ++i) {
-    const GlobalOrdinal row = myGIDs[i];
-    for (size_t j=0; j<ncol; ++j) {
-      for (LocalOrdinal k=0; k<VectorSize; ++k) {
-        BaseScalar v =
-          generate_multi_vector_coefficient<BaseScalar,size_t>(
-            nrow, ncol, VectorSize, row, j, k);
-        val1.fastAccessCoeff(k) = v;
-        val2.fastAccessCoeff(k) = 0.12345 * v;
+  {
+    auto x1_view = x1->getLocalViewHost(Tpetra::Access::OverwriteAll);
+    auto x2_view = x2->getLocalViewHost(Tpetra::Access::OverwriteAll);
+    Scalar val1(VectorSize, BaseScalar(0.0)), val2(VectorSize, BaseScalar(0.0));
+    for (size_t i=0; i<num_my_row; ++i) {
+      const GlobalOrdinal row = myGIDs[i];
+      for (size_t j=0; j<ncol; ++j) {
+        for (LocalOrdinal k=0; k<VectorSize; ++k) {
+          BaseScalar v =
+            generate_multi_vector_coefficient<BaseScalar,size_t>(
+              nrow, ncol, VectorSize, row, j, k);
+          val1.fastAccessCoeff(k) = v;
+          val2.fastAccessCoeff(k) = 0.12345 * v;
+        }
+        x1_view(i,j) = val1;
+        x2_view(i,j) = val2;
       }
-      x1_view[j][i] = val1;
-      x2_view[j][i] = val2;
     }
   }
-  x1_view = Teuchos::null;
-  x2_view = Teuchos::null;
 
   // Dot product
   Array<dot_type> dots(ncol);
@@ -540,25 +540,25 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
   size_t ncol = 5;
   RCP<Tpetra_MultiVector> x1 = Tpetra::createMultiVector<Scalar>(map, ncol);
   RCP<Tpetra_MultiVector> x2 = Tpetra::createMultiVector<Scalar>(map, ncol);
-  ArrayRCP< ArrayRCP<Scalar> > x1_view = x1->get2dViewNonConst();
-  ArrayRCP< ArrayRCP<Scalar> > x2_view = x2->get2dViewNonConst();
-  Scalar val1(VectorSize, BaseScalar(0.0)), val2(VectorSize, BaseScalar(0.0));
-  for (size_t i=0; i<num_my_row; ++i) {
-    const GlobalOrdinal row = myGIDs[i];
-    for (size_t j=0; j<ncol; ++j) {
-      for (LocalOrdinal k=0; k<VectorSize; ++k) {
-        BaseScalar v =
-          generate_multi_vector_coefficient<BaseScalar,size_t>(
-            nrow, ncol, VectorSize, row, j, k);
-        val1.fastAccessCoeff(k) = v;
-        val2.fastAccessCoeff(k) = 0.12345 * v;
+  {
+    auto x1_view = x1->getLocalViewHost(Tpetra::Access::OverwriteAll);
+    auto x2_view = x2->getLocalViewHost(Tpetra::Access::OverwriteAll);
+    Scalar val1(VectorSize, BaseScalar(0.0)), val2(VectorSize, BaseScalar(0.0));
+    for (size_t i=0; i<num_my_row; ++i) {
+      const GlobalOrdinal row = myGIDs[i];
+      for (size_t j=0; j<ncol; ++j) {
+        for (LocalOrdinal k=0; k<VectorSize; ++k) {
+          BaseScalar v =
+            generate_multi_vector_coefficient<BaseScalar,size_t>(
+              nrow, ncol, VectorSize, row, j, k);
+          val1.fastAccessCoeff(k) = v;
+          val2.fastAccessCoeff(k) = 0.12345 * v;
+        }
+        x1_view(i,j) = val1;
+        x2_view(i,j) = val2;
       }
-      x1_view[j][i] = val1;
-      x2_view[j][i] = val2;
     }
   }
-  x1_view = Teuchos::null;
-  x2_view = Teuchos::null;
 
   // Get subviews
   size_t ncol_sub = 2;
@@ -700,15 +700,16 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
 
   // Fill vector
   RCP<Tpetra_Vector> x = Tpetra::createVector<Scalar>(map);
-  ArrayRCP<Scalar> x_view = x->get1dViewNonConst();
-  for (size_t i=0; i<num_my_row; ++i) {
-    const GlobalOrdinal row = myGIDs[i];
-    for (LocalOrdinal j=0; j<VectorSize; ++j)
-      val.fastAccessCoeff(j) = generate_vector_coefficient<BaseScalar,size_t>(
-        nrow, VectorSize, row, j);
-    x_view[i] = val;
+  {
+    auto x_view = x->getLocalViewHost(Tpetra::Access::OverwriteAll);
+    for (size_t i=0; i<num_my_row; ++i) {
+      const GlobalOrdinal row = myGIDs[i];
+      for (LocalOrdinal j=0; j<VectorSize; ++j)
+        val.fastAccessCoeff(j) = generate_vector_coefficient<BaseScalar,size_t>(
+          nrow, VectorSize, row, j);
+      x_view(i,0) = val;
+    }
   }
-  x_view = Teuchos::null;
 
   // matrix->describe(*(Teuchos::fancyOStream(rcp(&std::cout,false))),
   //                  Teuchos::VERB_EXTREME);
@@ -724,7 +725,7 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
   //             Teuchos::VERB_EXTREME);
 
   // Check
-  ArrayRCP<Scalar> y_view = y->get1dViewNonConst();
+  auto y_view = y->getLocalViewHost(Tpetra::Access::ReadOnly);
   BaseScalar tol = 1.0e-14;
   for (size_t i=0; i<num_my_row; ++i) {
     const GlobalOrdinal row = myGIDs[i];
@@ -740,9 +741,9 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
         val.fastAccessCoeff(j) += v*v;
       }
     }
-    TEST_EQUALITY( y_view[i].size(), VectorSize );
+    TEST_EQUALITY( y_view(i,0).size(), VectorSize );
     for (LocalOrdinal j=0; j<VectorSize; ++j)
-      TEST_FLOATING_EQUALITY( y_view[i].fastAccessCoeff(j), val.fastAccessCoeff(j), tol );
+      TEST_FLOATING_EQUALITY( y_view(i,0).fastAccessCoeff(j), val.fastAccessCoeff(j), tol );
   }
 }
 
@@ -822,20 +823,21 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
   // Fill multi-vector
   size_t ncol = 5;
   RCP<Tpetra_MultiVector> x = Tpetra::createMultiVector<Scalar>(map, ncol);
-  ArrayRCP< ArrayRCP<Scalar> > x_view = x->get2dViewNonConst();
-  for (size_t i=0; i<num_my_row; ++i) {
-    const GlobalOrdinal row = myGIDs[i];
-    for (size_t j=0; j<ncol; ++j) {
-      for (LocalOrdinal k=0; k<VectorSize; ++k) {
-        BaseScalar v =
-          generate_multi_vector_coefficient<BaseScalar,size_t>(
-            nrow, ncol, VectorSize, row, j, k);
-        val.fastAccessCoeff(k) = v;
+  {
+    auto x_view = x->getLocalViewHost(Tpetra::Access::OverwriteAll);
+    for (size_t i=0; i<num_my_row; ++i) {
+      const GlobalOrdinal row = myGIDs[i];
+      for (size_t j=0; j<ncol; ++j) {
+        for (LocalOrdinal k=0; k<VectorSize; ++k) {
+          BaseScalar v =
+            generate_multi_vector_coefficient<BaseScalar,size_t>(
+              nrow, ncol, VectorSize, row, j, k);
+          val.fastAccessCoeff(k) = v;
+        }
+        x_view(i,j) = val;
       }
-      x_view[j][i] = val;
     }
   }
-  x_view = Teuchos::null;
 
   // matrix->describe(*(Teuchos::fancyOStream(rcp(&std::cout,false))),
   //                  Teuchos::VERB_EXTREME);
@@ -851,7 +853,7 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
   //             Teuchos::VERB_EXTREME);
 
   // Check
-  ArrayRCP< ArrayRCP<Scalar> > y_view = y->get2dViewNonConst();
+  auto y_view = y->getLocalViewHost(Tpetra::Access::ReadOnly);
   BaseScalar tol = 1.0e-14;
   for (size_t i=0; i<num_my_row; ++i) {
     const GlobalOrdinal row = myGIDs[i];
@@ -872,9 +874,9 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
           val.fastAccessCoeff(k) += v1*v2;
         }
       }
-      TEST_EQUALITY( y_view[j][i].size(), VectorSize );
+      TEST_EQUALITY( y_view(i,j).size(), VectorSize );
       for (LocalOrdinal k=0; k<VectorSize; ++k)
-        TEST_FLOATING_EQUALITY( y_view[j][i].fastAccessCoeff(k),
+        TEST_FLOATING_EQUALITY( y_view(i,j).fastAccessCoeff(k),
                                 val.fastAccessCoeff(k), tol );
     }
   }
@@ -958,13 +960,15 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
 
   // Fill vector
   RCP<Tpetra_Vector> x = Tpetra::createVector<Scalar>(map);
-  ArrayRCP<Scalar> x_view = x->get1dViewNonConst();
-  for (size_t i=0; i<num_my_row; ++i) {
-    const GlobalOrdinal row = myGIDs[i];
-    for (LocalOrdinal j=0; j<VectorSize; ++j)
-      val.fastAccessCoeff(j) = generate_vector_coefficient<BaseScalar,size_t>(
-        nrow, VectorSize, row, j);
-    x_view[i] = val;
+  {
+    auto x_view = x->getLocalViewHost(Tpetra::Access::OverwriteAll);
+    for (size_t i=0; i<num_my_row; ++i) {
+      const GlobalOrdinal row = myGIDs[i];
+      for (LocalOrdinal j=0; j<VectorSize; ++j)
+        val.fastAccessCoeff(j) = generate_vector_coefficient<BaseScalar,size_t>(
+          nrow, VectorSize, row, j);
+        x_view(i,0) = val;
+    }
   }
 
   // Multiply
@@ -980,25 +984,27 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
 
   // Multiply with flattened matix
   RCP<Tpetra_Vector> y2 = Tpetra::createVector<Scalar>(map);
-  RCP<Flat_Tpetra_Vector> flat_x =
-    Stokhos::create_flat_vector_view(*x, flat_x_map);
-  RCP<Flat_Tpetra_Vector> flat_y =
-    Stokhos::create_flat_vector_view(*y2, flat_y_map);
-  flat_matrix->apply(*flat_x, *flat_y);
+  {
+    RCP<Flat_Tpetra_Vector> flat_x =
+      Stokhos::create_flat_vector_view(*x, flat_x_map);
+    RCP<Flat_Tpetra_Vector> flat_y =
+      Stokhos::create_flat_vector_view(*y2, flat_y_map);
+    flat_matrix->apply(*flat_x, *flat_y);
+  }
 
   // flat_y->describe(*(Teuchos::fancyOStream(rcp(&std::cout,false))),
   //                  Teuchos::VERB_EXTREME);
 
   // Check
   BaseScalar tol = 1.0e-14;
-  ArrayRCP<Scalar> y_view = y->get1dViewNonConst();
-  ArrayRCP<Scalar> y2_view = y2->get1dViewNonConst();
+  auto y_view  = y-> getLocalViewHost(Tpetra::Access::ReadOnly);
+  auto y2_view = y2->getLocalViewHost(Tpetra::Access::ReadOnly);
   for (size_t i=0; i<num_my_row; ++i) {
-    TEST_EQUALITY( y_view[i].size(), VectorSize );
-    TEST_EQUALITY( y2_view[i].size(), VectorSize );
+    TEST_EQUALITY( y_view( i,0).size(), VectorSize );
+    TEST_EQUALITY( y2_view(i,0).size(), VectorSize );
     for (LocalOrdinal j=0; j<VectorSize; ++j)
-      TEST_FLOATING_EQUALITY( y_view[i].fastAccessCoeff(j),
-                              y2_view[i].fastAccessCoeff(j), tol );
+      TEST_FLOATING_EQUALITY( y_view( i,0).fastAccessCoeff(j),
+                              y2_view(i,0).fastAccessCoeff(j), tol );
   }
 }
 
@@ -1127,18 +1133,21 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
 
   // Fill RHS vector
   RCP<Tpetra_Vector> b = Tpetra::createVector<Scalar>(map);
-  ArrayRCP<Scalar> b_view = b->get1dViewNonConst();
-  Scalar b_val(VectorSize, BaseScalar(0.0));
-  for (LocalOrdinal j=0; j<VectorSize; ++j) {
-    b_val.fastAccessCoeff(j) =
-      BaseScalar(-1.0) + BaseScalar(j) / BaseScalar(VectorSize);
-  }
-  for (size_t i=0; i<num_my_row; ++i) {
-    const GlobalOrdinal row = myGIDs[i];
-    if (row == 0 || row == nrow-1)
-      b_view[i] = Scalar(0.0);
-    else
-      b_view[i] = -Scalar(b_val * h * h);
+  Scalar b_val;
+  {
+    auto b_view = b->getLocalViewHost(Tpetra::Access::OverwriteAll);
+    b_val = Scalar(VectorSize, BaseScalar(0.0));
+    for (LocalOrdinal j=0; j<VectorSize; ++j) {
+      b_val.fastAccessCoeff(j) =
+        BaseScalar(-1.0) + BaseScalar(j) / BaseScalar(VectorSize);
+    }
+    for (size_t i=0; i<num_my_row; ++i) {
+      const GlobalOrdinal row = myGIDs[i];
+      if (row == 0 || row == nrow-1)
+        b_view(i,0) = Scalar(0.0);
+      else
+        b_view(i,0) = -Scalar(b_val * h * h);
+    }
   }
 
   // Solve
@@ -1158,7 +1167,7 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
 
   // Check -- For a*y'' = b, correct answer is y = 0.5 *(b/a) * x * (x-1)
   btol = 1000*btol;
-  ArrayRCP<Scalar> x_view = x->get1dViewNonConst();
+  auto x_view = x->getLocalViewHost(Tpetra::Access::ReadOnly);
   Scalar val(VectorSize, BaseScalar(0.0));
   for (size_t i=0; i<num_my_row; ++i) {
     const GlobalOrdinal row = myGIDs[i];
@@ -1167,10 +1176,10 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
       val.fastAccessCoeff(j) =
         BaseScalar(0.5) * (b_val.coeff(j)/a_val.coeff(j)) * xx * (xx - BaseScalar(1.0));
     }
-    TEST_EQUALITY( x_view[i].size(), VectorSize );
+    TEST_EQUALITY( x_view(i,0).size(), VectorSize );
 
     // Set small values to zero
-    Scalar v = x_view[i];
+    Scalar v = x_view(i,0);
     for (LocalOrdinal j=0; j<VectorSize; ++j) {
       if (BST::abs(v.coeff(j)) < btol)
         v.fastAccessCoeff(j) = BaseScalar(0.0);
@@ -1269,18 +1278,21 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
 
   // Fill RHS vector
   RCP<Tpetra_Vector> b = Tpetra::createVector<Scalar>(map);
-  ArrayRCP<Scalar> b_view = b->get1dViewNonConst();
-  Scalar b_val(VectorSize, BaseScalar(0.0));
-  for (LocalOrdinal j=0; j<VectorSize; ++j) {
-    b_val.fastAccessCoeff(j) =
-      BaseScalar(-1.0) + BaseScalar(j) / BaseScalar(VectorSize);
-  }
-  for (size_t i=0; i<num_my_row; ++i) {
-    const GlobalOrdinal row = myGIDs[i];
-    if (row == 0 || row == nrow-1)
-      b_view[i] = Scalar(0.0);
-    else
-      b_view[i] = -Scalar(b_val * h * h);
+  Scalar b_val;
+  {
+    auto b_view = b->getLocalViewHost(Tpetra::Access::OverwriteAll);
+    b_val = Scalar(VectorSize, BaseScalar(0.0));
+    for (LocalOrdinal j=0; j<VectorSize; ++j) {
+      b_val.fastAccessCoeff(j) =
+        BaseScalar(-1.0) + BaseScalar(j) / BaseScalar(VectorSize);
+    }
+    for (size_t i=0; i<num_my_row; ++i) {
+      const GlobalOrdinal row = myGIDs[i];
+      if (row == 0 || row == nrow-1)
+        b_view(i,0) = Scalar(0.0);
+      else
+        b_view(i,0) = -Scalar(b_val * h * h);
+    }
   }
 
   // Create preconditioner
@@ -1308,7 +1320,7 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
 
   // Check -- For a*y'' = b, correct answer is y = 0.5 *(b/a) * x * (x-1)
   btol = 1000*btol;
-  ArrayRCP<Scalar> x_view = x->get1dViewNonConst();
+  auto x_view = x->getLocalViewHost(Tpetra::Access::ReadOnly);
   Scalar val(VectorSize, BaseScalar(0.0));
   for (size_t i=0; i<num_my_row; ++i) {
     const GlobalOrdinal row = myGIDs[i];
@@ -1317,10 +1329,10 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
       val.fastAccessCoeff(j) =
         BaseScalar(0.5) * (b_val.coeff(j)/a_val.coeff(j)) * xx * (xx - BaseScalar(1.0));
     }
-    TEST_EQUALITY( x_view[i].size(), VectorSize );
+    TEST_EQUALITY( x_view(i,0).size(), VectorSize );
 
     // Set small values to zero
-    Scalar v = x_view[i];
+    Scalar v = x_view(i,0);
     for (LocalOrdinal j=0; j<VectorSize; ++j) {
       if (BST::magnitude(v.coeff(j)) < btol)
         v.fastAccessCoeff(j) = BaseScalar(0.0);
@@ -1417,9 +1429,11 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
 
   // Fill RHS vector
   RCP<Tpetra_Vector> b = Tpetra::createVector<Scalar>(map);
-  ArrayRCP<Scalar> b_view = b->get1dViewNonConst();
-  for (size_t i=0; i<num_my_row; ++i) {
-    b_view[i] = Scalar(1.0);
+  {
+    auto b_view = b->getLocalViewHost(Tpetra::Access::OverwriteAll);
+    for (size_t i=0; i<num_my_row; ++i) {
+      b_view(i,0) = Scalar(1.0);
+    }
   }
 
   // Solve
@@ -1460,7 +1474,7 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
   //     [ 1, 1/2, ..., 1/VectorSize ]
   //     ....
   tol = 1000*tol;
-  ArrayRCP<Scalar> x_view = x->get1dViewNonConst();
+  auto x_view = x->getLocalViewHost(Tpetra::Access::ReadOnly);
   for (size_t i=0; i<num_my_row; ++i) {
     const GlobalOrdinal row = myGIDs[i];
     if (row % 2) {
@@ -1470,10 +1484,10 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
     }
     else
       val = Scalar(VectorSize, BaseScalar(0.0));
-    TEST_EQUALITY( x_view[i].size(), VectorSize );
+    TEST_EQUALITY( x_view(i,0).size(), VectorSize );
 
     // Set small values to zero
-    Scalar v = x_view[i];
+    Scalar v = x_view(i,0);
     for (LocalOrdinal j=0; j<VectorSize; ++j) {
       if (ST::magnitude(v.coeff(j)) < tol)
         v.fastAccessCoeff(j) = BaseScalar(0.0);
@@ -1548,15 +1562,17 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
   //     ...
 
   RCP<Tpetra_Vector> b = Tpetra::createVector<Scalar>(map);
-  ArrayRCP<Scalar> b_view = b->get1dViewNonConst();
-  for (size_t i=0; i<num_my_row; ++i) {
-    const GlobalOrdinal row = myGIDs[i];
-    b_view[i] = Scalar(0.0);
-    for (LocalOrdinal j=0; j<VectorSize; ++j)
-      if (int(j+2+row-VectorSize) > 0)
-        b_view[i].fastAccessCoeff(j) = BaseScalar(row+1);
+  {
+    auto b_view = b->getLocalViewHost(Tpetra::Access::OverwriteAll);
+    for (size_t i=0; i<num_my_row; ++i) {
+      const GlobalOrdinal row = myGIDs[i];
+      b_view(i,0) = Scalar(0.0);
+      for (LocalOrdinal j=0; j<VectorSize; ++j)
+        if (int(j+2+row-VectorSize) > 0)
+          b_view(i,0).fastAccessCoeff(j) = BaseScalar(row+1);
+    }
   }
-
+  
   // Solve
   typedef Teuchos::ScalarTraits<BaseScalar> ST;
 #ifdef HAVE_STOKHOS_ENSEMBLE_REDUCT
@@ -1617,10 +1633,10 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
   //     [ 0, 0, ..., 1, 1, 1]
   //     ...
   tol = 1000*tol;
-  ArrayRCP<Scalar> x_view = x->get1dViewNonConst();
+  auto x_view = x->getLocalViewHost(Tpetra::Access::ReadOnly);
   for (size_t i=0; i<num_my_row; ++i) {
     const GlobalOrdinal row = myGIDs[i];
-    Scalar v = x_view[i];
+    Scalar v = x_view(i,0);
 
     for (LocalOrdinal j=0; j<VectorSize; ++j) {
       if (ST::magnitude(v.coeff(j)) < tol)
@@ -1702,13 +1718,15 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
   //     ...
 
   RCP<Tpetra_Vector> b = Tpetra::createVector<Scalar>(map);
-  ArrayRCP<Scalar> b_view = b->get1dViewNonConst();
-  for (size_t i=0; i<num_my_row; ++i) {
-    const GlobalOrdinal row = myGIDs[i];
-    b_view[i] = Scalar(0.0);
-    for (LocalOrdinal j=0; j<VectorSize; ++j)
-      if (int(j+2+row-VectorSize) > 0)
-        b_view[i].fastAccessCoeff(j) = BaseScalar(row+1);
+  {
+    auto b_view = b->getLocalViewHost(Tpetra::Access::OverwriteAll);
+    for (size_t i=0; i<num_my_row; ++i) {
+      const GlobalOrdinal row = myGIDs[i];
+      b_view(i,0) = Scalar(0.0);
+      for (LocalOrdinal j=0; j<VectorSize; ++j)
+        if (int(j+2+row-VectorSize) > 0)
+          b_view(i,0).fastAccessCoeff(j) = BaseScalar(row+1);
+    }
   }
 
   // Solve
@@ -1771,10 +1789,10 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
   //     [ 0, 0, ..., 1, 1, 1]
   //     ...
   tol = 1000*tol;
-  ArrayRCP<Scalar> x_view = x->get1dViewNonConst();
+  auto x_view = x->getLocalViewHost(Tpetra::Access::ReadOnly);
   for (size_t i=0; i<num_my_row; ++i) {
     const GlobalOrdinal row = myGIDs[i];
-    Scalar v = x_view[i];
+    Scalar v = x_view(i,0);
 
     for (LocalOrdinal j=0; j<VectorSize; ++j) {
       if (ST::magnitude(v.coeff(j)) < tol)
@@ -1856,13 +1874,15 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
   //     ...
 
   RCP<Tpetra_Vector> b = Tpetra::createVector<Scalar>(map);
-  ArrayRCP<Scalar> b_view = b->get1dViewNonConst();
-  for (size_t i=0; i<num_my_row; ++i) {
-    const GlobalOrdinal row = myGIDs[i];
-    b_view[i] = Scalar(0.0);
-    for (LocalOrdinal j=0; j<VectorSize; ++j)
-      if (int(j+2+row-VectorSize) > 0)
-        b_view[i].fastAccessCoeff(j) = BaseScalar(row+1);
+  {
+    auto b_view = b->getLocalViewHost(Tpetra::Access::OverwriteAll);
+    for (size_t i=0; i<num_my_row; ++i) {
+      const GlobalOrdinal row = myGIDs[i];
+      b_view(i,0) = Scalar(0.0);
+      for (LocalOrdinal j=0; j<VectorSize; ++j)
+        if (int(j+2+row-VectorSize) > 0)
+          b_view(i,0).fastAccessCoeff(j) = BaseScalar(row+1);
+    }
   }
 
   // Solve
@@ -1925,10 +1945,10 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
   //     [ 0, 0, ..., 1, 1, 1]
   //     ...
   tol = 1000*tol;
-  ArrayRCP<Scalar> x_view = x->get1dViewNonConst();
+  auto x_view = x->getLocalViewHost(Tpetra::Access::ReadOnly);
   for (size_t i=0; i<num_my_row; ++i) {
     const GlobalOrdinal row = myGIDs[i];
-    Scalar v = x_view[i];
+    Scalar v = x_view(i,0);
 
     for (LocalOrdinal j=0; j<VectorSize; ++j) {
       if (ST::magnitude(v.coeff(j)) < tol)
@@ -2043,9 +2063,11 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
 
   // Fill RHS vector
   RCP<Tpetra_Vector> b = Tpetra::createVector<Scalar>(map);
-  ArrayRCP<Scalar> b_view = b->get1dViewNonConst();
-  for (size_t i=0; i<num_my_row; ++i) {
-    b_view[i] = Scalar(1.0);
+  {
+    auto b_view = b->getLocalViewHost(Tpetra::Access::OverwriteAll);
+    for (size_t i=0; i<num_my_row; ++i) {
+      b_view(i,0) = Scalar(1.0);
+    }
   }
 
   // Create preconditioner
@@ -2095,7 +2117,7 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
   //     [ 1, 1/2, ..., 1/VectorSize ]
   //     ....
   tol = 1000*tol;
-  ArrayRCP<Scalar> x_view = x->get1dViewNonConst();
+  auto x_view = x->getLocalViewHost(Tpetra::Access::ReadOnly);
   for (size_t i=0; i<num_my_row; ++i) {
     const GlobalOrdinal row = myGIDs[i];
     if (row % 2) {
@@ -2105,10 +2127,10 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
     }
     else
       val = Scalar(VectorSize, BaseScalar(0.0));
-    TEST_EQUALITY( x_view[i].size(), VectorSize );
+    TEST_EQUALITY( x_view(i,0).size(), VectorSize );
 
     // Set small values to zero
-    Scalar v = x_view[i];
+    Scalar v = x_view(i,0);
     for (LocalOrdinal j=0; j<VectorSize; ++j) {
       if (ST::magnitude(v.coeff(j)) < tol)
         v.fastAccessCoeff(j) = BaseScalar(0.0);
@@ -2212,18 +2234,21 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
 
   // Fill RHS vector
   RCP<Tpetra_Vector> b = Tpetra::createVector<Scalar>(map);
-  ArrayRCP<Scalar> b_view = b->get1dViewNonConst();
-  Scalar b_val(VectorSize, BaseScalar(0.0));
-  for (LocalOrdinal j=0; j<VectorSize; ++j) {
-    b_val.fastAccessCoeff(j) =
-      BaseScalar(-1.0) + BaseScalar(j) / BaseScalar(VectorSize);
-  }
-  for (size_t i=0; i<num_my_row; ++i) {
-    const GlobalOrdinal row = myGIDs[i];
-    if (row == 0 || row == nrow-1)
-      b_view[i] = Scalar(0.0);
-    else
-      b_view[i] = -Scalar(b_val * h * h);
+  Scalar b_val;
+  {
+    auto b_view = b->getLocalViewHost(Tpetra::Access::OverwriteAll);
+    b_val = Scalar(VectorSize, BaseScalar(0.0));
+    for (LocalOrdinal j=0; j<VectorSize; ++j) {
+      b_val.fastAccessCoeff(j) =
+        BaseScalar(-1.0) + BaseScalar(j) / BaseScalar(VectorSize);
+    }
+    for (size_t i=0; i<num_my_row; ++i) {
+      const GlobalOrdinal row = myGIDs[i];
+      if (row == 0 || row == nrow-1)
+        b_view(i,0) = Scalar(0.0);
+      else
+        b_view(i,0) = -Scalar(b_val * h * h);
+    }
   }
 
   // Create preconditioner
@@ -2280,7 +2305,7 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
 
   // Check -- For a*y'' = b, correct answer is y = 0.5 *(b/a) * x * (x-1)
   tol = 1000*tol;
-  ArrayRCP<Scalar> x_view = x->get1dViewNonConst();
+  auto x_view = x->getLocalViewHost(Tpetra::Access::ReadOnly);
   Scalar val(VectorSize, BaseScalar(0.0));
   for (size_t i=0; i<num_my_row; ++i) {
     const GlobalOrdinal row = myGIDs[i];
@@ -2289,10 +2314,10 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
       val.fastAccessCoeff(j) =
         BaseScalar(0.5) * (b_val.coeff(j)/a_val.coeff(j)) * xx * (xx - BaseScalar(1.0));
     }
-    TEST_EQUALITY( x_view[i].size(), VectorSize );
+    TEST_EQUALITY( x_view(i,0).size(), VectorSize );
 
     // Set small values to zero
-    Scalar v = x_view[i];
+    Scalar v = x_view(i,0);
     for (LocalOrdinal j=0; j<VectorSize; ++j) {
       if (ST::magnitude(v.coeff(j)) < tol)
         v.fastAccessCoeff(j) = BaseScalar(0.0);
@@ -2398,18 +2423,21 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
 
   // Fill RHS vector
   RCP<Tpetra_Vector> b = Tpetra::createVector<Scalar>(map);
-  ArrayRCP<Scalar> b_view = b->get1dViewNonConst();
-  Scalar b_val(VectorSize, BaseScalar(0.0));
-  for (LocalOrdinal j=0; j<VectorSize; ++j) {
-    b_val.fastAccessCoeff(j) =
-      BaseScalar(-1.0) + BaseScalar(j) / BaseScalar(VectorSize);
-  }
-  for (size_t i=0; i<num_my_row; ++i) {
-    const GlobalOrdinal row = myGIDs[i];
-    if (row == 0 || row == nrow-1)
-      b_view[i] = Scalar(0.0);
-    else
-      b_view[i] = -Scalar(b_val * h * h);
+  Scalar b_val;
+  {
+    auto b_view = b->getLocalViewHost(Tpetra::Access::OverwriteAll);
+    b_val = Scalar(VectorSize, BaseScalar(0.0));
+    for (LocalOrdinal j=0; j<VectorSize; ++j) {
+      b_val.fastAccessCoeff(j) =
+        BaseScalar(-1.0) + BaseScalar(j) / BaseScalar(VectorSize);
+    }
+    for (size_t i=0; i<num_my_row; ++i) {
+      const GlobalOrdinal row = myGIDs[i];
+      if (row == 0 || row == nrow-1)
+        b_view(i,0) = Scalar(0.0);
+      else
+        b_view(i,0) = -Scalar(b_val * h * h);
+    }
   }
 
   // Solve
@@ -2448,7 +2476,7 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
   // Check -- For a*y'' = b, correct answer is y = 0.5 *(b/a) * x * (x-1)
   typedef Teuchos::ScalarTraits<BaseScalar> ST;
   typename ST::magnitudeType tol = 1e-9;
-  ArrayRCP<Scalar> x_view = x->get1dViewNonConst();
+  auto x_view = x->getLocalViewHost(Tpetra::Access::ReadOnly);
   Scalar val(VectorSize, BaseScalar(0.0));
   for (size_t i=0; i<num_my_row; ++i) {
     const GlobalOrdinal row = myGIDs[i];
@@ -2457,10 +2485,10 @@ TEUCHOS_UNIT_TEST_TEMPLATE_4_DECL(
       val.fastAccessCoeff(j) =
         BaseScalar(0.5) * (b_val.coeff(j)/a_val.coeff(j)) * xx * (xx - BaseScalar(1.0));
     }
-    TEST_EQUALITY( x_view[i].size(), VectorSize );
+    TEST_EQUALITY( x_view(i,0).size(), VectorSize );
 
     // Set small values to zero
-    Scalar v = x_view[i];
+    Scalar v = x_view(i,0);
     for (LocalOrdinal j=0; j<VectorSize; ++j) {
       if (ST::magnitude(v.coeff(j)) < tol)
         v.fastAccessCoeff(j) = BaseScalar(0.0);