FadlunEtAlSolver.cu

Go to the documentation of this file.

#include <sys/stat.h>

#include "FadlunEtAlSolver.h"
#include "solvers/kernels/generateQT.h"


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


template <typename memoryType>
void FadlunEtAlSolver<memoryType>::generateQT()
{
    NavierStokesSolver<memoryType>::generateQT();
    G = NavierStokesSolver<memoryType>::Q;
    updateG();
} // generateQT


template <>
void FadlunEtAlSolver<device_memory>::updateG()
{
    const int blocksize = 256;
    
    int  nx = domInfo->nx,
         ny = domInfo->ny;
    
    int  GSize = 4*nx*ny-2*(nx+ny);
    
    int  *GRows = thrust::raw_pointer_cast(&(G.row_indices[0])),
         *GCols = thrust::raw_pointer_cast(&(G.column_indices[0]));
    int  *tags_r = thrust::raw_pointer_cast(&(tagsD[0]));

    real *GVals = thrust::raw_pointer_cast(&(G.values[0]));
    
    dim3 dimGrid( int((GSize-0.5)/blocksize) + 1, 1);
    dim3 dimBlock(blocksize, 1);
    
    kernels::updateQ <<<dimGrid, dimBlock>>> (GRows, GCols, GVals, GSize, tags_r);
} // updateG


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


template <typename memoryType>
void FadlunEtAlSolver<memoryType>::generateC()
{
    NavierStokesSolver<memoryType>::generateC();
} // generateC


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


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


template <typename memoryType>
void FEAModifiedSolver<memoryType>::generateQT()
{
    NavierStokesSolver<memoryType>::generateQT();
} // generateQT


template <typename memoryType>
void FEAModifiedSolver<memoryType>::generateC()
{
    NavierStokesSolver<memoryType>::generateC();
} // generateC


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