src/utils.c
@@ -11,6 +11,7 @@ { int i,j; for(i =0 ; i < M; ++i){ printf("%d ", i+1); for(j = 0; j < N; ++j){ printf("%10.6f, ", A[i*N+j]); } @@ -261,6 +262,16 @@ } } float mag_array(float *a, int n) { int i; float sum = 0; for(i = 0; i < n; ++i){ sum += a[i]*a[i]; } return sqrt(sum); } void scale_array(float *a, int n, float s) { int i;