search

Pressio / pressio

core C++ library
Other
45 stars 7 forks source link

ops: level3 tpetra #474

Closed fnrizzi closed 1 year ago

fnrizzi commented 1 year ago

fix this 

-  const auto numVecsA = A.getNumVectors();
-  const auto numVecsB = B.getNumVectors();
-  assert( (std::size_t)::pressio::ops::extent(A,0) == (std::size_t)::pressio::ops::extent(B,0));
-  assert( (std::size_t)::pressio::ops::extent(C,0) == (std::size_t)numVecsA );
-  assert( (std::size_t)::pressio::ops::extent(C,1) == (std::size_t)numVecsB );
+  // const auto numVecsA = A.getNumVectors();
+  // const auto numVecsB = B.getNumVectors();
+  // assert( (std::size_t)::pressio::ops::extent(A,0) == (std::size_t)::pressio::ops::extent(B,0));
+  // assert( (std::size_t)::pressio::ops::extent(C,0) == (std::size_t)numVecsA );
+  // assert( (std::size_t)::pressio::ops::extent(C,1) == (std::size_t)numVecsB );

   // for col-major, we should have inner stride == 1 and outerstride = num of rows
   // we should actually check at compile time somehow if C can be viewed in a kokkos view