mappedfevaluesbase.cpp

Go to the documentation of this file.

#include "mappedfevaluesbase.h"


void MappedFEValuesBase::setPoints(const ArrayOfVecDouble &v,const VecDouble &weights) 
{
  _pts=v;
  _weights=weights;
  unsigned n_pts=v.size();
  _J.reinit(n_pts);
  _JxW.reinit(n_pts);
  //clear the data
  _Values.reinit(n_pts,m_n_dofs);
  _Grad.resize(m_n_dofs*n_pts);
  for (unsigned i=0;i<_Grad.size();i++)
  {
    _Grad[i].reinit(3);
  }

  _Mapped_Points.resize(n_pts);
  for (unsigned i=0;i<_Mapped_Points.size();i++)
  {
    _Mapped_Points[i].reinit(3);
  }


}



void MappedFEValuesBase::compute_points(MappingOrthoMesh &map) 
{
  VecDoubleRef pt(3);

  for(unsigned i=0;i<_pts.size();i++)
  {
    //get pt i
    _pts.getVecValues(i,&pt); //Get the point

    //get GradT
    map.T(pt,_Mapped_Points[i]);
  }
  
}


void MappedFEValuesBase::print(std::ostream &out) 
{
  Matrix &M = _pts;

  out << "Points: \n";
  M.print_formatted(out,5,true,7,"0");
  out << "Jacobians: ";
  _J.print(out,5,true,true);
  out << "\nJacobians x W: ";
  _JxW.print(out,5,true,true);
  out<<std::endl;
  out << "\nPoints x Function Values:\n";
  _Values.print_formatted(out,5,true,7,"0");
  
  
    
}

---