/usr/bin/tetgen2geo is in rheolef 5.93-2.
This file is owned by root:root, with mode 0o755.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 | #!/bin/sh
#
# This file is part of Rheolef.
#
# Copyright (C) 2000-2009 Pierre Saramito
#
# Rheolef is free software; you can redistribute it and/or modify
# it under the terms of the GNU General Public License as published by
# the Free Software Foundation; either version 2 of the License, or
# (at your option) any later version.
#
# Rheolef is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
# GNU General Public License for more details.
#
# You should have received a copy of the GNU General Public License
# along with Rheolef; if not, write to the Free Software
# Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
#
# -------------------------------------------------------------------------
PACKAGE=@PACKAGE@
VERSION=@VERSION@
in_srcdir=true
top_srcdir=../..
#Prog:tetgen2geo
#NAME: tetgen2geo - convert tetgen mesh in geo format
#@cindex mesh
#@pindex tetgen2geo
#@pindex geo
#@fiindex @file{.node} tetgen mesh nodes
#@fiindex @file{.ele} tetgen mesh elements
#@fiindex @file{.face} tetgen mesh faces
#@fiindex @file{.dmn} domain names
#@fiindex @file{.geo} mesh
#@toindex @code{tetgen}
#SYNOPSIS:
#@example
# tetgen2geo @var{options} @var{input}[.node] @var{input}[.ele] @var{input}[.face] @var{input}[.dmn]
#@end example
#
#DESCRIPTION:
#Convert a tetgen mesh into @file{.geo} one.
#The output goes to standart output.
#The @file{.dmn} file specifies the domain names,
#since @code{tetgen} mesh generator uses numbers as domain labels.
#
#EXAMPLE:
#@example
# tetgen toto.smesh
# tetgen2geo toto.1.node toto.1.ele toto.1.face toto.dmn > toto.geo
#@end example
#
#DOMAIN NAME FILE:
#This auxilliary @file{.dmn} file defines the boundary domain names
#as used by Rheolef, since @code{tetgen} uses numeric
#labels for domains.
#OPTIONS:
#@table @code
# @itemx -upgrade
# @itemx -noupgrade
# Default is to output a version 2 @file{.geo} file format. @xref{geo command}.
# With the @code{-noupgrade}, a version 1 file format is assumed.
#@end table
#End:
#============================================================
# Usage
#============================================================
PROG=`expr $0 : '.*/\(.*\)' \| $0`
USAGE="usage: $PROG [-[no]upgrade] input[.node] input[.ele] input[.face] input[.dmn]"
#============================================================
# Arguments
#============================================================
if test $# = 0; then
echo $USAGE 1>&2; exit 0
fi
tmpdir=${TMPDIR-"/tmp"}
do_dom=false
set_dom=false
upgrade="-upgrade"
node_file=""
ele_file=""
face_file=""
dom_list=""
geobin="geo"
while test $# -ne 0; do
case $1 in
-help) echo $USAGE 1>&2; exit 0;;
-noupgrade) upgrade="";;
-upgrade) upgrade="-upgrade";;
*.node) node_file=$1;;
*.ele) ele_file=$1;;
*.face) face_file=$1;;
*.dmn) dmn_file=$1;;
-dmn) if test x"$2" = x""; then
echo "$PROG: no argument to -dmn option found." 1>&2
echo $USAGE 1>&2; exit 1
fi
dmn_file=$2
shift
;;
-geobin) if test x"$2" = x""; then
echo "$PROG: no argument to -geobin option found." 1>&2
echo $USAGE 1>&2; exit 1
fi
geobin=$2
shift
;;
*)
echo "$PROG: invalid $1 option." 1>&2
echo $USAGE 1>&2; exit 1
;;
esac
shift
done
if test x"${node_file}" = x""; then
echo "tetgen node file .node missing" 1>&2
echo $USAGE 1>&2
exit 1
fi
if test x"${ele_file}" = x""; then
echo "tetgen element file .ele missing" 1>&2
echo $USAGE 1>&2
exit 1
fi
if test x"${face_file}" = x""; then
echo "tetgen face file .face missing" 1>&2
echo $USAGE 1>&2
exit 1
fi
#============================================================
# buid the `.dmn' domain name file
#============================================================
if test x"${dom_list}" != x""; then
dmn_file=${TMPDIR}/tetgen2geo-$$.dmn
/bin/rm -f ${dmn_file}
echo "DomainNames" > ${dmn_file}
n=0
for dom in ${dom_list}; do
n=`expr ${n} + 1`
done
echo ${n} >> ${dmn_file}
for dom in ${dom_list}; do
echo $dom >> ${dmn_file}
done
tmpdmn_file=${dmn_file}
elif test x"${dmn_file}" = x""; then
echo "$PROG: undefined domain names" 1>&2
echo $USAGE 1>&2; exit 1
elif test ! -f ${dmn_file}; then
echo "$PROG: file \"${dmn_file}\" not found" 1>&2
echo $USAGE 1>&2; exit 1
fi
#============================================================
# buid the `.geo' mesh file
#============================================================
cat ${node_file} ${ele_file} ${face_file} ${dmn_file} | \
sed -e 's/#.*$//' | \
${geobin} -name toto -input-tetgen -upgrade -geo -noverbose -
status=$?
if test $status -ne 0; then
echo "$PROG(fatal): ${geobin} exit with $status status." 1>&2
exit $status
fi
|