From 12ba0f935031ac21235af125682e34ed7eb9860c Mon Sep 17 00:00:00 2001
From: fred <fred>
Date: Fri, 4 Oct 1996 19:59:40 +0000
Subject: [PATCH] lilypond-0.0.1

---
 vector.hh | 106 ++++++++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 106 insertions(+)
 create mode 100644 vector.hh

diff --git a/vector.hh b/vector.hh
new file mode 100644
index 0000000000..5be792c54b
--- /dev/null
+++ b/vector.hh
@@ -0,0 +1,106 @@
+#ifndef VECTOR_HH
+#define VECTOR_HH
+
+#include "glob.hh"
+#include "vray.hh"
+
+/// a row of numbers
+class Vector  {
+    svec<Real> dat;
+public:
+    void OK() const { dat.OK();}
+    int dim() const { return dat.sz(); }
+    Vector() { }
+    Vector(const Vector&n);
+    Vector(int n) {
+	dat.set_size(n);
+	fill(0);
+    }
+    void insert(Real v, int i) {
+	dat.insert(v,i);
+    }
+    void del(int i) { dat.del(i); }
+    operator String() const;
+    void fill(Real r) {
+	for (int i=0; i < dim(); i++)
+	    dat[i] =r;
+    }
+
+    void operator +=(Vector v) {
+	assert(v.dim() == dim());
+	for (int i=0; i < dim(); i++)
+	    dat[i] += v.dat[i];
+    }
+    
+    void operator /=(Real a) {
+	(*this) *= 1/a;
+    }
+
+    void operator *=(Real a) {
+	for (int i=0; i < dim(); i++)
+	    dat[i] *= a;
+    }
+
+    void operator -=(Vector v) {
+	assert(v.dim() == dim());
+	for (int i=0; i < dim(); i++)
+	    dat[i] -= v(i);	
+    }
+
+    Real &operator()(int i) { return dat[i]; }
+    Real operator()(int i) const { return dat[i]; }
+    Real elem(int i) { return dat[i]; }
+    Real operator *(Vector v) const {
+	Real ip=0;
+	assert(v.dim() == dim());
+	for (int i=0; i < dim(); i++)
+	    ip += dat[i] *v(i);
+	return ip;
+    }
+    Vector operator-() const;
+    Real norm() {
+	return sqrt(norm_sq() );
+    }
+    Real norm_sq() {
+	return ((*this) * (*this));
+    }
+    operator svec<Real> () { return dat; }
+    void print() const;
+    /// set to j-th element of unit-base
+    void set_unit(int j) ;
+};
+/**
+    a vector. Storage is handled in svec, Vector only does the mathematics.
+ */
+
+inline Vector
+operator+(Vector a, Vector const &b) {
+    a += b;
+    return a;
+}
+
+inline Vector
+operator-(Vector a, Vector const &b) {
+    a -= b;
+    return a;
+}
+
+inline Vector
+operator*(Vector v, Real a) {
+    v *= a;
+    return v;
+}
+
+inline Vector
+operator*( Real a,Vector v) {
+    v *= a;
+    return v;
+}
+
+inline Vector
+operator/(Vector v,Real a) {
+    v *= 1/a;
+    return v;
+}
+
+#endif
-- 
2.39.5