#include <fpcinvfracture.h>
Public Member Functions | |
FPcInvFracture (double gamma, double theta) | |
virtual double | operator() (double x, unsigned cmp) const |
virtual | ~FPcInvFracture () |
Protected Attributes | |
double | _gamma |
double | _theta |
Implement the Pc inverse function as described in the paper "A parallelizable method for two-phase flows in naturally-fractured reservoirs"
Definition at line 11 of file fpcinvfracture.h.
FPcInvFracture::FPcInvFracture | ( | double | gamma, | |
double | theta | |||
) |
Definition at line 5 of file fpcinvfracture.cpp.
virtual FPcInvFracture::~FPcInvFracture | ( | ) | [inline, virtual] |
Definition at line 22 of file fpcinvfracture.h.
double FPcInvFracture::operator() | ( | double | x, | |
unsigned | cmp | |||
) | const [virtual] |
Implements Function1D.
Definition at line 11 of file fpcinvfracture.cpp.
00012 { 00013 double s; 00014 assert(cmp==0); 00015 if( x > 0.0) 00016 { 00017 s = x-_theta+2.0*_gamma; 00018 s = s-pow((x - _theta)*(x - _theta)+4.0*x*_gamma,0.5); 00019 s = 0.5*s/(_gamma-_theta); 00020 return s; 00021 } 00022 else if( x <= 0.0) 00023 { 00024 return 0.99999; 00025 } 00026 else 00027 return 1.0e-12; 00028 }
double FPcInvFracture::_gamma [protected] |
Definition at line 15 of file fpcinvfracture.h.
double FPcInvFracture::_theta [protected] |
Definition at line 15 of file fpcinvfracture.h.