/* See {salamic_closer.h}. */
/* Last edited on 2015-10-01 00:46:42 by stolfilocal */

#define _GNU_SOURCE
#include <stdio.h>
#include <stdint.h>
#include <string.h>
#include <stdlib.h>
#include <assert.h>

#include <bool.h>
#include <vec.h>
#include <jsfile.h>
#include <affirm.h>

#include <salamic_mesh.h>
#include <salamic_closer.h>
#include <salamic_closer_trivial.h>
#include <salamic_closer_Minetto.h>
/* #include <salamic_closer_Bentley.h> */

void salamic_closer_close
  ( salamic_mesh_t *mesh,
    int32_t pZ,
    int32_t mf,
    int32_t ixFace[], 
    char *closer, 
    char *outPrefix, 
    salamic_stats_t *st
  )
  {
    char *fileName = NULL;
    asprintf(&fileName, "%s%+011d.txt", outPrefix, pZ);
    FILE *wr = open_write(fileName, TRUE);
    
    if (strcmp(closer, "TRIVIAL") == 0)
      { salamic_closer_trivial_close(mesh, pZ, mf, ixFace, wr, st); }
    else if (strcmp(closer, "MINETTO") == 0)
      { salamic_closer_Minetto_close(mesh, pZ, mf, ixFace, wr, st); }
    /* 
    else if (strcmp(closer, "BENTLEY") == 0)
      { salamic_closer_Bentley_close(mesh, pZ, mf, ixFace, wr, st); }
    */
    else 
      { demand(FALSE, "unknown or unimplemented loop closing algorithm"); }
      
    fclose(wr);
    free(fileName);
  }