hypar.h

Go to the documentation of this file.

#ifndef _HYPAR_H_
#define _HYPAR_H_

#include <basic.h>

typedef struct main_parameters {

  int     ndims;

  int     nvars;

  int     *dim_global;

  int     *dim_local;

  int     npoints_global;

  int     npoints_local, 
          npoints_local_wghosts; 
  int ndof_nodes; 
  int ndof_cells_wghosts; 
  int size_x;     
  int     ghosts;

  int     n_iter;

  int     restart_iter;

  int     my_idx;

  int     nsims;

  double  dt;

  int     *dim_global_ex;

  char    time_scheme         [_MAX_STRING_SIZE_];

  char    time_scheme_type    [_MAX_STRING_SIZE_];

  char    spatial_scheme_hyp  [_MAX_STRING_SIZE_];

  char    interp_type         [_MAX_STRING_SIZE_];

  char    SplitHyperbolicFlux [_MAX_STRING_SIZE_];

  char    spatial_type_par    [_MAX_STRING_SIZE_];

  char    spatial_scheme_par  [_MAX_STRING_SIZE_];

  int    *index;

  double *x;

  double *dxinv;

  double *u;

  double *hyp;

  double *par;

  double *source;

  double *fluxC;

  double *uC;

  double *fluxI;

  double *uL, 
         *uR, 
         *fL, 
         *fR; 
  double  *Deriv1, 
          *Deriv2; 
  int   nBoundaryZones;
  void  *boundary;
  int   *isPeriodic;

  void *time_integrator;

  int screen_op_iter;

  int file_op_iter;

  int write_residual;

  char input_mode    [_MAX_STRING_SIZE_];

  char ip_file_type  [_MAX_STRING_SIZE_];

  char output_mode   [_MAX_STRING_SIZE_];

  char op_file_format[_MAX_STRING_SIZE_];

  char op_overwrite  [_MAX_STRING_SIZE_];

  char plot_solution[_MAX_STRING_SIZE_];

  char *filename_index;
  int  index_length;
  char solnfilename_extn[_MAX_STRING_SIZE_];
  char plotfilename_extn[_MAX_STRING_SIZE_]; 

  char op_fname_root[_MAX_STRING_SIZE_];
  char aux_op_fname_root[_MAX_STRING_SIZE_];

  int (*WriteOutput)              (int,int,int*,double*,double*,char*,int*);

  int (*ApplyBoundaryConditions)  (void*,void*,double*,double*,double);

  int (*ApplyIBConditions)  (void*,void*,double*,double);

  int (*TimeIntegrate)            (void*);

  int (*InterpolateInterfacesHyp) (double*,double*,double*,double*,int,int,void*,void*,int);

  int (*NonlinearInterp)    (double*,void*,void*,double,
                             int(*)(double*,double*,int,void*,double));

  int (*SetInterpLimiterVar)      (double*,double*,double*,int,void*,void*);

  int (*InterpolateInterfacesPar) (double*,double*,int,void*,void*);

  int (*FirstDerivativePar)       (double*,double*,int,int,void*,void*);

  int (*SecondDerivativePar)      (double*,double*,int,void*,void*);

  int (*HyperbolicFunction) (double*,double*,void*,void*,double,int,
                             int(*)(double*,double*,int,void*,double),
                             int(*)(double*,double*,double*,double*,double*,
                                    double*,int,void*,double));

  int (*ParabolicFunction)  (double*,double*,void*,void*,double);

  int (*SourceFunction)     (double*,double*,void*,void*,double);

  char model[_MAX_STRING_SIZE_];

  void *physics;

  double (*ComputeCFL)         (void*,void*,double,double);

  double (*ComputeDiffNumber)  (void*,void*,double,double);

  int    (*FFunction)          (double*,double*,int,void*,double);
  int    (*dFFunction)         (double*,double*,int,void*,double);
  int    (*FdFFunction)         (double*,double*,int,void*,double);
  int flag_fdf_specified;

  int    (*Upwind)             (double*,double*,double*,double*,double*,double*,int,void*,double);
  int    (*UpwinddF)           (double*,double*,double*,double*,double*,double*,int,void*,double);
  int    (*UpwindFdF)          (double*,double*,double*,double*,double*,double*,int,void*,double);

  int    (*GFunction)          (double*,double*,int,void*,double);

  int    (*HFunction)          (double*,double*,int,int,void*,double);

  int    (*SFunction)          (double*,double*,void*,void*,double);

  int    (*UFunction)          (double*,double*,int,void*,void*,double);

  int    (*JFunction)          (double*,double*,void*,int,int,int);

  int    (*KFunction)          (double*,double*,void*,int,int);

  int    (*PreStage)           (int,double**,void*,void*,double);
  int    (*PostStage)          (double*,void*,void*,double);
  int    (*PreStep)            (double*,void*,void*,double);
  int    (*PostStep)           (double*,void*,void*,double,int);
  int    (*PrintStep)          (void*,void*,double);
  int    (*PhysicsOutput)      (void*,void*,double);

  int    (*PhysicsInput)       (void*, void*,int,int,int*);

  int   (*AveragingFunction)   (double*,double*,double*,void*);

  int   (*GetLeftEigenvectors) (double*,double*,void*,int);
  int   (*GetRightEigenvectors)(double*,double*,void*,int);

  void *interp;
  void *compact;
  void *msti;
  void *lusolver;

  double error[3];
  double *ConservationError;
  char   ConservationCheck[_MAX_STRING_SIZE_];
  double *VolumeIntegral;
  double *VolumeIntegralInitial;
  double *StageBoundaryIntegral;
  double *StepBoundaryIntegral;
  double *TotalBoundaryIntegral;
  int    (*VolumeIntegralFunction)    (double*,double*,void*,void*);
  int    (*BoundaryIntegralFunction)  (void*,void*);
  int    (*CalculateConservationError)(void*,void*);

#ifdef with_petsc
  int     use_petscTS;  
  double  *u0;          
  double  *uref;        
  double  *rhsref;      
  double  *rhs;         
#endif

#ifdef with_librom
  double* u_rom_predicted; 
  double rom_diff_norms[3]; 
  char op_rom_fname_root[_MAX_STRING_SIZE_];
#endif

  int flag_nonlinearinterp;

  int *stride_with_ghosts;
  int *stride_without_ghosts;

  int count_hyp, 
      count_par, 
      count_sou; 
#ifdef with_petsc
  int count_RHSFunction,    
      count_IFunction,      
      count_IJacobian,      
      count_RHSJacobian,    
      count_IJacFunction,   
      count_RHSJacFunction; 
#endif

  double *iblank;

  char ib_filename[_MAX_STRING_SIZE_];
  int flag_ib;
  void *ib;

  int (*IBFunction) (void*,void*,double*,double);

#if defined(HAVE_CUDA)
  int use_gpu;       
  int gpu_device_no; 
  int gpu_npoints_boundary_offset[3]; 
  int gpu_npoints_boundary[3]; 
  int *gpu_dim_local; 
  double *gpu_iblank; 
  double *gpu_x;      
  double *gpu_dxinv;  
  double *gpu_u;      
  int StageBoundaryBuffer_size; 
  double *StageBoundaryBuffer; 
#endif

#ifdef with_python
  void* py_plt_func; 
  void* py_plt_func_args; 
#endif

} HyPar;

/* Some definitions - types of discretizations available
   for the parabolic (2nd derivative) term  */
#define _NC_1STAGE_   "nonconservative-1stage"

#define _NC_2STAGE_   "nonconservative-2stage"

#define _NC_1_5STAGE_ "nonconservative-1.5stage"

#define _CONS_1STAGE_ "conservative-1stage"

#endif