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-"8"}
26
SBINDIR="${TOP_SRCDIR}/nfem/sbin"
27
BINDIR="${TOP_SRCDIR}/nfem/pbin"
28
. "${TOP_SRCDIR}/config/loop_mpirun.sh"
32
#echo " skipped (not yet)"
36
run "${SBINDIR}/mkgeo_grid_2d -v4 11 11 -a -1 -b 1 -c -1 -d 1 -boundary 2>/dev/null | ${BINDIR}/geo -upgrade - > mesh-2d.geo 2>/dev/null"
37
if test $? -ne 0; then exit 1; fi
39
progs="p_laplacian_newton"
41
pvalues="1.5 2.5" ;# 1.5 => singular matrix for P1 with mumps(np>1) and P2, P3 with mumps and pastix
42
for prog in $progs; do
43
for approx in $approxs; do
45
# Pk approx diverges for p < 2 && k >= 2
46
if test $p = "1.5" && test $approx = "P2" -o $approx = "P3"; then continue; fi
47
loop_mpirun "./$prog mesh-2d.geo $approx $p 1e-10 2>/dev/null | ./p_laplacian_post - -check >/dev/null 2>/dev/null"
48
if test $? -ne 0; then status=1; fi
52
# --------------------------------------------------------------------
54
# --------------------------------------------------------------------
55
run "/bin/rm -f mesh-2d.geo.geo"