Read/store ngspice simulation output

1 files changed, 132 insertions, 0 deletions

diff --git a/src/read_crossbar.c b/src/read_crossbar.c
new file mode 100644
index 0000000..3666d53
--- /dev/null
+++ b/src/read_crossbar.c
@@ -0,0 +1,132 @@
+#include "read_crossbar.h"
+
+#include <ctype.h>
+#include <errno.h> //maybe I could use this
+#include <fcntl.h>
+#include <stdio.h>
+#include <stdlib.h>
+#include <string.h>
+#include <sys/types.h>
+#include <sys/wait.h>
+#include <unistd.h>
+
+int run_ngspice(const char *crossbar_path) {
+    if (crossbar_path == NULL) {
+        fprintf(stderr, "One of the file paths returned null when running ngspice");
+        return -1;
+    }
+
+    char command[4096]; 
+    int written = snprintf(command, sizeof(command), 
+            "ngspice -b \"%s\" > /dev/null 2>&1", crossbar_path);
+   
+    if (written < 0) {
+        fprintf(stderr, "Failed to write spice command");
+        return -1;
+    }
+
+    int status = system(command);
+    if (status < 0) {
+        fprintf(stderr, "Failed to run ngspice command");
+    }
+
+    return 0;
+}
+
+int read_crossbar(const char *data_path, size_t num_outputs,
+        Crossbar_Output_Matrix *result) {
+    
+    FILE *file = fopen(data_path, "r"); 
+    if (!file){
+        fprintf(stderr, "Faile to open data file");
+        return -1;
+    }
+
+    result->num_samples = 0;
+    result->num_outputs = num_outputs;
+    result->time = NULL;
+    result->voltages = NULL;
+
+    size_t capacity = 100000;
+
+    result->time = malloc(capacity * sizeof(*result->time));
+    result->voltages = malloc(capacity * num_outputs * sizeof(*result->voltages));
+    
+    if (result->time == NULL || result->voltages == NULL) {
+        fprintf(stderr, "Failed to allocate memory for results");
+        fclose(file);
+        free_crossbar_output_matrix(result);
+        return -1;
+    }
+
+    size_t expected_fields = num_outputs * 2;
+    double *fields = malloc(expected_fields * sizeof(*fields));
+    if (fields == NULL) {
+        fprintf(stderr, "failed to allocate memory for fields");
+        fclose(file);
+        free_crossbar_output_matrix(result);
+        return -1;
+    }
+
+    char line[16384];
+    // size_t line_number = 0;
+
+    while(fgets(line, sizeof(line), file) != NULL) {
+        /* Raises capacity if needed, not working rn tho */
+        // if (result->num_samples == capacity) {
+        //     capacity *= 2;
+        //
+        //     result->time = realloc(result->time, capacity * sizeof(double));
+        //     result->voltages = realloc(result->voltages, capacity * sizeof(double));
+        //
+        //     if (result->time == NULL || result->voltages == NULL) {
+        //         fprintf(stderr, "Failed to reallocate memory");
+        //         fclose(file);
+        //         return -1;
+        //     }
+        // }
+
+        char *position = line;
+        double sample_time = 0.0;
+
+        //each loop reads one time-voltage pair, saves it
+        //then moves to next loop (next time-voltage pair)
+        for (size_t output = 0; output < num_outputs; output++) {
+            double time;
+            double voltage;
+            int char_count;
+
+            if(sscanf(position, "%lf %lf %n", &time, &voltage, &char_count) != 2) {
+                fprintf(stderr, "invalide data line");
+                fclose(file);
+                return -1;
+            }
+
+            if (output == 0) {
+                sample_time = time;
+            }
+
+            result->voltages[result->num_samples * num_outputs + output] = voltage;
+            position += char_count;
+
+        }
+        result->time[result->num_samples] = sample_time;
+        result->num_samples++;
+    }
+    fclose(file);
+    return 0;
+}
+
+void free_crossbar_output_matrix(Crossbar_Output_Matrix *result) {
+    if (result == NULL) {
+        return;
+    }
+
+    free(result->time);
+    free(result->voltages);
+
+    result->num_samples = 0;
+    result->num_outputs = 0;
+    result->time = NULL;
+    result->voltages = NULL;
+}