eigendunetransformations.hh

Go to the documentation of this file.

// SPDX-FileCopyrightText: 2021-2024 The Ikarus Developers mueller@ibb.uni-stuttgart.de
// SPDX-License-Identifier: LGPL-3.0-or-later
 
#pragma once
 
#include <assert.h>
 
#include <dune/common/diagonalmatrix.hh>
#include <dune/common/fmatrix.hh>
#include <dune/common/fvector.hh>
#include <dune/istl/scaledidmatrix.hh>
 
#include <Eigen/Core>
 
namespace Ikarus {
 
template <typename ST, int size>
Dune::FieldVector<ST, size> toDune(const Eigen::Vector<ST, size>& vec) {
  Dune::FieldVector<ST, size> fieldVector;
  for (int i = 0; i < size; ++i)
    fieldVector[i] = vec[i];
  return fieldVector;
}
 
template <typename ST, int rows>
Dune::FieldVector<ST, rows> toDune(const Eigen::Matrix<ST, rows, 0>& vec) {
  assert(vec.cols() == 1 && "The passed matrix needs to have a single column.");
  Dune::FieldVector<ST, rows> fieldVector{0.0};
 
  for (int i = 0; i < vec.rows(); ++i)
    fieldVector[i] = vec(i, 0);
  return fieldVector;
}
 
template <typename ST, int rows, int cols>
Dune::FieldMatrix<ST, rows, cols> toDune(const Eigen::Matrix<ST, rows, cols>& mat) {
  Dune::FieldMatrix<ST, rows, cols> fieldMatrix{0.0};
 
  for (int i = 0; i < mat.rows(); ++i)
    for (int j = 0; j < mat.cols(); ++j)
      fieldMatrix[i][j] = mat(i, j);
  return fieldMatrix;
}
 
template <typename ST, int size>
Eigen::Map<const Eigen::Vector<ST, size>> toEigenMap(const Dune::FieldVector<ST, size>& vec) {
  return {vec.data(), size};
}
 
template <typename ST, int size>
Eigen::Map<Eigen::Vector<ST, size>> toEigenMap(Dune::FieldVector<ST, size>& vec) {
  return {vec.data(), size};
}
 
} // namespace Ikarus