OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-01-16 01:41:56 -06:00
Code Issues Packages Projects Releases Wiki Activity

Added check for single-cell columns.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen 2012-05-14 11:04:00 +02:00
parent 67a6674495
commit 088c61ce00
1 changed files with 11 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
opm/polymer/GravityColumnSolverPolymer_impl.hpp
View File

@ -39,6 +39,7 @@
#include <opm/polymer/GravityColumnSolverPolymer.hpp> #include <opm/polymer/GravityColumnSolverPolymer.hpp>
#include <opm/core/linalg/blas_lapack.h> #include <opm/core/linalg/blas_lapack.h>
#include <opm/core/utility/ErrorMacros.hpp> #include <opm/core/utility/ErrorMacros.hpp>
#include <iterator>
namespace Opm namespace Opm
{ {
@ -178,6 +179,13 @@ namespace Opm
// This is written only to work with SinglePointUpwindTwoPhase, // This is written only to work with SinglePointUpwindTwoPhase,
// not with arbitrary problem models. // not with arbitrary problem models.
int col_size = column_cells.size(); int col_size = column_cells.size();
if (col_size == 1) {
sol_vec[2*column_cells[0] + 0] = 0.0;
sol_vec[2*column_cells[0] + 1] = 0.0;
return;
}
StateWithZeroFlux state(s, c, cmax); // This holds s by reference. StateWithZeroFlux state(s, c, cmax); // This holds s by reference.
// Assemble. // Assemble.
@ -249,6 +257,9 @@ namespace Opm
// Solution will be written to rhs. // Solution will be written to rhs.
dgbsv_(&N, &kl, &ku, &num_rhs, &hm[0], &nrow, &ipiv[0], &rhs[0], &N, &info); dgbsv_(&N, &kl, &ku, &num_rhs, &hm[0], &nrow, &ipiv[0], &rhs[0], &N, &info);
if (info != 0) { if (info != 0) {
std::cerr << "Failed column cells: ";
std::copy(column_cells.begin(), column_cells.end(), std::ostream_iterator<int>(std::cerr, " "));
std::cerr << "\n";
THROW("Lapack reported error in dgtsv: " << info); THROW("Lapack reported error in dgtsv: " << info);
} }
for (int ci = 0; ci < col_size; ++ci) { for (int ci = 0; ci < col_size; ++ci) {