createSolver.cu

Go to the documentation of this file.

#include "createSolver.h"


template <typename memoryType>
NavierStokesSolver<memoryType>* createSolver(parameterDB &paramDB, domain &domInfo)
{
    ibmScheme ibm = paramDB["simulation"]["ibmScheme"].get<ibmScheme>();
    NavierStokesSolver<memoryType> *solver = 0;
    switch(ibm)
    {
        case SAIKI_BIRINGEN:
            break;
        case TAIRA_COLONIUS:
            solver = new TairaColoniusSolver<memoryType>(&paramDB, &domInfo);
            break;
        case NAVIER_STOKES:
            solver = new NavierStokesSolver<memoryType>(&paramDB, &domInfo);
            break;
        case DIRECT_FORCING:
            solver = new DirectForcingSolver<memoryType>(&paramDB, &domInfo);
            break;
        case FADLUN_ET_AL:
            solver = new FadlunEtAlSolver<memoryType>(&paramDB, &domInfo);
            break;
        case DIFFUSION:
            solver = new DiffusionSolver<memoryType>(&paramDB, &domInfo);
            break;
        case DF_MODIFIED:
            solver = new DFModifiedSolver<memoryType>(&paramDB, &domInfo);
            break;
        case FEA_MODIFIED:
            solver = new FEAModifiedSolver<memoryType>(&paramDB, &domInfo);
            break;
        case DF_IMPROVED:
            solver = new DFImprovedSolver<memoryType>(&paramDB, &domInfo);
            break;
    }
    std::cout << "\nImmersed Boundary Method" << std::endl;
    std::cout <<   "------------------------" << std::endl;
    std::cout << solver->name() << '\n' << std::endl;
    return solver;
} // createSolver


// specializations of the template function
template NavierStokesSolver<device_memory>* createSolver(parameterDB &paramDB, domain &domInfo);