Question

这是一个晦涩难懂的问题，我真的不希望任何人回答，但我有这种方法需要（并返回）一个Eigen :: SparseMatrix。我想把它放到deal.ii库中，有没有办法从deal.ii / Eigen复制/转换SparseMatrix？我知道你可以将deal.ii复制到Trilinos SparseMatrix，例如：

#action-circle

Eigen :: SparseMatrix有类似的方式吗？我想Eigen在deal.ii中确实没有这种支持。所以也许有一些“蛮力”。类型方法，就像这种代码的尝试一样，显然不起作用：

`

  `SparseMatrix<double> matrix(sparsity);
...//fill matrix

  Epetra_Map map(TrilinosWrappers::types::int_type(5),
                 TrilinosWrappers::types::int_type(5),
                 0,
                 Utilities::Trilinos::comm_world());

  TrilinosWrappers::SparseMatrix tmatrix;
  tmatrix.reinit (map, map, matrix, 0, false);`

`

好的，所以我想出了如何转换来自dealii :: SparseMatrix - ＆gt;征::稀疏矩阵。

Eigen::SparseMatrix<double> ConvertToEigenMatrix(SparseMatrix<double> data)
{
    Eigen::SparseMatrix<double> eMatrix(data.m(), data.n());
    for (int i = 0; i < data.m(); ++i)
        eMatrix.row(i) =  Eigen::SparseMatrix<double> ::Map(&data[i][0], data.n());
    return eMatrix;

不，我只需要反过来说明。

Answer 1

这个问题已经很久了，但也许我仍然可以提供帮助。我是Deal.II开发者之一，我不记得在邮件列表上看到这个（这些问题比SO更活跃）。

deal.II中的SparseMatrix不存储自己的稀疏模式：相反，它存储指向SparsityPattern对象的指针。您需要在特征矩阵上循环两次：一次设置SparsityPattern，第二次复制矩阵值。类似下面的内容似乎有效：

#include <deal.II/lac/dynamic_sparsity_pattern.h>
#include <deal.II/lac/sparsity_pattern.h>
#include <deal.II/lac/sparse_matrix.h>

#include <eigen3/Eigen/Sparse>

#include <iostream>

int main()
{
  const std::size_t shape = 3;
  Eigen::SparseMatrix<double> matrix(shape, shape);
  matrix.insert(0, 0) = 1.0;
  matrix.insert(0, 1) = 2.0;
  matrix.insert(0, 2) = 1.0;
  matrix.insert(2, 2) = 2.0;
  matrix.makeCompressed();

  {
    dealii::SparsityPattern sparsity_pattern(matrix.rows(), matrix.cols());
    dealii::DynamicSparsityPattern dynamic_sparsity_pattern(matrix.rows(), matrix.cols());

    for (decltype(matrix.outerSize()) row_n = 0; row_n < matrix.outerSize(); ++row_n)
      for (Eigen::SparseMatrix<double>::InnerIterator it(matrix, row_n); it; ++it)
        dynamic_sparsity_pattern.add(it.row(), it.col());

    sparsity_pattern.copy_from(dynamic_sparsity_pattern);
    dealii::SparseMatrix<double> matrix2(sparsity_pattern);

    for (decltype(matrix.outerSize()) row_n = 0; row_n < matrix.outerSize(); ++row_n)
      for (Eigen::SparseMatrix<double>::InnerIterator it(matrix, row_n); it; ++it)
        matrix2.set(it.row(), it.col(), it.value());

    matrix2.print(std::cout); // prints the right matrix
  }
}

您还必须管理SparsityPattern对象的生命周期。

deal.II不使用CSR或CSC：它使用自己的类似CSR的格式，其中主对角线上的条目首先存储在包含该行的矩阵条目的数组中，因此我们确实需要复制迭代器接口。

将Eigen :: SparseMatrix <double>转换为deal.ii :: SparseMatrix <double>？

1 个答案: