aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorCalvin Morrison <mutantturkey@gmail.com>2014-03-09 17:03:30 -0400
committerCalvin Morrison <mutantturkey@gmail.com>2014-03-09 17:03:30 -0400
commitc4116bde303e76a56b59300fecbc21c1e4826efd (patch)
tree896e05d33654ea4bdf2ec705268023c9f67f4a91
parent4944915889527e96217a292324d00c7e8666a2db (diff)
formatting
-rw-r--r--src/c/quikr_functions.c83
1 files changed, 41 insertions, 42 deletions
diff --git a/src/c/quikr_functions.c b/src/c/quikr_functions.c
index 89f0698..595d2dc 100644
--- a/src/c/quikr_functions.c
+++ b/src/c/quikr_functions.c
@@ -13,10 +13,10 @@
void check_malloc(void *ptr, char *error) {
if (ptr == NULL) {
if(error != NULL) {
- fprintf(stderr,"Error: %s\n", error);
- }
+ fprintf(stderr,"Error: %s\n", error);
+ }
else { fprintf(stderr, "Error: Could not allocate enough memory - %s\n", strerror(errno));
- }
+ }
exit(EXIT_FAILURE);
}
}
@@ -46,42 +46,42 @@ void debug_arrays(double *count_matrix, struct matrix *sensing_matrix) {
void normalize_matrix(double *matrix, unsigned long long height, unsigned long long width) {
- unsigned long long x = 0;
- unsigned long long y = 0;
+ unsigned long long x = 0;
+ unsigned long long y = 0;
- for(x = 0; x < height; x++) {
+ for(x = 0; x < height; x++) {
- double row_sum = 0;
+ double row_sum = 0;
- for(y = 0; y < (width); y++)
- row_sum = row_sum + matrix[width * x + y];
- for(y = 0; y < (width); y++)
- matrix[width * x + y] = matrix[width * x + y] / row_sum;
- }
+ for(y = 0; y < (width); y++)
+ row_sum = row_sum + matrix[width * x + y];
+ for(y = 0; y < (width); y++)
+ matrix[width * x + y] = matrix[width * x + y] / row_sum;
+ }
}
-unsigned long long count_sequences(const char *filename) {
- char *line = NULL;
- size_t len = 0;
- ssize_t read;
+unsigned long long count_sequences(const char *filename) {
+ char *line = NULL;
+ size_t len = 0;
+ ssize_t read;
- unsigned long long sequences = 0;
+ unsigned long long sequences = 0;
- FILE *fh = fopen(filename, "r");
- if(fh == NULL) {
- fprintf(stderr, "could not open \"%s\"\n", filename );
+ FILE *fh = fopen(filename, "r");
+ if(fh == NULL) {
+ fprintf(stderr, "could not open \"%s\"\n", filename );
return 0;
- }
+ }
- while ((read = getline(&line, &len, fh)) != -1) {
- if(line[0] == '>')
- sequences++;
- }
+ while ((read = getline(&line, &len, fh)) != -1) {
+ if(line[0] == '>')
+ sequences++;
+ }
- free(line);
- fclose(fh);
+ free(line);
+ fclose(fh);
- return sequences;
+ return sequences;
}
@@ -94,20 +94,20 @@ struct matrix *load_sensing_matrix(const char *filename, unsigned int target_kme
unsigned int kmer = 0;
- unsigned long long i = 0;
+ unsigned long long i = 0;
unsigned long long *row = NULL;
unsigned long long sequences = 0;
unsigned long long width = 0;
- struct matrix *ret = NULL;
+ struct matrix *ret = NULL;
- gzFile fh = NULL;
+ gzFile fh = NULL;
- fh = gzopen(filename, "r");
- if(fh == NULL) {
- fprintf(stderr, "could not open %s", filename);
- exit(EXIT_FAILURE);
- }
+ fh = gzopen(filename, "r");
+ if(fh == NULL) {
+ fprintf(stderr, "could not open %s", filename);
+ exit(EXIT_FAILURE);
+ }
line = malloc(1024 * sizeof(char));
check_malloc(line, NULL);
@@ -150,18 +150,18 @@ struct matrix *load_sensing_matrix(const char *filename, unsigned int target_kme
width = pow_four(kmer);
// allocate a +1 size for the extra row
- matrix = malloc(sequences * (width) * sizeof(double));
+ matrix = malloc(sequences * (width) * sizeof(double));
check_malloc(matrix, NULL);
row = malloc((width) * sizeof(unsigned long long));
check_malloc(row, NULL);
- headers = malloc(sequences * sizeof(char *));
+ headers = malloc(sequences * sizeof(char *));
check_malloc(headers, NULL);
for(i = 0; i < sequences; i++) {
unsigned long long j = 0;
- // get header and add it to headers array
+ // get header and add it to headers array
char *header = malloc(256 * sizeof(char));
check_malloc(header, NULL);
gzgets(fh, header, 256);
@@ -189,17 +189,16 @@ struct matrix *load_sensing_matrix(const char *filename, unsigned int target_kme
}
}
- for(j = 0; j < width; j++) {
+ for(j = 0; j < width; j++) {
matrix[i*(width) + j] = ((double)row[j]);
}
}
// load the matrix of counts
- gzclose(fh);
+ gzclose(fh);
free(line);
free(row);
-
ret = malloc(sizeof(struct matrix));
(*ret).kmer = kmer;
(*ret).sequences = sequences;