From 37da2899f40661e3e9631e497da8dc59b971cbd0 Mon Sep 17 00:00:00 2001
From: "Charles.Forsyth" <devnull@localhost>
Date: Fri, 22 Dec 2006 17:07:39 +0000
Subject: 20060303a

---
 libmath/blas.c | 61 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 61 insertions(+)
 create mode 100644 libmath/blas.c

(limited to 'libmath/blas.c')

diff --git a/libmath/blas.c b/libmath/blas.c
new file mode 100644
index 00000000..a5cc9e0b
--- /dev/null
+++ b/libmath/blas.c
@@ -0,0 +1,61 @@
+#include "lib9.h"
+#include "mathi.h"
+
+double
+dot(int n, double *x, double *y)
+{
+	double	sum = 0;
+	if (n <= 0) 
+		return 0;
+	while (n--) {
+		sum += *x++ * *y++;
+	}
+	return sum;
+}
+
+
+int
+iamax(int n, double *x)
+{
+	int	i, m;
+	double	xm, a;
+	if (n <= 0) 
+		return 0;
+	m = 0;
+	xm = fabs(*x);
+	for (i = 1; i < n; i++) {
+		a = fabs(*++x);
+		if (xm < a) {
+			m = i;
+			xm = a;
+		}
+	}
+	return m;
+}
+
+
+double
+norm1(int n, double *x)
+{
+	double	sum = 0;
+	if (n <= 0) 
+		return 0;
+	while (n--) {
+		sum += fabs(*x++);
+	}
+	return sum;
+}
+
+
+double
+norm2(int n, double *x)
+{
+	double	sum = 0;
+	if (n <= 0) 
+		return 0;
+	while (n--) {
+		sum += *x * *x;
+		x++;
+	}
+	return sum;
+}
-- 
cgit v1.2.3