Geometry fitting example

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:

Evenly distributed over 0 <= phi <= 90
Evenly distributed over 0 <= theta <= 360
Evenly distributed over 1.9 <= radius <= 2.1

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:

	Status	Tested	Real time (s)
i686-linux
cm	Success	Sun Mar 6 00:01:41 2016	1
cm-debug	Success	Sat Mar 5 00:12:45 2016	5
mips-irix
cm	Success	Sun Aug 19 01:47:34 2007	16
cm-debug	Success	Wed Aug 15 01:47:53 2007	52
cm-debug-clear-malloc	Success	Sat Aug 18 01:57:28 2007	56
cm-debug-clear-malloc7	Success	Mon Aug 20 01:52:43 2007	56
cm64	Success	Sun Aug 19 01:47:52 2007	16
cm64-debug	Success	Tue Aug 21 01:50:19 2007	54
cm64-debug-clear-malloc	Success	Thu Apr 1 10:12:25 2004	23
rs6000-aix
cm	Success	Wed Mar 4 01:08:50 2009	2
cm-debug	Success	Mon Mar 2 01:11:28 2009	13
cm64	Success	Wed Mar 4 01:08:50 2009	2
cm64-debug	Success	Tue Mar 3 01:16:00 2009	13
x86_64-linux
cm	Success	Sun Mar 6 00:01:01 2016	1
cm-debug	Success	Sat Mar 5 00:01:47 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.

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