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:

StatusTestedReal time (s)
i686-linux
cmSuccessSun Mar 6 00:03:17 20162
cm-debugSuccessSat Mar 5 00:15:45 20166
mips-irix
cmSuccessSun Aug 19 01:57:38 200728
cm-debugSuccessWed Aug 15 02:17:10 2007115
cm-debug-clear-mallocSuccessSat Aug 18 02:32:49 2007122
cm-debug-clear-malloc7SuccessMon Aug 20 02:31:35 2007121
cm64SuccessSun Aug 19 02:01:48 200731
cm64-debugSuccessTue Aug 21 02:20:13 2007127
rs6000-aix
cmSuccessWed Mar 4 01:09:53 20093
cm-debugSuccessMon Mar 2 01:16:14 200926
cm64SuccessWed Mar 4 01:10:02 20093
cm64-debugSuccessTue Mar 3 01:20:09 200925
x86_64-linux
cmSuccessSun Mar 6 00:01:07 20161
cm-debugSuccessSat Mar 5 00:01:59 20164

Testing status by file:


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