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
DATADIR=${TOP_SRCDIR}/nfem/ptst
25
BINDIR="../../nfem/bin"
26
SBINDIR="../../nfem/sbin"
27
NPROC_MAX=${NPROC_MAX-"7"}
28
. "${TOP_SRCDIR}/config/loop_mpirun.sh"
32
# --------------------------------------------------------
33
# using direct surface mesh
34
# --------------------------------------------------------
36
# geo err_P1 err_P2 err_P3
38
circle_s-40-fix 2e-3 5e-7 7e-8
39
sphere_s-10-fix 3e-2 7e-4 6e-5
40
sphere_s_q-10-fix 3e-2 7e-4 1e-4
42
while test "$L" != ""; do
43
geo=`echo $L | awk '{print $1}'`
44
L=`echo $L | awk '{for (i=2; i <= NF; i++) print $i}'`
45
for Pk in P1 P2 P3; do
46
err=`echo $L | awk '{print $1}'`
47
L=`echo $L | awk '{for (i=2; i <= NF; i++) print $i}'`
48
loop_mpirun "./helmholtz_s_sphere $DATADIR/$geo-$Pk $Pk 2>/dev/null | RHEOPATH=$DATADIR ./helmholtz_s_sphere_error $err >/dev/null 2>/dev/null"
49
if test $? -ne 0; then status=1; fi