#!/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


