Two BEM one GRID Laplace equation test

Example b213b: Two BEM one GRID Laplace equation test

This example demonstrates how to couple multiple boundary element regions by direct assembly and tests this against an analytic solution to Laplace's equation.

The comfile run by this example is as follows:

#Example b213b
# This example demonstrates how to couple multiple ventriclar
#  regions by direct assembly.
#

# Problem regions
$GRIDHT=1 
$LVENTL=2 
$RVENTL=3 
$BEMHRT=4 

# Problem classes
$GRDCLS=1 
$LVECLS=2 
$RVECLS=3 

# Define problem geometry
fem define para;r;minimal;example
fem define regi;r;torsoslice;example
fem define node;r;LV3;example reg $LVENTL
fem define node;r;RV3;example reg $RVENTL
fem define node;r;bemhrt3;example reg $BEMHRT
fem define node;r;gridhtr;example reg $GRIDHT
fem define base;r;linear;example
fem define elem;r;LV3;example reg $LVENTL
fem define elem;r;RV3;example reg $RVENTL
fem define elem;r;bemhrt3;example reg $BEMHRT
fem define elem;r;gridhtr;example reg $GRIDHT
fem define fibr;r;gridhtr;example reg $GRIDHT
fem define elem;r;gridhtr;example reg $GRIDHT fibre

# Define finite difference grid
fem define grid;r;gridhtr;example reg $GRIDHT
fem update grid geomet reg $GRIDHT
fem update grid metric reg $GRIDHT
fem group grid xi2 0 elem 11..20,27,28 as LVendo reg $GRIDHT
fem group grid xi2 1 elem 11,13,19,20 as RVendo1 reg $GRIDHT
fem group grid xi2 0 elem 1,3,9,10,21,22,25,26 as RVendo2 reg $GRIDHT
fem group grid grid RVendo1,RVendo2 as RVendo reg $GRIDHT
fem group grid grid LVendo,RVendo as endocardium reg $GRIDHT
fem group grid xi2 1 elem 1..10,21..26 as epicardium  reg $GRIDHT
fem group grid external as boundary reg $GRIDHT
fem group grid line 1845 xidirn 2 as stimulus reg $GRIDHT

# Create node to grid mapping
fem update node grid grid_reg $GRIDHT bem_reg $BEMHRT exclude node 77,78
fem update node geomet from grid reg $BEMHRT except node 77,78

# Define equations
fem define equa;r;solve8;example reg $GRIDHT class $GRDCLS lock
fem define equa;r;LV;example reg $LVENTL class $LVECLS lock
fem define equa;r;RV;example reg $RVENTL class $RVECLS lock

# Define material parameters
fem define mate;r;LV;example reg $LVENTL class $LVECLS
fem define mate;d reg $GRIDHT class $GRDCLS
fem update grid materi reg $GRIDHT class $GRDCLS
fem define mate;r;RV;example reg $RVENTL class $RVECLS

# Define initial conditions
fem define init;d reg $GRIDHT class $GRDCLS
fem define init;r;LV;example reg $LVENTL class $LVECLS
fem define init;r;RV;example reg $RVENTL class $RVECLS

# Define right ventricle normal reversals
fem define norm;r;ventcl;example

# Define solution parameters
fem define solv;r;solve8;example reg $GRIDHT class $GRDCLS
fem define solv;r;LV;example reg $LVENTL class $LVECLS
fem define solv;r;RV;example reg $RVENTL class $RVECLS

# Solve the passive torso problems
fem solve reg $LVENTL class $LVECLS
fem solve reg $RVENTL class $RVECLS

# Solve the bidomain problem for 5ms
fem solve activ reg $GRIDHT,$LVENTL,$RVENTL class $GRDCLS,$LVECLS,$RVECLS grid endocardium cptype 4

Additional testing commands:

#Testing example b213b
fem list node sol reg 2 cl 2
fem list node sol reg 3 cl 3

Files used by this example are:

Name               Modified     Size

example_b213b.com  06-Jun-2001  2.7k
LV.ipequa          26-May-2003  1.2k
LV.ipinit          06-Jun-2001  1.2k
LV.ipmate          06-Jun-2001  606 
LV.ipsolv          13-Apr-2007  1.2k
LV3.ipelem         06-Jun-2001   24k
LV3.ipnode         06-Jun-2001   28k
RV.ipequa          26-May-2003  1.2k
RV.ipinit          06-Jun-2001  1.2k
RV.ipmate          06-Jun-2001  606 
RV.ipsolv          13-Apr-2007  1.2k
RV3.ipelem         06-Jun-2001   24k
RV3.ipnode         06-Jun-2001   28k
bemhrt3.ipelem     06-Jun-2001   80k
bemhrt3.ipnode     06-Jun-2001   94k
gridhtr.ipelem     06-Jun-2001  7.5k
gridhtr.ipelfb     06-Jun-2001  4.3k
gridhtr.ipfibr     06-Jun-2001  7.0k
gridhtr.ipgrid     06-Mar-2003  561 
gridhtr.ipnode     06-Jun-2001   14k
linear.ipbase      06-Jun-2001  2.2k
minimal.ippara     12-Nov-2002  5.9k
solve8.ipequa      26-May-2003  1.2k
solve8.ipsolv      13-Apr-2007  1.2k
test_output.com    06-Jun-2001   82 
torsoslice.ipregi  06-Jun-2001   93 
ventcl.ipnorm      06-Jun-2001  418

Download the entire example:

Name                                 Modified     Size

examples_b_b2_b21_b213_b213b.tar.gz  14-Apr-2007   64k

Testing status by version:

	Status	Tested	Real time (s)
i686-linux
cm	Success	Sun Mar 6 00:02:26 2016	1
cm-debug	Success	Sat Mar 5 00:06:22 2016	2
mips-irix
cm	Success	Sun Aug 19 01:40:45 2007	11
cm-debug	Success	Wed Aug 15 01:36:49 2007	26
cm-debug-clear-malloc	Success	Sat Aug 18 01:43:32 2007	31
cm-debug-clear-malloc7	Success	Mon Aug 20 01:40:32 2007	32
cm64	Success	Sun Aug 19 01:41:46 2007	11
cm64-debug	Success	Tue Aug 21 01:35:05 2007	28
cm64-debug-clear-malloc	Success	Thu Apr 1 11:19:50 2004	13
rs6000-aix
cm	Success	Wed Mar 4 01:09:01 2009	2
cm-debug	Success	Mon Mar 2 01:08:58 2009	7
cm64	Success	Wed Mar 4 01:08:36 2009	2
cm64-debug	Success	Tue Mar 3 01:14:13 2009	6
x86_64-linux
cm	Success	Sun Mar 6 00:01:06 2016	1
cm-debug	Success	Sat Mar 5 00:01:43 2016	1

Testing status by file:

cmiss_test.log

i686-linux
Success	cm:	cmiss_test.log.retain.
Success	cm-debug:	cmiss_test.log.retain.
mips-irix
Success	cm:	cmiss_test.log.retain.
Success	cm-debug:	cmiss_test.log.retain.
Success	cm-debug-clear-malloc:	cmiss_test.log.retain.
Success	cm-debug-clear-malloc7:	cmiss_test.log.retain.
Success	cm64:	cmiss_test.log.retain.
Success	cm64-debug:	cmiss_test.log.retain.
Success	cm64-debug-clear-malloc:	cmiss_test.log.retain.
rs6000-aix
Success	cm:	cmiss_test.log.retain.
Success	cm-debug:	cmiss_test.log.retain.
Success	cm64:	cmiss_test.log.retain.
Success	cm64-debug:	cmiss_test.log.retain.
x86_64-linux
Success	cm:	cmiss_test.log.retain.
Success	cm-debug:	cmiss_test.log.retain.

cmiss_test.out

i686-linux
Success	cm:	ndiff test: no significant differences with generic answer.
Success	cm-debug:	ndiff test: no significant differences with generic answer.
mips-irix
Success	cm:	ndiff test: no significant differences with generic answer.
Success	cm-debug:	ndiff test: no significant differences with generic answer.
Success	cm-debug-clear-malloc:	ndiff test: no significant differences with generic answer.
Success	cm-debug-clear-malloc7:	ndiff test: no significant differences with generic answer.
Success	cm64:	ndiff test: no significant differences with generic answer.
Success	cm64-debug:	ndiff test: no significant differences with generic answer.
Success	cm64-debug-clear-malloc:	ndiff test: no significant differences with generic answer.
rs6000-aix
Success	cm:	ndiff test: no significant differences with generic answer.
Success	cm-debug:	ndiff test: no significant differences with generic answer.
Success	cm64:	ndiff test: no significant differences with generic answer.
Success	cm64-debug:	ndiff test: no significant differences with generic answer.
x86_64-linux
Success	cm:	ndiff test: no significant differences with generic answer.
Success	cm-debug:	ndiff test: no significant differences with generic answer.

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

Input last modified: Fri Apr 13 10:39:17 2007

CMISS Help / Examples / b / b2 / b21 / b213 / b213b