This example solves an activation problem using the Beeler-Reuter model on a "knife" type element. There are two ways of ordering the collapse in this type of element, which reverse the order of xi collapsing. See b21527 to compare. The order used in this example does not work for activation and there is a floating point error due to numerical instability at time 3.7
# This is the example of a knife type element that gives incorrect
# result.
# Set up for bidomain.
$VM=1
$PHI_E=2
fem define parameters;r;knife
fem define coordinates;r;knife
fem define node;r;knife
fem define base;r;knife
fem define elem;r;knife
fem define grid;r;knife
fem update grid geometry
fem update grid metric
fem ex data;normals normals as normals
fem group grid external as boundary
fem group grid xi2 high as stim_point oneoff
fem def equa;r;knife class $VM,$PHI_E
fem def mate;r;knife class $VM,$PHI_E
fem def cell;r;knife class $VM,$PHI_E
fem update grid material class $VM,$PHI_E
fem def init;r;knife class $VM,$PHI_E
fem def solve;r;knife class $VM
fem def solve;r;knife class $PHI_E
fem solve to 0 class $VM,$PHI_E
# Export the grid elements and nodes
fem export node;Domain as Domain reg 1
fem export elem;Domain grid_numbers as Domain reg 1
fem export elem;"MemPot_000" field as Solution class $VM
for $M ( 0..10 )
{
for $L ( 0..9 )
{
#Solve for new time
print "Time step: $M.$L\n";
fem solve restart to "${M}.$L" class $VM,$PHI_E
# Export solution
fem export elem;"MemPot_${M}.$L" field as Solution class $VM
}
}
Name Modified Size
example_b21528.com 29-Apr-2003 1.2k Domain.exelem 29-Apr-2003 74k Domain.exnode 29-Apr-2003 1.4k MemPot_0.0.exelem 29-Apr-2003 95k MemPot_0.1.exelem 29-Apr-2003 95k MemPot_0.2.exelem 29-Apr-2003 95k MemPot_0.3.exelem 29-Apr-2003 95k MemPot_0.4.exelem 29-Apr-2003 95k MemPot_0.5.exelem 29-Apr-2003 95k MemPot_0.6.exelem 29-Apr-2003 95k MemPot_0.7.exelem 29-Apr-2003 95k MemPot_0.8.exelem 29-Apr-2003 95k MemPot_0.9.exelem 29-Apr-2003 95k MemPot_000.exelem 29-Apr-2003 95k MemPot_1.0.exelem 29-Apr-2003 95k MemPot_1.1.exelem 29-Apr-2003 95k MemPot_1.2.exelem 29-Apr-2003 95k MemPot_1.3.exelem 29-Apr-2003 95k MemPot_1.4.exelem 29-Apr-2003 95k MemPot_1.5.exelem 29-Apr-2003 95k MemPot_1.6.exelem 29-Apr-2003 95k MemPot_1.7.exelem 29-Apr-2003 95k MemPot_1.8.exelem 29-Apr-2003 95k MemPot_1.9.exelem 29-Apr-2003 95k MemPot_2.0.exelem 29-Apr-2003 95k MemPot_2.1.exelem 29-Apr-2003 95k MemPot_2.2.exelem 29-Apr-2003 95k MemPot_2.3.exelem 29-Apr-2003 95k MemPot_2.4.exelem 29-Apr-2003 95k MemPot_2.5.exelem 29-Apr-2003 95k MemPot_2.6.exelem 29-Apr-2003 95k MemPot_2.7.exelem 29-Apr-2003 95k MemPot_2.8.exelem 29-Apr-2003 95k MemPot_2.9.exelem 29-Apr-2003 95k MemPot_3.0.exelem 29-Apr-2003 95k MemPot_3.1.exelem 29-Apr-2003 95k MemPot_3.2.exelem 29-Apr-2003 95k MemPot_3.3.exelem 29-Apr-2003 95k MemPot_3.4.exelem 29-Apr-2003 95k MemPot_3.5.exelem 29-Apr-2003 95k MemPot_3.6.exelem 29-Apr-2003 95k knife.ipbase 29-Apr-2003 2.4k knife.ipcell 29-Apr-2003 2.0k knife.ipcoor 29-Apr-2003 570 knife.ipelem 29-Apr-2003 721 knife.ipequa 26-May-2003 1.6k knife.ipgrid 29-Apr-2003 606 knife.ipinit 29-Apr-2003 903 knife.ipmate 29-Apr-2003 1.8k knife.ipnode 29-Apr-2003 2.4k knife.ippara 29-Apr-2003 5.9k knife.ipsolv 13-Apr-2007 1.8k normals.exdata 29-Apr-2003 190k
Name Modified Size
examples_b_b2_b21_b215_b2152_b21528.tar.gz 14-Apr-2007 401k
Html last generated: Sun Mar 6 05:51:15 2016
Input last modified: Fri Apr 13 10:39:55 2007