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

Add MATLAB-compatible output of ADB::M values.

Browse Source 
This is mostly intended as an interim aid in development.  Not enabled
if 'NDEBUG' is defined.
This commit is contained in:
Bård Skaflestad 2013-05-03 23:45:13 +02:00
parent fc5a3f181e
commit 8b12d1506d
1 changed files with 42 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

42
sim_simple.cpp
View File

@ -112,7 +112,49 @@ struct HelperOps
}
};
#if !defined(NDEBUG)
#include <cstdio>
#endif // !defined(NDEBUG)
namespace {
#if !defined(NDEBUG)
void
printSparseMatrix(const Eigen::SparseMatrix<double>& A,
std::FILE* fp)
{
typedef Eigen::SparseMatrix<double>::Index Index;
const Index* const p = A.outerIndexPtr();
const Index* const i = A.innerIndexPtr();
const double* const x = A.valuePtr();
const Index cols = A.outerSize();
assert (A.innerSize() == cols);
for (Index j = 0; j < cols; j++) {
for (Index k = p[j]; k < p[j + 1]; k++) {
std::fprintf(fp, "%lu %lu %26.18e\n",
static_cast<unsigned long>(i[k] + 1),
static_cast<unsigned long>(j + 1), x[k]);
}
}
}
void
printSparseMatrix(const Eigen::SparseMatrix<double>& A ,
const char* const fn)
{
std::FILE* fp;
fp = std::fopen(fn, "w");
if (fp != 0) {
printSparseMatrix(A, fp);
}
std::fclose(fp);
}
#endif // !defined(NDEBUG)
template <typename Scalar>
class UpwindSelector {
public: