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

Undo changes to NewtonIterationsBlackoilInterleaved.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen 2015-10-08 15:11:50 +02:00
parent 9e3d360b5b
commit 9be05f4206
1 changed files with 16 additions and 18 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

34
opm/autodiff/NewtonIterationBlackoilInterleaved.cpp
View File

@ -207,7 +207,7 @@ namespace Opm
*/ */
template<typename Scalar> struct PointOneOp { template<typename Scalar> struct PointOneOp {
EIGEN_EMPTY_STRUCT_CTOR(PointOneOp) EIGEN_EMPTY_STRUCT_CTOR(PointOneOp)
Scalar operator()(const Scalar&, const Scalar&) const { return 0.1; } Scalar operator()(const Scalar& a, const Scalar& b) const { return 0.1; }
}; };
} }
@ -234,25 +234,23 @@ namespace Opm
assert(size == row_major.cols()); assert(size == row_major.cols());
// Create ISTL matrix with interleaved rows and columns (block structured). // Create ISTL matrix with interleaved rows and columns (block structured).
{ assert(np == 3);
assert(np == 3); istlA.setSize(row_major.rows(), row_major.cols(), row_major.nonZeros());
istlA.setSize(row_major.rows(), row_major.cols(), row_major.nonZeros()); istlA.setBuildMode(Mat::row_wise);
istlA.setBuildMode(Mat::row_wise); const int* ia = row_major.outerIndexPtr();
const int* ia = row_major.outerIndexPtr(); const int* ja = row_major.innerIndexPtr();
const int* ja = row_major.innerIndexPtr(); for (Mat::CreateIterator row = istlA.createbegin(); row != istlA.createend(); ++row) {
for (Mat::CreateIterator row = istlA.createbegin(); row != istlA.createend(); ++row) { const int ri = row.index();
const int ri = row.index(); for (int i = ia[ri]; i < ia[ri + 1]; ++i) {
for (int i = ia[ri]; i < ia[ri + 1]; ++i) { row.insert(ja[i]);
row.insert(ja[i]);
}
} }
}
// Set all blocks to zero. // Set all blocks to zero.
for (int row = 0; row < size; ++row) { for (int row = 0; row < size; ++row) {
for (int col_ix = ia[row]; col_ix < ia[row + 1]; ++col_ix) { for (int col_ix = ia[row]; col_ix < ia[row + 1]; ++col_ix) {
const int col = ja[col_ix]; const int col = ja[col_ix];
istlA[row][col] = 0.0; istlA[row][col] = 0.0;
}
} }
} }