Example e84: Fitting Activation Times to a Bilinear Field with Bicubic Geometry

Fits activation times measured experimentally on the epicardial surface. The heart mesh (region 6 of a 6 region mesh) has a bicubic geometry and the field is a bilinear field. Only the lower portion of the epicardium (corresponding to the ventricles) are displayed.

Front Activation - Right Ventricle  on the left side Back Activation - Right Ventricle on the right side


The comfile run by this example is as follows:

# Example e84
# -----------
#


# Inputs
#

$MODEL = "pig06_cust2_heart"
$NRHEART = 6


# Controls
#

$USEFIT = 1
$USEEXPORTCMGUI = 1

fem def para;r;fitting;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 $NRHEART

fem def elem;r;${MODEL}_lin;example reg $NRHEART


# Grouping specific elements 
#
fem group elem 3265..4192 as HEARTELEM reg $NRHEART
fem group elem 3265..3904 as EPIELEM   reg $NRHEART
fem group elem 3905..4192 as VENDOELEM reg $NRHEART


# Grouping specific elements corresponding to ventricular epicardium
#
fem group elem 3265..3289,3305..3329,3345..3369,3385..3409,3425..3449,3465..3489 as EPI_BOTTOM1 reg $NRHEART
fem group elem 3505..3529,3545..3569,3585..3609,3625..3649,3665..3684,3705..3724 as EPI_BOTTOM2 reg $NRHEART
fem group elem 3745..3769,3785..3809,3825..3844,3865..3884 as EPI_BOTTOM3 reg $NRHEART


if( $USEEXPORTCMGUI==1 )
#
#----------------------------#
# Export geometry
#----------------------------#
#

{
fem export node;pig as thenodes reg all
fem export elem;heart as heart    reg $NRHEART elem EPIELEM
} #USEEXPORTCMGUI

Additional testing commands:

#
#----------------------------#
# Define fields
#----------------------------#
#
fem def fiel;r;heart;example reg $NRHEART
fem def elem;r;heart;example  reg $NRHEART field


fem write iod;pig06_cust2


#
#----------------------------#
# Define the data
#----------------------------#
#
fem def data;r;epi_elec_field;example field reg $NRHEART
fem def xi;c close elem EPIELEM             reg $NRHEART
fem def xi;w;elect                          reg $NRHEART
fem export data;epi_elec as epi_elec error  reg $NRHEART


if( $USEFIT==1 )


# Command updated by fixcom.sh on Wed Aug 23 17:50:30 NZT 2000
# Old command: fem def fit;r;heart_lin;example field reg $NRHEART
{
fem def fit;r;heart_lin;example field reg $NRHEART

fem fit reg $NRHEART



# Command updated by fixcom.sh on Wed Aug 23 17:50:30 NZT 2000
# Old command: endif #USEFIT
} #USEFIT



if( $USEEXPORTCMGUI==1 )
#
#----------------------------#
# Export gemetry
#----------------------------#
#

# Command updated by fixcom.sh on Wed Aug 23 17:50:30 NZT 2000
# Old command: fem export node;activ       as thenodes2   reg $NRHEART offset      10000
{
fem export node;activ as thenodes2 reg $NRHEART offset 10000

fem export elem;activ_epi   as activ_epi   reg $NRHEART offset_elem 10000 elem EPI_BOTTOM1,EPI_BOTTOM2,EPI_BOTTOM3
fem export elem;activ_vendo as activ_vendo reg $NRHEART offset_elem 10000 elem VENDOELEM



# Command updated by fixcom.sh on Wed Aug 23 17:50:30 NZT 2000
# Old command: endif #USEEXPORTCMGUI
} #USEEXPORTCMGUI








Files used by this example are:

Name                          Modified     Size

example_e84.com 14-May-2001 1.1k activation.data 10-Apr-2000 1.1k bicubic_bilinear.ipbase 10-Apr-2000 7.8k epi_elec_field.ipdata 10-Apr-2000 9.7k fitted.binhis 10-Apr-2000 28k fitting.ippara 12-Nov-2002 5.9k heart.ipelfd 10-Apr-2000 138k heart.ipfiel 10-Apr-2000 78k heart_lin.ipfit 13-Apr-2007 1.6k pig06_cust2.irnode 10-Apr-2000 3.5M pig06_cust2_heart.ipnode 14-May-2001 791k pig06_cust2_heart_lin.ipelem 14-May-2001 368k pig06_cust2_lin.irelem 10-Apr-2000 1.6M six.ipregi 10-Apr-2000 93 test_output.com 20-Nov-2001 1.5k

Download the entire example:

Name                      Modified     Size

examples_e_e8_e84.tar.gz 12-Mar-2008 1.4M

Testing status by version:

StatusTestedReal time (s)
i686-linux
cmSuccessSun Mar 6 00:03:28 20163
cm-debugSuccessSat Mar 5 00:15:05 20166
mips-irix
cmSuccessSun Aug 19 02:00:30 200729
cm-debugSuccessWed Aug 15 01:55:56 200772
cm-debug-clear-mallocSuccessSat Aug 18 02:18:47 200783
cm-debug-clear-malloc7SuccessMon Aug 20 02:11:40 200784
cm64SuccessSun Aug 19 02:01:00 200731
cm64-debugSuccessTue Aug 21 01:56:47 200777
cm64-debug-clear-mallocSuccessThu Apr 1 12:11:17 200435
rs6000-aix
cmSuccessWed Mar 4 01:10:04 20093
cm-debugSuccessMon Mar 2 01:12:06 200918
cm64SuccessWed Mar 4 01:11:01 20093
cm64-debugSuccessTue Mar 3 01:17:34 200919
x86_64-linux
cmSuccessSun Mar 6 00:01:08 20162
cm-debugSuccessSat Mar 5 00:01:59 20163

Testing status by file:


Graphical output from this problem is given here.


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

Input last modified: Tue Mar 11 16:32:22 2008


CMISS Help / Examples / e / e8 / e84