CO2INJECTION: dfpcfracture.cpp Source File

00001 #include "dfpcfracture.h"
00002 //#include "numericmethods.h"
00003 
00004 
00005 DFPcFracture::DFPcFracture(double gamma,double theta) 
00006 {
00007   _gamma = gamma;
00008   _theta = theta;
00009 }
00010 
00011  double DFPcFracture::operator()(double x,unsigned cmp) const 
00012 {
00013   double sr=1.0e-12;
00014   assert(cmp==0);
00015   if(x > 0 && x <= 1.0)
00016     return _gamma - _gamma/pow(x,2) - _theta;
00017   else if (x<=0)
00018     return -1.0e+10; //_gamma - _gamma/pow(sr,2) - _theta;
00019   else
00020     return _theta;
00021 }