OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2026-07-30 04:57:56 -05:00
Code Issues Packages Projects Releases Wiki Activity

Use constant intermediate quantities.

Browse Source
This commit is contained in:
Bård Skaflestad
2013-05-12 22:48:15 +02:00
parent b0cf134479
commit 61f184cdc9
1 changed files with 11 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files
ImpesTPFAAD.hpp
+11 -7
View File
@@ -106,10 +106,12 @@ namespace Opm {
assert (0 <= phase);
assert (phase < np_ );
typename ADB::V A = A_ .block(0, phase * (np_ + 1), nc_, 1);
typename ADB::V dA = dA_.block(0, phase * (np_ + 1), nc_, 1);
typedef typename ADB::V V;
std::vector<typename ADB::M> jac(1, spdiag(dA));
const V A = A_ .block(0, phase * (np_ + 1), nc_, 1);
const V dA = dA_.block(0, phase * (np_ + 1), nc_, 1);
const std::vector<typename ADB::M> jac(1, spdiag(dA));
return one_ / ADB::function(A, jac);
}
@@ -117,7 +119,7 @@ namespace Opm {
typename ADB::V
phaseRelPerm(const int phase) const
{
typename ADB::V kr = kr_.block(0, phase, nc_, 1);
const typename ADB::V kr = kr_.block(0, phase, nc_, 1);
return kr;
}
@@ -128,10 +130,12 @@ namespace Opm {
assert (0 <= phase);
assert (phase < np_ );
typename ADB::V mu = mu_ .block(0, phase, nc_, 1);
typename ADB::V dmu = dmu_.block(0, phase, nc_, 1);
typedef typename ADB::V V;
std::vector<typename ADB::M> jac(1, spdiag(dmu));
const V mu = mu_ .block(0, phase, nc_, 1);
const V dmu = dmu_.block(0, phase, nc_, 1);
const std::vector<typename ADB::M> jac(1, spdiag(dmu));
return ADB::function(mu, jac);
}