Example 135: Generating a 3 generation lung mesh

This example uses the 'define mesh' command to create a 3 generation lung mesh. This is some of Merryn's boring stuff (sorry, Merryn!).

The comfile run by this example is as follows:

#Example_135  Generating a simple lung mesh

# If the example path is not set, default to current directory
if (! defined $example) {
    $example = "./";
}
# Drop off the trailing / in the example path
$chopped = chop $example;
if ($chopped ne "/") {
    $example .= $chopped;
}

$RADIUS = 1;

fem de param;r;$example/lungmesh;
fem de coor;r;$example/lungmesh;     #define 3d rectangular cartesian
fem define base;r;$example/lungmesh; #define a 1D linear Lagrange basis func

#CREATE SYMMETRIC CONDUCTING MODEL AND MULTI-BRANCHING RESPIRATORY MODEL
#     read a file that defines a symmetric airway tree (number of 
#     generations, number of elements and their length in each generation).
#     the elements are labelled as 'conducting'. 
fem de mesh;r;$example/lungmesh as conducting field;
#     add a multi-branching respiratory airway mesh to the peripheral 
#     elements in the conducting airway tree. the basis function must 
#     be specified (1d). the elements are labelled as 'respiratory'.
fem de;add mesh;c airway multi_hbw basis_function 1 as respiratory field radius_field 1 alveolar_field 2; 

#CALCULATE RADII AND ALVEOLAR CROSS-SECTIONAL AREA
#     update the conducting airway radii using anatomical data from Weibel (1963)
fem update mesh geometry elem conducting radius_field $RADIUS anatomical weibel;
#     update the respiratory airway radii using anatomical data from 
#     Haefeli-Bleuer & Weibel (1988)
fem update mesh geometry elem respiratory radius_field $RADIUS anatomical hbw;

Additional testing commands:

FEM list node;

Files used by this example are:

Name             Modified     Size

example_135.com  07-Apr-2004  1.5k
lungmesh.ipbase  10-Apr-2000  893 
lungmesh.ipcoor  10-Apr-2000  570 
lungmesh.ipmesh  31-Mar-2004  2.4k
lungmesh.ippara  31-Mar-2004  5.7k
test_output.com  31-Mar-2004   16 

Download the entire example:

Name                      Modified     Size

examples_1_13_135.tar.gz  18-Aug-2006   24k

Testing status by version:

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
  Failure	cm:	ndiff test: significant differences with generic answer.
  Failure	cm-debug:	ndiff test: 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
  Failure	cm:	ndiff test: significant differences with generic answer.
  Failure	cm-debug:	ndiff test: significant differences with generic answer.

  
  

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

Input last modified: Wed Jul 7 17:34:46 2004