Initialization.H

Go to the documentation of this file.

#include "PlasCom2.H"
#include "PCPPCommUtil.H"
#include "PCPPIntervalUtils.H"
#include "PCPPReport.H"
#include "EulerUtil.H"
#include "ViscidUtil.H"
#include "MaxwellUtil.H"

template<typename GridType,typename StateType,typename BoundaryType>
int InitializeSolution(const GridType &inGrid,StateType &inState,
                       StateType &inParam,std::vector<BoundaryType> &inBoundaries,
                       const std::string &initName,const std::vector<double> &inParams,
                       const std::vector<int> &inFlags,int threadId,
                       std::ostream *messageStream = NULL)
{
  if(!messageStream)
    messageStream = &std::cout;

  if(initName == "DensityPulse2"){
    return(euler::util::InitializeDensityPulse(inGrid,inState,inParam,inParams,inFlags,threadId,messageStream));
  } else if(initName == "Uniform"){
    return(euler::util::InitializeUniformFlow(inGrid,inState,inParam,inParams,threadId,messageStream));
  } else if(initName == "DensityPulse"){
    return(euler::util::InitializeDensityPulse(inGrid,inState,inParam,inParams,threadId,messageStream));
  } else if(initName == "AcousticPulse"){
    return(euler::util::InitializeAcousticPulse(inGrid,inState,inParam,inParams,inFlags,threadId,messageStream));
  } else if(initName == "Shock1D"){
    return(euler::util::InitializeShock1D(inGrid,inState,inParam,inParams,threadId,messageStream));
  } else if(initName == "Poiseuille"){
    return(viscid::util::InitializePoiseuille(inGrid,inState,inParam,inParams,threadId,messageStream));
  } else if(initName == "ScalarAdvection"){
    return(euler::util::InitializeGaussianScalar(inGrid,inState,inParam,inParams,threadId,messageStream));
  } else if(initName == "AdvectionDiffusion"){
    return(euler::util::InitializeAdvectionDiffusion(inGrid,inState,inParam,inParams,threadId,messageStream));
  } else if(initName == "ProtoY4Test1"){
    return(euler::util::InitializeProtoY4Test1(inGrid,inState,inParam,
                                               inBoundaries,inParams,
                                               inFlags,threadId));
  } else if(initName == "ConvectingVortex") {
    return(euler::util::InitializeConvectingVortex(inGrid, inState, inParam, inParams, inFlags, threadId));
  } else if(initName == "Riemann1D") {
    return(euler::util::InitializeRiemann1D(inGrid, inState, inParam, inParams, inFlags, threadId));
  } else if(initName == "Shocktube") {
    return(euler::util::InitializeShocktube(inGrid, inState, inParam, inParams, inFlags, threadId));
  } else if(initName == "GaussianPulse1DXDir") {
    return(Maxwell::util::InitializeGaussianPulse1DXDir(inGrid, inState, inParam, inParams, inFlags, threadId));
  }
  return(0);
};