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:

StatusTestedReal time (s)
i686-linux
cmSuccessSun Mar 6 00:03:07 20162
cm-debugSuccessSat Mar 5 00:12:32 20164
mips-irix
cmSuccessSun Aug 19 01:55:21 200725
cm-debugSuccessWed Aug 15 02:07:38 2007104
cm-debug-clear-mallocSuccessSat Aug 18 02:27:27 2007116
cm-debug-clear-malloc7SuccessMon Aug 20 02:29:27 2007114
cm64SuccessSun Aug 19 01:54:43 200725
cm64-debugSuccessTue Aug 21 02:11:47 2007110
cm64-debug-clear-mallocSuccessThu Apr 1 11:17:39 200454
rs6000-aix
cmSuccessWed Mar 4 01:09:49 20092
cm-debugSuccessMon Mar 2 01:15:24 200920
cm64SuccessWed Mar 4 01:09:49 20092
cm64-debugSuccessTue Mar 3 01:18:29 200919
x86_64-linux
cmSuccessSun Mar 6 00:01:05 20160
cm-debugSuccessSat Mar 5 00:01:57 20163

Testing status by file:


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