Changeset 6eaf89ea in sasview for src/sas/sascalc/file_converter/c_ext

Timestamp:

Aug 5, 2016 3:49:33 AM (8 years ago)

Author:

lewis

Branches:

master, ESS_GUI, ESS_GUI_Docs, ESS_GUI_batch_fitting, ESS_GUI_bumps_abstraction, ESS_GUI_iss1116, ESS_GUI_iss879, ESS_GUI_iss959, ESS_GUI_opencl, ESS_GUI_ordering, ESS_GUI_sync_sascalc, costrafo411, magnetic_scatt, release-4.1.1, release-4.1.2, release-4.2.2, ticket-1009, ticket-1094-headless, ticket-1242-2d-resolution, ticket-1243, ticket-1249, ticket885, unittest-saveload

Children:

c3f0114

Parents:

e5308b6

Message:

Implement reading BSL/OTOKO header file

Location:

src/sas/sascalc/file_converter/c_ext

Files:

• bsl_loader.c (modified) (6 diffs)
• bsl_loader.h (modified) (1 diff)

src/sas/sascalc/file_converter/c_ext/bsl_loader.c

@@ -24 +24 @@
     const int n_pixels;
     const int n_rasters;
+    const int swap_bytes;
     if (self != NULL) {
-        if (!PyArg_ParseTuple(args, "siii", &filename, &frame, &n_pixels, &n_rasters))
+        if (!PyArg_ParseTuple(args, "siiii", &filename, &frame, &n_pixels, &n_rasters, &swap_bytes))
             Py_RETURN_NONE;
         if (!(self->params.filename = malloc(strlen(filename) + 1)))
@@ -33 +34 @@
         self->params.n_pixels = n_pixels;
         self->params.n_rasters = n_rasters;
+        self->params.swap_bytes = swap_bytes;
     }
 
@@ -46 +48 @@
     char str[100];
     sprintf(str,
-        "Filename: %s\nframe: %d\nn_pixels: %d\nn_rasters: %d\n",
+        "Filename: %s\nframe: %d\nn_pixels: %d\nn_rasters: %d\nswap_bytes: %d",
         self->params.filename,
         self->params.frame,
         self->params.n_pixels,
-        self->params.n_rasters);
+        self->params.n_rasters,
+        self->params.swap_bytes);
     return Py_BuildValue("s", str);
 }
@@ -104 +107 @@
 
     return Py_BuildValue("i", self->params.n_rasters);
+}
+
+static PyObject *get_swap_bytes(CLoader *self, PyObject *args) {
+    return Py_BuildValue("i", self->params.swap_bytes);
+}
+
+static PyObject *set_swap_bytes(CLoader *self, PyObject *args) {
+    int new_swap;
+    if (!PyArg_ParseTuple(args, "i", &new_swap))
+        return NULL;
+    self->params.swap_bytes = new_swap;
+
+    return Py_BuildValue("i", self->params.swap_bytes);
+}
+
+float reverse_float(const float in_float){
+    float retval;
+    char *to_convert = (char *)&in_float;
+    char *return_float = (char *)&retval;
+
+    return_float[0] = to_convert[3];
+    return_float[1] = to_convert[2];
+    return_float[2] = to_convert[1];
+    return_float[3] = to_convert[0];
+
+    return retval;
 }
 
@@ -138 +167 @@
         for (pixel = 0; pixel < self->params.n_pixels; pixel++) {
             fread(&cur_val, sizeof(float), 1, input_file);
+            if (self->params.swap_bytes == 0)
+                cur_val = reverse_float(cur_val);
             PyArray_SETITEM(data, PyArray_GETPTR2(data, raster, pixel), PyFloat_FromDouble(cur_val));
             read_val = PyArray_GETITEM(data, PyArray_GETPTR2(data, raster, pixel));
@@ -161 +192 @@
     { "get_n_rasters", (PyCFunction)get_n_rasters, METH_VARARGS, "Get n_rasters" },
     { "set_n_rasters", (PyCFunction)set_n_rasters, METH_VARARGS, "Set n_rasters" },
+    { "get_swap_bytes", (PyCFunction)get_swap_bytes, METH_VARARGS, "Get swap_bytes" },
+    { "set_swap_bytes", (PyCFunction)set_swap_bytes, METH_VARARGS, "Set swap_bytes" },
     { "load_data", (PyCFunction)load_data, METH_VARARGS, "Load the data into a numpy array" },
     {NULL}

src/sas/sascalc/file_converter/c_ext/bsl_loader.h

@@ -11 +11 @@
     // Number of rasters in the file
     int n_rasters;
+    // Whether or not the bytes are in reverse order
+    int swap_bytes;
 } CLoader_params;