OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-12-27 01:31:00 -06:00
Code Issues Packages Projects Releases Wiki Activity

Work in progress.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen 2014-11-27 14:33:31 +01:00
parent 5d81306007
commit 654e06de00
1 changed files with 5 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
opm/core/tof/AnisotropicEikonal.cpp
View File

@ -250,9 +250,11 @@ namespace Opm
const double* g;
double operator()(const double theta) const
{
double a[2] = { x[0] - (1-theta)*x1[0] + theta*x2[0], x[1] - (1-theta)*x1[1] + theta*x2[1] };
double b[2] = { x1[0] - x2[0], x1[1] - x2[1] };
return u2 - u1 + 2*(a[0]*b[0]*g[0] + a[0]*b[1]*g[1] + a[1]*b[0]*g[2] + a[1]*b[1]*g[3]);
const double a[2] = { x[0] - (1-theta)*x1[0] + theta*x2[0], x[1] - (1-theta)*x1[1] + theta*x2[1] };
const double b[2] = { x1[0] - x2[0], x1[1] - x2[1] };
const double val = u2 - u1 + 2*(a[0]*b[0]*g[0] + a[0]*b[1]*g[1] + a[1]*b[0]*g[2] + a[1]*b[1]*g[3]);
std::cout << theta << " " << val << std::endl;
return val;
}
};