eqvart.cpp

Go to the documentation of this file.

#include "eqvart.h"

double EqVarT::getSolValue(Vector<double> &sol) 
{
  unsigned i;
  unsigned n_dofs_per_cell = m_fev->get_fe().dofs_per_cell;
  double result=0.0;
  for (i=0;i<n_dofs_per_cell;i++)
  {
    result+=sol(m_dofIndices[i])*m_fev->shape_value(i,m_qPoint);
  }
  return result;
}

double EqVarT::getCoefSolValue(unsigned index,VecDouble &sol) 
{
  return sol(m_dofIndices[index]);
}
 
double EqVarT::getSolGradComponent(Vector<double> &sol,int component) 
{
  unsigned i;
  unsigned n_dofs_per_cell = m_fev->get_fe().dofs_per_cell;
  double result=0.0;
  for (i=0;i<n_dofs_per_cell;i++)
  {
    result+=sol(m_dofIndices[i])*m_fev->shape_grad(i,m_qPoint)[component];
  }
  return result;
}


Tensor<1,DIM> EqVarT::getSolGrad(VecDouble &sol) 
{
  unsigned i;
  unsigned n_dofs_per_cell = m_fev->get_fe().dofs_per_cell;
  Tensor<1,DIM> result;
  result.clear();
  for (i=0;i<n_dofs_per_cell;i++)
  {
    result+=sol(m_dofIndices[i])*m_fev->shape_grad(i,m_qPoint);
  }
  return result;
}

/***************************************************
Objetivo:  Obter o valor de um vertor  baseado no indice
           da celula.

Parametros: sol = vetor contendo o valor de todas as celulas
********************************************************/
double EqVarT::getValueByCellIndex(const VecDouble &sol) 
{
  return sol(getCell()->index());
}  




/***************************************************
Objetivo: Configurar o membro m_face_values;  

Parametros: 

Retorno: 
********************************************************/
void EqVarT::setFaceValues(FEFaceValues<DIM> *face_values) 
{
  m_face_values = face_values;
}

void EqVarT::setFeValue(FEValues<DIM>  *fe) 
{
 m_fev = fe;
}

/***************************************************
Objetivo:  Dada o numero da funcao teste, contruir
           a condicao de Newmann relativo a linha
           da matriz local que a funcao pertence.

Parametros: i = Numero da equacao.

Retorno: 
********************************************************/
double EqVarT::localNeumannBoundary(unsigned i) 
{
  i=0;
  assert(0);   //Se chamado nessa classe, retorne erro imediatamente.
  return NAN;
}


/***************************************************
Objetivo:  

Parametros: 

Retorno: 
********************************************************/

double EqVarT::localEqVar(unsigned i,unsigned j) 
{
  printf("EqvarT::localEqVar = This function shoud never be called. The childrens must override this functions\nAbortting...\n");
  abort();
  return NAN;
}

double EqVarT::localRhs(unsigned i) 
{
  printf("EqvarT::localRhs = This function shoud never be called. The childrens must override this functions\nAbortting...\n");
  abort();
  return NAN;

}

/***************************************************
Objetivo:  

Parametros: 

Retorno: 
********************************************************/

unsigned EqVarT::component(unsigned i) 
{
  //Verifica se component nao esta sendo chamado dentro de addNeummanCondition.
  assert(validFeValues());
  return m_fev->get_fe().system_to_component_index(i).first;
}

/***************************************************
Objetivo:  

Parametros: 

Retorno: 
********************************************************/

double EqVarT::shape_div(unsigned i) 
{
  assert(validFeValues());
  assert(!m_fev->get_fe().is_primitive());
  return  m_fev->shape_grad_component(i,getQPoint(),0)[0] + m_fev->shape_grad_component(i,getQPoint(),1)[1];
}

---