Skip to content

Commit

Permalink
Merge Pull Request #9864 from brian-kelley/Trilinos/FixRitzinger
Browse files Browse the repository at this point in the history
Automatically Merged using Trilinos Pull Request AutoTester
PR Title: MueLu: fix ReitzingerPFactory for non-UVM
PR Author: brian-kelley
  • Loading branch information
trilinos-autotester authored Oct 27, 2021
2 parents 66f7efc + b155965 commit bc4e9f6
Showing 1 changed file with 5 additions and 5 deletions.
Original file line number Diff line number Diff line change
Expand Up @@ -388,13 +388,13 @@ namespace MueLu {
ArrayView<const LO> columns;
ArrayView<const SC> values;
Pe->getLocalRowView(i,columns,values);
// FIXME: This won't work on fancy nodes
ArrayView<SC> values_nc = Teuchos::av_const_cast<SC>(values);
Teuchos::Array<SC> newValues(values.size());
for (LO j=0; j<(LO)values.size(); j++)
if ((values_nc[j] == one || values_nc[j] == neg_one))
values_nc[j] = zero;
if ((values[j] == one || values[j] == neg_one))
newValues[j] = zero;
else
values_nc[j] /= two;
newValues[j] = values[j] / two;
Pe->replaceLocalValues(i,columns,newValues());
}//end for i < Ne
Pe->fillComplete(Pe->getDomainMap(),Pe->getRangeMap());
}
Expand Down

0 comments on commit bc4e9f6

Please sign in to comment.