BCInitialize.c

Go to the documentation of this file.

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <boundaryconditions.h>

int BCInitialize( void *b, 
                  int flag_gpu   )
{
  DomainBoundary *boundary = (DomainBoundary*) b;

#if defined(HAVE_CUDA)
  if (flag_gpu) {

    if      (!strcmp(boundary->bctype,_PERIODIC_)) {
      boundary->BCFunctionU = gpuBCPeriodicU;
    } else if (!strcmp(boundary->bctype,_SLIP_WALL_)) {
      boundary->BCFunctionU = gpuBCSlipWallU;
    } else {
      fprintf(stderr,"[GPU] Error in BCInitialize(): \"%s\" is not a supported boundary condition.\n",
              boundary->bctype);
      return 1;
    }

  } else {
#endif

    if      (!strcmp(boundary->bctype,_PERIODIC_                    )) boundary->BCFunctionU = BCPeriodicU;
    else if (!strcmp(boundary->bctype,_EXTRAPOLATE_                 )) boundary->BCFunctionU = BCExtrapolateU;
    else if (!strcmp(boundary->bctype,_DIRICHLET_                   )) boundary->BCFunctionU = BCDirichletU;  
    else if (!strcmp(boundary->bctype,_REFLECT_                     )) boundary->BCFunctionU = BCReflectU;    
    else if (!strcmp(boundary->bctype,_SPONGE_                      )) boundary->BCFunctionU = BCSpongeUDummy;    
    else if (!strcmp(boundary->bctype,_NOSLIP_WALL_                 )) boundary->BCFunctionU = BCNoslipWallU;    
    else if (!strcmp(boundary->bctype,_SLIP_WALL_                   )) boundary->BCFunctionU = BCSlipWallU;    
    else if (!strcmp(boundary->bctype,_THERMAL_SLIP_WALL_           )) boundary->BCFunctionU = BCThermalSlipWallU;    
    else if (!strcmp(boundary->bctype,_THERMAL_NOSLIP_WALL_         )) boundary->BCFunctionU = BCThermalNoslipWallU;    
    else if (!strcmp(boundary->bctype,_SW_SLIP_WALL_                )) boundary->BCFunctionU = BCSWSlipWallU;    
    else if (!strcmp(boundary->bctype,_SUBSONIC_OUTFLOW_            )) boundary->BCFunctionU = BCSubsonicOutflowU;    
    else if (!strcmp(boundary->bctype,_SUBSONIC_INFLOW_             )) boundary->BCFunctionU = BCSubsonicInflowU;    
    else if (!strcmp(boundary->bctype,_SUBSONIC_AMBIVALENT_         )) boundary->BCFunctionU = BCSubsonicAmbivalentU;    
    else if (!strcmp(boundary->bctype,_SUPERSONIC_OUTFLOW_          )) boundary->BCFunctionU = BCSupersonicOutflowU;    
    else if (!strcmp(boundary->bctype,_SUPERSONIC_INFLOW_           )) boundary->BCFunctionU = BCSupersonicInflowU;    
    else if (!strcmp(boundary->bctype,_TURBULENT_SUPERSONIC_INFLOW_ )) boundary->BCFunctionU = BCTurbulentSupersonicInflowU;    
    else if (!strcmp(boundary->bctype,_NO_FLUX_BC_                  )) boundary->BCFunctionU = BCNoFluxU;    
    else {
      fprintf(stderr,"Error in BCInitialize(): \"%s\" is not a supported boundary condition.\n",
              boundary->bctype);
      return(1);
    }

#if defined(HAVE_CUDA)
  }
#endif

  return 0;
}