timeintegration_struct.h

Go to the documentation of this file.

#ifndef _TIME_INTEGRATION_STRUCT_H_
#define _TIME_INTEGRATION_STRUCT_H_

#include <sys/time.h>
#include <basic.h>

/* definitions */
#define _FORWARD_EULER_ "euler"

#define _RK_            "rk"

#define _GLM_GEE_       "glm-gee"

typedef struct time_integration_variables {
  int     iter;
  int     n_iter;
  int     restart_iter;
  double  waqt;
  double  dt;
  double  norm;
  double  max_cfl;
  double  max_diff;

  void    *simulation;
  int     nsims;

  long    *u_offsets;

  long    *u_sizes;

  int     *bf_offsets;

  int     *bf_sizes;

  double  *u;

  double  *rhs;

  long    u_size_total;

  long    bf_size_total;

  double  **U;
  double  **Udot;

  int     rank;
  int     nproc;

  double **BoundaryFlux;

  void *ResidualFile;

  int (*TimeIntegrate) (void*);
  int (*RHSFunction)   (double*,double*,void*,void*,double);

  struct timeval iter_start_time;
  struct timeval iter_end_time;
  double iter_wctime;
  double iter_wctime_total;

#if defined(HAVE_CUDA)

  double *gpu_U;
  double *gpu_Udot;
  double *gpu_BoundaryFlux;
#endif

} TimeIntegration;

/* Explicit Runge-Kutta Methods */
#define _RK_1FE_        "1fe"

#define _RK_22_         "22"

#define _RK_33_         "33"

#define _RK_44_         "44"

#define _RK_SSP3_       "ssprk3"
#define _RK_TVD3_       "tvdrk3"  
typedef struct _explicit_rungekutta_time_integration_ {
  int nstages; 
  double *A, 
         *b, 
         *c; 
} ExplicitRKParameters;

/* General Linear Methods with Global Error Estimate */
#define _GLM_GEE_YYT_     "yyt"

#define _GLM_GEE_YEPS_    "yeps"

#define _GLM_GEE_23_      "23"

#define _GLM_GEE_24_      "24"

#define _GLM_GEE_25I_     "25i"

#define _GLM_GEE_35_      "35"

#define _GLM_GEE_EXRK2A_  "exrk2a"

#define _GLM_GEE_RK32G1_  "rk32g1"

#define _GLM_GEE_RK285EX_ "rk285ex"

typedef struct _glm_gee_time_integration_ {
  int nstages, 
      r;       
  char ee_mode[_MAX_STRING_SIZE_]; 
  double  *A_yyt,   
          *B_yyt ,  
          *C_yyt ,  
          *D_yyt ,  
          *c_yyt;   
  double  *A_yeps,  
          *B_yeps,  
          *C_yeps,  
          *D_yeps,  
          *c_yeps;  
  double *A, 
         *B, 
         *C, 
         *D, 
         *c; 
  double gamma; 
} GLMGEEParameters;

#endif