Bicubic Signal Fitting with No Cross Derivatives

Example e83: Bicubic Signal Fitting with No Cross Derivatives

Fitting descrete field to a bicubic field with cross derivatives set to 0.

The comfile run by this example is as follows:

#
# Example e83 - Signal fitting with no cross derivatives
#

$FORMAT = "binary" #set a fileformat

fem def para;r;full_pig;example         #define parameters
fem realloc
fem def coor 3,1                        #rc coordinates
fem def regi;r;full_pig;example         #one region
fem def base;r;full_pig;example         #bicubic Herm surfaces with sectors
fem def node;r;torso_cust;example       #p111 customised surface
fem def elem;r;torso_cust;example

#
#
# Export surface to CMGUI
#
fem export node;porky as porky
fem export elem;porky as porky


fem def data;r;p111_trans;example        #define data
fem def xi;r;p111;example
fem export data;p111 as elec error offset 10500
#                                        !outputs the data positions and
#                                        !error vectors denoted by error

#
#   Fit the signal 
#      - between 0.1 and 0.2 skipping every 5 time steps
#
fem def field;r;torso_surface;example      # bicubic Hermite field
fem def;add base;r;no_cross;example        # bicubic Herm surfaces with sectors
fem def elem;r;torso_surface;example field   # element fields
fem update scale_factors estimate          # need scale factors for new bases
fem def fit;r;torso_surface;example signal class 2
# link the binsig file to the current directory for fit sig
symlink $example.'green_up2.binsig', 'green_up2.binsig';
fem fit sig class 2 $FORMAT tstart 0.10 tend 0.20

#
#
#   Export the fitted (nodally based) signals to UNEMAP
#
fem eva elec;fitted hist fitted fr node elec 13..213 $FORMAT us fit cl 2
fem conv sig from bin infile fitted  outfile fitted
fem def export;r;unemap;example
fem export sig;fitted elec signal fitted

#
#
#   Export reconstructed signals (from the input data pts) 
#    for a comparison
#
fem def data;r;p111_surf;example
fem def xi;r;p111;example
fem eva elec;fitted_data hist fitted fr data $FORMAT us fit cl 2
fem conv sig from bin infile fitted_data  outfile fitted_data

#
#
# Export the nodal signals to UNEMAP
#
fem eva elec;p111 hist fitted fr node elec 13..213 $FORMAT us fit cl 2
fem def export;r;unemap;example
fem export sig;p111 elec signal p111 

#
# Export signals to CMGUI
#

fem export node;solutions history fitted as solutions using fit cl 2 $FORMAT
fem export elem;solutions field as solutions using fit cl 2

Additional testing commands:

#
# start of testing
#

#
#
#   Comparisons between INPUT (data based) and FITTED (data based) values
#
fem comp sig mast green_up2 comp fitted_data skip 5 start 0.10 end 0.20 $FORMAT
fem comp sig mast green_up2 comp fitted_data skip 5 start 0.13 end 0.17 $FORMAT

Files used by this example are:

Name                  Modified     Size

example_e83.com       25-Nov-2005  2.2k
full_pig.ipbase       16-Nov-2000  4.5k
full_pig.ippara       12-Nov-2002  5.9k
full_pig.ipregi       10-Apr-2000  104 
green_up2.binsig      10-Apr-2000  546k
no_cross.ipbase       16-Nov-2000  4.5k
p111.ipxi             10-Apr-2000   15k
p111_surf.ipdata      10-Apr-2000   24k
p111_trans.ipdata     10-Apr-2000   27k
test_output.com       16-Nov-2000  264 
torso_cust.ipelem     10-Apr-2000   86k
torso_cust.ipnode     10-Apr-2000  226k
torso_surface.ipelfd  16-Nov-2000   36k
torso_surface.ipfiel  10-Apr-2000   53k
torso_surface.ipfit   13-Apr-2007  2.0k
unemap.ipexpo         10-Apr-2000  660

Download the entire example:

Name                      Modified     Size

examples_e_e8_e83.tar.gz  14-Apr-2007  572k

Testing status by version:

	Status	Tested	Real time (s)
i686-linux
cm	Success	Sun Mar 6 00:03:17 2016	2
cm-debug	Success	Sat Mar 5 00:15:45 2016	6
mips-irix
cm	Success	Sun Aug 19 01:57:38 2007	28
cm-debug	Success	Wed Aug 15 02:17:10 2007	115
cm-debug-clear-malloc	Success	Sat Aug 18 02:32:49 2007	122
cm-debug-clear-malloc7	Success	Mon Aug 20 02:31:35 2007	121
cm64	Success	Sun Aug 19 02:01:48 2007	31
cm64-debug	Success	Tue Aug 21 02:20:13 2007	127
rs6000-aix
cm	Success	Wed Mar 4 01:09:53 2009	3
cm-debug	Success	Mon Mar 2 01:16:14 2009	26
cm64	Success	Wed Mar 4 01:10:02 2009	3
cm64-debug	Success	Tue Mar 3 01:20:09 2009	25
x86_64-linux
cm	Success	Sun Mar 6 00:01:07 2016	1
cm-debug	Success	Sat Mar 5 00:01:59 2016	4

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.
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.
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.

solutions_1001.exnode

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.
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:28 2016

Input last modified: Fri Apr 13 10:42:08 2007

CMISS Help / Examples / e / e8 / e83