only compute field fluid in place.

This commit is contained in:
Liu Ming 2016-07-29 11:07:21 +08:00 committed by Atgeirr Flø Rasmussen
parent 056708574c
commit 90086f690e

View File

@ -2347,18 +2347,16 @@ namespace detail {
fip[4] = state.rv.value() * fip[pg];
}
const int dims = *std::max_element(fipnum.begin(), fipnum.end()) + 1;
const int dims = 1;
std::vector<V> values(dims, V::Zero(5));
for (int d = 0; d < dims; ++d) {
for (int c = 0; c < nc; ++c) {
if (fipnum[c] == d) {
for (int i = 0; i < 5; ++i) {
values[d][i] += fip[i][c];
}
}
}
}
return values;