Example e34: Solving for Potentials from an Activation Map

Creates potentials from an activation map fitted from example e33. The activation sequence from 0 to 50 ms is used to construct potentials on the torso surface with the use of a double layer transfer matrix.

The images below show front and back view of the torso and the displaced epicardial and ventricular surfaces. The black line on the heart shows the location of the wavefront. The activation sequence on the heart and ventricles (0-50ms) are shown with the resulting potential field on the torso surface. The left image is a snapshot at 10 ms and the right image at 35 ms.

10ms 35ms


The comfile run by this example is as follows:

# Example e34
# -----------
#


# Times to use
#

$STARTTIME = 0.0
$ENDTIME = 50.0


# Inputs
#

$VERSION 	= "01"
$MODEL 		= "pig$VERSION"
$OUTPUT 	= "forward_$VERSION"

 

# Controls
# 

$CALC_TRANSFER 	= 1
$EXPORT_CMGUI 	= 1
$EXPORT_UNEMAP 	= 1
$FORMAT 	= "ascii"



fem def para;r;activation;example 
fem def regi;r;six;example
fem def coor 3,1
fem def base;r;bicubic_bilinear;example
fem def node;r;$MODEL;example reg all
fem def elem;r;${MODEL}_lin;example reg all

 
 
# Grouping specific nodes
#
fem group node 225      as FIXEDNODE
fem group node 186,439  as SURFENDS
fem group node 187..438 as SURFNORM
 

# Grouping specific elements
#
fem group elem 697..736  as EPI    reg 6
fem group elem 737..808  as VENDO  reg 6
fem group elem 201..464  as SKIN   reg 1


if( $EXPORT_CMGUI==1 )
{
#
#----------------------------#
# Export geometry
#----------------------------#
#

fem export node;pig as thenodes reg all

fem export elem;skin   as skin      reg 1     elem SKIN
fem export elem;lungs  as lungs     reg 2,3
fem export elem;lvendo as lvendo    reg 4
fem export elem;rvendo as rvendo    reg 5
fem export elem;heart  as heart     reg 6     #elem EPI
 
} #EXPORT_CMGUI

 

#
#----------------------------#
# Setup equations
#----------------------------#
#
fem def equa;r;forward;example   reg all 
fem def mate;r;forward;example reg all
fem def init;r;forward_$VERSION;example reg all

fem def coup;r;forward;example
fem def solv;r;forward;example coup reg all

Additional testing commands:


fem def tran;r;forward_$VERSION;example
fem eval transfer svd 
fem write matrix;forward matric T_BH FORMAT
fem def init;r;forward_$VERSION;example reg all

fem app tran activ outf $OUTPUT $FORMAT tsta $STARTTIME tend $ENDTIME
fem eval elec;$OUTPUT hist $OUTPUT from node $FORMAT


fem li sign sign $OUTPUT $FORMAT 
fem conv sign from ascii infile $OUTPUT outfile $OUTPUT
fem def export;r;unemap_439;example
fem export sign;$OUTPUT sign $OUTPUT elec


if( $EXPORT_CMGUI==1 )
{
#
#----------------------------#
# Export field
#----------------------------#
#

fem export node;$OUTPUT history $OUTPUT $FORMAT as soln

fem export elem;$OUTPUT field as soln elem SKIN
fem exp node;activ       as thenodes2   fiel name activ reg 6
fem exp elem;activ_epi   as activ_epi   fiel name activ reg 6 elem EPI
fem exp elem;activ_vendo as activ_vendo fiel name activ reg 6 elem VENDO

} #EXPORT_CMGUI


Files used by this example are:

Name                     Modified     Size

example_e34.com 01-Nov-2001 1.4k activation.ippara 12-Nov-2002 5.9k bicubic_bilinear.ipbase 10-Apr-2000 7.8k draw.com 09-Nov-2001 3.1k forward.ipcoup 21-Aug-2002 437 forward.ipcoup.old 10-Apr-2000 411 forward.irequa 26-May-2003 7.2k forward.irmate 10-Apr-2000 4.7k forward.irsolv 13-Apr-2007 1.6k forward_01.binsig 10-Apr-2000 106k forward_01.iptrsf 01-Nov-2001 1.1k forward_01.irinit 10-Apr-2000 13k pig01.irnode 10-Apr-2000 673k pig01_lin.irelem 10-Apr-2000 330k six.ipregi 10-Apr-2000 93 test_output.com 01-Nov-2001 888 unemap_439.ipexpo 10-Apr-2000 661

Download the entire example:

Name                      Modified     Size

examples_e_e3_e34.tar.gz 14-Apr-2007 468k

Testing status by version:

StatusTestedReal time (s)
i686-linux
cmSuccessSun Mar 6 00:07:54 20169
cm-debugSuccessSat Mar 5 00:44:14 201633
mips-irix
cmSuccessSun Aug 19 03:22:23 2007125
cm-debugSuccessWed Aug 15 05:15:47 2007449
cm-debug-clear-mallocSuccessSat Aug 18 05:49:54 2007459
cm-debug-clear-malloc7SuccessMon Aug 20 06:02:06 2007476
cm64SuccessSun Aug 19 03:26:34 2007124
cm64-debugSuccessTue Aug 21 03:51:08 2007476
cm64-debug-clear-mallocSuccessThu Apr 1 12:09:29 2004185
rs6000-aix
cmSuccessWed Mar 4 01:15:49 200912
cm-debugSuccessMon Mar 2 01:30:38 2009109
cm64SuccessWed Mar 4 01:16:03 200914
cm64-debugSuccessTue Mar 3 01:36:55 2009127
x86_64-linux
cmSuccessSun Mar 6 00:01:21 20164
cm-debugSuccessSat Mar 5 00:04:16 201617

Testing status by file:


Graphical output from this problem is given here.


Html last generated: Sun Mar 6 05:51:27 2016

Input last modified: Fri Apr 13 10:41:52 2007


CMISS Help / Examples / e / e3 / e34