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

Question

这是一个令人费解的问题，我真的不希望有人会回答，但我有这个方法，它接受(并返回)一个Eigen::SparseMatrix。我想把它放到deal.ii库中，有没有办法从deal.ii/Eigen复制/转换SparseMatrix？我知道您可以将deal.ii复制到Trilinos SparseMatrix中，如下所示：

  `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);`    

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

`

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; 

`

好的，我想出了如何从dealii::SparseMatrix转换成->：：SparseMatrix。

  SparseMatrix<double>::iterator smi = matrix.begin();
  SparseMatrix<double>::iterator smi_end = matrix.end();

  unsigned int row,col;
  double val;
  for (; smi!=smi_end; ++smi)
  {
       row = smi->row();
       col = smi->column();
       val = smi->value();

       spMat.insert(row, col) = val; 
       std::cout << val << std::endl; 
  }

不，我只需要找出相反的方向。