bench_utils.h

Go to the documentation of this file.

// ----------------------------------------------------------------------------
// openCARP is an open cardiac electrophysiology simulator.
//
// Copyright (C) 2020 openCARP project
//
// This program is licensed under the openCARP Academic Public License (APL)
// v1.0: You can use and redistribute it and/or modify it in non-commercial
// academic environments under the terms of APL as published by the openCARP
// project v1.0, or (at your option) any later version. Commercial use requires
// a commercial license (info@opencarp.org).
//
// This program is distributed without any warranty; see the openCARP APL for
// more details.
//
// You should have received a copy of the openCARP APL along with this program
// and can find it online: http://www.opencarp.org/license
// ----------------------------------------------------------------------------
 
#ifndef _BENCH_UTILS_H
#define _BENCH_UTILS_H
 
#include "cmdline.h"
#include "restitute.h"
 
class Bisect {
    public:
 
        Bisect( float min, float max, float tol=0.01 ):_min(min),_max(max),_tol(tol){}
        bool update( float &val, bool success ){
          if( success ){
            _curr = (val+_min)/2.;
            _max = val;
          } else {
            _curr = _max==INFINITY ? val*2 : (val+_max)/2.;
            _min = val;
          }
          _found = fabs((val-_curr)/val) <= _tol;
          val = _curr;
          return _found;
        }
        float current(){ return _curr; }
    private:
        float _min;           
        float _max;           
        float _curr;          
        float _tol;           
        bool  _found = false;
};
 
namespace limpet {
 
enum TimerIDs {
  CON_TM_IDX,   // IO timer for console output
  SVD_TM_IDX,   // IO timer for state variable output
  STA_TM_IDX,   // IO timer for state saves (checkpointing)
  SSV_TM_IDX,   // IO timer for single state vector saves
  STM_TM_IDX,   // timer for stimulus trigger
  LIGHT_TM_IDX, // timer for illumination trigger
  TRACE_TM_IDX,  // timer for trace output
  CLAMP_TM_IDX,  // timer for trace output
  DOPPLE_TM_IDX,
  RES_SAVE_TM_IDX,  // timer for saving state vectors during restitution
  N_TIMERS      // number of IO Timers we are going to use
};
 
#define TIME_OUT     1.     // output period for data
#define NUMSEG       1          // number of nodes or nodes
 
// STATE VARIABLE DUMPS
#define DT_SV_DUMP 1.           // output interval for sv dumps
 
// typedefs
struct IOCtrl {
  char w2file;    
  char wbin;      
  char wsplt;     
  char w2stdout;  
  char first;     
};
 
struct GVEC_DUMP {
  opencarp::FILE_SPEC hdls[NUM_IMP_DATA_TYPES+1];   
  char*     fn[NUM_IMP_DATA_TYPES+1];   
  int       dtype[NUM_IMP_DATA_TYPES+1]; 
  int       n_dumps;    
};
 
// record statistics on the duration of events
#define us_RESOLUTION 1
struct event_timing {
  bool   init;   
  double mn;     
  double mx;     
  double avg;    
  double tot;    
  int    count;  
};
 
enum timingIDs {
  SETUP_IDX,        
  INIT_IDX,         
  LOOP_IDX,         
  ODE_IDX,          
  N_TIMINGS         
};
 
typedef enum {
  dtype_Gatetype,
  dtype_Real,
  dtype_Float,
  dtype_Char,
  dtype_Integer,
  dtype_Double,
  dtype_Short,
  dtype_GlobalData_t,
  dtype_long,
  dtype_bool
} DataType;
 
static const char * const data_type_names[] = {
  "Gatetype",
  "Real",
  "float",
  "char",
  "int",
  "double",
  "short",
  "GlobalData_t",
  "long",
  "bool"
};
 
// Gatetype is defined in LIMPET only and typically defined as float
// if defined otherwise, the size has to be overruled within LIMPET
// in this case the const qualifier has to be removed.
#define GATETYPE_SIZE      4
#define GLOBALDATA_T_SIZE  8
static const int data_type_sizes[] = {
  GATETYPE_SIZE,
  sizeof(limpet::Real),
  sizeof(float),
  sizeof(char),
  sizeof(int),
  sizeof(double),
  sizeof(short),
  GLOBALDATA_T_SIZE,
  sizeof(long),
  sizeof(bool)
};
 
 
int  write_dump_header(GVEC_DUMP *gvd, SV_DUMP *svd, const char *ExpID);
void open_globalvec_dump(FILE **fhdls, GVEC_DUMP *gvd, MULTI_IF *pMIIF, char *base_name, IOCtrl *io);
void globalvec_dump(FILE **fhdls, GVEC_DUMP *gvd, MULTI_IF *pMIIF, opencarp::timer_manager *tm, IOCtrl *io, int numNode);
void close_globalvec_dump(FILE **fhdls, GVEC_DUMP *gvd, IOCtrl *io);
void dump_all(MULTI_IF *, int, char *, char *, double, double, char *);
void update_timing(event_timing *t, double event_duration);
void initialize_timings(event_timing *t);
double getCellVal(opencarp::sf_vec* v, int ind);
void initial_SVs(MULTI_IF *miif, char *SVs, char *imp, char *plgins, int num);
void print_param_help(IonType* im, IonTypeList& plugs);
float determine_duration(struct gengetopt_args_info *p, TrgList *stim_lst);
void determine_stim_list(char *stl, TrgList *trg, bool DIAs);
 
}  // namespace limpet
 
#endif