Phie Recovery

Module: devtests.bidomain.phie_recovery.run

Section author: Anton Prassl <anton.prassl@medunigraz.at>

Test phie recovery in some nodes above a simple slab.

Problem Setup

This example defines a thin cuboid domain (units in microns):

-1500  \leq x \leq 1500

-2500 \leq y \leq  2500

- 500 \leq z \leq   500

Run a modomain simulation in a cuboid portion of tissue and recover the phi_e signals above the corners of the geometry. Run the simulation with the same settings in postprocessing mode to recover phi_e based on the previous simulation. Investigate the difference of the two simulation modes (i.e. experiment “0” [normal run] and experiment “4” [postprocess only]).

../../_images/phie_recovery_setup.png

As with other examples, add the --visualize option to automatically load the results in meshalyzer:

./run.py --visualize

Tests

serial_recv_meth_1

Phie recovery during simulation and as postprocessing step on single core (using FE mass and stiffness matrices)

Tags: FAST SERIAL

Checks:

  • Compare against stored reference: max_error(phie_recovery.igb)
  • Compare against stored reference: max_error(phie_recovery_pp.igb)

serial_recv_meth_2

Phie recovery during simulation and as postprocessing step on single core (by integration over element-centered Im using integral solution of Poisson PDE)

Tags: FAST SERIAL

Checks:

  • Compare against stored reference: max_error(phie_recovery.igb)
  • Compare against stored reference: max_error(phie_recovery_pp.igb)

parallel_recv_meth_1

Phie recovery during simulation and as postprocessing step on multiple cores (using FE mass and stiffness

Tags: FAST PARALLEL

Checks:

  • Compare against stored reference: max_error(phie_recovery.igb)
  • Compare against stored reference: max_error(phie_recovery_pp.igb)

parallel_recv_meth_2

Phie recovery during simulation and as postprocessing step on multiple cores (by integration over element-centered Im using integral solution of Poisson PDE)

Tags: FAST PARALLEL

Checks:

  • Compare against stored reference: max_error(phie_recovery.igb)
  • Compare against stored reference: max_error(phie_recovery_pp.igb)