Example b2111: 2D activation using collocation

2D activation using the collocation method


In this example, the collocation method is used to solve for 2D activation. A 9x9 collocation point basis is used, and on a 16x16 element mesh, this totals 16641 collocation points. A constant zero degree fibre direction is defined, and a 2:1 anisotropic conductivity is used.

A non-reetrant wave is generated using initial conditions specified, and after solving the problem, a field is fitted to the computed activation times. WARNING: This may take several hours to solve on a fast machine!

The comfile run by this example is as follows:

#Example_b2111 2D activation using the collocation method
#  See Description.html for further explanation.
fem de param;r;example_b1;example           #Parameter settings

fem de bas;r;colloc-9x9;example 		#9x9 grid pt collocation basis
fem de node;r;square16x16;example 		#16x16 element square mesh
fem de elem;r;square16x16;example 		#Total 128mmx128mm
fem de fibr;r;sq16x16-00;example 		#Constant 0 degree fibre angle
fem de elem;r;;example fibre
fem de grid;r;9x9;example 				#Create 16641 grid pts
fem update grid geom
fem update grid metric
fem de equ;r;vcdc-is;example lock 		#Modified VCD model w/ ischemic
#					!  behaviour (permits reentry)
#  !The next lines define groups of collocation pts which may be used
#  !  as initial stimulus points.
fem group grid grid 7872,7873,7864,7865,7922,7930,8848,8849,8906 as center
#fem group grid xi1=0 oneoff as lhs-line
#fem group grid xi1=1 oneoff as rhs-line
#fem group grid xi2=0 oneoff as bottom-line
#fem group grid xi2=1 oneoff as top-line
#fem group grid grid 7873,7865,7857,7849,7841,7833,7825,7817,6841,6833,6825,6817,6809,6801,6793,6785 as cline1
#fem group grid grid 5809,5801,5793,5785,5777,5769,5761,5753,4777,4769,4761,4753,4745,4737,4729,4721 as cline2
#fem group grid grid 3745,3737,3729,3721,3713,3705,3697,3689,2713,2705,2697,2689,2681,2673,2665,2657 as cline3
#fem group grid grid 1681,1673,1665,1657,1649,1641,1633,1625,585,577,569,561,553,545,537,529 as cline4
fem de mat;r;vcd-21;example 		#2:1 conductivity ratio
fem de cell;r;vcd-21;example 		#2:1 conductivity ratio
fem up grid material 			#Update local conductivities
fem de ini;r;sq16x16-center;example 	#Stimulus in center of square
fem de sol;r;sq16x16;example 		#Solution parameters
#This section commented to allow example to run overnight.
#fem solve 				!Run solution
#fem up gauss activation			!Store activation times at gauss
#					!  points
#fem de win 				!Create output window
fem de field;r;sq16;example 		#Read field info
fem de elem;r;;example field   #Read field element info
fem de fit;r;sq16;example field gauss lock 	#And fitting parameters
#fem fit 				!Fit field from gauss points
#fem dr field field 1  			!Draw field
#fem dr cont field label off value 0,20,40,60,80,100,120,140,160,180,200
#					!Draw contour lines
#fem dr axes x=32,64,96 y=32,64,96 	!Draw axes
#fem dr fibr dxi 1,1,0.4 rgb=black 	!Draw fibre angles

Additional testing commands:

fem list base
fem list grid yq index v grid 200..250
fem list fit
fem list solve

Files used by this example are:

Name                        Modified     Size

example_b2111.com 20-Nov-2001 2.3k 9x9.ipgrid 06-Mar-2003 22k bak.ippara 12-Nov-2002 5.9k bak.ippara_old 10-Apr-2000 4.3k colloc-9x9.ipbase 10-Apr-2000 1.1k example_b1.ippara 12-Nov-2002 5.9k example_b1.ippara_old 10-Apr-2000 4.4k example_b1_full.ippara 12-Nov-2002 5.8k example_b1_full.ippara_old 10-Apr-2000 4.3k sq16.ipelfd 10-Apr-2000 36k sq16.ipfiel 10-Apr-2000 25k sq16.ipfit 13-Apr-2007 1.6k sq16x16-00.ipelfb 10-Apr-2000 35k sq16x16-00.ipfibr 30-Jan-2001 22k sq16x16-center.ipinit 10-Apr-2000 487 sq16x16.ipsolv 03-Oct-2001 987 sq16x16.ipsolv.old 10-Apr-2000 941 square16x16.ipelem 10-Apr-2000 65k square16x16.ipnode 10-Apr-2000 44k test_output.com 23-Apr-2001 82 vcd-21.ipcell 10-Apr-2000 2.3k vcd-21.ipmate 10-Apr-2000 1.7k vcd-21.ipmate_cell_old 10-Apr-2000 4.0k vcdc-is.ipequa 26-May-2003 1.5k

Download the entire example:

Name                                 Modified     Size

examples_b_b2_b21_b211_b2111.tar.gz 14-Apr-2007 31k

Testing status by version:

StatusTestedReal time (s)
cmSuccessSun Mar 6 00:02:11 20161
cm-debugSuccessSat Mar 5 00:09:19 20162
cmSuccessSun Aug 19 01:39:55 200711
cm-debugSuccessWed Aug 15 01:37:20 200729
cm-debug-clear-mallocSuccessSat Aug 18 01:45:40 200737
cm-debug-clear-malloc7SuccessMon Aug 20 01:44:02 200737
cm64SuccessSun Aug 19 01:40:07 200711
cm64-debugSuccessTue Aug 21 01:36:39 200730
cm64-debug-clear-mallocSuccessThu Apr 1 11:10:29 200415
cmSuccessWed Mar 4 01:08:07 20092
cm-debugSuccessMon Mar 2 01:09:43 20098
cm64SuccessWed Mar 4 01:08:07 20092
cm64-debugSuccessTue Mar 3 01:14:58 20099
cmSuccessSun Mar 6 00:01:04 20160
cm-debugSuccessSat Mar 5 00:01:42 20162

Testing status by file:

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

Input last modified: Fri Apr 13 10:37:40 2007

CMISS Help / Examples / b / b2 / b21 / b211 / b2111