Example b2111: 2D activation using collocation

2D activation using the collocation method

Example_B2111

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:

	Status	Tested	Real time (s)
i686-linux
cm	Success	Sun Mar 6 00:02:11 2016	1
cm-debug	Success	Sat Mar 5 00:09:19 2016	2
mips-irix
cm	Success	Sun Aug 19 01:39:55 2007	11
cm-debug	Success	Wed Aug 15 01:37:20 2007	29
cm-debug-clear-malloc	Success	Sat Aug 18 01:45:40 2007	37
cm-debug-clear-malloc7	Success	Mon Aug 20 01:44:02 2007	37
cm64	Success	Sun Aug 19 01:40:07 2007	11
cm64-debug	Success	Tue Aug 21 01:36:39 2007	30
cm64-debug-clear-malloc	Success	Thu Apr 1 11:10:29 2004	15
rs6000-aix
cm	Success	Wed Mar 4 01:08:07 2009	2
cm-debug	Success	Mon Mar 2 01:09:43 2009	8
cm64	Success	Wed Mar 4 01:08:07 2009	2
cm64-debug	Success	Tue Mar 3 01:14:58 2009	9
x86_64-linux
cm	Success	Sun Mar 6 00:01:04 2016	0
cm-debug	Success	Sat Mar 5 00:01:42 2016	2

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

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

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