这是一个晦涩的问题,我并不真正希望有人回答,但是我有这种方法(并返回)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-> eigen :: 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;
}
不,我只需要找出反面。
这个问题是旧的,但也许我仍然可以帮忙。我是Deal.ii开发人员,我不记得在邮件列表中看到此内容(这对这些类型的问题比这要活跃得多)。
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的格式,其中主角上的条目首先存储在包含该行的矩阵条目的数组中,因此我们确实需要使用迭代器接口。