types.h

Go to the documentation of this file.

#pragma once

#include <fstream>
#include <iostream>
#include <cmath>
#include <cstdio>

#include <cusp/csr_matrix.h>
#include <cusp/array1d.h>
#include <cusp/print.h>
#include <cusp/transpose.h>
#include <cusp/blas/blas.h>
#include <cusp/krylov/cg.h>
#include <cusp/krylov/bicgstab.h>
#include <cusp/krylov/gmres.h>
#include <cusp/multiply.h>
#include <cusp/elementwise.h>
#include <thrust/device_ptr.h>

#include "logger.h"


enum bcType
{
    DIRICHLET,  
    NEUMANN,    
    CONVECTIVE, 
    PERIODIC,   
    SPECIAL
};


enum boundary
{
    XMINUS, 
    XPLUS,  
    YMINUS, 
    YPLUS   
};


enum timeScheme
{
    EULER_EXPLICIT,    
    EULER_IMPLICIT,    
    ADAMS_BASHFORTH_2, 
    RUNGE_KUTTA_3,     
    CRANK_NICOLSON     
};


enum ibmScheme
{
    NAVIER_STOKES,  
    SAIKI_BIRINGEN, 
    DIRECT_FORCING, 
    FADLUN_ET_AL,   
    TAIRA_COLONIUS, 
    DIFFUSION,      
    DF_MODIFIED,    
    FEA_MODIFIED,   
    DF_IMPROVED     
};


enum interpolationType
{
    CONSTANT,  
    LINEAR,    
    QUADRATIC  
};


enum preconditionerType
{
    NONE,                 
    DIAGONAL,             
    SMOOTHED_AGGREGATION, 
    AINV                  
};

typedef double real;


using cusp::device_memory;
using cusp::host_memory;

typedef cusp::coo_matrix<int, real, host_memory>   cooH;


typedef cusp::coo_matrix<int, real, device_memory> cooD;


typedef cusp::csr_matrix<int, real, host_memory>   csrH;


typedef cusp::csr_matrix<int, real, device_memory> csrD;


typedef cusp::array1d<real, host_memory>   vecH;


typedef cusp::array1d<real, device_memory> vecD;