fixedvaluecondition.cpp

Go to the documentation of this file.

#include "fixedvaluecondition.h"


void FixedValueCondition::applyConditions(VecDouble &v) 
{
  
  for (unsigned i=0;i<data.size();i++)
  {
    v(data[i].index)=data[i].values(0);
  }
}

void FixedValueCondition::applyConditions(ArrayOfVecDouble &v) 
{
  for (unsigned i=0;i<data.size();i++)
  {
    assert(v.vecs_size() == data[i].values.size());
    for (unsigned j=0;j<v.vecs_size();j++)
    {
      v(data[i].index,j)=data[i].values(j);
    }
  }
  
}




void FixedValueCondition::addFixedCondition(OrthoMesh &mesh, Point3D &p,double value) 
{
  OrthoMesh::Cell_It cell = mesh.getCellAt(p);
  Node node;
  node.index = cell->index();
  node.values.reinit(1);
  node.values(0)= value;
  data.push_back(node);
  
}


void FixedValueCondition::addFixedCondition(unsigned index, double value) 
{
  Node node;
  node.index = index;
  node.values.reinit(1);
  node.values(0)= value;
  data.push_back(node);
}


 FixedValueCondition::FixedValueCondition() 
{

}


void FixedValueCondition::addTransportFixedCondition(OrthoMesh &mesh,const VecWellInfo &wells) 
{
  OrthoMesh::Cell_It cell = mesh.begin_cell();
  OrthoMesh::Cell_It endc = mesh.end_cell();
  Node node;
  for(;cell!=endc;cell++)
  {
    Point3D p;
    cell->barycenter(p);
    for (unsigned i=0;i<wells.size();i++)
    {
      if (wells[i].isPointInWell(p))
      {
          node.index=cell->index();
          node.values=wells[i].getTransportBC();
      }
    }
  }
}  

---