testfem.cc

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) NumeriCor GmbH
// SPDX-License-Identifier: LicenseRef-APL-1.1
 
#include <cstdio>
#include <csignal>
#include <vector>
 
// we set PrMGLOBAL ONLY ONCE. this will initialize the global variables
// in only one compilation unit
#define PrMGLOBAL
#include "simulator.h"
#include "numerics.h"
#include "physics.h"
#include "fem.h"
 
 
// globals
namespace opencarp {
namespace user_globals {
  SF::scatter_registry<mesh_int_t> scatter_reg;
  std::map<mesh_t, sf_mesh> mesh_reg;
  std::map<SF::quadruple<int>, SF::index_mapping<mesh_int_t> > map_reg;
  std::map<physic_t, Basic_physic*> physics_reg;
  timer_manager* tm_manager;
  std::map<datavec_t, sf_vec*> datavec_reg;
  FILE* petsc_error_fd = NULL;
  bool using_legacy_stimuli = false;
}  // namespace user_globals
}  // namespace opencarp
 
using namespace opencarp;
 
#undef  __FUNCT__
#define __FUNCT__ "main"
int main(int argc, char** argv)
{
  char *DBfile   = NULL;  // could be used to pass in ~/.petscrc
  char *help_msg = NULL;
  parse_params_cpy(argc, argv);
  initialize_PETSc(&argc, argv, DBfile, help_msg);
 
  setup_petsc_err_log();
  check_and_convert_params();
 
  // set up meshes
  std::map<mesh_t, sf_mesh> & mesh_registry = user_globals::mesh_reg;
  mesh_registry[reference_msh] = sf_mesh();
  mesh_registry[reference_msh].comm = PETSC_COMM_WORLD;
  mesh_registry[extra_elec_msh] = sf_mesh();
 
  setup_meshes();
 
  MaterialType mtype[2];
  FILE_SPEC logger = f_open("electrics.log", param_globals::experiment != 4 ? "w" : "r");
 
  set_elec_tissue_properties(mtype, Electrics::extra_grid, logger);
  region_mask(intra_elec_msh, mtype[Electrics::extra_grid].regions,
              mtype[Electrics::extra_grid].regionIDs, true, "gregion_e");
 
  SF::vector<stimulus> stimuli;
  elliptic_solver ellip_solver;
  ellip_solver.init();
  ellip_solver.rebuild_matrices(mtype, stimuli, logger);
 
  cleanup_and_exit();
}