lib/field/derivative.cc
void derivative::calcDerivative2xx(blitz::Array<real, 3> outputMat) {
    tempMat(blitz::Range(0, F.F.ubound(0) - 1), fullRange, fullRange) = central12n(F.F, 0);
    tempMat *= invDelx;

    outputMat(blitz::Range(0, F.F.ubound(0) - 1), fullRange, fullRange) = central22n(F.F, 0);
    outputMat *= invDelx*invDelx;

    if (inputParams.iScheme == 1) {
        outputMat = xxMetric(i)*tempMat(i, j, k) + 0.5*x2Metric(i)*outputMat(i, j, k);
    } else {
        outputMat = xxMetric(i)*tempMat(i, j, k) + x2Metric(i)*outputMat(i, j, k);
    }
    
}
/**
 ********************************************************************************************************************************************
 * \brief   Function to compute the second derivative of the field with respect to x [second order accurate central difference]
 *
 *          This function must be called using an output array whose shape and size should be same as that of the field
 *          
 ********************************************************************************************************************************************
 */
\brief   Function to compute the second derivative of the field with respect to x [second order accurate central difference]
This function must be called using an output array whose shape and size should be same as that of the field
\ brief Function to compute the second derivative of the field with respect to x [ second order accurate central difference ] This function must be called using an output array whose shape and size should be same as that of the field
void derivative::calcDerivative2xx(blitz::Array<real, 3> outputMat) {
    tempMat(blitz::Range(0, F.F.ubound(0) - 1), fullRange, fullRange) = central12n(F.F, 0);
    tempMat *= invDelx;
    outputMat(blitz::Range(0, F.F.ubound(0) - 1), fullRange, fullRange) = central22n(F.F, 0);
    outputMat *= invDelx*invDelx;
    if (inputParams.iScheme == 1) {
        outputMat = xxMetric(i)*tempMat(i, j, k) + 0.5*x2Metric(i)*outputMat(i, j, k);
    } else {
        outputMat = xxMetric(i)*tempMat(i, j, k) + x2Metric(i)*outputMat(i, j, k);
    }
}
void derivative :: calcDerivative2xx ( blitz :: Array < real , 3 > outputMat ) { tempMat ( blitz :: Range ( 0 , F . F . ubound ( 0 ) - 1 ) , fullRange , fullRange ) = central12n ( F . F , 0 ) ; tempMat *= invDelx ; outputMat ( blitz :: Range ( 0 , F . F . ubound ( 0 ) - 1 ) , fullRange , fullRange ) = central22n ( F . F , 0 ) ; outputMat *= invDelx * invDelx ; if ( inputParams . iScheme == 1 ) { outputMat = xxMetric ( i ) * tempMat ( i , j , k ) + 0.5 * x2Metric ( i ) * outputMat ( i , j , k ) ; } else { outputMat = xxMetric ( i ) * tempMat ( i , j , k ) + x2Metric ( i ) * outputMat ( i , j , k ) ; } }
\brief   Function to compute the second derivative of the field with respect to x [second order accurate central difference]
          This function must be called using an output array whose shape and size should be same as that of the field
\ brief Function to compute the second derivative of the field with respect to x [ second order accurate central difference ] This function must be called using an output array whose shape and size should be same as that of the field