When a DiaMatrix or EllMatrix is constructed with Kokkos::LayoutRight returns invalid results on Morpheus::Cuda and Morpheus::HIP for multiply() and update_diagonal() algorithms, due to improper construction of the linearized index in each of the following kernels:
When a
DiaMatrixorEllMatrixis constructed withKokkos::LayoutRightreturns invalid results onMorpheus::CudaandMorpheus::HIPformultiply()andupdate_diagonal()algorithms, due to improper construction of the linearized index in each of the following kernels:update_dia_diagonal_kernel()incore/src/impl/Dia/Kernels/Morpheus_MatrixOperations_Impl.hppspmv_dia_kernel()incore/src/impl/Dia/Kernels/Morpheus_Multiply_Impl.hppupdate_ell_diagonal_kernel()incore/src/impl/Ell/Kernels/Morpheus_MatrixOperations_Impl.hppspmv_ell_kernel()incore/src/impl/Ell/Kernels/Morpheus_MatrixOperations_Impl.hpp