Dia/Morpheus_Copy_Impl.hpp

 
#ifndef MORPHEUS_DIA_COPY_IMPL_HPP
#define MORPHEUS_DIA_COPY_IMPL_HPP
 
#include <Morpheus_FormatTags.hpp>
 
#include <impl/Morpheus_Utils.hpp>
#include <impl/DenseVector/Morpheus_Copy_Impl.hpp>
#include <impl/DenseMatrix/Morpheus_Copy_Impl.hpp>
 
namespace Morpheus {
namespace Impl {
 
template <typename SourceType, typename DestinationType>
void copy(const SourceType& src, DestinationType& dst,
          typename std::enable_if_t<
              Morpheus::is_dia_matrix_format_container_v<SourceType> &&
              Morpheus::is_dia_matrix_format_container_v<DestinationType>>* =
              nullptr) {
  MORPHEUS_ASSERT(
      (dst.nrows() == src.nrows()) && (dst.ncols() == src.ncols()),
      "Destination matrix must have equal shape to the source matrix");
  MORPHEUS_ASSERT(dst.nnnz() == src.nnnz(),
                  "Destination matrix must have equal number of non-zeros to "
                  "the source matrix");
  MORPHEUS_ASSERT(
      dst.diagonal_offsets().size() == src.cdiagonal_offsets().size(),
      "Destination matrix must have equal number of diagonals to the source "
      "matrix");
 
  Morpheus::Impl::copy(src.cdiagonal_offsets(), dst.diagonal_offsets());
  Morpheus::Impl::copy(src.cvalues(), dst.values());
}
 
}  // namespace Impl
}  // namespace Morpheus
 
#endif  // MORPHEUS_DIA_COPY_IMPL_HPP