2d circles with equivalent dipole sources

Example b2125: 2d circles with equivalent dipole sources

Example b2125 - 2d circles with equivalent dipole sources This example solves the VCD ionic equtions over a finite difference grid generated from an annulus of finite elements. An equivalent dipole source is calculated in each of the finite elements and these sources are used as the stimulus for a coupled fem-bem generalised Laplace problem.

The comfile run by this example is as follows:

$OUTBEM = 1
$FEMBEM = 2
$CENBEM = 3
$OUTFEM = 4
$GRDFEM = 5

$LAPLACE = 1
$TMEMBRN = 2
$EXTCLLR = 3


fem define para;r;fembem;example

fem define regi;r;fembem;example
fem define node;r;test;example reg all
fem define base;r;fembem;example
fem define elem;r;test;example reg all
fem define fibr;r;test;example reg $GRDFEM
fem define elem;r;test;example fibre reg $GRDFEM

fem define grid;r;heart;example reg $GRDFEM
fem update grid geometry reg $GRDFEM
fem update grid metric reg $GRDFEM
fem group grid external as boundary reg $GRDFEM

fem define equa;r;bem;example reg $OUTBEM,$FEMBEM,$CENBEM,$OUTFEM class $LAPLACE lock
fem define equa;r;heart;example reg $GRDFEM class $TMEMBRN,$EXTCLLR
fem define mate;r;bem;example reg $OUTBEM,$FEMBEM,$CENBEM,$OUTFEM class $LAPLACE
fem define mate;r;heart;example reg $GRDFEM class $TMEMBRN,$EXTCLLR
fem define cell;r;heart;example reg $GRDFEM class $TMEMBRN,$EXTCLLR
fem update grid material reg $GRDFEM class $TMEMBRN,$EXTCLLR
fem define init;r;bem;example reg $OUTBEM,$FEMBEM,$CENBEM,$OUTFEM class $LAPLACE
fem define init;r;heart;example reg $GRDFEM class $TMEMBRN,$EXTCLLR
fem define coup;r;bem;example class $LAPLACE
fem define solv;r;heart;example reg $GRDFEM class $TMEMBRN
fem define solv;r;heart2;example reg $GRDFEM class $EXTCLLR
fem define solv;r;bem;example coupled reg $OUTBEM,$FEMBEM,$CENBEM,$OUTFEM class $LAPLACE

fem solve to 0 reg $GRDFEM class $TMEMBRN,$EXTCLLR
fem define sour;c grid diregion $FEMBEM grregion $GRDFEM time 0
fem update source reg $FEMBEM class $LAPLACE
fem solve coupled reg $OUTBEM,$FEMBEM,$CENBEM,$OUTFEM class $LAPLACE

fem solve restart to 5 reg $GRDFEM class $TMEMBRN,$EXTCLLR
fem define sour;c grid diregion $FEMBEM grregion $GRDFEM time 5
fem update source reg $FEMBEM class $LAPLACE
fem solve coupled reg $OUTBEM,$FEMBEM,$CENBEM,$OUTFEM class $LAPLACE

fem solve restart to 10 reg $GRDFEM class $TMEMBRN,$EXTCLLR
fem define sour;c grid diregion $FEMBEM grregion $GRDFEM time 10
fem update source reg $FEMBEM class $LAPLACE
fem solve coupled reg $OUTBEM,$FEMBEM,$CENBEM,$OUTFEM class $LAPLACE

fem define wind;c
fem draw fiel grid 1 zmin=-50 zmax=2 noscale reg $GRDFEM class $EXTCLLR
fem draw fiel reg $OUTFEM noscale zmax=0.07
fem draw dipo reg $FEMBEM scal 10
fem draw line reg $GRDFEM

fem define sour;w;dipoles grid diregion $FEMBEM grregion $GRDFEM

Additional testing commands:

fem list sour reg 2
fem list node solution reg 1

Files used by this example are:

Name                   Modified     Size

example_b2125.com      12-Jun-2001  2.3k
bem.ipcoup             21-Aug-2002  437 
bem.ipcoup.old         10-Apr-2000  411 
bem.irequa             26-May-2003  4.8k
bem.irinit             10-Apr-2000  1.4k
bem.irmate             10-Apr-2000  2.3k
bem.irsolv             13-Apr-2007  1.8k
fembem.ipbase          10-Apr-2000  2.5k
fembem.ippara          14-Mar-2003  5.9k
fembem.ipregi          10-Apr-2000   93 
heart.ipcell           20-Nov-2001  2.3k
heart.ipequa           26-May-2003  1.5k
heart.ipgrid           06-Mar-2003  1.1k
heart.ipinit           20-Nov-2001  5.2k
heart.ipmate           10-Apr-2000  1.7k
heart.ipmate_cell_old  10-Apr-2000  4.0k
heart.ipsolv           13-Apr-2007  1.8k
heart.ipsolv.old       10-Apr-2000  1.5k
heart2.ipsolv          13-Apr-2007  1.8k
heart2.ipsolv.old      10-Apr-2000  1.5k
test.ipelfb            10-Apr-2000  1.2k
test.ipfibr            12-Jun-2001  2.6k
test.irelem            10-Apr-2000   14k
test.irnode            10-Apr-2000   22k
test_output.com        10-Apr-2000   49

Download the entire example:

Name                                 Modified     Size

examples_b_b2_b21_b212_b2125.tar.gz  14-Apr-2007   13k

Testing status by version:

	Status	Tested	Real time (s)
i686-linux
cm	Success	Sun Mar 6 00:03:07 2016	2
cm-debug	Success	Sat Mar 5 00:12:32 2016	4
mips-irix
cm	Success	Sun Aug 19 01:55:21 2007	25
cm-debug	Success	Wed Aug 15 02:07:38 2007	104
cm-debug-clear-malloc	Success	Sat Aug 18 02:27:27 2007	116
cm-debug-clear-malloc7	Success	Mon Aug 20 02:29:27 2007	114
cm64	Success	Sun Aug 19 01:54:43 2007	25
cm64-debug	Success	Tue Aug 21 02:11:47 2007	110
cm64-debug-clear-malloc	Success	Thu Apr 1 11:17:39 2004	54
rs6000-aix
cm	Success	Wed Mar 4 01:09:49 2009	2
cm-debug	Success	Mon Mar 2 01:15:24 2009	20
cm64	Success	Wed Mar 4 01:09:49 2009	2
cm64-debug	Success	Tue Mar 3 01:18:29 2009	19
x86_64-linux
cm	Success	Sun Mar 6 00:01:05 2016	0
cm-debug	Success	Sat Mar 5 00:01:57 2016	3

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:10 2016

Input last modified: Fri Apr 13 10:38:32 2007

CMISS Help / Examples / b / b2 / b21 / b212 / b2125