1 files changed, 228 insertions, 0 deletions
diff --git a/funtools/util/swap.c b/funtools/util/swap.c
new file mode 100644
index 0000000..320b241
--- /dev/null
+++ b/funtools/util/swap.c
@@ -0,0 +1,228 @@
+/*
+ *	Copyright (c) 1999-2003 Smithsonian Astrophysical Observatory
+ */
+
+/*
+ *
+ * swap.c -- swap data between native and IEEE format
+ *
+ *
+ */
+
+#include <stdio.h>
+#include <prsetup.h>
+#include <swap.h>
+
+/*
+ *
+ * swap.c -- routines to swap parts of data
+ *
+ */
+
+#define _idx(x, n)    ((unsigned char *) x)[n]
+
+/* from Harbison&Steele by way of GNU cinfigure ...
+   returns 1 for bigendian, 0 for littleendian */
+#ifdef ANSI_FUNC
+int
+is_bigendian (void)
+#else
+int is_bigendian()
+#endif
+{
+  union
+  {
+    long l;
+    char c[sizeof (long)];
+  } u;
+  u.l = 1;
+  return(u.c[sizeof (long) - 1] == 1);
+}
+
+#ifdef ANSI_FUNC
+void
+swap_short (short *buf, int n)
+#else
+void swap_short(buf, n)
+     short *buf;
+     int n;
+#endif
+{
+  int i;
+  register short *ptr;
+  short val;
+
+  ptr = buf;
+  for(i=0; i<n; i++){
+    _idx(&val, 0) = _idx(ptr, 1);
+    _idx(&val, 1) = _idx(ptr, 0);
+    *ptr++ = val;
+  }
+}
+
+#ifdef ANSI_FUNC
+void
+swap_ushort (unsigned short *buf, int n)
+#else
+void swap_ushort(buf, n)
+     unsigned short *buf;
+     int n;
+#endif
+{
+  int i;
+  register unsigned short *ptr;
+  unsigned short val;
+
+  ptr = buf;
+  for(i=0; i<n; i++){
+    _idx(&val, 0) = _idx(ptr, 1);
+    _idx(&val, 1) = _idx(ptr, 0);
+    *ptr++ = val;
+  }
+}
+
+#ifdef ANSI_FUNC
+void
+swap_int (int *buf, int n)
+#else
+void swap_int(buf, n)
+     int *buf;
+     int n;
+#endif
+{
+  int i;
+  register int *ptr;
+  int val;
+
+  ptr = buf;
+  for(i=0; i<n; i++){
+    _idx(&val, 0) = _idx(ptr, 3);
+    _idx(&val, 1) = _idx(ptr, 2);
+    _idx(&val, 2) = _idx(ptr, 1);
+    _idx(&val, 3) = _idx(ptr, 0);
+    *ptr++ = val;
+  }
+}
+
+#ifdef ANSI_FUNC
+void
+swap_uint (unsigned int *buf, int n)
+#else
+void swap_uint(buf, n)
+     unsigned int *buf;
+     int n;
+#endif
+{
+  int i;
+  register unsigned int *ptr;
+  unsigned int val;
+
+  ptr = buf;
+  for(i=0; i<n; i++){
+    _idx(&val, 0) = _idx(ptr, 3);
+    _idx(&val, 1) = _idx(ptr, 2);
+    _idx(&val, 2) = _idx(ptr, 1);
+    _idx(&val, 3) = _idx(ptr, 0);
+    *ptr++ = val;
+  }
+}
+
+#ifdef ANSI_FUNC
+void
+swap_float (float *buf, int n)
+#else
+void swap_float(buf, n)
+     float *buf;
+     int n;
+#endif
+{
+  int i;
+  register float *ptr;
+  float val;
+
+  ptr = buf;
+  for(i=0; i<n; i++){
+    _idx(&val, 0) = _idx(ptr, 3);
+    _idx(&val, 1) = _idx(ptr, 2);
+    _idx(&val, 2) = _idx(ptr, 1);
+    _idx(&val, 3) = _idx(ptr, 0);
+    *ptr++ = val;
+  }
+}
+
+#ifdef ANSI_FUNC
+void
+swap_double (double *buf, int n)
+#else
+void swap_double(buf, n)
+     double *buf;
+     int n;
+#endif
+{
+  int i;
+  register double *ptr;
+  double val;
+
+  ptr = buf;
+  for(i=0; i<n; i++){
+    _idx(&val, 0) = _idx(ptr, 7);
+    _idx(&val, 1) = _idx(ptr, 6);
+    _idx(&val, 2) = _idx(ptr, 5);
+    _idx(&val, 3) = _idx(ptr, 4);
+    _idx(&val, 4) = _idx(ptr, 3);
+    _idx(&val, 5) = _idx(ptr, 2);
+    _idx(&val, 6) = _idx(ptr, 1);
+    _idx(&val, 7) = _idx(ptr, 0);
+    *ptr++ = val;
+  }
+}
+
+/*
+ *
+
+ *
+ */
+#ifdef ANSI_FUNC
+void
+swap_data(void *buf, int len, int dtype)
+#else
+void swap_data(buf, len, dtype)
+     void *buf;
+     int len;
+     int dtype;
+#endif
+{
+  /* set up input data pointer */
+  switch(dtype){
+  case TY_CHAR:
+  case 'B':
+    break;
+  case TY_SHORT:
+  case 'I':
+    swap_short((short *)buf, len);
+    break;
+  case TY_USHORT:
+  case 'U':
+    swap_ushort((unsigned short *)buf, len);
+    break;
+  case TY_INT:
+  case 'J':
+    swap_int((int *)buf, len);
+    break;
+  case 'V':
+    swap_uint((unsigned int *)buf, len);
+    break;
+  case TY_FLOAT:
+  case 'E':
+    swap_float((float *)buf, len);
+    break;
+  case TY_DOUBLE:
+  case 'F':
+    swap_double((double *)buf, len);
+    break;
+  default:
+    fprintf(stderr, "ERROR: unknown input data type %d\n", dtype);
+    break;
+  }
+}
+