/*----------------------------------------------------------*/
/*                                                                                                                      */
/*                                              LIBMESH V 5.46                                          */
/*                                                                                                                      */
/*----------------------------------------------------------*/
/*                                                                                                                      */
/*      Description:            handle .meshb file format I/O           */
/*      Author:                         Loic MARECHAL                                           */
/*      Creation date:          feb 16 2007                                                     */
/*      Last modification:      apr 03 2012                                                     */
/*                                                                                                                      */
/*----------------------------------------------------------*/


/*----------------------------------------------------------*/
/* Includes                                                                                                     */
/*----------------------------------------------------------*/

#include <stdio.h>
#include <stdlib.h>
#include <stdarg.h>
#include <string.h>
#include <float.h>
#include <math.h>
#include <ctype.h>
#include "libmesh5.h"


/*----------------------------------------------------------*/
/* Defines                                                                                                      */
/*----------------------------------------------------------*/

#define Asc 1
#define Bin 2
#define MshFil 4
#define SolFil 8
#define MaxMsh 100
#define InfKwd 1
#define RegKwd 2
#define SolKwd 3
#define WrdSiz 4
#define BufSiz 10000


/*----------------------------------------------------------*/
/* Structures                                                                                           */
/*----------------------------------------------------------*/

typedef struct
{
        int typ, SolSiz, NmbWrd, NmbLin, NmbTyp, TypTab[ GmfMaxTyp ];
        long pos;
        char fmt[ GmfMaxTyp*9 ];
}KwdSct;

typedef struct
{
        int dim, ver, mod, typ, cod, pos;
        long NexKwdPos, siz;
        KwdSct KwdTab[ GmfMaxKwd + 1 ];
        FILE *hdl;
        int *IntBuf;
        float *FltBuf;
        unsigned char *buf;
        char FilNam[ GmfStrSiz ];
        double DblBuf[1000/8];
        unsigned char blk[ BufSiz + 1000 ];
}GmfMshSct;


/*----------------------------------------------------------*/
/* Global variables                                                                                     */
/*----------------------------------------------------------*/

static int GmfIniFlg=0;
static GmfMshSct *GmfMshTab[ MaxMsh + 1 ];
static const char *GmfKwdFmt[ GmfMaxKwd + 1 ][4] = 
{       {"Reserved", "", "", ""},
        {"MeshVersionFormatted", "", "", "i"},
        {"Reserved", "", "", ""},
        {"Dimension", "", "", "i"},
        {"Vertices", "Vertex", "i", "dri"},
        {"Edges", "Edge", "i", "iii"},
        {"Triangles", "Triangle", "i", "iiii"},
        {"Quadrilaterals", "Quadrilateral", "i", "iiiii"},
        {"Tetrahedra", "Tetrahedron", "i", "iiiii"},
        {"Prisms", "Prism", "i", "iiiiiii"},
        {"Hexahedra", "Hexahedron", "i", "iiiiiiiii"},
        {"IterationsAll", "IterationAll","","i"},
        {"TimesAll", "TimeAll","","r"},                                 
        {"Corners", "Corner", "i", "i"},
        {"Ridges", "Ridge", "i", "i"},
        {"RequiredVertices", "RequiredVertex", "i", "i"},
        {"RequiredEdges", "RequiredEdge", "i", "i"},
        {"RequiredTriangles", "RequiredTriangle", "i", "i"},
        {"RequiredQuadrilaterals", "RequiredQuadrilateral", "i", "i"},
        {"TangentAtEdgeVertices", "TangentAtEdgeVertex", "i", "iii"},
        {"NormalAtVertices", "NormalAtVertex", "i", "ii"},
        {"NormalAtTriangleVertices", "NormalAtTriangleVertex", "i", "iii"},
        {"NormalAtQuadrilateralVertices", "NormalAtQuadrilateralVertex", "i", "iiii"},
        {"AngleOfCornerBound", "", "", "r"},
        {"TrianglesP2", "TriangleP2", "i", "iiiiiii"},
        {"EdgesP2", "EdgeP2", "i", "iiii"},
        {"SolAtPyramids", "SolAtPyramid", "i", "sr"},
        {"QuadrilateralsQ2", "QuadrilateralQ2", "i", "iiiiiiiiii"},
        {"ISolAtPyramids", "ISolAtPyramid", "i", "iiiii"},
        {"SubDomainFromGeom", "SubDomainFromGeom", "i", "iiii"},
        {"TetrahedraP2", "TetrahedronP2", "i", "iiiiiiiiiii"},
        {"Fault_NearTri", "Fault_NearTri", "i", "i"},
        {"Fault_Inter", "Fault_Inter", "i", "i"},
        {"HexahedraQ2", "HexahedronQ2", "i", "iiiiiiiiiiiiiiiiiiiiiiiiiiii"},
        {"ExtraVerticesAtEdges", "ExtraVerticesAtEdge", "i", "in"},
        {"ExtraVerticesAtTriangles", "ExtraVerticesAtTriangle", "i", "in"},
        {"ExtraVerticesAtQuadrilaterals", "ExtraVerticesAtQuadrilateral", "i", "in"},
        {"ExtraVerticesAtTetrahedra", "ExtraVerticesAtTetrahedron", "i", "in"},
        {"ExtraVerticesAtPrisms", "ExtraVerticesAtPrism", "i", "in"},
        {"ExtraVerticesAtHexahedra", "ExtraVerticesAtHexahedron", "i", "in"},
        {"VerticesOnGeometricVertices", "VertexOnGeometricVertex", "i", "iir"},
        {"VerticesOnGeometricEdges", "VertexOnGeometricEdge", "i", "iirr"},
        {"VerticesOnGeometricTriangles", "VertexOnGeometricTriangle", "i", "iirrr"},
        {"VerticesOnGeometricQuadrilaterals", "VertexOnGeometricQuadrilateral", "i", "iirrr"},
        {"EdgesOnGeometricEdges", "EdgeOnGeometricEdge", "i", "iir"},
        {"Fault_FreeEdge", "Fault_FreeEdge", "i", "i"},
        {"Polyhedra", "Polyhedron", "i", "iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii"},
        {"Polygons", "Polygon", "", "iiiiiiiii"},
        {"Fault_Overlap", "Fault_Overlap", "i", "i"},
        {"Pyramids", "Pyramid", "i", "iiiiii"},
        {"BoundingBox", "", "", "drdr"},
        {"Body","i", "drdrdrdr"},
        {"PrivateTable", "PrivateTable", "i", "i"},
        {"Fault_BadShape", "Fault_BadShape", "i", "i"},
        {"End", "", "", ""},
        {"TrianglesOnGeometricTriangles", "TriangleOnGeometricTriangle", "i", "iir"},
        {"TrianglesOnGeometricQuadrilaterals", "TriangleOnGeometricQuadrilateral", "i", "iir"},
        {"QuadrilateralsOnGeometricTriangles", "QuadrilateralOnGeometricTriangle", "i", "iir"},
        {"QuadrilateralsOnGeometricQuadrilaterals", "QuadrilateralOnGeometricQuadrilateral", "i", "iir"},
        {"Tangents", "Tangent", "i", "dr"},
        {"Normals", "Normal", "i", "dr"},
        {"TangentAtVertices", "TangentAtVertex", "i", "ii"},
        {"SolAtVertices", "SolAtVertex", "i", "sr"},
        {"SolAtEdges", "SolAtEdge", "i", "sr"},
        {"SolAtTriangles", "SolAtTriangle", "i", "sr"},
        {"SolAtQuadrilaterals", "SolAtQuadrilateral", "i", "sr"},
        {"SolAtTetrahedra", "SolAtTetrahedron", "i", "sr"},
        {"SolAtPrisms", "SolAtPrism", "i", "sr"},
        {"SolAtHexahedra", "SolAtHexahedron", "i", "sr"},
        {"DSolAtVertices", "DSolAtVertex", "i", "sr"},
        {"ISolAtVertices", "ISolAtVertex", "i", "i"},
        {"ISolAtEdges", "ISolAtEdge", "i", "ii"},
        {"ISolAtTriangles", "ISolAtTriangle", "i", "iii"},
        {"ISolAtQuadrilaterals", "ISolAtQuadrilateral", "i", "iiii"},
        {"ISolAtTetrahedra", "ISolAtTetrahedron", "i", "iiii"},
        {"ISolAtPrisms", "ISolAtPrism", "i", "iiiiii"},
        {"ISolAtHexahedra", "ISolAtHexahedron", "i", "iiiiiiii"},
        {"Iterations", "","","i"},
        {"Time", "","","r"},
        {"Fault_SmallTri", "Fault_SmallTri","i","i"},
        {"CoarseHexahedra", "CoarseHexahedron", "i", "i"}
 };


/*----------------------------------------------------------*/
/* Prototypes of local procedures                                                       */
/*----------------------------------------------------------*/

static void ScaWrd(GmfMshSct *, unsigned char *);
static void ScaDblWrd(GmfMshSct *, unsigned char *);
static void ScaBlk(GmfMshSct *, unsigned char *, int);
static long GetPos(GmfMshSct *);
static void RecWrd(GmfMshSct *, unsigned char *);
static void RecDblWrd(GmfMshSct *, unsigned char *);
static void RecBlk(GmfMshSct *, unsigned char *, int);
static void SetPos(GmfMshSct *, long);
static int ScaKwdTab(GmfMshSct *);
static void ExpFmt(GmfMshSct *, int);
static void ScaKwdHdr(GmfMshSct *, int);


/*----------------------------------------------------------*/
/* Open a mesh file in read or write mod                                        */
/*----------------------------------------------------------*/

int GmfOpenMesh(const char *FilNam, int mod, ...)
{
        int i, KwdCod, res, *PtrVer, *PtrDim, MshIdx=0;
        char str[ GmfStrSiz ];
        va_list VarArg;
        GmfMshSct *msh;
        char *ptr;
        int k;

        if(!GmfIniFlg)
        {
                for(i=0;i<=MaxMsh;i++)
                        GmfMshTab[i] = NULL;

                GmfIniFlg = 1;
        }

        /*---------------------*/
        /* MESH STRUCTURE INIT */
        /*---------------------*/

        for(i=1;i<=MaxMsh;i++)
                if(!GmfMshTab[i])
                {
                        MshIdx = i;
                        break;
                }

        if( !MshIdx || !(msh = calloc(1, sizeof(GmfMshSct))) )
                return(0);

        /* Copy the FilNam into the structure */

        if(strlen(FilNam) + 7 >= GmfStrSiz)
        {
                free (msh);
                return(0);
        }

        strcpy(msh->FilNam, FilNam);

        /* Store the opening mod (read or write) and guess the filetype (binary or ascii) depending on the extension */

        msh->mod = mod;
        msh->buf = (unsigned char *)msh->DblBuf;
        msh->FltBuf = (float *)msh->DblBuf;
        msh->IntBuf = (int *)msh->DblBuf;

        k = strlen(msh->FilNam) - 6;
        if(k < 0)
                k = 0;
        ptr = msh->FilNam+k;
        if(strstr(ptr, ".meshb"))
                msh->typ |= (Bin | MshFil);
        else if(strstr(ptr, ".mesh"))
                msh->typ |= (Asc | MshFil);
        else if(strstr(ptr, ".solb"))
                msh->typ |= (Bin | SolFil);
        else if(strstr(ptr, ".sol"))
                msh->typ |= (Asc | SolFil);
        else {
                free (msh);
                return(0);
        }

        /* Open the file in the required mod and initialyse the mesh structure */

        if(msh->mod == GmfRead)
        {

                /*-----------------------*/
                /* OPEN FILE FOR READING */
                /*-----------------------*/

                va_start(VarArg, mod);
                PtrVer = va_arg(VarArg, int *);
                PtrDim = va_arg(VarArg, int *);
                va_end(VarArg);

                /* Create the name string and open the file */

                if(!(msh->hdl = fopen(msh->FilNam, "rb")))
                {
                        free (msh);
                        return(0);
                }

                /* Read the endian coding tag, the mesh version and the mesh dimension (mandatory kwd) */

                if(msh->typ & Bin)
                {
                        fread((unsigned char *)&msh->cod, WrdSiz, 1, msh->hdl);

                        if( (msh->cod != 1) && (msh->cod != 16777216) )
                        {
                                free (msh);
                                return(0);
                        }

                        ScaWrd(msh, (unsigned char *)&msh->ver);

                        if( (msh->ver < 1) || (msh->ver > 3) )
                        {
                                free (msh);
                                return(0);
                        }

                        if( (msh->ver == 3) && (sizeof(long) == 4) )
                        {
                                free (msh);
                                return(0);
                        }

                        ScaWrd(msh, (unsigned char *)&KwdCod);

                        if(KwdCod != GmfDimension)
                        {
                                free (msh);
                                return(0);
                        }

                        GetPos(msh);
                        ScaWrd(msh, (unsigned char *)&msh->dim);
                }
                else
                {
                        do
                        {
                                res = fscanf(msh->hdl, "%s", str);
                        }while( (res != EOF) && strcmp(str, "MeshVersionFormatted") );

                        if(res == EOF)
                        {
                                free (msh);
                                return(0);
                        }

                        fscanf(msh->hdl, "%d", &msh->ver);

                        if( (msh->ver < 1) || (msh->ver > 3) )
                        {
                                free (msh);
                                return(0);
                        }

                        do
                        {
                                res = fscanf(msh->hdl, "%s", str);
                        }while( (res != EOF) && strcmp(str, "Dimension") );

                        if(res == EOF)
                        {
                                free (msh);
                                return(0);
                        }

                        fscanf(msh->hdl, "%d", &msh->dim);
                }

                if( (msh->dim != 2) && (msh->dim != 3) )
                {
                        free (msh);
                        return(0);
                }

                (*PtrVer) = msh->ver;
                (*PtrDim) = msh->dim;

                /*------------*/
                /* KW READING */
                /*------------*/

                /* Read the list of kw present in the file */

                if(!ScaKwdTab(msh))
                {
                        free (msh);
                        return(0);
                }

                GmfMshTab[ MshIdx ] = msh;

                return(MshIdx);
        }
        else if(msh->mod == GmfWrite)
        {

                /*-----------------------*/
                /* OPEN FILE FOR WRITING */
                /*-----------------------*/

                msh->cod = 1;

                /* Check if the user provided a valid version number and dimension */

                va_start(VarArg, mod);
                msh->ver = va_arg(VarArg, int);
                msh->dim = va_arg(VarArg, int);
                va_end(VarArg);

                if( (msh->ver < 1) || (msh->ver > 3) )
                {
                        free (msh);
                        return(0);
                }

                if( (msh->ver == 3) && (sizeof(long) == 4) )
                {
                        free (msh);
                        return(0);
                }

                if( (msh->dim != 2) && (msh->dim != 3) )
                {
                        free (msh);
                        return(0);
                }

                /* Create the mesh file */

                if(!(msh->hdl = fopen(msh->FilNam, "wb")))
                {
                        free (msh);
                        return(0);
                }

                GmfMshTab[ MshIdx ] = msh;


                /*------------*/
                /* KW WRITING */
                /*------------*/

                /* Write the mesh version and dimension */

                if(msh->typ & Asc)
                {
                        fprintf(msh->hdl, "%s %d\n\n", GmfKwdFmt[ GmfVersionFormatted ][0], msh->ver);
                        fprintf(msh->hdl, "%s %d\n", GmfKwdFmt[ GmfDimension ][0], msh->dim);
                }
                else
                {
                        RecWrd(msh, (unsigned char *)&msh->cod);
                        RecWrd(msh, (unsigned char *)&msh->ver);
                        GmfSetKwd(MshIdx, GmfDimension, 0);
                        RecWrd(msh, (unsigned char *)&msh->dim);
                }

                return(MshIdx);
        }
        else
        {
                free (msh);
                return(0);
        }
}


/*----------------------------------------------------------*/
/* Close a meshfile in the right way                                            */
/*----------------------------------------------------------*/

int GmfCloseMesh(int MshIdx)
{
        int res = 1;
        GmfMshSct *msh;

        if( (MshIdx < 1) || (MshIdx > MaxMsh) )
                return(0);

        msh = GmfMshTab[ MshIdx ];
        RecBlk(msh, msh->buf, 0);

        /* In write down the "End" kw in write mode */

        if(msh->mod == GmfWrite){
                if(msh->typ & Asc)
                        fprintf(msh->hdl, "\n%s\n", GmfKwdFmt[ GmfEnd ][0]);
                else
                        GmfSetKwd(MshIdx, GmfEnd, 0);
        }
        /* Close the file and free the mesh structure */

        if(fclose(msh->hdl))
                res = 0;

        free(msh);
        GmfMshTab[ MshIdx ] = NULL;

        return(res);
}


/*----------------------------------------------------------*/
/* Read the number of lines and set the position to this kwd*/
/*----------------------------------------------------------*/

int GmfStatKwd(int MshIdx, int KwdCod, ...)
{
        int i, *PtrNmbTyp, *PtrSolSiz, *TypTab;
        GmfMshSct *msh;
        KwdSct *kwd;
        va_list VarArg;

        if( (MshIdx < 1) || (MshIdx > MaxMsh) )
                return(0);

        msh = GmfMshTab[ MshIdx ];

        if( (KwdCod < 1) || (KwdCod > GmfMaxKwd) )
                return(0);

        kwd = &msh->KwdTab[ KwdCod ];

        if(!kwd->NmbLin)
                return(0);

        /* Read further arguments if this kw is a sol */

        if(kwd->typ == SolKwd)
        {
                va_start(VarArg, KwdCod);

                PtrNmbTyp = va_arg(VarArg, int *);
                *PtrNmbTyp = kwd->NmbTyp;

                PtrSolSiz = va_arg(VarArg, int *);
                *PtrSolSiz = kwd->SolSiz;

                TypTab = va_arg(VarArg, int *);

                for(i=0;i<kwd->NmbTyp;i++)
                        TypTab[i] = kwd->TypTab[i];

                va_end(VarArg);
        }

        return(kwd->NmbLin);
}


/*----------------------------------------------------------*/
/* Set the current file position to a given kwd                         */
/*----------------------------------------------------------*/

int GmfGotoKwd(int MshIdx, int KwdCod)
{
        GmfMshSct *msh;
        KwdSct *kwd;

        if( (MshIdx < 1) || (MshIdx > MaxMsh) )
                return(0);

        msh = GmfMshTab[ MshIdx ];

        if( (KwdCod < 1) || (KwdCod > GmfMaxKwd) )
                return(0);

        kwd = &msh->KwdTab[ KwdCod ];

        if(!kwd->NmbLin)
                return(0);

        return(fseek(msh->hdl, kwd->pos, SEEK_SET));
}


/*----------------------------------------------------------*/
/* Write the kwd and set the number of lines                            */
/*----------------------------------------------------------*/

int GmfSetKwd(int MshIdx, int KwdCod, ...)
{
        int i, NmbLin=0, *TypTab;
        long CurPos;
        va_list VarArg;
        GmfMshSct *msh;
        KwdSct *kwd;

        if( (MshIdx < 1) || (MshIdx > MaxMsh) )
                return(0);

        msh = GmfMshTab[ MshIdx ];
        RecBlk(msh, msh->buf, 0);

        if( (KwdCod < 1) || (KwdCod > GmfMaxKwd) )
                return(0);

        kwd = &msh->KwdTab[ KwdCod ];

        /* Read further arguments if this kw has a header */

        if(strlen(GmfKwdFmt[ KwdCod ][2]))
        {
                va_start(VarArg, KwdCod);
                NmbLin = va_arg(VarArg, int);

                if(!strcmp(GmfKwdFmt[ KwdCod ][3], "sr"))
                {
                        kwd->NmbTyp = va_arg(VarArg, int);
                        TypTab = va_arg(VarArg, int *);

                        for(i=0;i<kwd->NmbTyp;i++)
                                kwd->TypTab[i] = TypTab[i];
                }

                va_end(VarArg);
        }

        /* Setup the kwd info */

        ExpFmt(msh, KwdCod);

        if(!kwd->typ)
                return(0);
        else if(kwd->typ == InfKwd)
                kwd->NmbLin = 1;
        else
                kwd->NmbLin = NmbLin;

        /* Store the next kwd position in binary file */

        if( (msh->typ & Bin) && msh->NexKwdPos )
        {
                CurPos = ftell(msh->hdl);
                fseek(msh->hdl, msh->NexKwdPos, SEEK_SET);
                SetPos(msh, CurPos);
                fseek(msh->hdl, CurPos, SEEK_SET);
        }

        /* Write the header */

        if(msh->typ & Asc)
        {
                fprintf(msh->hdl, "\n%s\n", GmfKwdFmt[ KwdCod ][0]);

                if(kwd->typ != InfKwd)
                        fprintf(msh->hdl, "%d\n", kwd->NmbLin);

                /* In case of solution field, write the extended header */

                if(kwd->typ == SolKwd)
                {
                        fprintf(msh->hdl, "%d ", kwd->NmbTyp);

                        for(i=0;i<kwd->NmbTyp;i++)
                                fprintf(msh->hdl, "%d ", kwd->TypTab[i]);

                        fprintf(msh->hdl, "\n\n");
                }
        }
        else
        {
                RecWrd(msh, (unsigned char *)&KwdCod);
                msh->NexKwdPos = ftell(msh->hdl);
                SetPos(msh, 0);

                if(kwd->typ != InfKwd)
                        RecWrd(msh, (unsigned char *)&kwd->NmbLin);

                /* In case of solution field, write the extended header at once */

                if(kwd->typ == SolKwd)
                {
                        RecWrd(msh, (unsigned char *)&kwd->NmbTyp);

                        for(i=0;i<kwd->NmbTyp;i++)
                                RecWrd(msh, (unsigned char *)&kwd->TypTab[i]);
                }
        }

        /* Reset write buffer position */
        msh->pos = 0;

        /* Estimate the total file size and check whether it crosses the 2GB threshold */

        msh->siz += kwd->NmbLin * kwd->NmbWrd * WrdSiz;

        if(msh->siz > 2E9)
                return(0);
        else
                return(kwd->NmbLin);
}


/*----------------------------------------------------------*/
/* Read a full line from the current kwd                                        */
/*----------------------------------------------------------*/

void GmfGetLin(int MshIdx, int KwdCod, ...)
{
        int i, j;
        float *FltSolTab;
        double *DblSolTab;
        va_list VarArg;
        GmfMshSct *msh = GmfMshTab[ MshIdx ];
        KwdSct *kwd = &msh->KwdTab[ KwdCod ];

        /* Start decoding the arguments */

        va_start(VarArg, KwdCod);

        if(kwd->typ != SolKwd)
        {
                int k, nb_repeat = 0;

                if(msh->ver == 1)
                {
                        if(msh->typ & Asc)
                        {
                                for(i=0;i<kwd->SolSiz;i++)
                                        if(kwd->fmt[i] == 'r')
                                                fscanf(msh->hdl, "%f", va_arg(VarArg, float *));
                                        else if(kwd->fmt[i] == 'n') {
                                                fscanf(msh->hdl, "%d", &nb_repeat);
                                                *(va_arg(VarArg,  int *)) = nb_repeat;
                                                for(k=0;k<nb_repeat;k++)
                                                        fscanf(msh->hdl, "%d", va_arg(VarArg, int *));
                                        }
                                        else
                                                fscanf(msh->hdl, "%d", va_arg(VarArg, int *));
                        }
                        else
                        {
                                for(i=0;i<kwd->SolSiz;i++)
                                        if(kwd->fmt[i] == 'r')
                                                ScaWrd(msh, (unsigned char *)va_arg(VarArg, float *));
                                        else if(kwd->fmt[i] == 'n') {
                                                ScaWrd(msh, (unsigned char *)&nb_repeat);
                                                *(va_arg(VarArg,  int *)) = nb_repeat;
                                                for(k=0;k<nb_repeat;k++)
                                                        ScaWrd(msh, (unsigned char *)va_arg(VarArg, int *));
                                        }
                                        else
                                                ScaWrd(msh, (unsigned char *)va_arg(VarArg, int *));
                        }
                }
                else
                {
                        if(msh->typ & Asc)
                        {
                                for(i=0;i<kwd->SolSiz;i++)
                                        if(kwd->fmt[i] == 'r')
                                                fscanf(msh->hdl, "%lf", va_arg(VarArg, double *));
                                        else if(kwd->fmt[i] == 'n') {
                                                fscanf(msh->hdl, "%d", &nb_repeat);
                                                *(va_arg(VarArg,  int *)) = nb_repeat;
                                                for(k=0;k<nb_repeat;k++)
                                                        fscanf(msh->hdl, "%d", va_arg(VarArg, int *));
                                        }
                                        else
                                                fscanf(msh->hdl, "%d", va_arg(VarArg, int *));
                        }
                        else
                                for(i=0;i<kwd->SolSiz;i++)
                                        if(kwd->fmt[i] == 'r')
                                                ScaDblWrd(msh, (unsigned char *)va_arg(VarArg, double *));
                                        else if(kwd->fmt[i] == 'n') {
                                                ScaWrd(msh, (unsigned char *)&nb_repeat);
                                                *(va_arg(VarArg,  int *)) = nb_repeat;
                                                for(k=0;k<nb_repeat;k++)
                                                        ScaWrd(msh, (unsigned char *)va_arg(VarArg, int *));
                                        }
                                        else
                                                ScaWrd(msh, (unsigned char *)va_arg(VarArg, int *));
                }
        }
        else
        {
                if(msh->ver == 1)
                {
                        FltSolTab = va_arg(VarArg, float *);

                        if(msh->typ & Asc)
                                for(j=0;j<kwd->SolSiz;j++)
                                        fscanf(msh->hdl, "%f", &FltSolTab[j]);
                        else
                                ScaBlk(msh, (unsigned char *)FltSolTab, kwd->NmbWrd);
                }
                else
                {
                        DblSolTab = va_arg(VarArg, double *);

                        if(msh->typ & Asc)
                                for(j=0;j<kwd->SolSiz;j++)
                                        fscanf(msh->hdl, "%lf", &DblSolTab[j]);
                        else
                                for(j=0;j<kwd->SolSiz;j++)
                                        ScaDblWrd(msh, (unsigned char *)&DblSolTab[j]);
                }
        }

        va_end(VarArg);
}


/*----------------------------------------------------------*/
/* Write a full line from the current kwd                                       */
/*----------------------------------------------------------*/

void GmfSetLin(int MshIdx, int KwdCod, ...)
{
        int i, j, pos, *IntBuf;
        float *FltSolTab;
        double *DblSolTab, *DblBuf;
        va_list VarArg;
        GmfMshSct *msh = GmfMshTab[ MshIdx ];
        KwdSct *kwd = &msh->KwdTab[ KwdCod ];

        /* Start decoding the arguments */

        va_start(VarArg, KwdCod);

        if(kwd->typ != SolKwd)
        {
                int k, nb_repeat = 0;

                if(msh->ver == 1)
                {
                        if(msh->typ & Asc)
                        {
                                for(i=0;i<kwd->SolSiz;i++)
                                        if(kwd->fmt[i] == 'r')
                                                fprintf(msh->hdl, "%g ", (float)va_arg(VarArg, double));
                                        else if(kwd->fmt[i] == 'n') {
                                                nb_repeat = va_arg(VarArg, int);
                                                fprintf(msh->hdl, "%d ", nb_repeat);
                                                for(k=0;k<nb_repeat;k++)
                                                        fprintf(msh->hdl, "%d ", va_arg(VarArg, int));
                                        }
                                        else
                                                fprintf(msh->hdl, "%d ", va_arg(VarArg, int));
                        }
                        else
                        {
                                int size_of_block = kwd->SolSiz;
                                for(i=0;i<kwd->SolSiz;i++)
                                        if(kwd->fmt[i] == 'r')
                                                msh->FltBuf[i] = va_arg(VarArg, double);
                                        else if(kwd->fmt[i] == 'n') {
                                                nb_repeat = va_arg(VarArg, int);
                                                msh->FltBuf[i] = nb_repeat;
                                                for(k=0;k<nb_repeat;k++) {
                                                        msh->IntBuf[i+1+k] = va_arg(VarArg, int);
                                                        size_of_block ++;
                                                }
                                        }
                                        else
                                                msh->IntBuf[i] = va_arg(VarArg, int);

                                RecBlk(msh, msh->buf, size_of_block);
                        }
                }
                else
                {
                        if(msh->typ & Asc)
                        {
                                for(i=0;i<kwd->SolSiz;i++)
                                        if(kwd->fmt[i] == 'r')
                                                fprintf(msh->hdl, "%.15lg ", va_arg(VarArg, double));
                                        else if(kwd->fmt[i] == 'n') {
                                                nb_repeat = va_arg(VarArg, int);
                                                fprintf(msh->hdl, "%d ", nb_repeat);
                                                for(k=0;k<nb_repeat;k++)
                                                        fprintf(msh->hdl, "%d ", va_arg(VarArg, int));
                                        }
                                        else
                                                fprintf(msh->hdl, "%d ", va_arg(VarArg, int));
                        }
                        else
                        {
                                pos = 0;

                                for(i=0;i<kwd->SolSiz;i++)
                                        if(kwd->fmt[i] == 'r')
                                        {
                                                DblBuf = (double *)&msh->buf[ pos ];
                                                *DblBuf = va_arg(VarArg, double);
                                                pos += 8;
                                        }
                                        else if(kwd->fmt[i] == 'n')
                                        {
                                                IntBuf = (int *)&msh->buf[ pos ];
                                                nb_repeat = va_arg(VarArg, int);
                                                *IntBuf = nb_repeat;
                                                pos += 4;
                                                for(k=0;k<nb_repeat;k++) {
                                                        IntBuf = (int *)&msh->buf[ pos ];
                                                        *IntBuf = va_arg(VarArg, int);
                                                        pos += 4;
                                                }
                                        }
                                        else
                                        {
                                                IntBuf = (int *)&msh->buf[ pos ];
                                                *IntBuf = va_arg(VarArg, int);
                                                pos += 4;
                                        }
                                RecBlk(msh, msh->buf, pos/4);
                        }
                }
        }
        else
        {
                if(msh->ver == 1)
                {
                        FltSolTab = va_arg(VarArg, float *);

                        if(msh->typ & Asc)
                                for(j=0;j<kwd->SolSiz;j++)
                                        fprintf(msh->hdl, "%g ", FltSolTab[j]);
                        else
                                RecBlk(msh, (unsigned char *)FltSolTab, kwd->NmbWrd);
                }
                else
                {
                        DblSolTab = va_arg(VarArg, double *);

                        if(msh->typ & Asc)
                                for(j=0;j<kwd->SolSiz;j++)
                                        fprintf(msh->hdl, "%.15lg ", DblSolTab[j]);
                        else
                                RecBlk(msh, (unsigned char *)DblSolTab, kwd->NmbWrd);
                }
        }

        va_end(VarArg);

        if(msh->typ & Asc)
                fprintf(msh->hdl, "\n");
}


/*----------------------------------------------------------*/
/* Private procedure for transmesh : copy a whole line          */
/*----------------------------------------------------------*/

void GmfCpyLin(int InpIdx, int OutIdx, int KwdCod)
{
        double d;
        float f;
        int i, a;
        GmfMshSct *InpMsh = GmfMshTab[ InpIdx ], *OutMsh = GmfMshTab[ OutIdx ];
        KwdSct *kwd = &InpMsh->KwdTab[ KwdCod ];

        for(i=0;i<kwd->SolSiz;i++)
        {
                if(kwd->fmt[i] == 'r')
                {
                        if(InpMsh->ver == 1)
                        {
                                if(InpMsh->typ & Asc)
                                        fscanf(InpMsh->hdl, "%f", &f);
                                else
                                        ScaWrd(InpMsh, (unsigned char *)&f);

                                d = f;
                        }
                        else
                        {
                                if(InpMsh->typ & Asc)
                                        fscanf(InpMsh->hdl, "%lf", &d);
                                else
                                        ScaDblWrd(InpMsh, (unsigned char *)&d);

                                f = (float)d;
                        }

                        if(OutMsh->ver == 1)
                                if(OutMsh->typ & Asc)
                                        fprintf(OutMsh->hdl, "%g ", f);
                                else
                                        RecWrd(OutMsh, (unsigned char *)&f);
                        else
                                if(OutMsh->typ & Asc)
                                        fprintf(OutMsh->hdl, "%.15g ", d);
                                else
                                        RecDblWrd(OutMsh, (unsigned char *)&d);
                }
                else if(kwd->fmt[i] == 'n')
                {
                        int k, nb_repeat = 0;

                        if(InpMsh->typ & Asc)
                                fscanf(InpMsh->hdl, "%d", &a);
                        else
                                ScaWrd(InpMsh, (unsigned char *)&a);

                        nb_repeat = a;

                        if(OutMsh->typ & Asc)
                                fprintf(OutMsh->hdl, "%d ", a);
                        else
                                RecWrd(OutMsh, (unsigned char *)&a);

                        for(k=0;k<nb_repeat;k++) {
                                if(InpMsh->typ & Asc)
                                        fscanf(InpMsh->hdl, "%d", &a);
                                else
                                        ScaWrd(InpMsh, (unsigned char *)&a);

                                if(OutMsh->typ & Asc)
                                        fprintf(OutMsh->hdl, "%d ", a);
                                else
                                        RecWrd(OutMsh, (unsigned char *)&a);
                        }
                }
                else
                {
                        if(InpMsh->typ & Asc)
                                fscanf(InpMsh->hdl, "%d", &a);
                        else
                                ScaWrd(InpMsh, (unsigned char *)&a);

                        if(OutMsh->typ & Asc)
                                fprintf(OutMsh->hdl, "%d ", a);
                        else
                                RecWrd(OutMsh, (unsigned char *)&a);
                }
        }

        if(OutMsh->typ & Asc)
                fprintf(OutMsh->hdl, "\n");
}


/*----------------------------------------------------------*/
/* Find every kw present in a meshfile                                          */
/*----------------------------------------------------------*/

static int ScaKwdTab(GmfMshSct *msh)
{
        int KwdCod;
        long  NexPos, CurPos, EndPos;
        char str[ GmfStrSiz ];

        if(msh->typ & Asc)
        {
                /* Scan each string in the file until the end */

                while(fscanf(msh->hdl, "%s", str) != EOF)
                {
                        /* Fast test in order to reject quickly the numeric values */

                        if(isalpha(str[0]))
                        {
                                /* Search which kwd code this string is associated with, 
                                        then get its header and save the curent position in file (just before the data) */

                                for(KwdCod=1; KwdCod<= GmfMaxKwd; KwdCod++)
                                        if(!strcmp(str, GmfKwdFmt[ KwdCod ][0]))
                                        {
                                                ScaKwdHdr(msh, KwdCod);
                                                break;
                                        }
                        }
                        else if(str[0] == '#')
                                while(fgetc(msh->hdl) != '\n');
                }
        }
        else
        {
                /* Get file size */

                CurPos = ftell(msh->hdl);
                fseek(msh->hdl, 0, SEEK_END);
                EndPos = ftell(msh->hdl);
                fseek(msh->hdl, CurPos, SEEK_SET);

                /* Jump through kwd positions in the file */

                do
                {
                        /* Get the kwd code and the next kwd position */

                        ScaWrd(msh, (unsigned char *)&KwdCod);
                        NexPos = GetPos(msh);

                        if(NexPos > EndPos)
                                return(0);

                        /* Check if this kwd belongs to this mesh version */

                        if( (KwdCod >= 1) && (KwdCod <= GmfMaxKwd) )
                                ScaKwdHdr(msh, KwdCod);

                        /* Go to the next kwd */

                        if(NexPos)
                                fseek(msh->hdl, NexPos, SEEK_SET);
                }while(NexPos && (KwdCod != GmfEnd));
        }

        return(1);
}


/*----------------------------------------------------------*/
/* Read and setup the keyword's header                                          */
/*----------------------------------------------------------*/

static void ScaKwdHdr(GmfMshSct *msh, int KwdCod)
{
        int i;
        KwdSct *kwd = &msh->KwdTab[ KwdCod ];

        if(!strcmp("i", GmfKwdFmt[ KwdCod ][2]))
        {
                if(msh->typ & Asc)
                        fscanf(msh->hdl, "%d", &kwd->NmbLin);
                else
                        ScaWrd(msh, (unsigned char *)&kwd->NmbLin);
        }
        else
                kwd->NmbLin = 1;

        if(!strcmp("sr", GmfKwdFmt[ KwdCod ][3]))
        {
                if(msh->typ & Asc)
                {
                        fscanf(msh->hdl, "%d", &kwd->NmbTyp);

                        for(i=0;i<kwd->NmbTyp;i++)
                                fscanf(msh->hdl, "%d", &kwd->TypTab[i]);
                }
                else
                {
                        ScaWrd(msh, (unsigned char *)&kwd->NmbTyp);

                        for(i=0;i<kwd->NmbTyp;i++)
                                ScaWrd(msh, (unsigned char *)&kwd->TypTab[i]);
                }
        }

        ExpFmt(msh, KwdCod);
        kwd->pos = ftell(msh->hdl);
}


/*----------------------------------------------------------*/
/* Expand the compacted format and compute the line size        */
/*----------------------------------------------------------*/

static void ExpFmt(GmfMshSct *msh, int KwdCod)
{
        int i, j, TmpSiz=0;
        char chr;
        const char *InpFmt = GmfKwdFmt[ KwdCod ][3];
        KwdSct *kwd = &msh->KwdTab[ KwdCod ];

        /* Set the kwd's type */

        if(!strlen(GmfKwdFmt[ KwdCod ][2]))
                kwd->typ = InfKwd;
        else if(!strcmp(InpFmt, "sr"))
                kwd->typ = SolKwd;
        else
                kwd->typ = RegKwd;

        /* Get the solution-field's size */

        if(kwd->typ == SolKwd)
                for(i=0;i<kwd->NmbTyp;i++)
                        switch(kwd->TypTab[i])
                        {
                                case GmfSca    : TmpSiz += 1; break;
                                case GmfVec    : TmpSiz += msh->dim; break;
                                case GmfSymMat : TmpSiz += (msh->dim * (msh->dim+1)) / 2; break;
                                case GmfMat    : TmpSiz += msh->dim * msh->dim; break;
                        }

        /* Scan each character from the format string */

        i = kwd->SolSiz = kwd->NmbWrd = 0;

        while(i < strlen(InpFmt))
        {
                chr = InpFmt[ i++ ];

                if(chr == 'd')
                {
                        chr = InpFmt[i++];

                        for(j=0;j<msh->dim;j++)
                                kwd->fmt[ kwd->SolSiz++ ] = chr;
                }
                else if(chr == 's')
                {
                        chr = InpFmt[i++];

                        for(j=0;j<TmpSiz;j++)
                                kwd->fmt[ kwd->SolSiz++ ] = chr;
                }
                else
                        kwd->fmt[ kwd->SolSiz++ ] = chr;
        }

        for(i=0;i<kwd->SolSiz;i++)
                if(kwd->fmt[i] == 'i')
                        kwd->NmbWrd++;
                else if(msh->ver >= 2)
                        kwd->NmbWrd += 2;
                else
                        kwd->NmbWrd++;
}


/*----------------------------------------------------------*/
/* Read a four bytes word from a mesh file                                      */
/*----------------------------------------------------------*/

static void ScaWrd(GmfMshSct *msh, unsigned char *wrd)
{
        unsigned char swp;

        fread(wrd, WrdSiz, 1, msh->hdl);

        if(msh->cod == 1)
                return;

        swp = wrd[3];
        wrd[3] = wrd[0];
        wrd[0] = swp;

        swp = wrd[2];
        wrd[2] = wrd[1];
        wrd[1] = swp;
}


/*----------------------------------------------------------*/
/* Read an eight bytes word from a mesh file                            */
/*----------------------------------------------------------*/

static void ScaDblWrd(GmfMshSct *msh, unsigned char *wrd)
{
        int i;
        unsigned char swp;

        fread(wrd, WrdSiz, 2, msh->hdl);

        if(msh->cod == 1)
                return;

        for(i=0;i<4;i++)
        {
                swp = wrd[7-i];
                wrd[7-i] = wrd[i];
                wrd[i] = swp;
        }
}


/*----------------------------------------------------------*/
/* Read ablock of four bytes word from a mesh file                      */
/*----------------------------------------------------------*/

static void ScaBlk(GmfMshSct *msh, unsigned char *blk, int siz)
{
        int i, j;
        unsigned char swp, *wrd;

        fread(blk, WrdSiz, siz, msh->hdl);

        if(msh->cod == 1)
                return;

        for(i=0;i<siz;i++)
        {
                wrd = &blk[ i * 4 ];

                for(j=0;j<2;j++)
                {
                        swp = wrd[ 3-j ];
                        wrd[ 3-j ] = wrd[j];
                        wrd[j] = swp;
                }
        }
}


/*----------------------------------------------------------*/
/* Read a 4 or 8 bytes position in mesh file                            */
/*----------------------------------------------------------*/

static long GetPos(GmfMshSct *msh)
{
        int IntVal;
        long pos;

        if(msh->ver >= 3)
                ScaDblWrd(msh, (unsigned char*)&pos);
        else
        {
                ScaWrd(msh, (unsigned char*)&IntVal);
                pos = IntVal;
        }

        return(pos);
}


/*----------------------------------------------------------*/
/* Write a four bytes word to a mesh file                                       */
/*----------------------------------------------------------*/

static void RecWrd(GmfMshSct *msh, unsigned char *wrd)
{
        fwrite(wrd, WrdSiz, 1, msh->hdl);
}


/*----------------------------------------------------------*/
/* Write an eight bytes word to a mesh file                                     */
/*----------------------------------------------------------*/

static void RecDblWrd(GmfMshSct *msh, unsigned char *wrd)
{
        fwrite(wrd, WrdSiz, 2, msh->hdl);
}


/*----------------------------------------------------------*/
/* Write a block of four bytes word to a mesh file                      */
/*----------------------------------------------------------*/

static void RecBlk(GmfMshSct *msh, unsigned char *blk, int siz)
{
        /* Copy this line-block into the main mesh buffer */

        if(siz)
        {
                memcpy(&msh->blk[ msh->pos ], blk, siz * WrdSiz);
                msh->pos += siz * WrdSiz;
        }

        /* When the buffer is full or this procedure is called with a 0 size, flush the cache on disk */

        if( (msh->pos > BufSiz) || (!siz && msh->pos) )
        {
                fwrite(msh->blk, 1, msh->pos, msh->hdl);
                msh->pos = 0;
        }
}


/*----------------------------------------------------------*/
/* Write a 4 or 8 bytes position in a mesh file                         */
/*----------------------------------------------------------*/

static void SetPos(GmfMshSct *msh, long pos)
{
        int IntVal;

        if(msh->ver >= 3)
                RecDblWrd(msh, (unsigned char*)&pos);
        else
        {
                IntVal = pos;
                RecWrd(msh, (unsigned char*)&IntVal);
        }
}