Vector.cs

Go to the documentation of this file.

using System;
using System.Collections;
using BenTools.Data;

namespace BenTools.Mathematics
{
    public class Vector : IEnumerable, IComparable
    {
        public static int Precision = 10;
        double[] data;
        public object Tag=null;
        public Vector(int dim)
        {
            data = new double[dim];
        }
        public Vector(params double[] X)
        {
            data = new double[X.Length];
            X.CopyTo(data,0);
        }
        public Vector(Vector O)
            : this(O.data)
        {}
        public Vector(string S)
        {
            if(S[0]!='(' || S[S.Length-1]!=')')
                throw new Exception("Formatfehler!");
            string[] P = MathTools.HighLevelSplit(S.Substring(1,S.Length-2),';');
            data = new double[P.Length];
            int i;
            for(i=0;i<data.Length;i++)
            {
                data[i] = double.Parse(P[i]);
            }
        }
        public double this [int i]
        {
            get
            {
                return data[i];
            }
            set
            {
                data[i] = Math.Round(value,Precision);
            }
        }

        public int Dim
        {
            get
            {
                return data.Length;
            }
        }

        public double SquaredLength
        {
            get
            {
                return this*this;
            }
        }

        public double ElementSum
        {
            get
            {
                int i;
                double E = 0;
                for(i=0;i<Dim;i++)
                    E+=data[i];
                return E;
            }
        }
        public void Randomize(double Min, double Max)
        {
            int i;
            for(i=0;i<data.Length;i++)
            {
                this[i] = Min + (Max-Min)*MathTools.R.NextDouble();
            }
        }
        public void Randomize(Vector[] MinMax)
        {
            int i;
            for(i=0;i<data.Length;i++)
            {
                this[i] = MinMax[0][i] + (MinMax[1][i]-MinMax[0][i])*MathTools.R.NextDouble();
            }
        }
        public void Multiply(double r)
        {
            int i;
            for(i=0;i<data.Length;i++)
            {
                this[i] *= r;
            }
        }
        public void Add(Vector V)
        {
            int i;
            for(i=0;i<data.Length;i++)
            {
                this[i] += V[i];
            }
        }
        public void Add(double d)
        {
            int i;
            for(i=0;i<data.Length;i++)
            {
                this[i] += d;
            }
        }
        IEnumerator IEnumerable.GetEnumerator()
        {
            return data.GetEnumerator();
        }

        public override string ToString()
        {
            string S = "(";
            int i;
            for(i=0;i<data.Length;i++)
            {
                S += data[i].ToString("G4");
                if(i<data.Length-1)
                    S+=";";
            }
            S+=")";
            return S;
        }

        public override bool Equals(object obj)
        {
            Vector B = obj as Vector;
            if(B==null || data.Length != B.data.Length)
                return false;
            int i;
            for(i=0;i<data.Length;i++)
            {
                if(Math.Abs(data[i]-B.data[i])>1e-10)
                    return false;
            }
            return true;
        }

        public override int GetHashCode()
        {
            int Erg = 0;
            foreach(double D in data)
                Erg = Erg ^ Math.Round(D,Precision).GetHashCode();
            return Erg;
        }

        public static Vector operator - (Vector A, Vector B)
        {
            if(A.Dim!=B.Dim)
                throw new Exception("Vectors of different dimension!");
            Vector Erg = new Vector(A.Dim);
            int i;
            for(i=0;i<A.Dim;i++)
                Erg[i]=A[i]-B[i];
            return Erg;
        }

        public static Vector operator + (Vector A, Vector B)
        {
            if(A.Dim!=B.Dim)
                throw new Exception("Vectors of different dimension!");
            Vector Erg = new Vector(A.Dim);
            int i;
            for(i=0;i<A.Dim;i++)
                Erg[i]=A[i]+B[i];
            return Erg;
        }

        public static double operator * (Vector A, Vector B)
        {
            if(A.Dim!=B.Dim)
                throw new Exception("Vectors of different dimension!");
            double Erg = 0;
            int i;
            for(i=0;i<A.Dim;i++)
                Erg+=A[i]*B[i];
            return Erg;
        }

        public static Vector operator * (Vector A, double B)
        {
            Vector Erg = new Vector(A.Dim);
            int i;
            for(i=0;i<A.Dim;i++)
                Erg[i]=A[i]*B;
            return Erg;
        }

        public static Vector operator * (double A, Vector B)
        {
            return B*A;
        }
        public static explicit operator double[](Vector A)
        {
            return A.data;
        }
        public static double Dist(Vector V1, Vector V2)
        {
            if(V1.Dim != V2.Dim)
                return -1;
            int i;
            double E = 0,D;
            for(i=0;i<V1.Dim;i++)
            {
                D=(V1[i]-V2[i]);
                E+=D*D;
            }
            return E;

        }

        public int CompareTo(object obj)
        {
            Vector A = this;
            Vector B = obj as Vector;
            if(A==null || B==null)
                return 0;
            double Al,Bl;
            Al = A.SquaredLength;
            Bl = B.SquaredLength;
            if(Al>Bl)
                return 1;
            if(Al<Bl)
                return -1;
            int i;
            for(i=0;i<A.Dim;i++)
            {
                if(A[i]>B[i])
                    return 1;
                if(A[i]<B[i])
                    return -1;
            }
            return 0;
        }
        public virtual Vector Clone()
        {
            return new Vector(data);
        }
    }

}