00001 #include "fwellboundarycondition.h" 00002 #include "numericmethods.h" 00003 00004 FWellBoundaryCondition::FWellBoundaryCondition(double value,Point3D &dx,VecHoles& holes) 00005 :m_value(value),m_holes(holes),m_DX(dx) 00006 { 00007 m_DX*=1.5; 00008 } 00009 00010 bool FWellBoundaryCondition::isInDomain(const Point3D &p,unsigned component) const 00011 { 00012 assert(component==0); 00013 for (unsigned i=0;i<m_holes.size();i++) 00014 { 00015 if (m_holes[i].isPointInViscinity(p,m_DX)) 00016 return true; 00017 } 00018 return false; 00019 } 00020 00021 double FWellBoundaryCondition::value (const Point3D &p, const unsigned int component) const 00022 { 00023 assert(component == 0); 00024 return m_value; 00025 } 00026 00027 00028 FWellBoundaryCondition::~FWellBoundaryCondition() 00029 { 00030 00031 }