properly create state.temperature in the FullyImplicitCompressiblePolymerSolver class

This commit is contained in:
Andreas Lauser
2016-02-22 15:14:09 +01:00
parent db86c367d5
commit 7c3020822a

View File

@@ -440,6 +440,10 @@ namespace {
int nextvar = 0;
state.pressure = vars[ nextvar++ ];
// Temperature. (this is always a constant so far)
const V T = Eigen::Map<const V>(& x.temperature()[0], nc, 1);
state.temperature = ADB::constant(T);
// Saturation.
const std::vector<int>& bpat = vars[0].blockPattern();
{