Fix when the matrix with no ghosts will be used.

This commit is contained in:
Markus Blatt 2020-05-12 11:34:38 +02:00
parent e8b79b2b63
commit 1ff1f24fbb

View File

@ -503,7 +503,7 @@ protected:
if( isParallel() ) if( isParallel() )
{ {
if ( ownersFirst_ && (!parameters_.linear_solver_use_amg_ || !parameters_.use_cpr_) ) { if ( ownersFirst_ && !parameters_.linear_solver_use_amg_ && !useFlexible_) {
typedef WellModelGhostLastMatrixAdapter< Matrix, Vector, Vector, WellModel, true > Operator; typedef WellModelGhostLastMatrixAdapter< Matrix, Vector, Vector, WellModel, true > Operator;
Operator opA(*matrix_, *matrix_, wellModel, interiorCellNum_); Operator opA(*matrix_, *matrix_, wellModel, interiorCellNum_);