PolySolve

This library contains a cross-platform Eigen wrapper for many different external linear solvers including (but not limited to):

• CHOLMOD
• Hypre
• AMGCL
• Pardiso

Example Usage

const std::string solver_name = "Hypre"
auto solver = LinearSolver::create(solver_name, "");

// Configuration parameters like iteration or accuracy for iterative solvers
// solver->setParameters(params);

// System sparse matrix
Eigen::SparseMatrix<double> A;

// Right-hand side
Eigen::VectorXd b;

// Solution
Eigen::VectorXd x(b.size());

solver->analyzePattern(A, A.rows());
solver->factorize(A);
solver->solve(b, x);

You can use LinearSolver::availableSolvers() to obtain the list of available solvers.

Parameters

Iterative solvers (AMGCL, Eigen Internal Solvers, HYPRE)

• max_iter controls the solver's iterations, default 1000
• conv_tol, tolerance controls the convergence tolerance, default 1e-10

Hypre Only

• pre_max_iter, number of pre iterations, default 1

AMGCL Only

The default parameters of the AMGCL solver are:

{
    "precond": {
        "relax": {
            "degree": 16,
            "type": "chebyshev",
            "power_iters": 100,
            "higher": 2,
            "lower": 0.008333333333,
            "scale": true
        },
        "class": "amg",
        "max_levels": 6,
        "direct_coarse": false,
        "ncycle": 2,
        "coarsening": {
            "type": "smoothed_aggregation",
            "estimate_spectral_radius": true,
            "relax": 1,
            "aggr": {
                "eps_strong": 0
            }
        }
    },
    "solver": {
        "tol": 1e-10,
        "maxiter": 1000,
        "type": "cg"
    }
}

For a more details and options refer to the AMGCL documentation.

Pardiso

mtype, sets the matrix type, default 11