Example 241: Geometry fitting example

This example tests the geometry fitting capabilities of CMISS. A half sphere initial mesh is defined using 2D cubic Hermite and Hermite simplex elements. The initial mesh has a radius of 1.8. 400 data points have been generated (via a 'c' program) about the half sphere. The data is generated according to the following scheme:

Boundary conditions on the problem are that nodes on the edge of the half sphere retain derivatives normal to the plane intersecting the sphere. The theoretical solution to this problem is a half sphere with radius 2.0. Due to the randomness of the data distribution, this will not be obtained, but the error is reduced to 0.03.

If it is necessary to view the results, run the CMISS comfile geom_sol and then the CMGUI comfile geom_draw.

Any problems or queries, see Glen Harris.


The comfile run by this example is as follows:

##Geoemtry fitting of 2D cubic Hermite and Hermite simplex elements
#to a data field defined over a half sphere.
FEM de coor;r;geometry_fit;example
FEM de bas;r;geometry_fit;example
FEM de no;r;geometry_fit;example
FEM de ele;r;geometry_fit;example
FEM de data;r;geometry_fit;example
FEM de xi;c orthog
FEM li data error
FEM export node;geom_init as geom_init
FEM export elem;geom_init as geom_init
FEM export data;geom_init_error error as geom_init_error
FEM de fit;l;geometry_fit;example opti
FEM de opti;r;geometry_fit;example
FEM eval resid wrt data
FEM list opti sum
optimise
FEM update scale_factor
#The following steps are necessary to view the results with CMGUI
FEM de no;w;geom_sol
FEM de ele;w;geom_sol

Additional testing commands:


fem write iod;sphere

FEM list node
FEM list element

Files used by this example are:

Name                    Modified     Size

example_241.com 20-Nov-2001 715 generate 10-Apr-2000 18k generate.c 10-Apr-2000 1.4k geom_draw.com 10-Apr-2000 1.2k geom_init.exelem 10-Apr-2000 7.3k geom_init.exnode 10-Apr-2000 3.2k geom_init_error.exnode 10-Apr-2000 38k geom_sol.com 10-Apr-2000 435 geom_sol.exelem 10-Apr-2000 7.1k geom_sol.exnode 10-Apr-2000 3.0k geom_sol_error.exnode 10-Apr-2000 38k geometry_fit.exelem 10-Apr-2000 16k geometry_fit.exnode 10-Apr-2000 7.1k geometry_fit.ipbase 10-Apr-2000 2.7k geometry_fit.ipcoor 10-Apr-2000 570 geometry_fit.ipdata 10-Apr-2000 23k geometry_fit.ipelem 10-Apr-2000 2.7k geometry_fit.ipfit 10-Apr-2000 4.0k geometry_fit.ipmesh 10-Apr-2000 894 geometry_fit.ipnode 10-Apr-2000 8.4k geometry_fit.ipopti 19-Mar-2002 1.8k geometry_fit.ippara 12-Nov-2002 5.0k makefile 10-Apr-2000 572 test_output.com 10-Apr-2000 54

Download the entire example:

Name                      Modified     Size

examples_2_24_241.tar.gz 19-Aug-2006 67k

Testing status by version:

StatusTestedReal time (s)
i686-linux
cmSuccessSun Mar 6 00:01:41 20161
cm-debugSuccessSat Mar 5 00:12:45 20165
mips-irix
cmSuccessSun Aug 19 01:47:34 200716
cm-debugSuccessWed Aug 15 01:47:53 200752
cm-debug-clear-mallocSuccessSat Aug 18 01:57:28 200756
cm-debug-clear-malloc7SuccessMon Aug 20 01:52:43 200756
cm64SuccessSun Aug 19 01:47:52 200716
cm64-debugSuccessTue Aug 21 01:50:19 200754
cm64-debug-clear-mallocSuccessThu Apr 1 10:12:25 200423
rs6000-aix
cmSuccessWed Mar 4 01:08:50 20092
cm-debugSuccessMon Mar 2 01:11:28 200913
cm64SuccessWed Mar 4 01:08:50 20092
cm64-debugSuccessTue Mar 3 01:16:00 200913
x86_64-linux
cmSuccessSun Mar 6 00:01:01 20161
cm-debugSuccessSat Mar 5 00:01:47 20162

Testing status by file:


Graphical output from this problem is given here.


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

Input last modified: Tue Nov 12 12:05:20 2002


CMISS Help / Examples / 2 / 24 / 241