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)

Last run: 2021-09-27 00:21:56.924340, revision {‘base’: ‘778659c7’}, dependency revisions {PT_C: 31642c1e,eikonal: b58be0de,elasticity: e4d37cfa}

Runtime: 0:00:01.599096
ALL PASSED
  PASS max_error(phie_recovery.igb): 0.0
  PASS max_error(phie_recovery_pp.igb): 0.0

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)

Last run: 2021-09-27 00:21:58.550145, revision {‘base’: ‘778659c7’}, dependency revisions {PT_C: 31642c1e,eikonal: b58be0de,elasticity: e4d37cfa}

Runtime: 0:00:01.249375
ALL PASSED
  PASS max_error(phie_recovery.igb): 0.0
  PASS max_error(phie_recovery_pp.igb): 0.0

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)

Last run: 2021-09-27 00:21:59.826697, revision {‘base’: ‘778659c7’}, dependency revisions {PT_C: 31642c1e,eikonal: b58be0de,elasticity: e4d37cfa}

Runtime: 0:00:04.142681
ALL PASSED
  PASS max_error(phie_recovery.igb): 1.862645149230957e-09
  PASS max_error(phie_recovery_pp.igb): 3.725290298461914e-09

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)

Last run: 2021-09-27 00:22:03.995535, revision {‘base’: ‘778659c7’}, dependency revisions {PT_C: 31642c1e,eikonal: b58be0de,elasticity: e4d37cfa}

Runtime: 0:00:04.234751
ALL PASSED
  PASS max_error(phie_recovery.igb): 6.332993507385254e-08
  PASS max_error(phie_recovery_pp.igb): 6.332993507385254e-08