Beispiel #1
0
        // numerics
        public MatrixNd IdentityMatrix(int r, int c)
        {
            MatrixNd mat = new MatrixNd(r, c);

            for (int i = 0; i < r; ++i)
            {
                mat[i, i] = 1.0;
            }
            return(mat);
        }
Beispiel #2
0
        static public MatrixNd operator *(double factor, MatrixNd m)
        {
            MatrixNd mat = new MatrixNd(m);

            for (int i = 0; i < m.length; ++i)
            {
                mat[i] *= factor;
            }
            return(mat);
        }
Beispiel #3
0
 public MatrixNd(MatrixNd mat)
 {
     M      = mat.Row;
     N      = mat.Col;
     length = M * N;
     arr    = new double[length];
     for (int i = 0; i < length; ++i)
     {
         arr[i] = mat.arr[i];
     }
 }
Beispiel #4
0
        public MatrixNd Transpose()
        {
            MatrixNd mat = new MatrixNd(N, M);

            for (int i = 0; i < M; ++i)
            {
                for (int j = 0; j < N; ++j)
                {
                    mat[j, i] = this[i, j];
                }
            }
            return(mat);
        }
Beispiel #5
0
        static public MatrixNd operator /(double factor, MatrixNd m)
        {
            MatrixNd mat = new MatrixNd(m);

            if (factor == 0)
            {
                throw new DivideByZeroException();
            }
            for (int i = 0; i < m.length; ++i)
            {
                mat[i] /= factor;
            }
            return(mat);
        }
Beispiel #6
0
        static public MatrixNd operator -(MatrixNd m1, MatrixNd m2)
        {
            if (m1.M != m2.M || m1.N != m2.N)
            {
                return(null);
            }
            MatrixNd m = new MatrixNd(m1);

            for (int i = 0; i < m1.length; ++i)
            {
                m[i] -= m2[i];
            }
            return(m);
        }
Beispiel #7
0
        static public MatrixNd operator *(MatrixNd m1, MatrixNd m2)
        {
            if (m1.Col != m2.Row)
            {
                return(null);
            }
            MatrixNd mat = new MatrixNd(m1.Row, m2.Col);

            for (int i = 0; i < m1.Row; ++i)
            {
                for (int j = 0; j < m2.Col; ++j)
                {
                    for (int k = 0; k < m1.Col; ++k)
                    {
                        mat[i, j] += m1[i, k] * m2[k, j];
                    }
                }
            }
            return(mat);
        }