vecincbc.cpp

Go to the documentation of this file.

#include "vecincbc.h"
#include "sfunctions.h"


void VecIncBC::reinit(unsigned nBC,unsigned vectors_size) 
{
  _vectors_size=vectors_size;
  m_v.clear();
  m_v.reserve(nBC*vectors_size);
  this->reserve(nBC);
}




void VecIncBC::add_bc(Index faceIndex,Index nearCell,const VecDouble &bc) 
{
  for (Index i=0;i<_vectors_size;i++)
    m_v.push_back(bc(i));
  this->push_back(faceIndex,VecIncBCData(&(m_v.back())-_vectors_size +1,nearCell));
}



void VecIncBC::set_bc_from_function(OrthoMesh &mesh,Function3D &f) 
{
  OrthoMesh::Face_It  face = mesh.begin_face();
  OrthoMesh::Face_It  end = mesh.end_face();
  Point3D P;
  VecDouble v(f.getImageDim());
  unsigned nBC=0;
  for (;face!=end;face++)
  {
    if (face->at_boundary())
    {
       face->barycenter(P);
        if (!f.isInDomain(P))
          continue;
        else
          nBC++;
    }
  }
  printf("Reinit BC to %u\n",nBC);
  reinit(nBC,f.getImageDim());

  for (face=mesh.begin_face();face!=end;face++)
  {
    if (face->at_boundary())
    {
       face->barycenter(P);
        if (!f.isInDomain(P))
          continue;
        else
        {
          unsigned dim = f.getImageDim();
          for (unsigned j=0;j<dim;j++)
            v(j)=f(P,j);
          unsigned c1,c2;
          face->getValidCellIndices(c1,c2);
          add_bc(face->index(),c1,v);
        }
    }
  }
}


void VecIncBC::rewind() 
{
  assert(!m_v.empty());
  it=begin();
}



bool VecIncBC::get_ref_bc(Index index, VecDoubleRef &vSw) 
{
  assert(vSw.size() == _vectors_size);
  VecIncBCData *data = it.inc_search(index);
  if (data)
  {
    vSw.setData(data->p);
    return true;
  }
  else
    return false;
  
}

bool VecIncBC::copy_bc(Index index, VecDouble &vSw) 
{
  assert(vSw.size() == _vectors_size);
  assert(vSw.is_owner());
  VecIncBCData *data = it.inc_search(index);
  if (data)
  {
    const double *pd=data->p;
    for (unsigned i=0;i<_vectors_size;i++)
      vSw(i)=pd[i];
    return true;
  }
  else
    return false;
  

}



bool VecIncBC::has_bc(Index index) 
{
    return it.inc_search(index)!=NULL;
}


std::ostream& operator<< (std::ostream& out,VecIncBC &vbc) 
{
  //VecIncBC::iterator it = vbc.begin();
  std::vector<VecIncBC::Node >::const_iterator it = vbc.m_list.begin();
  std::vector<VecIncBC::Node >::const_iterator itlast = vbc.m_list.end();
  itlast--;
  while(it != itlast)
  {
    out << it->index <<") " << "Near Cell" << it->data.nearCell << " Add: " <<  it->data.p << ": ";
    for (unsigned i=0;i<vbc._vectors_size;i++)
    {
      out << it->data.p[i] << " ";
    }
    out << std::endl;
    it++;
  }
  return out;
}


/*
  Get the pointer to the array of doubles
  containing the boundary condition for index
  bc. If such boundary condition does not exist  
  a NULL pointer is returned
 */
double* VecIncBC::get_bc_ptr(Index index) 
{
  VecIncBCData *data = it.inc_search(index);
  return (data) ? data->p : NULL;
}

---