DFModifiedSolver.cu

Go to the documentation of this file.

#include <sys/stat.h>

#include "DFModifiedSolver.h"


template <typename memoryType>
DFModifiedSolver<memoryType>::DFModifiedSolver(parameterDB *pDB, domain *dInfo)
{
    NavierStokesSolver<memoryType>::paramDB = pDB;
    NavierStokesSolver<memoryType>::domInfo = dInfo;
} // DFModifiedSolver


template <typename memoryType>
void DFModifiedSolver<memoryType>::calculateExplicitLambdaTerms()
{
    // temp_1 = Q.pressure
    cusp::multiply(NavierStokesSolver<memoryType>::Q, DirectForcingSolver<memoryType>::pressure, NavierStokesSolver<memoryType>::temp1);
    // r^n = r^n - temp1
    cusp::blas::axpy(NavierStokesSolver<memoryType>::temp1, NavierStokesSolver<memoryType>::rn, -1.0);
} // calculateExplicitLambdaTerms


template <typename memoryType>
void DFModifiedSolver<memoryType>::projectionStep()
{
    // u^{n+1} = u^{*} - B^N.Q.\phi 
    NavierStokesSolver<memoryType>::projectionStep();

    // p^{n+1} = p^{n} + \phi
    NavierStokesSolver<memoryType>::logger.startTimer("projectionStep");
    cusp::blas::axpy(NavierStokesSolver<memoryType>::lambda, DirectForcingSolver<memoryType>::pressure, 1.0);
    NavierStokesSolver<memoryType>::logger.stopTimer("projectionStep");
} // projectionStep


// specialization of the class
template class DFModifiedSolver<device_memory>;