Commit bf18ee17 authored by Marcel Koch's avatar Marcel Koch
Browse files

fix ilu issues

parent 645004c0
Pipeline #34003 failed with stages
in 9 minutes and 15 seconds
...@@ -6,8 +6,8 @@ Parameters ...@@ -6,8 +6,8 @@ Parameters
elementType p1; elementType p1;
minLevel 0; minLevel 3;
maxLevel 2; maxLevel 3;
// ########## // ##########
// # domain # // # domain #
...@@ -59,13 +59,13 @@ Parameters ...@@ -59,13 +59,13 @@ Parameters
solverType gmg_wjac; solverType gmg_wjac;
coarseGridSolverType cg_ginkgo; coarseGridSolverType cg_ginkgo;
cgHytegVerbose false; cgHytegVerbose true;
preSmooth 3; preSmooth 3;
postSmooth 3; postSmooth 3;
relativeResidualReduction 1e-06; relativeResidualReduction 1e-06;
maxCoarseGridSolverIterations 1000000000; maxCoarseGridSolverIterations 1000000000;
gkoExecutor reference; gkoExecutor cuda;
useAgglomeration true; useAgglomeration false;
// ################ // ################
// # app settings # // # app settings #
......
...@@ -404,10 +404,13 @@ class GinkgoCGSolver : public Solver< OperatorType > ...@@ -404,10 +404,13 @@ class GinkgoCGSolver : public Solver< OperatorType >
solver_ = solver_factory_->generate( matrix_ ); solver_ = solver_factory_->generate( matrix_ );
if ( matrix_->get_size()[0] > 0 && matrix_->get_size()[1] > 0 ) if ( matrix_->get_size()[0] > 0 && matrix_->get_size()[1] > 0 )
{ {
//auto par_ilu = gko::factorization::Ilu< valueType, int32_t >::build().on( solver_exec_ )->generate( matrix_ ); auto ilu = gko::preconditioner::Ilu<>::build()
auto ilu = gko::preconditioner::Ilu<>::build().on( solver_exec_ )->generate( matrix_ ); .with_factorization_factory(
auto jac = gko::preconditioner::Jacobi<valueType, int32_t>::build().on(solver_exec_)->generate(matrix_); gko::share( gko::factorization::Ilu< valueType, int32_t >::build().on( solver_exec_ ) ) )
solver_->set_preconditioner( gko::share( jac ) ); .on( solver_exec_ )
->generate( matrix_ );
auto jac = gko::preconditioner::Jacobi< valueType, int32_t >::build().on( solver_exec_ )->generate( matrix_ );
solver_->set_preconditioner( gko::share( ilu ) );
} }
x.getStorage()->getTimingTree()->stop( "Ginkgo CG Solver Set-Up" ); x.getStorage()->getTimingTree()->stop( "Ginkgo CG Solver Set-Up" );
} }
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment