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

fix dirichlet handling of initial guess

The initial guess of the original system is already used in the RHS creation. To stay consistent, i.e. using the solution as initial guess results in immediate convergence, the initial guess for the constraint system has to be set to zero.
parent 018f685d
...@@ -215,13 +215,7 @@ class ZeroRowsDirichletHandler : public DirichletHandlerBase ...@@ -215,13 +215,7 @@ class ZeroRowsDirichletHandler : public DirichletHandlerBase
orig_initial_guess_ = cur_initial_guess; orig_initial_guess_ = cur_initial_guess;
z_ = gko::clone( cur_initial_guess ); z_ = gko::clone( cur_initial_guess );
for ( auto idx : this->bcIndices_ ) z_->fill(gko::zero<valueType >());
{
if ( auto lidx = to_local_idx( idx, part, rank ); lidx )
{
z_->get_local()->at( *lidx ) = gko::zero< valueType >();
}
}
return z_.get(); return z_.get();
} }
......
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