OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

Allow VexCL to use runtime parameters

Browse Source
This commit is contained in:
Tong Dong Qiu 2021-07-05 10:42:25 +02:00
parent b4eb711c71
commit 47ffeaf24a
1 changed files with 3 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
opm/simulators/linalg/bda/amgclSolverBackend.cpp
View File

@ -245,16 +245,13 @@ void amgclSolverBackend<block_size>::solve_system(double *b, WellContributions &
vex::scoped_program_header h1(ctx, amgcl::backend::vexcl_static_matrix_declaration<double, block_size>());
typedef amgcl::backend::vexcl<dmat_type> Backend;
typedef amgcl::make_block_solver<
amgcl::relaxation::as_preconditioner<Backend, amgcl::relaxation::ilu0>,
amgcl::solver::bicgstab<Backend>
> Solver;
typedef amgcl::make_solver<amgcl::runtime::preconditioner<Backend>, amgcl::runtime::solver::wrapper<Backend> > Solver;
typename Backend::params bprm;
bprm.q = ctx; // set vexcl context
auto A = std::tie(N, A_rows, A_cols, A_vals);
auto Atmp = std::tie(N, A_rows, A_cols, A_vals);
auto A = amgcl::adapter::block_matrix<dmat_type>(Atmp);
Solver solve(A, prm, bprm);