OilfieldToolsNet/opm-core
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Honour exact interface of LAPACK routine DGTSV.

Browse Source 
The "MAT_SIZE_T" is not necessarily equivalent to 'int'.
This commit is contained in:
Bård Skaflestad
2012-10-05 00:53:47 +02:00
parent de3b7f002f
commit a6dc97a706
1 changed files with 3 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files
opm/core/transport/GravityColumnSolver_impl.hpp
+3 -3
View File
@@ -175,10 +175,10 @@ namespace Opm
}
// model_.sourceTerms(); // Not needed
// Solve.
const int num_rhs = 1;
int info = 0;
const MAT_SIZE_T num_rhs = 1, colSize = col_size;
MAT_SIZE_T info = 0;
// Solution will be written to rhs.
dgtsv_(&col_size, &num_rhs, DL, D, DU, &rhs[0], &col_size, &info);
dgtsv_(&colSize, &num_rhs, DL, D, DU, &rhs[0], &colSize, &info);
if (info != 0) {
THROW("Lapack reported error in dgtsv: " << info);
}