3
# This file is part of Rheolef.
5
# Copyright (C) 2000-2009 Pierre Saramito
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.
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.
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
21
# -------------------------------------------------------------------------
23
TOP_SRCDIR=${TOP_SRCDIR-"../.."}
24
GEODIR=${TOP_SRCDIR}/nfem/ptst
25
NPROC_MAX=${NPROC_MAX-"2"}
26
. "${TOP_SRCDIR}/config/loop_mpirun.sh"
28
#echo " not yet (skipped)"
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
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
38
# geo approx tolerance (err_linf max(u,p))
48
while test "$L" != ""; do
49
geo=`echo $L | awk '{print $1}'`
50
approx=`echo $L | awk '{print $2}'`
51
tol=`echo $L | awk '{print $3}'`
52
L=`echo $L | awk '{for (i=4; i <= NF; i++) print $i}'`
53
loop_mpirun "./stokes_taylor_dg $geo $approx 2>/dev/null | RHEOPATH=.:$GEODIR \$RUN ./stokes_taylor_error_dg $tol >/dev/null 2>/dev/null"
54
if test $? -ne 0; then status=1; fi
56
run "rm -f mesh-2d-t.geo mesh-2d-q.geo"