`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]).

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: 2023-03-25 00:57:33.976244, revision {‘base’: ‘6b57420d’}, dependency revisions {PT_C: 31642c1e,cvsys: b052f3a9,eikonal: b58be0de,elasticity: fadd9ded}

Runtime: 0:00:01.832196
ALL PASSED
  PASS max_error(phie_recovery.igb): 7.450580596923828e-09
  PASS max_error(phie_recovery_pp.igb): 7.450580596923828e-09

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: 2023-03-25 00:57:35.836998, revision {‘base’: ‘6b57420d’}, dependency revisions {PT_C: 31642c1e,cvsys: b052f3a9,eikonal: b58be0de,elasticity: fadd9ded}

Runtime: 0:00:01.558766
ALL PASSED
  PASS max_error(phie_recovery.igb): 7.450580596923828e-09
  PASS max_error(phie_recovery_pp.igb): 1.4901161193847656e-08

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: 2023-03-25 00:57:37.432870, revision {‘base’: ‘6b57420d’}, dependency revisions {PT_C: 31642c1e,cvsys: b052f3a9,eikonal: b58be0de,elasticity: fadd9ded}

Runtime: 0:00:01.370919
ALL PASSED
  PASS max_error(phie_recovery.igb): 7.450580596923828e-09
  PASS max_error(phie_recovery_pp.igb): 7.450580596923828e-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: 2023-03-25 00:57:38.835297, revision {‘base’: ‘6b57420d’}, dependency revisions {PT_C: 31642c1e,cvsys: b052f3a9,eikonal: b58be0de,elasticity: fadd9ded}

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