solve (const PETScWrappers::MatrixBase &A,
std::vector<double> &kr,
std::vector<OutputVector> &vr,
- const unsigned int n_eigenvectors);
+ const unsigned int n_eigenvectors = 1);
/**
* Same as above, but here a composite method for solving the
const PETScWrappers::MatrixBase &B,
std::vector<double> &kr,
std::vector<OutputVector> &vr,
- const unsigned int n_eigenvectors);
+ const unsigned int n_eigenvectors = 1);
/**
* Initialize solver for the linear system $Ax=\lambda x$. (Note:
SolverBase::solve (const PETScWrappers::MatrixBase &A,
std::vector<double> &kr,
std::vector<OutputVector> &vr,
- const unsigned int n_eigenvectors = 1)
+ const unsigned int n_eigenvectors)
{
// Panic if no eigenpairs are wanted.
AssertThrow (n_eigenvectors != 0, ExcSLEPcWrappersUsageError());
const PETScWrappers::MatrixBase &B,
std::vector<double> &kr,
std::vector<OutputVector> &vr,
- const unsigned int n_eigenvectors = 1)
+ const unsigned int n_eigenvectors)
{
// Panic if no eigenpairs are wanted.
AssertThrow (n_eigenvectors != 0, ExcSLEPcWrappersUsageError());
#endif
/*---------------------------- slepc_solver.h ---------------------------*/
-