Example b212d: Implicit 2D Cubic Ionic Activation (bidomain) using Grid-based FE


This example uses the Grid-based Finite Element Method. It solves the Cubic Ionic Activation model, using the bidomain equations over a two dimensional slice with fibres.

After running this example, the output from this problem may be viewed by running draw.com.


The comfile run by this example is as follows:

#Example b211t
# ------------
fem define para;r;box;example
fem define coor;r;box;example

fem define node;r;box;example
fem define base;r;box;example
fem define elem;r;box;example

fem define fibr;r;box;example
fem define elem;r;box;example fibre

fem define grid;r;box;example
fem update grid geometry
fem group grid external as boundary
fem group grid xi1=low as stimulus

fem define equa;r;box;example class 1,2
fem define mate;r;box;example class 1,2
fem define cell;r;box;example class 1,2
fem define init;r;box;example class 1,2
fem define solve;r;box;example class 1
fem define solve;r;box;example class 2

# cmgui output
fem export node;box as box
fem export elem;box grid_numbers as box

fem export elem;"gridfield_000" field

# loop through time from 0 to 10 msec
foreach $count ( '001'..'020' ) 
{
 my $time = $count / 2;
 if( $count == 1 ) {
   fem solve to $time class 1,2
 } else {
   fem solve to "$time" class 1,2 restart
 }
 print "time: $time\n";
 #fem export elem;"gridfield_$count" field class 1
}

#See draw.com to view results

Additional testing commands:

#Testing example b212d
fem li grid yq index v

Files used by this example are:

Name               Modified     Size

example_b212d.com 13-Apr-2007 1.0k box.exelem 07-Feb-2001 3.4k box.exnode 07-Feb-2001 896 box.ipbase 07-Feb-2001 1.1k box.ipcell 07-Feb-2001 1.7k box.ipcoor 07-Feb-2001 858 box.ipelem 07-Feb-2001 584 box.ipelfb 07-Feb-2001 358 box.ipequa 26-May-2003 1.5k box.ipfibr 07-Feb-2001 828 box.ipgrid 06-Mar-2003 537 box.ipinit 07-Feb-2001 852 box.ipmate 07-Feb-2001 1.7k box.ipnode 07-Feb-2001 1.0k box.ippara 12-Nov-2002 5.9k box.ipsolv 13-Apr-2007 2.1k box.ipsolv.old 05-Jun-2001 1.6k draw.com 09-Nov-2001 4.2k test_output.com 23-Apr-2001 46

Download the entire example:

Name                                 Modified     Size

examples_b_b2_b21_b212_b212d.tar.gz 14-Apr-2007 10k

Testing status by version:

StatusTestedReal time (s)
i686-linux
cmSuccessSun Mar 6 00:02:20 20160
cm-debugSuccessSat Mar 5 00:04:18 20160
mips-irix
cmSuccessSun Aug 19 01:28:39 20073
cm-debugSuccessWed Aug 15 01:25:54 20078
cm-debug-clear-mallocSuccessSat Aug 18 01:28:08 200710
cm-debug-clear-malloc7SuccessMon Aug 20 01:25:13 200711
cm64SuccessSun Aug 19 01:29:27 20073
cm64-debugSuccessTue Aug 21 01:23:26 20078
cm64-debug-clear-mallocSuccessThu Apr 1 11:18:13 20044
rs6000-aix
cmSuccessWed Mar 4 01:07:40 20091
cm-debugSuccessMon Mar 2 01:07:00 20091
cm64SuccessWed Mar 4 01:07:40 20091
cm64-debugSuccessTue Mar 3 01:12:16 20091
x86_64-linux
cmSuccessSun Mar 6 00:01:05 20160
cm-debugSuccessSat Mar 5 00:01:35 20160

Testing status by file:


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

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


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