Example 345: Curved 33x33 mesh f=0 k11=1 k22=0.1


The comfile run by this example is as follows:

#Example_3415 Multigrid colloc soln of div(kGradu)=d(u)/dt
#             on 16 elements using cubic_linear basis
 
fem de para;r;mg_curved16;example

fem de node;r;mg_curved16;example  #25 nodes
fem de base;r;;example             #cubic-lin + 9*9 pt extended basis
fem de elem;r;;example             #16 element
fem de fibr;r;;example
fem de elem;r;;example fibre
fem de grid;r;;example
fem update grid geom
fem update grid metric level 1
fem define equation;r;;example     #div(k.grad(u))=f
fem de window
fem dr lines
fem group grid external as bdry
fem define material;r;;example     #f=0, k11=1, k22=0.1
fem update grid material
fem de init;r;;example             #bdry pts =1
fem de solv;r;;example
fem dr grid level 5            #multigrid level 5 (coarsest grid)
fem dr grid level 4            #multigrid level 4
fem dr grid level 3            #multigrid level 3
fem dr grid level 2            #multigrid level 2
fem dr grid level 1            #multigrid level 1 (finest grid)
fem dr grid values             #shows initial field
sub solve
{
fem up grid time
fem solve
fem li grid YQ grid 719
fem dr grid values
}
solve                          #use anisotropic, dynamic, dt=.01
fem dr grid values             #shows computed field
fem dr scale

Additional testing commands:

fem list grid grid 1..100
fem list init
fem list solve


Files used by this example are:

Name                Modified     Size

example_345.com 20-Nov-2001 1.2k active.ipinit 10-Apr-2000 473 active.ipsolv 13-Apr-2007 1.1k mg_curved16.ipbase 10-Apr-2000 2.8k mg_curved16.ipelem 10-Apr-2000 5.0k mg_curved16.ipelfb 10-Apr-2000 2.3k mg_curved16.ipequa 26-May-2003 1.1k mg_curved16.ipfibr 30-Jan-2001 3.9k mg_curved16.ipgrid 06-Mar-2003 1.8k mg_curved16.ipinit 10-Apr-2000 376 mg_curved16.ipmate 10-Apr-2000 839 mg_curved16.ipnode 10-Apr-2000 7.6k mg_curved16.ippara 12-Nov-2002 5.9k mg_curved16.ipsolv 13-Apr-2007 1.2k test_output.com 20-Nov-2001 56

Download the entire example:

Name                      Modified     Size

examples_3_34_345.tar.gz 14-Apr-2007 17k

Testing status by version:

StatusTestedReal time (s)
i686-linux
cmSuccessSun Mar 6 00:01:53 20160
cm-debugSuccessSat Mar 5 00:02:09 20160
mips-irix
cmSuccessSun Aug 19 01:30:03 20074
cm-debugSuccessWed Aug 15 01:26:06 200710
cm-debug-clear-mallocSuccessSat Aug 18 01:34:28 200716
cm-debug-clear-malloc7SuccessMon Aug 20 01:31:51 200716
cm64SuccessSun Aug 19 01:29:06 20074
cm64-debugSuccessTue Aug 21 01:26:13 200711
cm64-debug-clear-mallocSuccessThu Apr 1 10:13:46 20046
rs6000-aix
cmSuccessWed Mar 4 01:07:11 20091
cm-debugSuccessMon Mar 2 01:07:34 20092
cm64SuccessWed Mar 4 01:07:11 20091
cm64-debugSuccessTue Mar 3 01:12:10 20092
x86_64-linux
cmSuccessSun Mar 6 00:01:02 20160
cm-debugSuccessSat Mar 5 00:01:13 20160

Testing status by file:


Html last generated: Sun Mar 6 05:50:19 2016

Input last modified: Fri Apr 13 10:32:56 2007


CMISS Help / Examples / 3 / 34 / 345