lib/field/vfield.cc
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
 *
 *          The diffusion term (grad-squared) is caulculated here.
 *          The second derivatives of each component field are calculated along x, y and z.
 *          These terms are added to the corresponding components of the given plain
 *          vector field (plainvf), which is usually the RHS of the PDE being solved.
 *
 * \param   H is a reference to the plainvf into which the output is written
 ********************************************************************************************************************************************
 */
\brief   Function to compute the diffusion term
The diffusion term (grad-squared) is caulculated here.
The second derivatives of each component field are calculated along x, y and z.
These terms are added to the corresponding components of the given plain
vector field (plainvf), which is usually the RHS of the PDE being solved.

\param   H is a reference to the plainvf into which the output is written
\ brief Function to compute the diffusion term The diffusion term ( grad - squared ) is caulculated here . The second derivatives of each component field are calculated along x y and z . These terms are added to the corresponding components of the given plain vector field ( plainvf ) which is usually the RHS of the PDE being solved . \ param H is a reference to the plainvf into which the output is 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
          The diffusion term (grad-squared) is caulculated here.
\ brief Function to compute the diffusion term The diffusion term ( grad - squared ) is caulculated here .