/*******************************************************************************
*									       *
* mcf_NTuIOFiles.c -- Utilities to manipulate files within the MCFIO Gen.      *
*        				Ntuple schema                          *
*									       *
*	P. Lebrun, September 1995.					       *
*									       *
*******************************************************************************/
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <sys/param.h>
#include <limits.h>
#include <time.h>
#include <rpc/types.h>
#include <sys/types.h>
#include <rpc/xdr.h>
#ifdef _HPUX_SOURCE
#include <unistd.h>
#endif
#include "mcf_nTupleDescript.h"
#include "mcf_xdr.h"
#include "mcfio_Dict.h"
#include "mcfio_Direct.h"
#include "mcf_NTuIOFiles.h"
#include "mcf_NTuIOUtils.h"
#include "mcf_ntubld_db.h"
#ifndef False
#define False 0
#endif
#ifndef True
#define True 1
#endif

extern char *VarTypesNamesF77[N_VAR_TYPES];
extern char *VarTypesNamesC[N_VAR_TYPES];

extern struct line_title_c line_title_c_;
extern struct header_c header_c_;
extern struct variable_c variable_c_;

/*
** Ntuple identifiers list, initialized here and in mcfio_Util1
*/
nTuDDL **NTuDDLList = NULL;
int NumOfNTuples = 0;
bool_t McfNTuPleSaveDecoding = True;

static char *makeStructName(char *title, int orgStyle);
static size_t nDatVariable(varGenNtuple *varTmp);
static size_t sizeVariable(varGenNtuple *varTmp);
static char *mcf_copyNtrim(char *fromString);

int mcfioC_DeclareNtuple(int uid, char *title, char *category, 
                                int stream, char *filename)
{
    nTuDDL *ddl, *ddlRef;
    int i, l, j, jstr, dejaVu, id, **ip, **ipo, **ipt;
    
    
    if ((stream < 1) || (stream > MCF_STREAM_NUM_MAX)) {
     fprintf(stderr,
  " mcfio_NtupleDDLRead: Illegal MCFIO stream number.\n"); 
     return -1;
    }
    jstr = stream-1;
    if (McfStreamPtrList[jstr] == NULL) { 
        fprintf(stderr,
   " mcfio_DeclareNtuple: First, declare the stream by calling mcfio_Open...\n"); 
     return -1;
    }
      
    if (McfStreamPtrList[jstr]->row != MCFIO_WRITE) {
        fprintf(stderr,
   " mcfio_DeclareNtuple: You must declare an Ntuple for an Output Stream\n"); 
     return -1;
    }
        
    if (!mcf_CheckValidCat(category, False)) return 0;
    
    /* Check that this item characterized by uid/Category has not already been 
         created. If so, do not create a new one.  If associated to the same
         stream, flag this as an error. */
    
    id = mcf_NTuId(uid, category);			
    if (id != -1) {
       ddl = mcf_GetNTuByPtrID(id);
       if (ddl->streamId == stream) {  
        fprintf(stderr,
        "Mcfio Declare Ntuple:  An item with this uid/Category already exists.\n");
        fprintf(stderr, "  uid = %d, Category = %s, ", uid, category);
        fprintf(stderr, "Ntuple not created.\n");
        return -1;
        }
     }
     /*
     ** May be this dbin template has already been digested. If so, refer
     ** to it, to avoid re-computing all the offsets.
    */
     ip = (int **) NTuDDLList;
     for (i=0, dejaVu=False; i< NumOfNTuples; i++, ip++) {
	 ddlRef = (nTuDDL *) *ip;
         if (strcmp(filename, ddlRef->dbinFileName) == 0) {
             dejaVu = True;
            /* Create a holder for this Ntuple Description */
             ddl = (nTuDDL * ) malloc(sizeof(nTuDDL));
             /*
             ** back up in the linked list if need be, until we 
             ** a fully documented descriptor.
             */
             while (ddlRef->descrNtu == NULL) ddlRef = ddlRef->reference;
             ddl->reference = ddlRef; ddl->descrNtu = NULL;
             break;
         }
     }
     if (dejaVu == False)  {
         ddl = mcf_GetFileNTuDDL(filename);
         if (ddl == NULL) { 
              fprintf(stderr,
                  " mcfio_NtupleDDLRead: Error reading %s\n", filename );
                 return -1; 
         }
         ddl->reference = NULL;
    }
    ddl->title = mcf_copyNtrim(mcf_ValidStr(title, NTU_MAX_TITLE_LENGTH, 
    				 "title"));
    if (category == NULL)
      ddl->category = 
         mcf_copyNtrim(mcf_ValidStr(category, NTU_MAX_CATEGORY_LENGTH, 
    				 "category"));
    else {
	ddl->category = mcf_copyNtrim(category);
    }
    AddNTuDDLtoList(ddl);
/*
** Now we compute the offssets. 
*/
    if (dejaVu == False)  {
        mcf_ComputeNTuOffsets(ddl);    
/*
** Now we compute the lengths..
*/
       mcf_ComputeNTuLengths(ddl);
    }
    ddl->uid = uid;
    ddl->streamId = stream;
    /*
    ** Set the sequential id for this particular stream
    */
    for (i=0, j=0; i<NumOfNTuples; i++) 
        if (NTuDDLList[i]->streamId == ddl->streamId) j++;
    ddl->seqNTuId = j;
    (McfStreamPtrList[jstr]->fhead->nNTuples)++;
    return ddl->id;
}
     
int mcfioC_EndDeclNTuples(int stream)
/*
** Routine to end theNtuple delcaration and rewrite the beginning of the 
** file.
*/
{
   int i, j, jstr, idtmp, ntot;
   u_int p1, p2;
   FILE *ff;
   mcfStream *str;
   
  if (McfStreamPtrList == NULL) { 
     fprintf(stderr,
  " mcfio_EndDeclNtuple: No stream open, No inialization.\n"); 
     return -1;
  }
  jstr = stream-1;
  if (McfStreamPtrList[jstr] == NULL) { 
     fprintf(stderr,
 " mcfio_EndDeclNtuple: First, declare the stream by calling mcfio_Open...\n"); 
     return -1;
  }
  str = McfStreamPtrList[jstr];
  if (str->row != MCFIO_WRITE) {
     fprintf(stderr,
 " mcfio_EndDeclNtuple: This routine is not applicable to Input streams...\n"); 
     return -1;
  }
  if (str->fhead->nNTuples < 1) {
     fprintf(stderr,
 " mcfio_EndDeclNtuple: No Ntuple declared for this stream...\n"); 
     return 0;
  }     
   /*
   ** Now we can try toto complete the file header. As it is now bigger, 
   ** and it is the first structure written, it is easier to start over.
   ** Destroy the XDR stream, close the file, and reopen it.
   */  
   xdr_destroy(str->xdr);
   fclose(str->filePtr);
   remove(str->filename);
   ff = fopen(str->filename, "w");
   if (ff == NULL) {
     fprintf(stderr,
  " mcfio_EndDeclNtuple: Problem re-opening file %s, message \n", 
        str->filename);
     return -1;
   }
   xdrstdio_create(str->xdr, ff, XDR_ENCODE);
   p1 = xdr_getpos(str->xdr);
   str->firstPos = p1;
   str->currentPos = p1;
   /*
   ** In the file header, we do not store the NTuple Ids, as they are 
   ** not necessarily valid in an other context, where we have different
   ** streams/NTuples combinations. The SeqNTuId are trivial, 
   ** within a stream, at the file header (1,2,3,..) 
   **  But, of course, we must provide an array for the event header..
   */
   str->ehead->dimNTuples = str->fhead->nNTuples;
   str->ehead->nNTuples = 0;
   str->ehead->nTupleIds = 
          (int *) malloc(sizeof(int) * str->fhead->nNTuples);
   
   str->ehead->ptrNTuples = 
          (u_int *) malloc(sizeof(u_int) * str->fhead->nNTuples);
   for (i=0; i<str->ehead->dimNTuples; i++) str->ehead->ptrNTuples[i]=0;       
   
   str->status = MCFIO_BOF;
   if (mcfioC_Wrtfhead(str, INITIATE) == FALSE){
       mcfioC_FreeStream(&McfStreamPtrList[jstr]);
       fclose(ff);
       return -1;
   }
   /*
   ** Write the first dummy table 
   */
   if (mcfioC_Wrttable(str, INITIATE) == FALSE) return -1;
   /*
   ** Write the first dummy event header
   */
   if (mcfioC_WrtEvt(str, INITIATE) == FALSE) return -1;
   str->ehead->evtnum = 0;
   str->status = MCFIO_RUNNING;
   return (str->fhead->nNTuples);
}

nTuDDL *mcf_GetFileNTuDDL(char*filename)
{
    nTuDDL *ddl;
    int i, l, j;
    char *text, *tc;
    varGenNtuple *varTmp;
    descrGenNtuple *dNTu;
    
    /* Create a holder for this Ntuple Description */
    ddl = (nTuDDL * ) malloc(sizeof(nTuDDL));
    ddl->dbinFileName = (char *) malloc(sizeof(char) * (strlen(filename) +1));
    strcpy(ddl->dbinFileName, filename);
    ddl->descrNtu = (descrGenNtuple *) malloc(sizeof(descrGenNtuple));
    dNTu = ddl->descrNtu;
    
    header_c_.n_obj_header = 0;
    line_title_c_.n_obj_line_title = 0;
    mcf_ntubldRead(filename);
    if ((line_title_c_.n_obj_line_title < 1)  ||
        (header_c_.n_obj_header != 1)) {
        fprintf(stderr, 
                " This file was not created by the ntuBuild aplication!");
        return NULL;
    }
    if (strcmp(line_title_c_.line_title[0].line,
               "ntuBuild Database, v1.0") != 0) {
        fprintf(stderr, 
                " This file was not created by a wrong version of ntuBuild!");
        return NULL;
    }
    /*
    ** There are 80 character per lines in dbin..
    */
    text = (char *)
            malloc(sizeof(char) * 80 * (line_title_c_.n_obj_line_title -1));
    for (i=1, tc=text; i<line_title_c_.n_obj_line_title; i++) {
        strcpy(tc, line_title_c_.line_title[i].line); 
        tc += strlen(line_title_c_.line_title[i].line);
        *tc = '\n'; tc++;
    }
    *tc = '\0';
    dNTu->description = text;
    
    l = strlen(header_c_.header[0].title);
    dNTu->title = (char *) malloc(sizeof(char) * (l+1));
    strcpy(dNTu->title, header_c_.header[0].title);
    
    strcpy(dNTu->version, header_c_.header[0].version);
    
    strcpy(dNTu->nameIndex, header_c_.header[0].namemaxindex);
    
    dNTu->maxMultiplicity = header_c_.header[0].maxmult;
    
    dNTu->orgStyle = header_c_.header[0].orgstyle;
    dNTu->numVariables = header_c_.header[0].nvar;
    dNTu->numAvailable = dNTu->numVariables;
    dNTu->variables =
     (varGenNtuple **) malloc(sizeof(varGenNtuple *) * dNTu->numVariables);
    /*
    ** Now the variables
    */
    for (i=0; i<variable_c_.n_obj_variable; i++) {
        dNTu->variables[i] =   
                (varGenNtuple *) malloc(sizeof(varGenNtuple));
        varTmp = dNTu->variables[i];
        varTmp->nameBlank = False;
        varTmp->name = (char *)
            malloc(sizeof(char) * (strlen(variable_c_.variable[i].name) + 1));
        strcpy(varTmp->name, variable_c_.variable[i].name);
        
        if ((strlen(variable_c_.variable[i].description) > 1) ||
                    variable_c_.variable[i].description[0] != ' ') { 
           varTmp->description = (char *) malloc(sizeof(char) * 
                (strlen(variable_c_.variable[i].description) + 1));
           strcpy(varTmp->description, variable_c_.variable[i].description);
        } else varTmp->description = NULL;  
        varTmp->type = variable_c_.variable[i].type;
        varTmp->isFixedSize = True;
        if (strncmp(variable_c_.variable[i].isfixedsize,"Yes",3))
            varTmp->isFixedSize = False;
        varTmp->numDim = variable_c_.variable[i].numdim;
        if (varTmp->numDim > 0) 
           for (j=0; j< varTmp->numDim; j++)
               varTmp->dimensions[j] = variable_c_.variable[i].dimensions[j];
                
    }
    /*
    ** Set the ordering. Trivial in this case, it has been ordered in 
    ** the save routine.
    */
    dNTu->varOrdering = (int *) malloc(sizeof(int) * dNTu->numAvailable);
    for (i=0; i<dNTu->numVariables; i++)
       dNTu->varOrdering[i] = i; 
    dNTu->subOffset = NULL;
    dNTu->subXDROffset = NULL;    
    return ddl;
    
}

/*
** Compute the offsets by writing a simple program, stand alone, that uses
** the d/s
*/
void    mcf_ComputeNTuOffsets(nTuDDL *ddl)    
{
     char tmpName[128], *tc, *tc1, *nameCom;
     int i, j, l, firstIndexed, nDat;
     char filenameInclude[128], filenameProgram[128], filenameExec[128];
     char filenameData[128], nameMaxIndex[32];
     char line[256];
     void **ptrBegVar;
     varGenNtuple *varTmp;
     descrGenNtuple *dNTu;
     FILE *Ffp;
     
     dNTu =ddl->descrNtu;  
     for (i=0; i< dNTu->numVariables; i++) { 
        varTmp = dNTu->variables[i];
        varTmp->offset = 0;
     }
        
     memset(tmpName, 0, 127);
     tc = tmpName;
     sprintf(tc, "tmp_%s_XXXXXX", ddl->descrNtu->title);
     tc1 = mktemp(tmpName);
     if (tc1 == NULL) {
         fprintf(stderr, 
        " Can not compose a tempoary name in mcf_ComputeOffsets!");
         return;
     }
     sprintf(filenameInclude, "%s.h", tc1);
     sprintf(filenameProgram, "%s.c", tc1);
     sprintf(filenameData, "%s.dat", tc1);
     strcpy(filenameExec, tc1);
     mcf_ComposeDoth(ddl->descrNtu, filenameInclude);
/*
** Compose a little moronic program that establishes the addresses of all 
** variables. There might be a better way, though.. However, this ought to be
** safe. 
*/
     Ffp =  fopen( filenameProgram, "w");
     fprintf(Ffp, "#include <stdio.h>\n");
     fprintf(Ffp, "#include \"%s\"\n",filenameInclude);
     if (dNTu->orgStyle == PARALLEL_ARRAY_NTU)
         fprintf(Ffp, "#define NUM_VAR %d\n", (dNTu->numVariables+3));
     else 
         fprintf(Ffp, "#define NUM_VAR %d\n",
                (dNTu->numVariables + 3 + dNTu->maxMultiplicity) );
         
     nameCom = makeStructName(dNTu->title, dNTu->orgStyle);
     
     fprintf(Ffp, "%s_struct tmpStruct; \n", nameCom);
     fprintf(Ffp, "main(int argc, char **argv)\n");
     fprintf(Ffp, "{\n");
     fprintf(Ffp, "    void *ptrBegVar[NUM_VAR];\n");
     fprintf(Ffp, "    FILE *Ffp;\n");
     fprintf(Ffp, "    int i;\n");
     fprintf(Ffp, "\n");
     fprintf(Ffp, "    ptrBegVar[0] = (void *) &tmpStruct.version[0];\n");
     fprintf(Ffp,
      "    ptrBegVar[1] = (void *) &tmpStruct.%s;\n",dNTu->nameIndex);
     for(i=0, firstIndexed=-1; i<dNTu->numVariables; i++) {
           if (dNTu->variables[i]->isFixedSize == False) {
                  firstIndexed = i; break;
           }       
     }
     dNTu->firstIndexed = firstIndexed;             
     if (dNTu->orgStyle == PARALLEL_ARRAY_NTU) {
         for(i=0; i<dNTu->numVariables; i++) {
           varTmp =  dNTu->variables[i];
           /*
           ** Assume that all the variables are properly 
           ** defined at this stage (e..g, coming from a valid DDL dbin file)
           ** and in order
           */
           tc = line;
           if ((varTmp->numDim == 0) && (varTmp->isFixedSize == True)) 
               sprintf(tc,
                "    ptrBegVar[%d] = (void *) &tmpStruct.%s%n",
                    (i+2), varTmp->name, &l);
           else           
               sprintf(tc,
                "    ptrBegVar[%d] = (void *) tmpStruct.%s%n",
                    (i+2), varTmp->name, &l);
           tc+=l;
           fprintf(Ffp, "%s;\n", line);    
         }
         fprintf(Ffp,
      "    ptrBegVar[%d] = (void *) tmpStruct.fence;\n",dNTu->numVariables+2);
     } else {
         for(i=0; i<dNTu->numVariables; i++) {
           varTmp =  dNTu->variables[i];
           tc = line;
           if (varTmp->isFixedSize == True) {
              if (varTmp->numDim == 0)
                sprintf(tc,
               "    ptrBegVar[%d] = (void *) &tmpStruct.%s%n",
                    (i+2), varTmp->name, &l);
               else 
                sprintf(tc,
               "    ptrBegVar[%d] = (void *) tmpStruct.%s%n",
                    (i+2), varTmp->name, &l);
           } else {
               if (varTmp->numDim == 0)         
                   sprintf(tc,
                   "    ptrBegVar[%d] = (void *) &tmpStruct.var[0].%s%n",
                    (i+2), varTmp->name, &l);
               else 
                   sprintf(tc,
                   "    ptrBegVar[%d] = (void *) tmpStruct.var[0].%s%n",
                    (i+2), varTmp->name, &l);
           }                                     
           fprintf(Ffp, "%s;\n", line); 
         }
         tc1 = dNTu->nameIndex;
         strcpy(nameMaxIndex, tc1);
         l = strlen(tc1); 
         if (l > 26) {
            strncpy(nameMaxIndex, tc1, 26);
            sprintf(&nameMaxIndex[26],"_max");
         } else
            sprintf(nameMaxIndex, "%s_max", tc1);
         fprintf(Ffp,"    for (i=0; i<%s; i++) \n", nameMaxIndex);
         tc = line;
         if (firstIndexed != -1) {
             varTmp = dNTu->variables[firstIndexed];
             sprintf(tc,
             "       ptrBegVar[i+%d] = (void *) &tmpStruct.var[i].%s%n",
              (2+dNTu->numVariables), varTmp->name, &l); tc+=l;
              if (varTmp->numDim > 0) for (j=0; j<varTmp->numDim; j++, tc+=l) 
                                             sprintf(tc, "[0]%n", &l);   
              fprintf(Ffp, "%s;\n", line);
         }      
         fprintf(Ffp,
      "    ptrBegVar[%d] = (void *) tmpStruct.fence;\n",
              dNTu->numVariables+2+dNTu->maxMultiplicity);
     }
     fprintf(Ffp, " ");
     fprintf(Ffp,"    Ffp = fopen(\"%s\",\"w\");\n",filenameData);
     fprintf(Ffp,"    fwrite((void *) ptrBegVar, sizeof(void *),\
(size_t) NUM_VAR, Ffp);\n");
     fprintf(Ffp,"    fclose(Ffp);\n");
     fprintf(Ffp,"}\n");
     fclose(Ffp);
     free(nameCom);
     /*
     ** Now compile, link and load this exec, read the result
     */
#ifdef _HPUX_SOURCE
     sprintf(line,"cc -Aa -D_HPUX_SOURCE -o %s %s", 
                            filenameExec, filenameProgram);
#else
     sprintf(line,"cc -o %s %s", filenameExec, filenameProgram);
#endif
     system(line);
     sprintf(line,"./%s", filenameExec);
     system(line);
     if (dNTu->orgStyle == PARALLEL_ARRAY_NTU) nDat = dNTu->numVariables+3;
        else nDat = dNTu->numVariables+3+dNTu->maxMultiplicity;
     if (firstIndexed == -1) nDat = dNTu->numVariables+3;   
     ptrBegVar = (void **) malloc (sizeof(void *) * (nDat));
     Ffp = fopen(filenameData, "r");
     fread((void *) ptrBegVar, sizeof(void *), (size_t) nDat, Ffp);
     fclose(Ffp);
     /*
     ** remove garbage files..
     */
     remove(filenameData); remove(filenameProgram); remove(filenameExec);
     remove(filenameInclude);
     /*
     ** Convert these addresses to offsets
     */
     dNTu->multOffset =  ((long) ptrBegVar[1] - (long) ptrBegVar[0]);
     if (dNTu->orgStyle == PARALLEL_ARRAY_NTU) {
         dNTu->fenceOffset =
          ((long) ptrBegVar[dNTu->numVariables+2] - (long) ptrBegVar[0]); 
         for (i=0; i< dNTu->numVariables; i++) 
            dNTu->variables[i]->offset = 
                   ((long) ptrBegVar[i+2] - (long) ptrBegVar[0]);
     } else {
         for (i=0; i< dNTu->numVariables; i++) {
            varTmp = dNTu->variables[i]; 
            if (varTmp->isFixedSize)
                varTmp->offset = 
                    ((long) ptrBegVar[i+2] - (long) ptrBegVar[0]);
            else 
                varTmp->offset = 
                   ((long) ptrBegVar[i+2] - (long)ptrBegVar[firstIndexed+2]);
         }
         if (dNTu->subOffset != NULL) free(dNTu->subOffset);
         dNTu->subOffset =
            (long *) malloc(sizeof(long) * dNTu->maxMultiplicity);
         if (firstIndexed != -1) {    
            for (i=0; i<dNTu->maxMultiplicity; i++) 
                dNTu->subOffset[i] =          
               ((long) ptrBegVar[i+2+dNTu->numVariables]  - 
               (long) ptrBegVar[0]);
         }          
         dNTu->fenceOffset =
            ((long) ptrBegVar[dNTu->numVariables+2+dNTu->maxMultiplicity]
                     - (long) ptrBegVar[0]); 
     }   
     free(ptrBegVar);
}

/*
** Compute the lengths for the XDR Array statements. It is assumed that the 
** NTUple descriptor is sorted, no blank variables.
*/
void    mcf_ComputeNTuLengths(nTuDDL *ddl)    
{
     int i, j, lastTmp, sameType;
     size_t nDat, sizeItem;
     varGenNtuple *var1, *var2;
     descrGenNtuple *dNTu;
     
     dNTu =ddl->descrNtu;
     if (dNTu->firstIndexed != -1) lastTmp = dNTu->firstIndexed;
         else lastTmp = dNTu->numVariables;
     /*
     ** fixed size first.. 
     */   
     for (i=0; i<lastTmp; i++)
         dNTu->variables[i]->lengthW = nDatVariable(dNTu->variables[i]); 
/*
** This, in principle, is the optimized version, where we collaps single
** fields of the same type into an array. However, this is machine 
** dependant.
*/                                    
     for (i=0; i<lastTmp; i++) {
         var1 = dNTu->variables[i];
         if (var1->lengthW != 0) {
            nDat = nDatVariable(var1);
            j=i+1;
            sizeItem = sizeVariable(var1); 
            sameType = True;
            while ((j<lastTmp) && (sameType)) { 
               var2 = dNTu->variables[j];
               if (var2->type != var1->type) sameType = False;
               if (sameType  && ((( var2->offset - 
                                    var1->offset)/sizeItem) ==
                                     nDat)) {
                  nDat += nDatVariable(var2);
                  var2->lengthW = 0; j++;
               }
            }
            var1->lengthW = nDat;
            var1->lengthB = nDat*sizeItem;
         }
     } 
     /*
     ** The variable size, similar code. This fill is very simple if the 
     ** if the organisation is parallel arrays, as we can not implmenent
     ** compaction 
     */           
     if (dNTu->firstIndexed == -1) return;
     if (dNTu->orgStyle == PARALLEL_ARRAY_NTU) {
         for (i=dNTu->firstIndexed; i<dNTu->numVariables; i++) { 
              dNTu->variables[i]->lengthW
                          =  nDatVariable(dNTu->variables[i]);
              dNTu->variables[i]->lengthB = dNTu->variables[i]->lengthW   
                         * sizeVariable(dNTu->variables[i]);
         }                
     } else {
         for (i=dNTu->firstIndexed; i<dNTu->numVariables; i++) 
            dNTu->variables[i]->lengthW =  nDatVariable(dNTu->variables[i]);
         for (i=dNTu->firstIndexed; i<dNTu->numVariables; i++) {
             var1 = dNTu->variables[i];
             if (var1->lengthW != 0) {
                 nDat = nDatVariable(var1);
                 j=i+1;
                 sizeItem = sizeVariable(var1); 
                 sameType = True;
                 while ((j<dNTu->numVariables) && (sameType)) { 
                    var2 = dNTu->variables[j];
                    if (var2->type != var1->type) sameType = False;
                    if (sameType  && (((var2->offset - 
                                   var1->offset)/sizeItem) ==
                                     nDat)) {
                     nDat += nDatVariable(var2);
                     var2->lengthW = 0; j++; 
                  }
                }
                var1->lengthW = nDat;
                var1->lengthB = nDat*sizeItem;
             }
         }
     }
     
}
/*
** Compute, in size_t units (bytes, I hope) the length of a particular
** variable.  Only the fixed size part, we will have to multiplity 
** by the multiplicty in the XDR filter.
*/

static size_t nDatVariable(varGenNtuple *var)
{
   size_t n;
   int i;
    
    n=1;
    for (i=0; i<var->numDim; i++) n = n * var->dimensions[i];
    return n;
}
static size_t sizeVariable(varGenNtuple *var)
{
   size_t n;
   
   switch (var->type) {
        case BYTE_NTU: case CHARACTER_NTU:
           n = sizeof(char);
           break;
        case INTEGER2_NTU:
           n = sizeof(short);
           break;
        case LOGICAL_NTU: case INTEGER_NTU:
           n = sizeof(int);
           break;
        case REAL_NTU:
           n = sizeof(float);
           break;
        case DBL_PRECISION_NTU:
           n = sizeof(double);
           break;
        case COMPLEX_NTU:
           n = 2 * sizeof(float);
           break;
        case DBL_COMPLEX_NTU:
           n = 2 * sizeof(double);
           break;
        case POINTER_NTU:
           n = sizeof(void *);
           break;
        default : 
           fprintf(stderr, " mcf_ComputNTuLength, internal error \n");
           n = 0;
           break;
    }
    return n;
}

/*
** Compose the .h file. Called from NTuBldMenu and this file. The structure 
** is assumed valid. 
*/
void    mcf_ComposeDoth(descrGenNtuple *dNTu, char *filename)
{
    char *nameCom, line[FILENAME_MAX+500], *tmp, *version, *text, *tc, *tc2;
    char nameMaxIndex[32], nameTmpIndex[32];
    char nullDescr[4], *descrTmp;
    int i, j, l, kmode, nc, ncTot, nl, iv;
    time_t clock;
    FILE *Ffp;
    varGenNtuple *var;
    
    nameCom = makeStructName(dNTu->title, dNTu->orgStyle);
    strcpy(nullDescr, "? ");
    strcpy(line, filename);
    tc = strchr(line, '.');
    if (tc == NULL) {
         l = strlen(filename);
         tc = line; tc+=l;
    }     
    strcpy(tc,".h");   
    Ffp = fopen(line, "w");
    fprintf(Ffp,"/* ntuBuild\n");
    time(&clock);
    tmp = line; sprintf(tmp,"** Creation Date : %n", &l); tmp += l;
    strncpy(tmp,ctime(&clock), 24); tmp += 24; *tmp='\n'; tmp++; *tmp = '\0';
    fprintf(Ffp,line);
    fprintf(Ffp,"**  User Comments\n");
    text = dNTu->description;
    tc = text; 
    if (*tc == '\0') 
       fprintf(Ffp,"** no user comments\n");
    else {
       ncTot = strlen(tc); nc =0;
       while (nc < ncTot) {
            tc2 = strchr(tc,'\n');
            nl = (int) (tc2-tc)/sizeof(char);
            if ((tc2 == NULL) || (nl > 75)) nl = 75;
            strncpy(line, tc, nl); line[nl] = '\0';
            fprintf (Ffp,"** %s\n", line);
            tc += nl; nc += nl;
            if (*tc == '\n') {
               tc++;
               nc++;
            }
       }
    }
    fprintf(Ffp,"*/    \n");
    version = dNTu->version;
    text = dNTu->nameIndex;
    strcpy(nameTmpIndex, text);
    l = strlen(text); 
    if (l > 26) {
            strncpy(nameMaxIndex, text, 26);
            sprintf(&nameMaxIndex[26],"_max");
    } else
            sprintf(nameMaxIndex, "%s_max", text);
    fprintf(Ffp,"#define %s %d\n", nameMaxIndex, dNTu->maxMultiplicity);        
    if (dNTu->orgStyle == PARALLEL_ARRAY_NTU) {
        fprintf(Ffp, "typedef struct _%s_struct {\n", nameCom);
        /*
        ** The first 64 bits contain the version token, as a char[8] string
        ** floowed by the multiplicty variable, followed by an integer pad
        */ 
        fprintf(Ffp,"    char version[8]; /* Version token */\n");
        fprintf(Ffp,
        "    int %s; /* Generalized Ntuple Multiplicity value */ \n",
                                              nameTmpIndex);
        fprintf(Ffp,
        "    int padding; /* Padding for 64 bit architecture */ \n");
        for (iv=0; iv< dNTu->numVariables; iv++) {
           for (j=0; j<dNTu->numAvailable; j++)  
                if (dNTu->varOrdering[j] == iv)  i = j; 
           var = dNTu->variables[i];
           kmode = 0; if (var->isFixedSize != True) kmode = 1;
           if (var->description == NULL) descrTmp = nullDescr;
                  else descrTmp = var->description;
           tc = line;
           if ((var->type != COMPLEX_NTU) &&
               (var->type != DBL_COMPLEX_NTU)) { 
                  sprintf(tc,"    %s %n", VarTypesNamesC[var->type], &l);
                  tc +=l;            
                  if ((var->numDim == 0) && (kmode ==0))
                            sprintf(tc," %s; /* %s */",
                                              var->name, descrTmp);
                  else if (var->numDim == 0) {
                        sprintf(tc," %s[%s]; /* %s */",
                        var->name, nameMaxIndex, descrTmp); 
                  } else { 
                       sprintf(tc," %s%n",var->name, &l); tc+=l;
                       if (kmode == 1) {
                           sprintf(tc, "[%s]%n", nameMaxIndex, &l);
                           tc +=l;
                   }    
                   for (j=var->numDim-1; j>-1; j--, tc+=l) 
                             sprintf(tc,"[%d]%n", var->dimensions[j], &l); 
                     
                        sprintf (tc,"; /* %s */", descrTmp);
                  }  
               } else { /* got to convert to float or dbl */
                    if (var->type == COMPLEX_NTU)
                       sprintf(tc,"    float %n", &l);
                       
                     else if (var->type == DBL_COMPLEX_NTU)
                       sprintf(tc,"    double %n", &l);
                       
                     tc +=l;            
                     if ((var->numDim == 0) && (kmode ==0)) 
                          sprintf(tc," %s[2]; /* %s */", var->name, descrTmp);
                     else if (var->numDim == 0) {
                       sprintf(tc," %s[%s][2]; /* %s */",
                        var->name, nameMaxIndex, descrTmp); 
                     } else { 
                       sprintf(tc," %s%n",var->name, &l); tc+=l;
                       if (kmode == 1) {
                          sprintf(tc, "[%s]%n", nameMaxIndex, &l);
                         tc +=l;
                       }    
                       for (j=var->numDim-1; j>-1; j--, tc+=l) 
                             sprintf(tc,"[%d]%n", var->dimensions[j], &l); 
                       sprintf (tc,"[2]; /* %s */", descrTmp);
                    }
               }
               fprintf(Ffp,"%s\n", line);
           }
        fprintf(Ffp,"    int fence[2]; \n");
        fprintf(Ffp,"} %s_struct; \n", nameCom);
     }else { 
     /*
     ** The other type of organisation, using structure
     */
        fprintf(Ffp, "typedef struct _%s_v_struct{\n", nameCom);
        for (iv=0; iv< dNTu->numVariables; iv++) {
           for (j=0; j<dNTu->numAvailable; j++)  
                if (dNTu->varOrdering[j] == iv)  i = j; 
           var = dNTu->variables[i];
           if (var->isFixedSize == False) {
               tc = line;
               if (var->type == COMPLEX_NTU)
                       sprintf(tc,"    float %n", &l);
               else if (var->type == DBL_COMPLEX_NTU)
                       sprintf(tc,"    double %n", &l);
               else 
                     sprintf(tc,"    %s %n", VarTypesNamesC[var->type], &l);
               tc +=l;            
               sprintf(tc," %s%n",var->name, &l); tc+=l;
               if (var->numDim != 0) {
                    for (j=var->numDim-1; j>-1; j--, tc+=l) 
                           sprintf(tc,"[%d]%n", var->dimensions[j], &l); 
               }
               if ((var->type == COMPLEX_NTU) ||
                   (var->type == DBL_COMPLEX_NTU)) {
                       sprintf (tc,"[2]%n",&l);
                       tc += l;
               }       
               if (var->description == NULL) descrTmp = nullDescr;
                    else descrTmp = var->description;
               sprintf(tc,"; /* %s */%n", descrTmp, &l); tc += l;
               fprintf(Ffp,"%s\n", line);
            }
        }
        fprintf(Ffp,"} %s_v_struct; \n", nameCom);
        fprintf(Ffp,"/* ----- */  \n", line);
        /*
        ** the mother structure now
        */
        fprintf(Ffp, "typedef struct _%s_struct{\n", nameCom);
        fprintf(Ffp,"    char version[8]; /* Version token */\n");
        fprintf(Ffp,
        "    int %s; /* Generalized Ntuple Multiplicity value */ \n",
                                              nameTmpIndex);
        fprintf(Ffp,
        "    int padding; /* Padding for 64 bit architecture */ \n");
        for (iv=0; iv< dNTu->numVariables; iv++) {
           for (j=0; j<dNTu->numAvailable; j++)  
                if (dNTu->varOrdering[j] == iv)  i = j; 
           var = dNTu->variables[i];
           if (var->isFixedSize == True) {
              tc = line;
               if (var->type == COMPLEX_NTU)
                       sprintf(tc,"    float %n", &l);
               else if (var->type == DBL_COMPLEX_NTU)
                       sprintf(tc,"    double %n", &l);
               else 
                     sprintf(tc,"    %s %n", VarTypesNamesC[var->type], &l);
               tc +=l;            
               sprintf(tc," %s%n",var->name, &l); tc+=l;
               if (var->numDim != 0) {
                    for (j=var->numDim-1; j>-1; j--, tc+=l) 
                           sprintf(tc,"[%d]%n", var->dimensions[j], &l); 
               }
               if ((var->type == COMPLEX_NTU) ||
                     (var->type == DBL_COMPLEX_NTU)) {
                       sprintf (tc,"[2]%n",&l);
                       tc += l;
               }       
               if (var->description == NULL) descrTmp = nullDescr;
                    else descrTmp = var->description;
               sprintf(tc,"; /* %s */%n", descrTmp, &l); tc += l;
               fprintf(Ffp,"%s\n", line);
            }
        }
        fprintf(Ffp,
        "    %s_v_struct var[%s]; /* The array of substructures */\n",  
                              nameCom, nameMaxIndex); 
        fprintf(Ffp,"    int fence[2]; \n");
        fprintf(Ffp,"} %s_struct; \n", nameCom);
    }    
    free(nameCom);
    fclose(Ffp);
    
}

void mcfioC_SetForSaveDecoding(int val)
{
     if(val != 0) McfNTuPleSaveDecoding = True;
     else McfNTuPleSaveDecoding = False;
}
     
static char *makeStructName(char *title, int orgStyle)
{
    char *out;
    int i, l, nMax;
    
    l = strlen(title);
    if (orgStyle  == PARALLEL_ARRAY_NTU) nMax = 23;  
    else nMax = 21;
    if (l > nMax) l = nMax;
    out = (char *) malloc(sizeof(char) * (l+1));
    strncpy(out, title, l); out[l]='\0';    
    for (i=0; i<l; i++) if (out[i] == ' ') out[i] == '_';
    return out;
} 
/*
** CopyNtrim -  Copy "fromString" to a malloc'd new string,
** 		trimming off leading and trailing spaces & tabs.
**		The newly malloc'd string is returned.
** 		If fromString is NULL, NULL is returned.
*/
static char *mcf_copyNtrim(char *fromString)
{
    char *c, *toString;
    int len, i;
    
    if (fromString == NULL)
    	return NULL;
    toString = (char *) malloc(strlen(fromString)+1);
    
    /* Find the first non-white character */
    for (c=fromString; *c == ' ' || *c == '\t'; c++);

    /* Copy the remainder of fromString to toString */
    strcpy(toString, c);
    
    /* Remove trailing spaces and tabs by converting to nulls */
    len = strlen(toString);
    if (len == 0)			/* special case for empty strings */
    	return toString;
    for (i = len-1; i >= 0; --i) {
    	if (isspace(toString[i]))
    	    toString[i] = '\0';
    	else
    	    break;
    }
    return toString;
}