#!/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:mesh2geo
#NAME: mesh2geo - convert mmg3d mesh in geo format
#@cindex mesh
#@pindex mesh2geo
#@pindex geo
#@fiindex @file{.mesh} mmg3d mesh
#@fiindex @file{.dmn} domain names
#@fiindex @file{.geo} mesh
#@toindex @code{mmg3d}
#SYNOPSIS:
#@example
#  mesh2geo @var{options} @var{input}[.mesh] @var{input}[.dmn]
#@end example
#
#DESCRIPTION:
#Convert a mesh @file{.mesh} into @file{.geo} one.
#The output goes to standart output.
#The @file{.dmn} file specifies the domain names,
#since @code{mmg3d} mesh generator uses numbers as domain labels.
#
#EXAMPLE:
#@example
#  mmg3d -O 1 -in toto.mesh -out toto-opt.mesh
#  mesh2geo toto-opt.mesh toto.dmn > toto-opt.geo
#@end example
#
#DOMAIN NAME FILE:
#This auxilliary @file{.dmn} file defines the boundary domain names
#as used by Rheolef, since @code{mmg3d} 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[.mesh] 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"
mesh_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";;
      	*.mesh)       mesh_file=$1;;
      	*.dmn)        dmn_file=$1;;
      	-mesh)        if test x"$2" = x""; then
			echo "$PROG: no argument to -mesh option found." 1>&2
		        echo $USAGE 1>&2; exit 1
                      fi
      		      mesh_file=$2
		      shift
		      ;;
      	-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"${mesh_file}" = x""; then
   echo $USAGE 1>&2
   exit 1
fi

#============================================================
# buid the `.dmn' domain name file
#============================================================

if test x"${dom_list}" != x""; then
    dmn_file=${TMPDIR}/mesh2geo-$$.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 ${mesh_file} ${dmn_file} | ${geobin} -name toto -input-mmg3d -upgrade -geo -noverbose -
status=$?
if test $status -ne 0; then
    echo "$PROG(fatal): ${geobin} exit with $status status." 1>&2
    exit $status
fi


