physics_types.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) NumeriCor GmbH
// SPDX-License-Identifier: LicenseRef-APL-1.1
 
#ifndef PHYSICS_TYPES_H
#define PHYSICS_TYPES_H
 
#include "basics.h"
 
namespace opencarp {
 
// unit conversion
#define UM2_to_CM2                1.0e-8              
#define UM_to_CM                  1.0e-4              
#define UM_to_MM                  1.0e-3              
#define MM_to_UM                  1.0e3               
 
 
enum physic_t {
  ion_phys=0, elec_phys, laplace_phys, eikonal_phys, emi_phys, mech_phys, diffK_phys, magnet_phys, opt_phys,
  fluid_phys, NUM_PHYSICS
};
 
class Basic_physic {
public:
  const char* name = NULL;
  FILE_SPEC logger = NULL;
  int timer_idx = -1;
 
  // For very heavy objects like the physics, we might not want to use the
  // constructor / destructor mechanisms, to avoid expensive operations when
  // someone (accidentaly) creates a local copy
  virtual void initialize() = 0;
  virtual void destroy()    = 0;
  virtual void compute_step() = 0;
  virtual void output_step() = 0;
 
  virtual inline void output_timings()
  {
    if(name) log_msg(0,0,0, "%s:", name);
    else     log_msg(0,0,0, "<unnamed>:");
    log_msg(0,0,0, "  Init:    %.2lf sec", initialize_time);
    log_msg(0,0,0, "  Compute: %.2lf sec", compute_time);
    log_msg(0,0,0, "  Output:  %.2lf sec", output_time);
    log_msg(0,0,0, "");
  }
 
  // we also store timings for the main Basic_physic API
  double initialize_time = 0.0;
  double compute_time    = 0.0;
  double output_time     = 0.0;
 
  virtual std::string timer_unit(const int timer_id) = 0;
  virtual double timer_val(const int timer_id) = 0;
};
 
 
enum datavec_t {
  vm_vec, iion_vec, vm_emi_itf, iion_emi_itf
};
 
}  // namespace opencarp
 
#endif