Fix norm2 so it adheres to specification: norm2(x) = sqrt(dot(x , x))

1 files changed, 1 insertions, 1 deletions

diff --git a/libmath/blas.c b/libmath/blas.c
index e40d2310..4e005674 100644
--- a/libmath/blas.c
+++ b/libmath/blas.c
@@ -58,5 +58,5 @@ norm2(int n, double *x)
 		sum += *x * *x;
 		x++;
 	}
-	return sum;
+	return sqrt(sum);
 }