From 4cf4a89f05c671041c217a724fc9b4b026eb9de9 Mon Sep 17 00:00:00 2001
From: AlexeyAB <alexeyab84@gmail.com>
Date: Fri, 01 Jun 2018 21:33:17 +0000
Subject: [PATCH] Minor fix

---
 src/matrix.c |   98 +++++++++++++++++++++++++++++++++++++------------
 1 files changed, 74 insertions(+), 24 deletions(-)

diff --git a/src/matrix.c b/src/matrix.c
index 96bd332..08e6109 100644
--- a/src/matrix.c
+++ b/src/matrix.c
@@ -13,16 +13,53 @@
     free(m.vals);
 }
 
-float matrix_accuracy(matrix truth, matrix guess)
+float matrix_topk_accuracy(matrix truth, matrix guess, int k)
 {
-    int k = truth.cols;
-    int i;
-    int count = 0;
+    int *indexes = calloc(k, sizeof(int));
+    int n = truth.cols;
+    int i,j;
+    int correct = 0;
     for(i = 0; i < truth.rows; ++i){
-        int class = max_index(guess.vals[i], k);
-        if(truth.vals[i][class]) ++count;
+        top_k(guess.vals[i], n, k, indexes);
+        for(j = 0; j < k; ++j){
+            int class_id = indexes[j];
+            if(truth.vals[i][class_id]){
+                ++correct;
+                break;
+            }
+        }
     }
-    return (float)count/truth.rows;
+    free(indexes);
+    return (float)correct/truth.rows;
+}
+
+void scale_matrix(matrix m, float scale)
+{
+    int i,j;
+    for(i = 0; i < m.rows; ++i){
+        for(j = 0; j < m.cols; ++j){
+            m.vals[i][j] *= scale;
+        }
+    }
+}
+
+matrix resize_matrix(matrix m, int size)
+{
+    int i;
+    if (m.rows == size) return m;
+    if (m.rows < size) {
+        m.vals = realloc(m.vals, size*sizeof(float*));
+        for (i = m.rows; i < size; ++i) {
+            m.vals[i] = calloc(m.cols, sizeof(float));
+        }
+    } else if (m.rows > size) {
+        for (i = size; i < m.rows; ++i) {
+            free(m.vals[i]);
+        }
+        m.vals = realloc(m.vals, size*sizeof(float*));
+    }
+    m.rows = size;
+    return m;
 }
 
 void matrix_add_matrix(matrix from, matrix to)
@@ -80,30 +117,43 @@
 
 matrix csv_to_matrix(char *filename)
 {
-	FILE *fp = fopen(filename, "r");
-	if(!fp) file_error(filename);
+    FILE *fp = fopen(filename, "r");
+    if(!fp) file_error(filename);
 
     matrix m;
     m.cols = -1;
 
-	char *line;
+    char *line;
 
-	int n = 0;
-	int size = 1024;
-	m.vals = calloc(size, sizeof(float*));
-	while((line = fgetl(fp))){
+    int n = 0;
+    int size = 1024;
+    m.vals = calloc(size, sizeof(float*));
+    while((line = fgetl(fp))){
         if(m.cols == -1) m.cols = count_fields(line);
-		if(n == size){
-			size *= 2;
-			m.vals = realloc(m.vals, size*sizeof(float*));
-		}
-		m.vals[n] = parse_fields(line, m.cols);
-		free(line);
-		++n;
-	}
-	m.vals = realloc(m.vals, n*sizeof(float*));
+        if(n == size){
+            size *= 2;
+            m.vals = realloc(m.vals, size*sizeof(float*));
+        }
+        m.vals[n] = parse_fields(line, m.cols);
+        free(line);
+        ++n;
+    }
+    m.vals = realloc(m.vals, n*sizeof(float*));
     m.rows = n;
-	return m;
+    return m;
+}
+
+void matrix_to_csv(matrix m)
+{
+    int i, j;
+
+    for(i = 0; i < m.rows; ++i){
+        for(j = 0; j < m.cols; ++j){
+            if(j > 0) printf(",");
+            printf("%.17g", m.vals[i][j]);
+        }
+        printf("\n");
+    }
 }
 
 void print_matrix(matrix m)

--
Gitblit v1.10.0