← Back to branch summary

~ubuntu-branches/ubuntu/wily/rheolef/wily

« back to all changes in this revision

Viewing changes to doc/pexamples/navier_stokes_taylor_cks_dg.sh

  • Committer: Package Import Robot
  • Author(s): Pierre Saramito, Pierre Saramito, Sylvestre Ledru
  • Date: 2013-04-30 11:29:51 UTC
  • mfrom: (1.2.1) (11.1.1 experimental)
  • Revision ID: package-import@ubuntu.com-20130430112951-g69nuqin77t75dek
Tags: 6.4-1
[ Pierre Saramito ]
* New upstream release 6.4 (major changes):
  - new powerful c++ expressions for FEM specifications
  - configure script improved (portability increased)
  - minor bug fixes
* control:
  - all architectures are considered (armel & armhf are back)
  - no more circular-dependency librheolef-dev <--> rheolef
* rules: add security compiler flags (hardening and -W,-l,relro)
* watch: include an updated version, thanks to B. Martens <bartm@debian.org>
* copyright: various file format fixes

[ Sylvestre Ledru ]
*  Standards-Version updated to 3.9.4

Show diffs side-by-side

added added

removed removed

Lines of Context:
doc/pexamples/navier_stokes_taylor_cks_dg.sh
 
1
#!/bin/sh
 
2
#
 
3
# This file is part of Rheolef.
 
4
#
 
5
# Copyright (C) 2000-2009 Pierre Saramito 
 
6
#
 
7
# Rheolef is free software; you can redistribute it and/or modify
 
8
# it under the terms of the GNU General Public License as published by
 
9
# the Free Software Foundation; either version 2 of the License, or
 
10
# (at your option) any later version.
 
11
#
 
12
# Rheolef is distributed in the hope that it will be useful,
 
13
# but WITHOUT ANY WARRANTY; without even the implied warranty of
 
14
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 
15
# GNU General Public License for more details.
 
16
#
 
17
# You should have received a copy of the GNU General Public License
 
18
# along with Rheolef; if not, write to the Free Software
 
19
# Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
 
20
#
 
21
# -------------------------------------------------------------------------
 
22
SRCDIR=${SRCDIR-"."}
 
23
TOP_SRCDIR=${TOP_SRCDIR-"../.."}
 
24
GEODIR=${TOP_SRCDIR}/nfem/ptst
 
25
NPROC_MAX=${NPROC_MAX-"2"}
 
26
. "${TOP_SRCDIR}/config/loop_mpirun.sh"
 
27
 
 
28
#echo "      not yet (skipped)"
 
29
#exit 0
 
30
 
 
31
status=0
 
32
run "$TOP_SRCDIR/nfem/sbin/mkgeo_grid_2d -t 10 -v4 2>/dev/null | $TOP_SRCDIR/nfem/pbin/geo -upgrade - > mesh-2d-t.geo 2>/dev/null"
 
33
if test $? -ne 0; then status=1; fi
 
34
 
 
35
run "$TOP_SRCDIR/nfem/sbin/mkgeo_grid_2d -q 10 -v4 2>/dev/null | $TOP_SRCDIR/nfem/pbin/geo -upgrade - > mesh-2d-q.geo 2>/dev/null"
 
36
if test $? -ne 0; then status=1; fi
 
37
 
 
38
Re=10
 
39
max_iter=3
 
40
not_yet=1e38
 
41
# geo                   approx  err_u_linf      err_p_linf
 
42
L="
 
43
mesh-2d-t               P1d     3e-2            4e-1
 
44
mesh-2d-t               P2d     2e-3            4e-2
 
45
mesh-2d-t               P3d     6e-5            4e-3
 
46
mesh-2d-q               P1d     3e-2            2e-0
 
47
mesh-2d-q               P2d     6e-4            5e-1
 
48
mesh-2d-q               P3d     2e-5            6e-3
 
49
"
 
50
 
 
51
while test "$L" != ""; do
 
52
  geo=`echo $L | awk '{print $1}'`
 
53
  approx=`echo $L | awk '{print $2}'`
 
54
  tol_u=`echo $L | awk '{print $3}'`
 
55
  tol_p=`echo $L | awk '{print $4}'`
 
56
  L=`echo $L | awk '{for (i=5; i <= NF; i++) print $i}'`
 
57
  loop_mpirun "./navier_stokes_taylor_cks_dg $geo $approx $Re $max_iter 2>/dev/null | RHEOPATH=.:$GEODIR \$RUN ./navier_stokes_taylor_error_dg $tol_u $tol_p -have-kinetic-energy >/dev/null 2>/dev/null"
 
58
  if test $? -ne 0; then status=1; fi
 
59
done
 
60
run "rm -f mesh-2d-t.geo mesh-2d-q.geo"
 
61
exit $status