1 files changed, 133 insertions, 0 deletions
diff --git a/Ganlib/src/kdi_c.c b/Ganlib/src/kdi_c.c
new file mode 100644
index 0000000..8c551bb
--- /dev/null
+++ b/Ganlib/src/kdi_c.c
@@ -0,0 +1,133 @@
+
+/**********************************/
+/* C API for kdi file support     */
+/* author: A. Hebert (30/04/2002) */
+/**********************************/
+
+/*
+Copyright (C) 2002 Ecole Polytechnique de Montreal
+
+This library is free software; you can redistribute it and/or
+modify it under the terms of the GNU Lesser General Public
+License as published by the Free Software Foundation; either
+version 2.1 of the License, or (at your option) any later version.
+*/
+
+#if defined(CRAY)
+#define lnword 8
+#else
+#define lnword 4
+#endif
+
+#if !defined(MSDOS)
+#include <unistd.h>
+#endif
+#include <fcntl.h>
+#include <stdlib.h>
+#include <string.h>
+#include "kdi.h"
+
+long long offset;
+
+kdi_file * kdiop_c(char *nomC,int_32 iactio)
+{
+  kdi_file *my_file;
+  my_file = (kdi_file *) malloc(sizeof(*my_file));
+  strcpy(my_file->nom,nomC);
+  my_file->fd = NULL;
+  if (iactio == 0) {
+     FILE *file ;
+     long fd;
+     if ( ( file = fopen(nomC,"rb") ) != NULL ) {
+        fclose(file) ;
+        perror ("open error 0 in kdiop_c ");
+        return NULL;
+     }
+     fd = creat(nomC,0600);
+     close(fd);
+     my_file->fd = fopen(nomC,"r+b");
+  } else if (iactio == 1) {
+     FILE *file ;
+     if ( ( file = fopen(nomC,"rb") ) == NULL ) {
+        perror ("open error 1 in kdiop_c ");
+        return NULL;
+     }
+     fclose(file) ;
+     my_file->fd = fopen(nomC,"r+b");
+  } else if (iactio == 2) {
+     FILE *file ;
+     if ( ( file = fopen(nomC,"rb") ) == NULL ) {
+        perror ("open error 2 in kdiop_c ");
+        return NULL;
+     }
+     my_file->fd = file;
+  } else {
+     return NULL;
+  }
+  if ( my_file->fd == NULL ) {
+     perror ("open error 3 in kdiop_c ");
+     return NULL;
+  }
+  return my_file;
+}
+
+int_32 kdiput_c(kdi_file *my_file,int_32 *data,int_32 iofset,int_32 length)
+{
+  int_32 irc=0;
+  offset=(long long)iofset*lnword;
+  if (my_file == NULL) {
+     irc = -1;
+  } else if (fseek(my_file->fd,offset,0) >= 0) {
+     long long n, iof=0;
+     while ((n = fwrite(&data[iof],lnword,length,my_file->fd)) < length-iof) {
+        if (n < 0) return n-1;
+        iof+=n;
+     }
+  } else {
+     irc = -3;
+  }
+  return irc;
+}
+
+int_32 kdiget_c(kdi_file *my_file,int_32 *data,int_32 iofset,int_32 length)
+{
+  int_32 irc=0;
+  offset=(long long)iofset*lnword;
+  if (my_file == NULL) {
+     irc = -1;
+  } else if (fseek(my_file->fd,offset,0) >= 0) {
+     long long n, iof=0;
+     while ((n = fread(&data[iof],lnword,length,my_file->fd)) < length-iof) {
+        if (n == 0) return -4;
+        if (n < 0) return n-1;
+        iof+=n;
+     }
+  } else {
+     irc = -3;
+  }
+  return irc;
+}
+
+int_32 kdicl_c(kdi_file *my_file,int_32 istatu)
+{
+  long irc;
+  if (my_file == NULL) {
+     irc = -1;
+  } else if (istatu == 1) {
+     irc = fclose(my_file->fd);
+     free(my_file);
+  } else if (istatu == 2) {
+     irc = fclose(my_file->fd);
+     if (irc != 0 ) {
+        perror ("close error 1 in kdicl_c ");
+        return irc;
+     }
+     irc = remove(my_file->nom);
+     free(my_file);
+  } else {
+     irc = -999;
+  }
+  my_file = NULL;
+  if (irc != 0) perror ("close error 2 in kdicl_c ");
+  return irc;
+}