text: #ifndef MONOALG3D_MODEL_DIFRANCESCO_1985_H
  #define MONOALG3D_MODEL_DIFRANCESCO_1985_H
  
  #include "../model_common.h"
  #include <stdint.h>
  
  #define NEQ 16
  #define INITIAL_V (-87.0)
  
  #ifdef __CUDACC__
  
  __constant__  size_t pitch;
  size_t pitch_h;
  
  __global__ void kernel_set_model_inital_conditions(real *sv, int num_volumes);
  
  __global__ void solve_gpu(real dt, real *sv, real* stim_currents,
                            uint32_t *cells_to_solve, uint32_t num_cells_to_solve,
                            int num_steps);
  
  inline __device__ void RHS_gpu(real *sv_, real *rDY_, real stim_current, int threadID_);
  
  #endif
  
  void solve_model_ode_cpu(real dt, real *sv, real stim_current);
  void RHS_cpu(const real *sv, real *rDY_, real stim_current);
  
  #endif //MONOALG3D_FHN_MOD_H
  
  