Added selected test models to be used as regression tests during development

git-svn-id: http://svn.sintef.no/trondheim/IFEM/trunk@903 e10b68d5-8a6e-419e-a041-bce267b0401d

Apps/FiniteDefElasticity/Test/CanTS-p1.inp

@@ -0,0 +1,31 @@
+# $Id$
+# Cantilever beam with a tip shear load.
+# Isotropic linear elastic material.
+# Nonlinear 3D analysis with linear elements.
+
+PATCHES 1
+slenderBeam.g2
+
+REFINE 1
+# patch ru rv rw
+  1     39 0  0
+
+CONSTRAINTS 1
+# patch face code
+  1     1    123
+
+PRESSURE 1
+# patch face dir pval     linear
+  1     2    3   18223.95 1.0
+
+ISOTROPIC 1
+# code E     nu  rho
+  0    1.0e8 0.0 0.0
+
+TIME_STEPPING
+# start stop dt
+  0.0   1.0  0.1
+
+NONLINEAR_SOLVER
+# maxit convTol divgLim
+  20    1.0e-16 1.0e5

Apps/FiniteDefElasticity/Test/CanTS-p2.inp

@@ -0,0 +1,35 @@
+# $Id$
+# Cantilever beam with a tip shear load.
+# Isotropic linear elastic material.
+# Nonlinear 3D analysis with quadratic elements.
+
+PATCHES 1
+slenderBeam.g2
+
+RAISEORDER 1
+# patch ru rv rw
+  1     1  1  1
+
+REFINE 1
+# patch ru rv rw
+  1     9  0  0
+
+CONSTRAINTS 1
+# patch face code
+  1     1    123
+
+PRESSURE 1
+# patch face dir pval     linear
+  1     2    3   18223.95 1.0
+
+ISOTROPIC 1
+# code E     nu  rho
+  0    1.0e8 0.0 0.0
+
+TIME_STEPPING
+# start stop dt
+  0.0   1.0  0.1
+
+NONLINEAR_SOLVER
+# maxit convTol divgLim
+  20    1.0e-16 1.0e5

Apps/FiniteDefElasticity/Test/CanTS-p3.inp

@@ -0,0 +1,35 @@
+# $Id$
+# Cantilever beam with a tip shear load.
+# Isotropic linear elastic material.
+# Nonlinear 3D analysis with cubic elements.
+
+PATCHES 1
+slenderBeam.g2
+
+RAISEORDER 1
+# patch ru rv rw
+  1     2  2  2 
+
+REFINE 1
+# patch ru rv rw
+  1     1  0  0
+
+CONSTRAINTS 1
+# patch face code
+  1     1    123
+
+PRESSURE 1
+# patch face dir pval     linear
+  1     2    3   18223.95 1.0
+
+ISOTROPIC 1
+# code E     nu  rho
+  0    1.0e8 0.0 0.0
+
+TIME_STEPPING
+# start stop dt
+  0.0   1.0  0.1
+
+NONLINEAR_SOLVER
+# maxit convTol divgLim
+  20    1.0e-16 1.0e5

Apps/FiniteDefElasticity/Test/FBlock-h8x2-Q4P3.inp

@@ -0,0 +1 @@
+FBlock-h8x2-p4.inp

Apps/FiniteDefElasticity/Test/FBlock-h8x2-Q4Q3.inp

@@ -0,0 +1 @@
+FBlock-h8x2-p4.inp

Apps/FiniteDefElasticity/Test/FBlock-h8x2-p4.inp

@@ -0,0 +1,36 @@
+# $Id$
+# Compression of a 2D rectangular rubber block.
+# Isotropic nearly-incompressible hyperelastic material.
+# Nonlinear analysis with quartic elements.
+
+PATCHFILE block.g2
+
+REFINE 1
+# patch ru rv
+  1     7  1
+
+RAISEORDER 1
+# patch ru rv
+  1     3  3
+
+CONSTRAINTS 3
+# patch edge code
+  1     1    1
+  1     2    1
+  1     3    2
+
+PRESSURE 1
+# patch edge dir pval      x0  x1  linear
+  1     4    2  -2.0 StepX 1.0 3.0 1.0
+
+ISOTROPIC 1
+# code kappa    mu      rho mVER
+  0    400942.0 80.1938 0.0 11
+
+TIME_STEPPING
+# start stop  dt
+  0.0   200.0 1.0
+
+NONLINEAR_SOLVER
+# maxit convTol
+  10    1.0e-16

Apps/FiniteDefElasticity/Test/FBlock-h8x3-Q3P2.inp

@@ -0,0 +1 @@
+FBlock-h8x3-p3.inp

Apps/FiniteDefElasticity/Test/FBlock-h8x3-Q3Q2.inp

@@ -0,0 +1 @@
+FBlock-h8x3-p3.inp

Apps/FiniteDefElasticity/Test/FBlock-h8x3-p3.inp

@@ -0,0 +1,36 @@
+# $Id$
+# Compression of a 2D rectangular rubber block.
+# Isotropic nearly-incompressible hyperelastic material.
+# Nonlinear analysis with cubic elements.
+
+PATCHFILE block.g2
+
+REFINE 1
+# patch ru rv
+  1     7  2
+
+RAISEORDER 1
+# patch ru rv
+  1     2  2
+
+CONSTRAINTS 3
+# patch edge code
+  1     1    1
+  1     2    1
+  1     3    2
+
+PRESSURE 1
+# patch edge dir pval      x0  x1  linear
+  1     4    2  -2.0 StepX 1.0 3.0 1.0
+
+ISOTROPIC 1
+# code kappa    mu      rho mVER
+  0    400942.0 80.1938 0.0 11
+
+TIME_STEPPING
+# start stop  dt
+  0.0   200.0 1.0
+
+NONLINEAR_SOLVER
+# maxit convTol
+  10    1.0e-16

Apps/FiniteDefElasticity/Test/FBlock-h9x5-Q2P1.inp

@@ -0,0 +1 @@
+FBlock-h9x5-p2.inp

Apps/FiniteDefElasticity/Test/FBlock-h9x5-Q2Q1.inp

@@ -0,0 +1 @@
+FBlock-h9x5-p2.inp

Apps/FiniteDefElasticity/Test/FBlock-h9x5-p2.inp

@@ -0,0 +1,38 @@
+# $Id$
+# Compression of a 2D rectangular rubber block.
+# Isotropic nearly-incompressible hyperelastic material.
+# Nonlinear analysis with quadratic elements.
+
+PATCHFILE block.g2
+
+REFINE 2
+# patch dir xi ...
+  1     1   0.125 0.25 0.35 0.45 0.55 0.65 0.75 0.875
+# patch ru rv
+  1     0  4
+
+RAISEORDER 1
+# patch ru rv
+  1     1  1
+
+CONSTRAINTS 3
+# patch edge code
+  1     1    1
+  1     2    1
+  1     3    2
+
+PRESSURE 1
+# patch edge dir pval      x0  x1  linear
+  1     4    2  -2.0 StepX 1.0 3.0 1.0
+
+ISOTROPIC 1
+# code kappa    mu      rho mVER
+  0    400942.0 80.1938 0.0 11
+
+TIME_STEPPING
+# start stop  dt
+  0.0   200.0 1.0
+
+NONLINEAR_SOLVER
+# maxit convTol
+  10    1.0e-16

Apps/FiniteDefElasticity/Test/block.g2

@@ -0,0 +1,10 @@
+200 1 0 0
+2 0
+2 2
+0 0 1 1
+2 2
+0 0 1 1
+0.0 0.0
+4.0 0.0
+0.0 1.0
+4.0 1.0

Apps/FiniteDefElasticity/Test/run.sh

@@ -0,0 +1,62 @@
+#!/bin/sh
+# $Id$
+##################################################################
+# This script defines a set of sample simulation cases that
+# is used to verify the integrity of the IFEM-based simulator.
+# Copy this file to a sub-folder Test in your App-directory and
+# insert the simulations that you want to use as regression tests.
+##################################################################
+
+# Define the name of the executable here
+mysim=NonLinEl
+
+run () {
+# This function runs a simulation with the specified options,
+# pipes the terminal output to a log-file, and compares it
+# with a previous simulation stored in the Reference folder.
+  inp=$1
+  log=`basename $1 .inp`.log
+  echo Running $inp ...
+  shift
+  time -f "real time %E, user time %U, system time %S" \
+  ../Release/bin/$mysim $inp $* > $log
+  if [ ! -e Reference/$log ]; then
+    mv $log Reference
+  elif cmp -s $log Reference/$log; then
+    echo Ok
+  else
+    echo "Warning: Discrepancies between current run and reference."
+    diff $log Reference
+  fi
+}
+
+if [ ! -d Reference ]; then mkdir Reference; fi
+
+#####################################
+# Enter the various cases below:
+# Format: run <inputfile> [<options>]
+#####################################
+
+# 3D cantilever beam, linear-elastic material
+for input in CanTS-p?.inp; do
+  ln -s $input TL-$input
+  run TL-$input -dense -vtf 1
+  ln -s $input UL-$input
+  run UL-$input -UL -dense -vtf 1
+  rm ?L-$input
+done
+
+# 3D Thick cylinder compression, linear-elastic material
+run Cyl-p2.inp -UL -vtf 1 -nGauss 3 -nviz 3
+run Cyl-p3.inp -UL -vtf 1 -nGauss 4 -nviz 3
+run Cyl-p4.inp -UL -vtf 1 -nGauss 4 -nviz 3
+
+# 2D Rubber block, hyperelastic Neo-Hooke material.
+# - mixed formulation with internal pressure modes.
+run FBlock-h9x5-Q2P1.inp -2Dpstrain -MX 1 -nGauss 3 -vtf 1 -lagrange
+run FBlock-h8x3-Q3P2.inp -2Dpstrain -MX 2 -nGauss 4 -vtf 1 -lagrange
+run FBlock-h8x2-Q4P3.inp -2Dpstrain -MX 3 -nGauss 5 -vtf 1 -lagrange
+# - mixed formulation with continuous pressure field.
+run FBlock-h9x5-Q2Q1.inp -2Dpstrain -mixed -nGauss 3 -vtf 1 -lagrange
+run FBlock-h8x3-Q3Q2.inp -2Dpstrain -mixed -nGauss 4 -vtf 1 -lagrange
+run FBlock-h8x2-Q4Q3.inp -2Dpstrain -mixed -nGauss 5 -vtf 1 -lagrange

Apps/FiniteDefElasticity/Test/slenderBeam.g2

@@ -0,0 +1,16 @@
+700 1 0 0
+3 0
+2 2
+0 0 1 1
+2 2
+0 0 1 1
+2 2
+0 0 1 1
+ 0.0 0.0    0.0
+10.0 0.0    0.0
+ 0.0 0.1478 0.0
+10.0 0.1478 0.0
+ 0.0 0.0    0.1
+10.0 0.0    0.1
+ 0.0 0.1478 0.1
+10.0 0.1478 0.1

Apps/LinearElasticity/Test/Cylinder-Lagrange.inp

@@ -0,0 +1 @@
+Cylinder-p4.inp

Apps/LinearElasticity/Test/Cylinder-NURBS.inp

@@ -0,0 +1 @@
+Cylinder-p4.inp

Apps/LinearElasticity/Test/Cylinder-Spectral.inp

@@ -0,0 +1 @@
+Cylinder-p4.inp

Apps/LinearElasticity/Test/Cylinder-p4.inp

@@ -0,0 +1,24 @@
+# $Id$
+# Thick cylinder with internal pressure.
+# Static linear-elastic analysis, quartic NURBS elements.
+
+PATCHES 1
+cylinder.g2
+
+CONSTRAINTS 4
+# patch face code
+  1     5    3
+  1     6    3
+# patch corner code
+  1     -1     123
+  1     -3      23
+
+PERIODIC 1
+# patch direction
+  1     1
+
+PRESSURE 1
+# patch face direction pressure
+  1     3    0         -1.0e6
+
+LOCAL_SYSTEM CYLINDRICZ

Apps/LinearElasticity/Test/Hole2D-Lagrange.inp

@@ -0,0 +1 @@
+Hole2D-p3.inp

Apps/LinearElasticity/Test/Hole2D-NURBS.inp

@@ -0,0 +1 @@
+Hole2D-p3.inp

Apps/LinearElasticity/Test/Hole2D-Spectral.inp

@@ -0,0 +1 @@
+Hole2D-p3.inp

Apps/LinearElasticity/Test/Hole2D-p3.inp

@@ -0,0 +1,28 @@
+# $Id$
+# Infinite membrane with a circular hole.
+# Neumann boundary conditions from analytical solution.
+# Static linear-elastic analysis, cubic NURBS elements.
+
+PATCHES 1
+hole2D.g2
+
+REFINE 1
+# patch ru rv
+  1     3  3
+
+CONSTRAINTS 2
+# patch edge code
+  1     1    1
+  1     2    2
+
+# Analytical solution
+# Specifier a   F0   nu
+ANASOL Hole 1.0 10.0 0.3
+
+PRESSURE 1
+# patch edge
+  1     4
+
+ISOTROPIC 1
+# code E      nu  rho
+  0    1000.0 0.3 0.0

Apps/LinearElasticity/Test/PipeJoint-Lagrange.inp

@@ -0,0 +1 @@
+PipeJoint.inp

Apps/LinearElasticity/Test/PipeJoint-NURBS.inp

@@ -0,0 +1 @@
+PipeJoint.inp

Apps/LinearElasticity/Test/PipeJoint-vibration.inp

@@ -0,0 +1 @@
+PipeJoint.inp

Apps/LinearElasticity/Test/PipeJoint.inp

@@ -0,0 +1,16 @@
+# $Id$
+# Pipe joint with shear-loaded brace.
+# Static linear-elastic analysis.
+# 10-patch model, cubic NURBS elements.
+
+PATCHFILE    pipe_bifurcation.g2
+NODEFILE     pipe_bifurcation.gno
+PROPERTYFILE pipe_bifurcation.prc
+
+DIRICHLET 1
+# code
+  123
+
+CONSTANT_PRESSURE 1
+# code dir pval
+  1001 1   1.0e8

Apps/LinearElasticity/Test/SquarePlate-Lagrange.inp

@@ -0,0 +1 @@
+SquarePlate-p2.inp

Apps/LinearElasticity/Test/SquarePlate-Spectral.inp

@@ -0,0 +1 @@
+SquarePlate-p2.inp

Apps/LinearElasticity/Test/SquarePlate-Splines.inp

@@ -0,0 +1 @@
+SquarePlate-p2.inp

Apps/LinearElasticity/Test/SquarePlate-p2.inp

@@ -0,0 +1,35 @@
+# $Id$
+# Square plate subjected to in-plane compression.
+# Linearized buckling analysis, quadratic splines.
+
+PATCHES 1
+plate_10x10x01.g2
+
+RAISEORDER 1
+# patch ru rv rw
+  1     1  1  1
+
+REFINE 1
+# patch ru rv rw
+  1     6  6  0
+
+CONSTRAINTS 4
+# patch line code par
+  1     11   3    0.5
+  1     12   3    0.5
+  1     23   3    0.5
+  1     24   3    0.5
+
+FIXPOINTS 2
+# patch rx  ry  rz  code
+  1     0.0 0.0 0.5 123
+  1     0.0 1.0 0.5 13
+
+PRESSURE 2
+# patch face dir pressure
+  1     1    1    10.0
+  1     2    1   -10.0
+
+ISOTROPIC 1
+# code E     nu  rho
+  0    10.92 0.3 1.0e-6

Apps/LinearElasticity/Test/hole2D.g2

@@ -0,0 +1,26 @@
+200 1 0 0
+2 1
+5 4
+0 0 0 0 1 2 2 2 2
+4 4
+0 0 0 0 1 1 1 1
+0 1 1
+0.235702 0.902369 0.902369
+0.638071 0.638071 0.804738
+0.902369 0.235702 0.902369
+1 0 1
+0 2 1
+0.471405 1.804738 0.902369
+1.276142 1.276142 0.804738
+1.804738 0.471405 0.902369
+2 0 1
+0 3 1
+0.707107 2.707107 0.902369
+1.914214 1.914214 0.804738
+2.707107 0.707107 0.902369
+3 0 1
+0 4 1
+4 4 1
+4 4 1
+4 4 1
+4 0 1

Apps/LinearElasticity/Test/pipe_bifurcation.g2

@@ -0,0 +1,362 @@
+700 1 0 0
+3 1
+4 4
+0 0 0 0 58 58 58 58
+5 3
+15.707963267949 15.707963267949 15.707963267949 31.4159265358979 31.4159265358979 47.1238898038469 47.1238898038469 47.1238898038469
+2 2
+0 0 1 1
+-58 6.12323399573677e-16 10 1
+-38.6666666666667 4.08215599715784e-16 10 1
+-19.3333333333333 2.04107799857892e-16 10 1
+0 0 10 1
+-41.0121933088198 -7.07106781186547 7.07106781186548 0.707106781186548
+-27.3414622058798 -7.07106781186548 7.07106781186548 0.707106781186548
+-13.6707311029399 -7.07106781186548 7.07106781186548 0.707106781186548
+0 -7.07106781186548 7.07106781186548 0.707106781186548
+-58 -10 1.22464679914735e-15 1
+-38.6666666666667 -10 1.22464679914735e-15 1
+-19.3333333333333 -10 1.22464679914735e-15 1
+0 -10 1.22464679914735e-15 1
+-41.0121933088198 -7.07106781186548 -7.07106781186547 0.707106781186548
+-27.3414622058798 -7.07106781186548 -7.07106781186547 0.707106781186548
+-13.6707311029399 -7.07106781186548 -7.07106781186547 0.707106781186548
+0 -7.07106781186548 -7.07106781186547 0.707106781186548
+-58 -1.83697019872103e-15 -10 1
+-38.6666666666667 -2.04107799857892e-15 -10 1
+-19.3333333333333 -2.24518579843681e-15 -10 1
+0 -2.44929359829471e-15 -10 1
+-58 5.51091059616309e-16 9 1
+-38.6666666666667 3.67394039744206e-16 9 1
+-19.3333333333333 1.83697019872103e-16 9 1
+0 0 9 1
+-41.0121933088198 -6.36396103067893 6.36396103067893 0.707106781186548
+-27.3414622058798 -6.36396103067893 6.36396103067893 0.707106781186548
+-13.6707311029399 -6.36396103067893 6.36396103067893 0.707106781186548
+0 -6.36396103067893 6.36396103067893 0.707106781186548
+-58 -9 1.10218211923262e-15 1
+-38.6666666666667 -9 1.10218211923262e-15 1
+-19.3333333333333 -9 1.10218211923262e-15 1
+0 -9 1.10218211923262e-15 1
+-41.0121933088198 -6.36396103067893 -6.36396103067893 0.707106781186548
+-27.3414622058798 -6.36396103067893 -6.36396103067893 0.707106781186548
+-13.6707311029399 -6.36396103067893 -6.36396103067893 0.707106781186548
+0 -6.36396103067893 -6.36396103067893 0.707106781186548
+-58 -1.65327317884893e-15 -9 1
+-38.6666666666667 -1.83697019872103e-15 -9 1
+-19.3333333333333 -2.02066721859313e-15 -9 1
+0 -2.20436423846524e-15 -9 1
+
+700 1 0 0
+3 1
+4 4
+0 0 0 0 58 58 58 58
+3 3
+0 0 0 15.707963267949 15.707963267949 15.707963267949
+2 2
+0 0 1 1
+-58 10 0 1
+-42 10 0 1
+-26 10 0 1
+-10 10 0 1
+-41.0121933088198 7.07106781186548 7.07106781186547 0.707106781186548
+-29.698484809835 7.07106781186548 7.07106781186547 0.707106781186548
+-18.3847763108502 7.07106781186548 7.07106781186547 0.707106781186548
+-7.07106781186548 7.07106781186547 7.07106781186547 0.707106781186548
+-58 6.12323399573677e-16 10 1
+-38.6666666666667 4.08215599715784e-16 10 1
+-19.3333333333333 2.04107799857892e-16 10 1
+0 0 10 1
+-58 9 0 1
+-41.6666666666667 9 0 1
+-25.3333333333333 9 0 1
+-9 9 0 1
+-41.0121933088198 6.36396103067893 6.36396103067893 0.707106781186548
+-29.4627825494395 6.36396103067893 6.36396103067893 0.707106781186548
+-17.9133717900592 6.36396103067893 6.36396103067893 0.707106781186548
+-6.36396103067893 6.36396103067893 6.36396103067893 0.707106781186548
+-58 5.51091059616309e-16 9 1
+-38.6666666666667 3.67394039744206e-16 9 1
+-19.3333333333333 1.83697019872103e-16 9 1
+0 0 9 1
+
+700 1 0 0
+3 1
+4 4
+0 0 0 0 58 58 58 58
+3 3
+47.1238898038469 47.1238898038469 47.1238898038469 62.8318530717959 62.8318530717959 62.8318530717959
+2 2
+0 0 1 1
+-58 -1.83697019872103e-15 -10 1
+-38.6666666666667 -2.04107799857892e-15 -10 1
+-19.3333333333333 -2.24518579843681e-15 -10 1
+0 -2.44929359829471e-15 -10 1
+-41.0121933088198 7.07106781186547 -7.07106781186548 0.707106781186548
+-29.698484809835 7.07106781186547 -7.07106781186548 0.707106781186548
+-18.3847763108502 7.07106781186547 -7.07106781186548 0.707106781186548
+-7.07106781186548 7.07106781186547 -7.07106781186548 0.707106781186548
+-58 10 -2.44929359829471e-15 1
+-42 10 -1.63286239886314e-15 1
+-26 10 -8.16431199431569e-16 1
+-10 10 0 1
+-58 -1.65327317884893e-15 -9 1
+-38.6666666666667 -1.83697019872103e-15 -9 1
+-19.3333333333333 -2.02066721859313e-15 -9 1
+0 -2.20436423846524e-15 -9 1
+-41.0121933088198 6.36396103067893 -6.36396103067893 0.707106781186548
+-29.4627825494395 6.36396103067893 -6.36396103067893 0.707106781186548
+-17.9133717900592 6.36396103067893 -6.36396103067893 0.707106781186548
+-6.36396103067893 6.36396103067892 -6.36396103067893 0.707106781186548
+-58 9 -2.20436423846524e-15 1
+-41.6666666666667 9 -1.46957615897682e-15 1
+-25.3333333333333 9 -7.34788079488412e-16 1
+-9 9 0 1
+
+700 1 0 0
+3 1
+4 4
+0 0 0 0 55 55 55 55
+3 3
+47.1238898038469 47.1238898038469 47.1238898038469 62.8318530717959 62.8318530717959 62.8318530717959
+2 2
+0 0 1 1
+55 -1.65327317884893e-15 -9 1
+36.6666666666667 -1.83697019872103e-15 -9 1
+18.3333333333333 -2.02066721859313e-15 -9 1
+0 -2.20436423846524e-15 -9 1
+38.8908729652601 6.36396103067893 -6.36396103067893 0.707106781186548
+28.0485689870664 6.36396103067893 -6.36396103067893 0.707106781186548
+17.2062650088727 6.36396103067893 -6.36396103067893 0.707106781186548
+6.36396103067893 6.36396103067892 -6.36396103067893 0.707106781186548
+55 9 -2.20436423846524e-15 1
+39.6666666666667 9 -1.46957615897682e-15 1
+24.3333333333333 9 -7.34788079488412e-16 1
+9 9 0 1
+55 -1.83697019872103e-15 -10 1
+36.6666666666667 -2.04107799857892e-15 -10 1
+18.3333333333333 -2.24518579843681e-15 -10 1
+0 -2.44929359829471e-15 -10 1
+38.8908729652601 7.07106781186547 -7.07106781186548 0.707106781186548
+28.2842712474619 7.07106781186547 -7.07106781186548 0.707106781186548
+17.6776695296637 7.07106781186547 -7.07106781186548 0.707106781186548
+7.07106781186548 7.07106781186547 -7.07106781186548 0.707106781186548
+55 10 -2.44929359829471e-15 1
+40 10 -1.63286239886314e-15 1
+25 10 -8.16431199431569e-16 1
+10 10 0 1
+
+700 1 0 0
+3 1
+4 4
+0 0 0 0 55 55 55 55
+3 3
+0 0 0 15.707963267949 15.707963267949 15.707963267949
+2 2
+0 0 1 1
+55 9 0 1
+39.6666666666667 9 0 1
+24.3333333333333 9 0 1
+9 9 0 1
+38.8908729652601 6.36396103067893 6.36396103067893 0.707106781186548
+28.0485689870664 6.36396103067893 6.36396103067893 0.707106781186548
+17.2062650088727 6.36396103067893 6.36396103067893 0.707106781186548
+6.36396103067893 6.36396103067893 6.36396103067893 0.707106781186548
+55 5.51091059616309e-16 9 1
+36.6666666666667 3.67394039744206e-16 9 1
+18.3333333333333 1.83697019872103e-16 9 1
+0 0 9 1
+55 10 0 1
+40 10 0 1
+25 10 0 1
+10 10 0 1
+38.8908729652601 7.07106781186548 7.07106781186547 0.707106781186548
+28.2842712474619 7.07106781186548 7.07106781186547 0.707106781186548
+17.6776695296637 7.07106781186548 7.07106781186547 0.707106781186548
+7.07106781186548 7.07106781186547 7.07106781186547 0.707106781186548
+55 6.12323399573677e-16 10 1
+36.6666666666667 4.08215599715784e-16 10 1
+18.3333333333333 2.04107799857892e-16 10 1
+0 0 10 1
+
+700 1 0 0
+3 1
+4 4
+0 0 0 0 55 55 55 55
+5 3
+15.707963267949 15.707963267949 15.707963267949 31.4159265358979 31.4159265358979 47.1238898038469 47.1238898038469 47.1238898038469
+2 2
+0 0 1 1
+55 5.51091059616309e-16 9 1
+36.6666666666667 3.67394039744206e-16 9 1
+18.3333333333333 1.83697019872103e-16 9 1
+0 0 9 1
+38.8908729652601 -6.36396103067893 6.36396103067893 0.707106781186548
+25.9272486435067 -6.36396103067893 6.36396103067893 0.707106781186548
+12.9636243217534 -6.36396103067893 6.36396103067893 0.707106781186548
+0 -6.36396103067893 6.36396103067893 0.707106781186548
+55 -9 1.10218211923262e-15 1
+36.6666666666667 -9 1.10218211923262e-15 1
+18.3333333333333 -9 1.10218211923262e-15 1
+0 -9 1.10218211923262e-15 1
+38.8908729652601 -6.36396103067893 -6.36396103067893 0.707106781186548
+25.9272486435067 -6.36396103067893 -6.36396103067893 0.707106781186548
+12.9636243217534 -6.36396103067893 -6.36396103067893 0.707106781186548
+0 -6.36396103067893 -6.36396103067893 0.707106781186548
+55 -1.65327317884893e-15 -9 1
+36.6666666666667 -1.83697019872103e-15 -9 1
+18.3333333333333 -2.02066721859313e-15 -9 1
+0 -2.20436423846524e-15 -9 1
+55 6.12323399573677e-16 10 1
+36.6666666666667 4.08215599715784e-16 10 1
+18.3333333333333 2.04107799857892e-16 10 1
+0 0 10 1
+38.8908729652601 -7.07106781186547 7.07106781186548 0.707106781186548
+25.9272486435067 -7.07106781186548 7.07106781186548 0.707106781186548
+12.9636243217534 -7.07106781186548 7.07106781186548 0.707106781186548
+0 -7.07106781186548 7.07106781186548 0.707106781186548
+55 -10 1.22464679914735e-15 1
+36.6666666666667 -10 1.22464679914735e-15 1
+18.3333333333333 -10 1.22464679914735e-15 1
+0 -10 1.22464679914735e-15 1
+38.8908729652601 -7.07106781186548 -7.07106781186547 0.707106781186548
+25.9272486435067 -7.07106781186548 -7.07106781186547 0.707106781186548
+12.9636243217534 -7.07106781186548 -7.07106781186547 0.707106781186548
+0 -7.07106781186548 -7.07106781186547 0.707106781186548
+55 -1.83697019872103e-15 -10 1
+36.6666666666667 -2.04107799857892e-15 -10 1
+18.3333333333333 -2.24518579843681e-15 -10 1
+0 -2.44929359829471e-15 -10 1
+
+700 1 0 0
+3 1
+4 4
+0 0 0 0 35 35 35 35
+3 3
+15.707963267949 15.707963267949 15.707963267949 31.4159265358979 31.4159265358979 31.4159265358979
+2 2
+0 0 1 1
+-2.52786151779387e-15 35 10 1
+-2.73196931765177e-15 23.3333333333333 10 1
+-2.93607711750966e-15 11.6666666666667 10 1
+-3.14018491736755e-15 3.14018491736755e-15 10 1
+-7.07106781186548 24.7487373415292 7.07106781186547 0.707106781186548
+-7.07106781186548 18.8561808316413 7.07106781186547 0.707106781186548
+-7.07106781186548 12.9636243217534 7.07106781186547 0.707106781186548
+-7.07106781186548 7.07106781186548 7.07106781186547 0.707106781186548
+-10 35 1.22464679914735e-15 1
+-10 26.6666666666667 1.22464679914735e-15 1
+-10 18.3333333333333 1.22464679914735e-15 1
+-10 10 1.22464679914735e-15 1
+-2.27507536601449e-15 35 9 1
+-2.45877238588659e-15 23.3333333333333 9 1
+-2.64246940575869e-15 11.6666666666667 9 1
+-2.8261664256308e-15 2.8261664256308e-15 9 1
+-6.36396103067893 24.7487373415292 6.36396103067893 0.707106781186548
+-6.36396103067893 18.6204785712458 6.36396103067893 0.707106781186548
+-6.36396103067893 12.4922198009623 6.36396103067893 0.707106781186548
+-6.36396103067893 6.36396103067893 6.36396103067893 0.707106781186548
+-9 35 1.10218211923262e-15 1
+-9 26.3333333333333 1.10218211923262e-15 1
+-9 17.6666666666667 1.10218211923262e-15 1
+-9 9 1.10218211923262e-15 1
+
+700 1 0 0
+3 1
+4 4
+0 0 0 0 35 35 35 35
+3 3
+47.1238898038469 47.1238898038469 47.1238898038469 62.8318530717959 62.8318530717959 62.8318530717959
+2 2
+0 0 1 1
+-1.83697019872103e-15 35 -10 1
+-2.04107799857892e-15 23.3333333333333 -10 1
+-2.24518579843681e-15 11.6666666666667 -10 1
+-2.44929359829471e-15 0 -10 1
+7.07106781186547 24.7487373415292 -7.07106781186548 0.707106781186548
+7.07106781186547 18.8561808316413 -7.07106781186548 0.707106781186548
+7.07106781186547 12.9636243217534 -7.07106781186548 0.707106781186548
+7.07106781186547 7.07106781186548 -7.07106781186548 0.707106781186548
+10 35 -2.44929359829471e-15 1
+10 26.6666666666667 -1.63286239886314e-15 1
+10 18.3333333333333 -8.16431199431569e-16 1
+10 10 0 1
+-1.65327317884893e-15 35 -9 1
+-1.83697019872103e-15 23.3333333333333 -9 1
+-2.02066721859313e-15 11.6666666666667 -9 1
+-2.20436423846524e-15 0 -9 1
+6.36396103067893 24.7487373415292 -6.36396103067893 0.707106781186548
+6.36396103067893 18.6204785712458 -6.36396103067893 0.707106781186548
+6.36396103067893 12.4922198009623 -6.36396103067893 0.707106781186548
+6.36396103067892 6.36396103067893 -6.36396103067893 0.707106781186548
+9 35 -2.20436423846524e-15 1
+9 26.3333333333333 -1.46957615897682e-15 1
+9 17.6666666666667 -7.34788079488412e-16 1
+9 9 0 1
+
+700 1 0 0
+3 1
+4 4
+0 0 0 0 35 35 35 35
+3 3
+31.4159265358979 31.4159265358979 31.4159265358979 47.1238898038469 47.1238898038469 47.1238898038469
+2 2
+0 0 1 1
+-10 35 1.22464679914735e-15 1
+-10 26.6666666666667 1.22464679914735e-15 1
+-10 18.3333333333333 1.22464679914735e-15 1
+-10 10 1.22464679914735e-15 1
+-7.07106781186548 24.7487373415292 -7.07106781186547 0.707106781186548
+-7.07106781186548 18.8561808316413 -7.07106781186547 0.707106781186548
+-7.07106781186548 12.9636243217534 -7.07106781186547 0.707106781186548
+-7.07106781186548 7.07106781186548 -7.07106781186547 0.707106781186548
+-1.83697019872103e-15 35 -10 1
+-2.04107799857892e-15 23.3333333333333 -10 1
+-2.24518579843681e-15 11.6666666666667 -10 1
+-2.44929359829471e-15 0 -10 1
+-9 35 1.10218211923262e-15 1
+-9 26.3333333333333 1.10218211923262e-15 1
+-9 17.6666666666667 1.10218211923262e-15 1
+-9 9 1.10218211923262e-15 1
+-6.36396103067893 24.7487373415292 -6.36396103067893 0.707106781186548
+-6.36396103067893 18.6204785712458 -6.36396103067893 0.707106781186548
+-6.36396103067893 12.4922198009623 -6.36396103067893 0.707106781186548
+-6.36396103067893 6.36396103067893 -6.36396103067893 0.707106781186548
+-1.65327317884893e-15 35 -9 1
+-1.83697019872103e-15 23.3333333333333 -9 1
+-2.02066721859313e-15 11.6666666666667 -9 1
+-2.20436423846524e-15 0 -9 1
+
+700 1 0 0
+3 1
+4 4
+0 0 0 0 35 35 35 35
+3 3
+0 0 0 15.707963267949 15.707963267949 15.707963267949
+2 2
+0 0 1 1
+10 35 0 1
+10 26.6666666666667 0 1
+10 18.3333333333333 0 1
+10 10 0 1
+7.07106781186548 24.7487373415292 7.07106781186547 0.707106781186548
+7.07106781186548 18.8561808316413 7.07106781186547 0.707106781186548
+7.07106781186548 12.9636243217534 7.07106781186547 0.707106781186548
+7.07106781186547 7.07106781186548 7.07106781186547 0.707106781186548
+6.12323399573677e-16 35 10 1
+4.08215599715784e-16 23.3333333333333 10 1
+2.04107799857892e-16 11.6666666666667 10 1
+0 0 10 1
+9 35 0 1
+9 26.3333333333333 0 1
+9 17.6666666666667 0 1
+9 9 0 1
+6.36396103067893 24.7487373415292 6.36396103067893 0.707106781186548
+6.36396103067893 18.6204785712458 6.36396103067893 0.707106781186548
+6.36396103067893 12.4922198009623 6.36396103067893 0.707106781186548
+6.36396103067893 6.36396103067893 6.36396103067893 0.707106781186548
+5.51091059616309e-16 35 9 1
+3.67394039744206e-16 23.3333333333333 9 1
+1.83697019872103e-16 11.6666666666667 9 1
+0 0 9 1
+

Apps/LinearElasticity/Test/pipe_bifurcation.gno

@@ -0,0 +1,210 @@
+0
+0 1 2 3 4 5 6 7 
+28 1
+30 1
+32 1
+34 1
+36 1
+39 1
+42 1
+45 1
+-1 0
+-1 0
+-1 0
+-1 0
+110 1 3
+110 1 3
+110 1 2
+110 1 2
+110 1 2
+116 1 2
+-1
+1
+8 9 0 1 10 11 4 5 
+48 1
+28 1
+50 1
+32 1
+52 1
+53 1
+54 1
+55 1
+-1 0
+-1 0
+-1 0
+-1 0
+122 1 1
+122 1 1
+122 1 2
+110 1 2
+122 1 2
+124 1 2
+-1
+2
+2 3 8 9 6 7 10 11 
+30 1
+48 1
+34 1
+50 1
+56 1
+57 1
+58 1
+59 1
+-1 0
+-1 0
+-1 0
+-1 0
+126 1 1
+126 1 1
+110 1 2
+122 1 2
+126 1 2
+128 1 2
+-1
+3
+12 7 13 14 15 3 16 17 
+60 1
+62 1
+64 1
+66 1
+68 1
+69 1
+70 1
+71 1
+-1 0
+-1 0
+-1 0
+-1 0
+130 1 1
+130 1 1
+130 1 2
+130 1 2
+130 1 2
+132 1 2
+-1
+4
+13 14 18 5 16 17 19 1 
+62 1
+72 1
+66 1
+74 1
+76 1
+77 1
+78 1
+79 1
+-1 0
+-1 0
+-1 0
+-1 0
+134 1 1
+134 1 1
+130 1 2
+134 1 2
+134 1 2
+136 1 2
+-1
+5
+18 5 12 7 19 1 15 3 
+72 1
+60 1
+74 1
+64 1
+80 1
+45 1
+83 1
+39 1
+-1 0
+-1 0
+-1 0
+-1 0
+138 1 3
+107 1 -3
+134 1 2
+130 1 2
+138 1 2
+144 1 2
+-1
+6
+20 1 21 9 22 5 23 11 
+86 1
+88 1
+90 1
+92 1
+94 1
+53 -1
+95 1
+55 -1
+-1 0
+-1 0
+-1 0
+-1 0
+150 1 1
+122 -1 1
+150 1 2
+150 1 2
+150 1 2
+152 1 2
+-1
+7
+24 3 25 17 26 7 27 14 
+96 1
+98 1
+100 1
+102 1
+104 1
+71 1
+105 1
+69 1
+-1 0
+-1 0
+-1 0
+-1 0
+154 1 1
+129 1 -1
+154 1 2
+154 1 2
+154 1 2
+156 1 2
+-1
+8
+21 9 24 3 23 11 26 7 
+88 1
+96 1
+92 1
+100 1
+106 1
+57 -1
+107 1
+59 -1
+-1 0
+-1 0
+-1 0
+-1 0
+158 1 1
+126 -1 1
+150 1 2
+154 1 2
+158 1 2
+160 1 2
+-1
+9
+25 17 20 1 27 14 22 5 
+98 1
+86 1
+102 1
+90 1
+108 1
+79 1
+109 1
+77 1
+-1 0
+-1 0
+-1 0
+-1 0
+162 1 1
+133 1 -1
+154 1 2
+150 1 2
+162 1 2
+164 1 2
+-1

Apps/LinearElasticity/Test/pipe_bifurcation.prc

@@ -0,0 +1,10 @@
+123 0 2 0
+123 1 2 0
+123 2 2 0
+123 3 2 0
+123 4 2 0
+123 5 2 0
+1001 6 2 0
+1001 7 2 0
+1001 8 2 0
+1001 9 2 0

Apps/LinearElasticity/Test/plate_10x10x01.g2

@@ -0,0 +1,16 @@
+700 1 0 0
+3 0
+2 2
+0 0 1 1
+2 2
+0 0 1 1
+2 2
+0 0 1 1
+ 0  0 0
+10  0 0
+ 0 10 0
+10 10 0
+ 0  0 0.1
+10  0 0.1
+ 0 10 0.1
+10 10 0.1

Apps/LinearElasticity/Test/run.sh

@@ -0,0 +1,55 @@
+#!/bin/sh
+# $Id$
+##################################################################
+# This script defines a set of sample simulation cases that
+# is used to verify the integrity of an IFEM-based simulator.
+# Copy this file to a sub-folder Test in your App-directory and
+# insert the simulations that you want to use as regression tests.
+##################################################################
+
+# Define the name of the executable here
+mysim=LinEl
+
+run () {
+# This function runs a simulation with the specified options,
+# pipes the terminal output to a log-file, and compares it
+# with a previous simulation stored in the Reference folder.
+  inp=$1
+  log=`basename $1 .inp`.log
+  echo Running $inp ...
+  shift
+  time -f "real time %E, user time %U, system time %S" \
+  ../../../Release/bin/$mysim $inp $* > $log
+  if [ ! -e Reference/$log ]; then
+    mv $log Reference
+  elif cmp -s $log Reference/$log; then
+    echo Ok
+  else
+    echo "Warning: Discrepancies between current run and reference."
+    diff $log Reference
+  fi
+}
+
+if [ ! -d Reference ]; then mkdir Reference; fi
+
+#####################################
+# Enter the various cases below:
+# Format: run <inputfile> [<options>]
+#####################################
+
+run Hole2D-NURBS.inp    -2D -vtf 1 -nviz 4
+run Hole2D-Lagrange.inp -2D -vtf 1 -lagrange
+run Hole2D-Spectral.inp -2D -vtf 1 -spectral
+
+run Cylinder-NURBS.inp    -checkRHS -vtf 1 -nviz 5
+run Cylinder-Lagrange.inp -checkRHS -vtf 1 -lagrange
+run Cylinder-Spectral.inp -checkRHS -vtf 1 -spectral
+
+eigOpt="-eig 5 -nev 20 -ncv 40 -shift 2.0e-4"
+run SquarePlate-Splines.inp  $eigOpt -nGauss 3 -vtf 1 -nviz 3
+run SquarePlate-Lagrange.inp $eigOpt -nGauss 3 -vtf 1 -lagrange
+run SquarePlate-Spectral.inp $eigOpt -nGauss 3 -vtf 1 -spectral
+
+run PipeJoint-vibration.inp -vtf 1 -nu 7 -nv 7 -free -eig 4 -nev 16 -ncv 32
+run PipeJoint-NURBS.inp     -vtf 1 -nu 7 -nv 7
+run PipeJoint-Lagrange.inp  -vtf 1 -lagrange

Apps/run.sh

@@ -0,0 +1,39 @@
+#!/bin/sh
+# $Id$
+##################################################################
+# This script defines a set of sample simulation cases that
+# is used to verify the integrity of an IFEM-based simulator.
+# Copy this file to a sub-folder Test in your App-directory and
+# insert the simulations that you want to use as regression tests.
+##################################################################
+
+# Define the name of the executable here
+mysim=
+
+run () {
+# This function runs a simulation with the specified options,
+# pipes the terminal output to a log-file, and compares it
+# with a previous simulation stored in the Reference folder.
+  inp=$1
+  log=`basename $1 .inp`.log
+  echo Running $inp ...
+  shift
+  time -f "real time %E, user time %U, system time %S" \
+  ../Release/bin/$mysim $inp $* > $log
+  if [ ! -e Reference/$log ]; then
+    mv $log Reference
+  elif cmp -s $log Reference/$log; then
+    echo Ok
+  else
+    echo "Warning: Discrepancies between current run and reference."
+    diff $log Reference
+  fi
+}
+
+if [ ! -d Reference ]; then mkdir Reference; fi
+
+#####################################
+# Enter the various cases below:
+# Format: run <inputfile> [<options>]
+#####################################
+