Example 231: Fit a field over two elements using orthogonal data projections

The comfile run by this example is as follows:

# Example_231 Fit a field over two elements using orthogonal data projections
fem                            #set the FEM environment
fem define node;r;example_231;example              #reads 4 nodal coordinates                                                                           
fem define base;r              #define bilinear basis                                                                               
fem define element;r           #reads 1 element                                                                                     
fem define window;c            #calculates plotting window coordinates                                                              
fem define data;r geometry     #reads 14 geometric data points from data                                                            
#                               file                                                                                                
fem def xi;c orthog	       #calculates xi projections onto linear elements
fem list data                  #shows data with weights and xi positions                                                            
fem define field;r             #reads a field for fitting (initially zero)                                                          
fem define elem;r field
fem define fit;r
fem list fit                   #shows fitting dof's                                                                                 
fem fit                #fits field coordinate increments to data                                                            
fem update node in 1           #transfers field increments to the geometric                                                         
#                               variable                                                                                            
fem define line;c              #to recalculate line information with new                                                            
#                               basis                                                                                               
fem cancel xi                  #to cancel current xi locations...                                                                   
fem define data;c xi orthog    #...and calculate new ones                                                                           
fem draw line                  #draw line segments                                                                                  
fem draw data                  #draw data segments                                                                                  
fem draw data proj             #draw data projections                                                                               
fem quit                                                                                                                            

Files used by this example are:

Name                Modified     Size

example_231.com 10-Apr-2000 2.9k example_231.ipbase 10-Apr-2000 1.1k example_231.ipcoor 10-Apr-2000 858 example_231.ipdata 10-Apr-2000 3.3k example_231.ipelem 10-Apr-2000 635 example_231.ipelfd 10-Apr-2000 368 example_231.ipfiel 10-Apr-2000 771 example_231.ipfit 13-Apr-2007 1.5k example_231.ipnode 10-Apr-2000 1.3k example_231.ipwind 10-Apr-2000 532 test.com 10-Apr-2000 2.9k test.ipbase 10-Apr-2000 1.5k test.ipcoor 10-Apr-2000 570 test.ipelem 10-Apr-2000 713 test.ipnode 10-Apr-2000 5.7k

Download the entire example:

Name                      Modified     Size

examples_2_23_231.tar.gz 14-Apr-2007 4.9k

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

Input last modified: Fri Apr 13 10:31:16 2007

CMISS Help / Examples / 2 / 23 / 231