Example e31: Male torso - no salu

Forward simulation on homogeneous (heart-torso BEM model) with tradition transfer matrix calculations. This problem does not use the Salu consistency criteria. An un-scaled dipole source is used.


The comfile run by this example is as follows:

# Example e31 
# -----------
#


#----------------------------------------------------------#
#
# 3D time dependent torso problem
# This example shows a heart in a torso with a moving 
# dipole inside the heart -ie "quasistatic" 3d.  The 
# dipole path was obtained from Green Lane Hospital VCG.
#
#----------------------------------------------------------#


fem def para;r;torsotimedep;example          #Sets up parameters for this prob.
fem def regi;r;torsotimedep;example          #Defines two regions.
fem def coor 3,1
fem def node;r;;example region all           #Defines 291 nodes.
fem def base;r;;example                      #Defines 3 types of basis function.
fem def elem;r;;example region all           #Defines 304 elements.


#
#----------------------------#
# Export geometry
#----------------------------#
#
fem export node;hearttorso as hearttorso
fem export elem;hearttorso as hearttorso


#
#----------------------------#
# Setup equations
#----------------------------#
#
fem def equa;r;torsotimedep;example region all
#                                            #Defines a quasistatic analysis of
#                                            # the standard Laplace equation.
fem def mate;r;;example region all
fem def init;r;;example region all
fem def sour;r;;example region 2             #Defines the dipole path.
fem list source region 2
fem def coup;r;;example
fem def solv;r;;example coupled region all
fem solve coupled

Additional testing commands:

fem list node solution
fem def tran;r;;example
fem evaluate transfer svd
fem def inve;r;;example

#fem evaluate inverse svd
#fem write matrix matrices T_BH,T_BH_INV
#fem apply inverse history infile 3dhearttorso outfile inverse_3dhearttorso
#fem check solution inverse infile1 3dhearttorso infile2 inverse_3dhearttorso

Files used by this example are:

Name                     Modified     Size

example_e31.com 20-Nov-2001 1.4k hearttorso.irequa 26-May-2003 3.5k test_output.com 20-Nov-2001 319 torsotimedep.ipbase 10-Apr-2000 4.5k torsotimedep.ipcoup 21-Aug-2002 464 torsotimedep.ipcoup.old 10-Apr-2000 413 torsotimedep.ipinve 06-Dec-2001 587 torsotimedep.ippara 12-Nov-2002 5.9k torsotimedep.ipregi 10-Apr-2000 120 torsotimedep.ipsour 10-Apr-2000 84k torsotimedep.iptrsf 01-Nov-2001 770 torsotimedep.irelem 10-Apr-2000 109k torsotimedep.irequa 26-May-2003 3.3k torsotimedep.irinit 10-Apr-2000 160k torsotimedep.irmate 10-Apr-2000 156 torsotimedep.irnode 10-Apr-2000 249k torsotimedep.irsolv 13-Apr-2007 2.0k

Download the entire example:

Name                      Modified     Size

examples_e_e3_e31.tar.gz 14-Apr-2007 118k

Testing status by version:

StatusTestedReal time (s)
i686-linux
cmSuccessSun Mar 6 00:08:05 201610
cm-debugSuccessSat Mar 5 00:40:01 201629
mips-irix
cmSuccessSun Aug 19 03:37:21 2007142
cm-debugSuccessWed Aug 15 06:13:33 2007536
cm-debug-clear-mallocSuccessSat Aug 11 04:12:34 2007544
cm-debug-clear-malloc7SuccessMon Aug 20 03:35:52 2007545
cm64SuccessSun Aug 19 03:42:23 2007146
cm64-debugSuccessFri Aug 17 04:28:43 2007614
cm64-debug-clear-mallocSuccessThu Apr 1 12:03:50 2004252
rs6000-aix
cmSuccessWed Mar 4 01:13:14 20098
cm-debugSuccessMon Mar 2 01:33:16 2009122
cm64SuccessWed Mar 4 01:13:35 20098
cm64-debugSuccessTue Mar 3 01:39:01 2009137
x86_64-linux
cmSuccessSun Mar 6 00:01:24 20165
cm-debugSuccessSat Mar 5 00:05:00 201624

Testing status by file:


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

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


CMISS Help / Examples / e / e3 / e31