Bugfix: size of vector for A should be np*np.

This commit is contained in:
Atgeirr Flø Rasmussen 2012-10-02 14:35:28 +02:00
parent 017663cc5b
commit 4f276a8870

View File

@ -297,7 +297,7 @@ namespace Opm
if (np != 2) {
THROW("computeTransportSource() requires a 2 phase case.");
}
std::vector<double> A(np, np);
std::vector<double> A(np*np);
for (int w = 0; w < nw; ++w) {
const double* comp_frac = wells->comp_frac + np*w;
for (int perf = wells->well_connpos[w]; perf < wells->well_connpos[w + 1]; ++perf) {