Morpheus_MatrixOperations.hpp

 
#ifndef MORPHEUS_MATRIXOPERATIONS_HPP
#define MORPHEUS_MATRIXOPERATIONS_HPP
 
#include <impl/Morpheus_MatrixOperations_Impl.hpp>
#include <impl/Dynamic/Morpheus_MatrixOperations_Impl.hpp>
 
namespace Morpheus {
template <typename ExecSpace, typename Matrix, typename Vector>
void update_diagonal(Matrix& A, const Vector& diagonal) {
  Impl::update_diagonal<ExecSpace>(A, diagonal);
}
 
template <typename ExecSpace, typename Matrix, typename Vector>
void get_diagonal(const Matrix& A, Vector& diagonal) {
  Impl::get_diagonal<ExecSpace>(A, diagonal);
}
 
template <typename ExecSpace, typename Matrix, typename SizeType,
          typename ValueType>
void set_value(Matrix& A, SizeType row, SizeType col, ValueType value) {
  Impl::set_value(A, row, col, value);
}
 
template <typename ExecSpace, typename Matrix, typename IndexVector,
          typename ValueVector>
void set_values(Matrix& A, typename IndexVector::value_type m,
                const IndexVector idxm, typename IndexVector::value_type n,
                const IndexVector idxn, ValueVector values) {
  Impl::set_values<ExecSpace>(A, m, idxm, n, idxn, values);
}
 
template <typename ExecSpace, typename Matrix, typename TransposeMatrix>
void transpose(const Matrix& A, TransposeMatrix& At) {
  Impl::transpose<ExecSpace>(A, At);
}
}  // namespace Morpheus
 
#endif  // MORPHEUS_MATRIXOPERATIONS_HPP