1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
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;
}
|
