lattice.c

#include "typedef.h"
#include "datei.h"

extern int INFO_LEVEL;

lattice_element **lattice(char *symb,int dim,int almost,int zclass,int *no,
                          int OPTION)
{
   lattice_element **res;

   FILE *infile;

   int i;

   /* saves space on the stack */
   static char filename[1000];

   /* get the appropiate filename */
   sprintf(filename,"%s%s%d/%s%s_%d_%d",TOPDIR,"/tables/lattices/dim",
             dim,"lattice_",symb,almost,zclass);

   infile = fopen(filename,"r");

   /* save cores */
   if (infile == NULL){
     fprintf(stderr,"Couldn't find my input file %s\n",filename);
     exit(4);
   }

   fscanf(infile,"#%d\n",no);

   if (INFO_LEVEL & 1){
      fprintf(stderr,"%d\n",no[0]);
   }

   /* now read the data for lattice */
   res = (lattice_element **) malloc(no[0] * sizeof(lattice_element *));

   /* read the subgroups */
   for (i=0;i<no[0];i++){

      /* get the element of lattice */
      res[i] = fget_lattice_element(infile,OPTION);

      /* output for debugging */
      if (INFO_LEVEL & 1){
         fput_lattice_element(res[i],NULL);
      }
   }

   fclose(infile);

   return res;
}