я использую Spectra C ++ библиотека, которая построена на вершине Eigen C ++ библиотека.
Мне нужно будет построить правильный SymGEigsSolver
объект на основе данного SelectionRule
а также GeigsMode
параметры в методе. Определение SymGEigsSolver
как следует:
template<typename Scalar, int SelectionRule, typename OpType, typename BOpType, int GEigsMode>
class Spectra::SymGEigsSolver< Scalar, SelectionRule, OpType, BOpType, GEigsMode >
Однако, когда я попытался скомпилировать следующую функцию:
typedef SparseMatrix<double, Eigen::RowMajor> SpMat;
typedef Spectra::SparseSymMatProd<double, Eigen::Upper, Eigen::RowMajor> SpSymMatProd;
typedef Spectra::SparseCholesky<double, Eigen::Upper> SpCholesky;
typedef Spectra::SparseRegularInverse<double, Eigen::Upper, Eigen::RowMajor> SpRegularInverse;template < typename Scalar,
int SelectionRule,
typename OpType,
typename BOpType,
int GEigsMode >
SymGEigsSolver<Scalar, SelectionRule, OpType, BOpType, GEigsMode>& CreateSolverEngine
(const SpMat& A, const SpMat& B, int selectRule, int gMode, int nMode, int ncv)
{
SpSymMatProd op(A);
if (gMode == Spectra::GEIGS_CHOLESKY)
{
SpCholesky Bop(B);
SymGEigsSolver<Scalar, SelectionRule, SpSymMatProd, SpCholesky, GEigsMode>
ges1(op, Bop, nMode, ncv);
return ges1;
}
if (gMode == Spectra::GEIGS_REGULAR_INVERSE)
{
SpRegularInverse Bop2(B);
SymGEigsSolver<Scalar, SelectionRule, SpSymMatProd, SpRegularInverse, GEigsMode>
ges(op, Bop2, nMode, ncv);
return ges;
}
throw std::out_of_range("out of range for "+gMode);}
Я получил ошибки компиляции, такие как:
Error C2664 'Spectra::SymGEigsSolver<Scalar,SelectionRule,OpType,BOpType,GEigsMode>::SymGEigsSolver(OpType *,BOpType *,int,int)'
: cannot convert parameter 1
from 'SpSymMatProd' to 'SpSymMatProd *'
А также
C2440 'return' : cannot convert from
'Spectra::SymGEigsSolver<Scalar,SelectionRule,OpType,BOpType,GEigsMode>'
to
'Spectra::SymGEigsSolver<Scalar,SelectionRule,OpType,BOpType,GEigsMode>&'
А также
C2664 Spectra::SymGEigsSolver<Scalar,SelectionRule,OpType,BOpType,GEigsMode>::SymGEigsSolver(OpType *,BOpType *,int,int)'
: cannot convert parameter 1 from 'SpSymMatProd' to 'SpSymMatProd *'
Я сам разобрался с ответом — я опечатка &
на имя метода. Это правильное решение:
typedef SparseMatrix<double, Eigen::RowMajor> SpMat;
typedef Spectra::SparseSymMatProd<double, Eigen::Upper, Eigen::RowMajor> SpSymMatProd;
typedef Spectra::SparseCholesky<double, Eigen::Upper> SpCholesky;
typedef Spectra::SparseRegularInverse<double, Eigen::Upper, Eigen::RowMajor> SpRegularInverse;template < typename Scalar,
int SelectionRule,
typename OpType,
typename BOpType,
int GEigsMode >
SymGEigsSolver<Scalar, SelectionRule, OpType, BOpType, GEigsMode> CreateSolverEngine
(const SpMat& A, const SpMat& B, int selectRule, int gMode, int nMode, int ncv)
{
SpSymMatProd op(A);
if (gMode == Spectra::GEIGS_CHOLESKY)
{
SpCholesky Bop(B);
SymGEigsSolver<Scalar, SelectionRule, SpSymMatProd, SpCholesky, GEigsMode>
ges1(op, Bop, nMode, ncv);
return ges1;
}
if (gMode == Spectra::GEIGS_REGULAR_INVERSE)
{
SpRegularInverse Bop2(B);
SymGEigsSolver<Scalar, SelectionRule, SpSymMatProd, SpRegularInverse, GEigsMode>
ges(op, Bop2, nMode, ncv);
return ges;
}
throw std::out_of_range("out of range for "+gMode);}
Других решений пока нет …