CO2INJECTION: fsquaremobt.cpp Source File

00001 
00002 #include "fsquaremobt.h"
00003 
00004  FSquareMobT::FSquareMobT(double v1,double v2, double Sr1,double Sr2) 
00005 {
00006   _v1 = v1;
00007   _v2 = v2;
00008   _Sr1 = Sr1;
00009   _Sr2 = Sr2;
00010   _srw = (1-_Sr1)*(1-_Sr1);
00011   _sro = 1-_Sr2; 
00012 }
00013 
00014 
00015  double FSquareMobT::operator ()(double x, unsigned cmp)const 
00016 {
00017   return (x-_Sr1)*(x-_Sr1)/(_v1*_srw) + (1- x/_sro)*(1- x/_sro)/_v2;
00018 }