lib/field/vfield.cc
/**
 ********************************************************************************************************************************************
 * \brief   Function to compute the diffusion term
 *
 *          It is assumed that the vector field is specified at face-centers, as required by the \ref sfield#computeNLin
 *          "computeNLin" function of sfield.
 *
 * \param   H is a pointer to a vector field (vfield) to which the output of the function is to be written
 ********************************************************************************************************************************************
 */
\brief   Function to compute the diffusion term
It is assumed that the vector field is specified at face-centers, as required by the \ref sfield#computeNLin
"computeNLin" function of sfield.

\param   H is a pointer to a vector field (vfield) to which the output of the function is to be written
\ brief Function to compute the diffusion term It is assumed that the vector field is specified at face - centers as required by the \ ref sfield#computeNLin " computeNLin " function of sfield . \ param H is a pointer to a vector field ( vfield ) to which the output of the function is to be written
void vfield::computeDiff(plainvf &H) {
    derivTempX = 0.0;
    derVx.calcDerivative2xx(derivTempX);
    H.Vx(Vx.fCore) += derivTempX(Vx.fCore);
#ifndef PLANAR
    derivTempX = 0.0;
    derVx.calcDerivative2yy(derivTempX);
    H.Vx(Vx.fCore) += derivTempX(Vx.fCore);
#endif
    derivTempX = 0.0;
    derVx.calcDerivative2zz(derivTempX);
    H.Vx(Vx.fCore) += derivTempX(Vx.fCore);
#ifndef PLANAR
    derivTempY = 0.0;
    derVy.calcDerivative2xx(derivTempY);
    H.Vy(Vy.fCore) += derivTempY(Vy.fCore);
    derivTempY = 0.0;
    derVy.calcDerivative2yy(derivTempY);
    H.Vy(Vy.fCore) += derivTempY(Vy.fCore);
    derivTempY = 0.0;
    derVy.calcDerivative2zz(derivTempY);
    H.Vy(Vy.fCore) += derivTempY(Vy.fCore);
#endif
    derivTempZ = 0.0;
    derVz.calcDerivative2xx(derivTempZ);
    H.Vz(Vz.fCore) += derivTempZ(Vz.fCore);
#ifndef PLANAR
    derivTempZ = 0.0;
    derVz.calcDerivative2yy(derivTempZ);
    H.Vz(Vz.fCore) += derivTempZ(Vz.fCore);
#endif
    derivTempZ = 0.0;
    derVz.calcDerivative2zz(derivTempZ);
    H.Vz(Vz.fCore) += derivTempZ(Vz.fCore);
}
void vfield :: computeDiff ( plainvf & H ) { derivTempX = 0.0 ; derVx . calcDerivative2xx ( derivTempX ) ; H . Vx ( Vx . fCore ) += derivTempX ( Vx . fCore ) ; #ifndef PLANAR derivTempX = 0.0 ; derVx . calcDerivative2yy ( derivTempX ) ; H . Vx ( Vx . fCore ) += derivTempX ( Vx . fCore ) ; #endif derivTempX = 0.0 ; derVx . calcDerivative2zz ( derivTempX ) ; H . Vx ( Vx . fCore ) += derivTempX ( Vx . fCore ) ; #ifndef PLANAR derivTempY = 0.0 ; derVy . calcDerivative2xx ( derivTempY ) ; H . Vy ( Vy . fCore ) += derivTempY ( Vy . fCore ) ; derivTempY = 0.0 ; derVy . calcDerivative2yy ( derivTempY ) ; H . Vy ( Vy . fCore ) += derivTempY ( Vy . fCore ) ; derivTempY = 0.0 ; derVy . calcDerivative2zz ( derivTempY ) ; H . Vy ( Vy . fCore ) += derivTempY ( Vy . fCore ) ; #endif derivTempZ = 0.0 ; derVz . calcDerivative2xx ( derivTempZ ) ; H . Vz ( Vz . fCore ) += derivTempZ ( Vz . fCore ) ; #ifndef PLANAR derivTempZ = 0.0 ; derVz . calcDerivative2yy ( derivTempZ ) ; H . Vz ( Vz . fCore ) += derivTempZ ( Vz . fCore ) ; #endif derivTempZ = 0.0 ; derVz . calcDerivative2zz ( derivTempZ ) ; H . Vz ( Vz . fCore ) += derivTempZ ( Vz . fCore ) ; }
\brief   Function to compute the diffusion term
          It is assumed that the vector field is specified at face-centers, as required by the \ref sfield#computeNLin
          "computeNLin" function of sfield.
\ brief Function to compute the diffusion term It is assumed that the vector field is specified at face - centers as required by the \ ref sfield#computeNLin " computeNLin " function of sfield .