//display specified matrix static void displayMatrix(GeneralMatrix displayMat) { for (int i = 0; i < displayMat.RowDimension; i++) { for (int j = 0; j < displayMat.ColumnDimension; j++) { Console.Write(displayMat.GetElement(i, j) + ","); } Console.WriteLine(); } }
public void Hessian() { double eps = 1e-5; double[] vector = {0,0,0}; double[][] result = {new double[]{10.0, 2.0, -2.0}, new double[]{ 2.0, 4.0, 2.0}, new double[]{-2.0, 2.0, 4.0}}; GeneralMatrix expectedMatrix = new GeneralMatrix(result); TestFunction f = new TestFunction(); GeneralMatrix hessian = f.CalculateHessian(vector); for (int i=0;i<3; i++) for (int j=0;j<3; j++) Assert.IsTrue(System.Math.Abs(hessian.GetElement(i,j)-expectedMatrix.GetElement(i,j))<eps); }
/// <summary>Set a submatrix.</summary> /// <param name="i0"> Initial row index /// </param> /// <param name="i1"> Final row index /// </param> /// <param name="c"> Array of column indices. /// </param> /// <param name="X"> A(i0:i1,c(:)) /// </param> /// <exception cref="System.IndexOutOfRangeException"> Submatrix indices /// </exception> public virtual void SetMatrix(int i0, int i1, int[] c, GeneralMatrix X) { try { for (int i = i0; i <= i1; i++) { for (int j = 0; j < c.Length; j++) { A[i][c[j]] = X.GetElement(i - i0, j); } } } catch (System.IndexOutOfRangeException e) { throw new System.IndexOutOfRangeException("Submatrix indices", e); } }
/// <summary>Set a submatrix.</summary> /// <param name="r"> Array of row indices. /// </param> /// <param name="j0"> Initial column index /// </param> /// <param name="j1"> Final column index /// </param> /// <param name="X"> A(r(:),j0:j1) /// </param> /// <exception cref="System.IndexOutOfRangeException"> Submatrix indices /// </exception> public virtual void SetMatrix(int[] r, int j0, int j1, GeneralMatrix X) { try { for (int i = 0; i < r.Length; i++) { for (int j = j0; j <= j1; j++) { A[r[i]][j] = X.GetElement(i, j - j0); } } } catch (System.IndexOutOfRangeException e) { throw new System.IndexOutOfRangeException("Submatrix indices", e); } }
/// <summary> /// Transforms a basis blade to another basis which is represented by the columns of M. /// </summary> /// <param name="a">The basis blade to transform according to <paramref name="M"/>.</param> /// <param name="M">The matrix to use to transform <paramref name="a"/>.</param> /// <returns>a list of BasisBlade whose sum represents <paramref name="a"/> on the new basis.</returns> public static ArrayList Transform(BasisBlade a, DotNetMatrix.GeneralMatrix M) { ArrayList A = new ArrayList(); A.Add(new BasisBlade(0, a.scale, a.symScale)); // start with just the scalar part; int dim = M.RowDimension; // for each 1 bit: convert to list of blades int i = 0; uint b = a.bitmap; while (b != 0) { if ((b & 1) != 0) { // take column 'i' out of the matrix, wedge it to 'A' ArrayList tmp = new ArrayList(); for (int j = 0; j < dim; j++) { double m = M.GetElement(j, i); if (m != 0.0) { for (int k = 0; k < A.Count; k++) { BasisBlade o = BasisBlade.op((BasisBlade)A[k], new BasisBlade((uint)(1 << j), m)); if (o.scale != 0.0) { tmp.Add(o); } } } } A = tmp; } b >>= 1; i++; } return(A); }
public void ZeroElementAccessThrowsAnIndexOutOfRangeExceptionWhenRowValueIsEqualToRowDimension() { var matrix = new GeneralMatrix(avals); matrix.GetElement(matrix.RowDimension, matrix.ColumnDimension-1); }
private double[] CalculateNextPoint(double[] pX, double[] pGrad, GeneralMatrix hessian) { int i=0; double xmin=0; double step = _step; GeneralMatrix alfaX = new GeneralMatrix(_nDim,1); GeneralMatrix prevX = new GeneralMatrix(pX,_nDim); GeneralMatrix prevGrad = new GeneralMatrix(pGrad,_nDim); double[] intermediate = new double[_nDim];; alfaX = hessian.Multiply(prevGrad); //doing a line search to minimize alpha OneDWrapper wrapper = new OneDWrapper(_f,prevX,alfaX); LineSearch search = new LineSearch(); double[] interval = new double[Constants.BRACKET_POINTS]; int it1 = search.FindMinInterval(wrapper,_alpha,step,50,ref interval); int it2 = search.FindMinimumViaBrent(wrapper,interval[0],interval[1],interval[2],50,_epsilon, ref xmin); for (i=0;i<_nDim; i++) intermediate[i] = prevX.GetElement(i,0) - xmin*alfaX.GetElement(i,0); _alpha = xmin; return intermediate; }
public void Substraction() { double[] columnwise = new double[]{1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0, 11.0, 12.0}; double[][] avals = {new double[]{1.0, 4.0, 7.0, 10.0}, new double[]{2.0, 5.0, 8.0, 11.0}, new double[]{3.0, 6.0, 9.0, 12.0}}; int validld = 3; /* leading dimension of intended test Matrices */ //one dimensional array of doubles packed by columns ala Fortran //is passed to the constructor //the integer value tell the constructor how to //breakup one domensional value into mulidimensional column GeneralMatrix A = new GeneralMatrix(columnwise, validld); GeneralMatrix B = new GeneralMatrix(avals); double tmp = B.GetElement(0, 0); avals[0][0] = 0.0; GeneralMatrix C = B.Subtract(A); avals[0][0] = tmp; B = GeneralMatrix.Create(avals); tmp = B.GetElement(0, 0); avals[0][0] = 0.0; Assert.IsTrue(tmp==B.GetElement(0, 0)); }
public void GetElement() { double[][] avals = {new double[]{1.0, 4.0, 7.0, 10.0}, new double[]{2.0, 5.0, 8.0, 11.0}, new double[]{3.0, 6.0, 9.0, 12.0}}; GeneralMatrix B = new GeneralMatrix(avals); Assert.AreEqual(B.GetElement(B.RowDimension-1,B.ColumnDimension-1),avals[B.RowDimension-1][B.ColumnDimension-1]); }
/// <summary> /// Average values of the priority matrix over sum of columns. /// set values of sum of averaged rows into a new matrix /// </summary> /// <param name="argMatrix"></param> /// <param name="selection"></param> public void PCalc(GeneralMatrix argMatrix, GeneralMatrix selection) { int n = argMatrix.ColumnDimension; GeneralMatrix sMatrix = new GeneralMatrix(argMatrix.ArrayCopy); double c=0.0; int i,j; for (i=0;i<sMatrix.ColumnDimension; i++) { c=0.0; for (j=0; j<sMatrix.RowDimension; j++) c+=sMatrix.GetElement(j,i); selection.SetElement(i,0,c); } for (i=0;i<sMatrix.ColumnDimension; i++) { for (j=0; j<sMatrix.RowDimension; j++) sMatrix.SetElement(j,i,sMatrix.GetElement(j,i)/selection.GetElement(i,0)); } for (i=0;i<sMatrix.RowDimension; i++) { c=0.0; for (j=0; j<sMatrix.ColumnDimension; j++) c+=sMatrix.GetElement(i,j); selection.SetElement(i,0,c/n); } }
/// <summary>Initializes this Metric object from metric matrix (called by constructor) </summary> /// <param name="m">The NxN metric matrix</param> private void init(DotNetMatrix.GeneralMatrix m) { if (!Util.IsSymmetric(m)) throw new Exception("The metric matrix must be symmetric"); m_matrix = m.Copy(); // System.Console.WriteLine("m_matrix: " + Util.ToString(m_matrix)); // compute eigen value decomposition m_eig = new DotNetMatrix.EigenvalueDecomposition(m_matrix); // System.Console.WriteLine("m_eig: " + Util.ToString(m_eig.GetV())); m_invEigMatrix = m_eig.GetV().Transpose(); m_eigenMetric = m_eig.RealEigenvalues; // { // DotNetMatrix.GeneralMatrix D = Util.Diagonal(m_eigenMetric); // DotNetMatrix.GeneralMatrix tmp = m_eig.GetV().Multiply(D).Multiply(m_invEigMatrix); // System.Console.WriteLine("input_matrix = " + Util.ToString(tmp)); // } m_isDiagonal = Util.IsDiagonal(m_matrix); if (!m_isDiagonal) { m_isEuclidean = m_isAntiEuclidean = false; } else { m_isEuclidean = m_isAntiEuclidean = true; for (int i = 0; i < m.RowDimension; i++) { if (m_matrix.GetElement(i, i) != 1.0) m_isEuclidean = false; if (m_matrix.GetElement(i, i) != -1.0) m_isAntiEuclidean = false; } } }
// Take SV columns from the U matrix private static GeneralMatrix GetSVs(GeneralMatrix matrix, int sv) { GeneralMatrix toPass = new GeneralMatrix(matrix.RowDimension, sv); for (int j = 0; j < sv; j++) { for (int i = 0; i < matrix.RowDimension; i++) { toPass.SetElement(i, j, matrix.GetElement(i, j)); } } return toPass; }
void PrintMatrix(GeneralMatrix matrix, StreamWriter output) { for (int j = 0; j < matrix.ColumnDimension; j++) { for (int i = 0; i < matrix.RowDimension; i++) { output.Write(matrix.GetElement(i, j) + ", "); } output.Write("\n"); } output.Write("\n"); }
protected RpropResult RPropLoop(double[] seed, bool precise) { //Console.WriteLine("RpropLoop"); InitialStepSize(); double[] curGradient; RpropResult ret = new RpropResult(); if (seed != null) { curGradient = InitialPointFromSeed(ret, seed); } else { curGradient = InitialPoint(ret); } double curUtil = ret.initialUtil; double oldUtil = curUtil; double[] formerGradient = new double[dim]; double[] curValue = new double[dim]; double[] testValue = new double[dim]; double lambda = 0.1; Tuple <double[], double> tup; Buffer.BlockCopy(ret.initialValue, 0, curValue, 0, sizeof(double) * dim); formerGradient = curGradient; //Buffer.BlockCopy(curGradient,0,formerGradient,0,sizeof(double)*dim); int itcounter = 0; int badcounter = 0; /*Console.WriteLine("Initial Sol:"); * for(int i=0; i<dim;i++) { * Console.Write("{0} ",curValue[i]); * } * Console.WriteLine(); * Console.WriteLine("Initial Util: {0}",curUtil); */ #if (GSOLVER_LOG) Log(curUtil, curValue); #endif int maxIter = 60; int maxBad = 30; double minStep = 1E-11; if (precise) { maxIter = 120; //110 maxBad = 60; //60 minStep = 1E-15; //15 } int convergendDims = 0; while (itcounter++ < maxIter && badcounter < maxBad) { convergendDims = 0; //First Order resp. approximated Second Order Gradient for (int i = 0; i < dim; i++) { if (curGradient[i] * formerGradient[i] > 0) { rpropStepWidth[i] *= 1.3; } else if (curGradient[i] * formerGradient[i] < 0) { rpropStepWidth[i] *= 0.5; } rpropStepWidth[i] = Math.Max(minStep, rpropStepWidth[i]); //rpropStepWidth[i] = Math.Max(0.000001,rpropStepWidth[i]); if (curUtil > 0.0) { //if (curGradient[i] > 0) curValue[i] += rpropStepWidth[i]; //else if (curGradient[i] < 0) curValue[i] -= rpropStepWidth[i]; } else { //linear assumption //curValue[i] += -curUtil/curGradient[i]; //quadratic assumption /*double ypSquare = 0; * for(int j=0; j<dim; j++) { * ypSquare += curGradient[j]*curGradient[j]; * } * double m=ypSquare/curUtil; * curValue[i] = -curGradient[i]/(2*m) + curValue[i];*/ } if (curValue[i] > limits[i, 1]) { curValue[i] = limits[i, 1]; } else if (curValue[i] < limits[i, 0]) { curValue[i] = limits[i, 0]; } if (rpropStepWidth[i] < rpropStepConvergenceThreshold[i]) { ++convergendDims; } } //Abort if all dimensions are converged if (!precise && convergendDims >= dim) { if (curUtil > ret.finalUtil) { ret.finalUtil = curUtil; Buffer.BlockCopy(curValue, 0, ret.finalValue, 0, sizeof(double) * dim); } return(ret); } // Conjugate Gradient if (curUtil < 0.5) { DotNetMatrix.GeneralMatrix X = new DotNetMatrix.GeneralMatrix(dim * 2 + 1, dim); DotNetMatrix.GeneralMatrix Y = new DotNetMatrix.GeneralMatrix(dim * 2 + 1, 1); for (int n = 0; n < dim; n++) { X.SetElement(0, n, curValue[n]); } Y.SetElement(0, 0, 0); for (int j = 0; j < dim * 2; j++) { for (int n = 0; n < dim; n++) { double rVal = rand.NextDouble() * rpropStepConvergenceThreshold[n] * 1000.0; testValue[n] = curValue[n] + rVal; } tup = term.Differentiate(testValue); for (int n = 0; n < dim; n++) { X.SetElement(j + 1, n, tup.Item1[n]); } Y.SetElement(j + 1, 0, tup.Item2 - curUtil); } DotNetMatrix.GeneralMatrix JJ = X.Transpose().Multiply(X); /*if(curUtil>oldUtil) lambda *= 10; * else lambda *= 0.1;*/ //DotNetMatrix.GeneralMatrix B = JJ.Add(GeneralMatrix.Identity(dim, dim).Multiply(lambda)).Inverse().Multiply(X.Transpose()).Multiply(Y); DotNetMatrix.GeneralMatrix B = JJ.Add(JJ.SVD().S.Multiply(lambda)).Inverse().Multiply(X.Transpose()).Multiply(Y); //DotNetMatrix.GeneralMatrix B = JJ.Inverse().Multiply(X.Transpose()).Multiply(Y); for (int j = 0; j < dim; j++) { curValue[j] += 0.01 * B.GetElement(j, 0); } Console.WriteLine(curUtil); Console.Write(B.Transpose()); Console.WriteLine(); } ///////////////////////// this.FEvals++; tup = term.Differentiate(curValue); bool allZero = true; for (int i = 0; i < dim; i++) { if (Double.IsNaN(tup.Item1[i])) { ret.aborted = true; #if (GSOLVER_LOG) LogStep(); #endif return(ret); } allZero &= (tup.Item1[i] == 0); } oldUtil = curUtil; curUtil = tup.Item2; formerGradient = curGradient; curGradient = tup.Item1; #if (GSOLVER_LOG) Log(curUtil, curValue); #endif //Console.WriteLine("CurUtil: {0} Final {1}",curUtil,ret.finalUtil); if (curUtil > ret.finalUtil) { badcounter = 0; //Math.Max(0,badcounter-1); //if (curUtil-ret.finalUtil < 0.00000000000001) { //Console.WriteLine("not better"); // badcounter++; //} else { //badcounter = 0; //} ret.finalUtil = curUtil; Buffer.BlockCopy(curValue, 0, ret.finalValue, 0, sizeof(double) * dim); //ret.finalValue = curValue; #if (ALWAYS_CHECK_THRESHOLD) if (curUtil > utilityThreshold) { return(ret); } #endif } else { //if (curUtil < ret.finalUtil || curUtil > 0) badcounter++; badcounter++; } if (allZero) { //Console.WriteLine("All Zero!"); /*Console.WriteLine("Util {0}",curUtil); * Console.Write("Vals: "); * for(int i=0; i < dim; i++) { * Console.Write("{0}\t",curValue[i]); * } * Console.WriteLine();*/ ret.aborted = false; #if (GSOLVER_LOG) LogStep(); #endif return(ret); } } #if (GSOLVER_LOG) LogStep(); #endif ret.aborted = false; return(ret); }
/// <summary> /// Set a submatrix. /// </summary> /// <param name = "r"> Array of row indices. /// </param> /// <param name = "j0"> Initial column index /// </param> /// <param name = "j1"> Final column index /// </param> /// <param name = "x"> A(r(:),j0:j1) /// </param> /// <exception cref = "System.IndexOutOfRangeException"> Submatrix indices /// </exception> public virtual void SetMatrix(int[] r, int j0, int j1, GeneralMatrix x) { try { for (int i = 0; i < r.Length; i++) { for (int j = j0; j <= j1; j++) { _a[r[i]][j] = x.GetElement(i, j - j0); } } } catch (IndexOutOfRangeException e) { throw new IndexOutOfRangeException("Submatrix indices", e); } }
public void TestChoiceMatrix() { double[][] mat1 = new double[][] { new double[] {1,2,3,4}, new double[] {0,1,2,5}, new double[] {0,0,1,5}, new double[] {0,0,0,1} }; GeneralMatrix gmat1 = new GeneralMatrix(mat1); gmat1 = AHPModel.ExpandUtility(gmat1); double[][] mat2 = new double[][] { new double[] {1,2,3,1}, new double[] {0,1,2,3}, new double[] {0,0,1,9}, new double[] {0,0,0,1} }; GeneralMatrix gmat2 = new GeneralMatrix(mat2); gmat2 = AHPModel.ExpandUtility(gmat2); double[][] mat3 = new double[][] { new double[] {1,4,2,2}, new double[] {0,1,2,4}, new double[] {0,0,1,4}, new double[] {0,0,0,1} }; GeneralMatrix gmat3 = new GeneralMatrix(mat3); gmat3 = AHPModel.ExpandUtility(gmat3); //3 criteria 4 choices AHPModel model = new AHPModel(3,4); model.AddCriterionRatedChoices(0,mat1); model.AddCriterionRatedChoices(1,mat2); model.AddCriterionRatedChoices(2,mat3); GeneralMatrix matrix = model.ChoiceMatrix; GeneralMatrix nmat1 = matrix.GetMatrix(0,3,0,3); GeneralMatrix nmat2 = matrix.GetMatrix(0,3,4,7); GeneralMatrix nmat3 = matrix.GetMatrix(0,3,8,11); bool alarm = false; int i,j; for (i=0;i<nmat1.ColumnDimension;i++) for (j=0;j<nmat1.RowDimension;j++) if (nmat1.GetElement(i,j)!=gmat1.GetElement(i,j)) alarm=true; for (i=0;i<nmat2.ColumnDimension;i++) for (j=0;j<nmat2.RowDimension;j++) if (nmat2.GetElement(i,j)!=gmat2.GetElement(i,j)) alarm=true; for (i=0;i<nmat3.ColumnDimension;i++) for (j=0;j<nmat3.RowDimension;j++) if (nmat3.GetElement(i,j)!=gmat3.GetElement(i,j)) alarm=true; Assert.IsFalse(alarm); }
public void ConstructionFromAnArrayAltersOriginalArrayValues() { double[][] array = avals; var b = new GeneralMatrix(array); //B still references avals underneath var temp = b.GetElement(0, 0); array[0][0] = 0.0; Assert.IsTrue((temp - b.GetElement(0, 0)) != 0.0); }
public void ConstructionFromAPackedArrayDoesNotAltersOriginalArrayValues() { double[] packedColumns = new[] { 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0, 11.0, 12.0 }; int numRows = 3; /* leading dimension of intended test Matrices */ var b = new GeneralMatrix(packedColumns, numRows); //B still references array underneath var temp = b.GetElement(0, 0); packedColumns[0] = 0.0; Assert.IsTrue((temp - b.GetElement(0, 0)) == 0.0); }
public void Negative_GetElementRow() { double[][] avals = {new double[]{1.0, 4.0, 7.0, 10.0}, new double[]{2.0, 5.0, 8.0, 11.0}, new double[]{3.0, 6.0, 9.0, 12.0}}; GeneralMatrix B = new GeneralMatrix(avals); tmp = B.GetElement(B.RowDimension, B.ColumnDimension - 1); }
public void SetElementCorrectlySetsAMatrixElement() { var matrix = new GeneralMatrix(avals); matrix.SetElement(ib,jb,0.0); Assert.AreEqual(0.0,matrix.GetElement(ib,jb)); }
public static void Main(System.String[] argv) { GeneralMatrix A, B, C, Z, O, I, R, S, X, SUB, M, T, SQ, DEF, SOL; int errorCount = 0; int warningCount = 0; double tmp; double[] columnwise = new double[]{1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0, 11.0, 12.0}; double[] rowwise = new double[]{1.0, 4.0, 7.0, 10.0, 2.0, 5.0, 8.0, 11.0, 3.0, 6.0, 9.0, 12.0}; double[][] avals = {new double[]{1.0, 4.0, 7.0, 10.0}, new double[]{2.0, 5.0, 8.0, 11.0}, new double[]{3.0, 6.0, 9.0, 12.0}}; double[][] rankdef = avals; double[][] tvals = {new double[]{1.0, 2.0, 3.0}, new double[]{4.0, 5.0, 6.0}, new double[]{7.0, 8.0, 9.0}, new double[]{10.0, 11.0, 12.0}}; double[][] subavals = {new double[]{5.0, 8.0, 11.0}, new double[]{6.0, 9.0, 12.0}}; double[][] rvals = {new double[]{1.0, 4.0, 7.0}, new double[]{2.0, 5.0, 8.0, 11.0}, new double[]{3.0, 6.0, 9.0, 12.0}}; double[][] pvals = {new double[]{1.0, 1.0, 1.0}, new double[]{1.0, 2.0, 3.0}, new double[]{1.0, 3.0, 6.0}}; double[][] ivals = {new double[]{1.0, 0.0, 0.0, 0.0}, new double[]{0.0, 1.0, 0.0, 0.0}, new double[]{0.0, 0.0, 1.0, 0.0}}; double[][] evals = {new double[]{0.0, 1.0, 0.0, 0.0}, new double[]{1.0, 0.0, 2e-7, 0.0}, new double[]{0.0, - 2e-7, 0.0, 1.0}, new double[]{0.0, 0.0, 1.0, 0.0}}; double[][] square = {new double[]{166.0, 188.0, 210.0}, new double[]{188.0, 214.0, 240.0}, new double[]{210.0, 240.0, 270.0}}; double[][] sqSolution = {new double[]{13.0}, new double[]{15.0}}; double[][] condmat = {new double[]{1.0, 3.0}, new double[]{7.0, 9.0}}; int rows = 3, cols = 4; int invalidld = 5; /* should trigger bad shape for construction with val */ int raggedr = 0; /* (raggedr,raggedc) should be out of bounds in ragged array */ int raggedc = 4; int validld = 3; /* leading dimension of intended test Matrices */ int nonconformld = 4; /* leading dimension which is valid, but nonconforming */ int ib = 1, ie = 2, jb = 1, je = 3; /* index ranges for sub GeneralMatrix */ int[] rowindexset = new int[]{1, 2}; int[] badrowindexset = new int[]{1, 3}; int[] columnindexset = new int[]{1, 2, 3}; int[] badcolumnindexset = new int[]{1, 2, 4}; double columnsummax = 33.0; double rowsummax = 30.0; double sumofdiagonals = 15; double sumofsquares = 650; /// <summary>Constructors and constructor-like methods: /// double[], int /// double[][] /// int, int /// int, int, double /// int, int, double[][] /// Create(double[][]) /// Random(int,int) /// Identity(int) /// /// </summary> print("\nTesting constructors and constructor-like methods...\n"); try { /// <summary>check that exception is thrown in packed constructor with invalid length *</summary> A = new GeneralMatrix(columnwise, invalidld); errorCount = try_failure(errorCount, "Catch invalid length in packed constructor... ", "exception not thrown for invalid input"); } catch (System.ArgumentException e) { try_success("Catch invalid length in packed constructor... ", e.Message); } try { /// <summary>check that exception is thrown in default constructor /// if input array is 'ragged' * /// </summary> A = new GeneralMatrix(rvals); tmp = A.GetElement(raggedr, raggedc); } catch (System.ArgumentException e) { try_success("Catch ragged input to default constructor... ", e.Message); } catch (System.IndexOutOfRangeException e) { errorCount = try_failure(errorCount, "Catch ragged input to constructor... ", "exception not thrown in construction...ArrayIndexOutOfBoundsException thrown later"); System.Console.Out.WriteLine(e.Message); } try { /// <summary>check that exception is thrown in Create /// if input array is 'ragged' * /// </summary> A = GeneralMatrix.Create(rvals); tmp = A.GetElement(raggedr, raggedc); } catch (System.ArgumentException e) { try_success("Catch ragged input to Create... ", e.Message); System.Console.Out.WriteLine(e.Message); } catch (System.IndexOutOfRangeException e) { errorCount = try_failure(errorCount, "Catch ragged input to Create... ", "exception not thrown in construction...ArrayIndexOutOfBoundsException thrown later"); System.Console.Out.WriteLine(e.Message); } A = new GeneralMatrix(columnwise, validld); B = new GeneralMatrix(avals); tmp = B.GetElement(0, 0); avals[0][0] = 0.0; C = B.Subtract(A); avals[0][0] = tmp; B = GeneralMatrix.Create(avals); tmp = B.GetElement(0, 0); avals[0][0] = 0.0; if ((tmp - B.GetElement(0, 0)) != 0.0) { /// <summary>check that Create behaves properly *</summary> errorCount = try_failure(errorCount, "Create... ", "Copy not effected... data visible outside"); } else { try_success("Create... ", ""); } avals[0][0] = columnwise[0]; I = new GeneralMatrix(ivals); try { check(I, GeneralMatrix.Identity(3, 4)); try_success("Identity... ", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "Identity... ", "Identity GeneralMatrix not successfully created"); System.Console.Out.WriteLine(e.Message); } /// <summary>Access Methods: /// getColumnDimension() /// getRowDimension() /// getArray() /// getArrayCopy() /// getColumnPackedCopy() /// getRowPackedCopy() /// get(int,int) /// GetMatrix(int,int,int,int) /// GetMatrix(int,int,int[]) /// GetMatrix(int[],int,int) /// GetMatrix(int[],int[]) /// set(int,int,double) /// SetMatrix(int,int,int,int,GeneralMatrix) /// SetMatrix(int,int,int[],GeneralMatrix) /// SetMatrix(int[],int,int,GeneralMatrix) /// SetMatrix(int[],int[],GeneralMatrix) /// /// </summary> print("\nTesting access methods...\n"); /// <summary>Various get methods: /// /// </summary> B = new GeneralMatrix(avals); if (B.RowDimension != rows) { errorCount = try_failure(errorCount, "getRowDimension... ", ""); } else { try_success("getRowDimension... ", ""); } if (B.ColumnDimension != cols) { errorCount = try_failure(errorCount, "getColumnDimension... ", ""); } else { try_success("getColumnDimension... ", ""); } B = new GeneralMatrix(avals); double[][] barray = B.Array; if (barray != avals) { errorCount = try_failure(errorCount, "getArray... ", ""); } else { try_success("getArray... ", ""); } barray = B.ArrayCopy; if (barray == avals) { errorCount = try_failure(errorCount, "getArrayCopy... ", "data not (deep) copied"); } try { check(barray, avals); try_success("getArrayCopy... ", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "getArrayCopy... ", "data not successfully (deep) copied"); System.Console.Out.WriteLine(e.Message); } double[] bpacked = B.ColumnPackedCopy; try { check(bpacked, columnwise); try_success("getColumnPackedCopy... ", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "getColumnPackedCopy... ", "data not successfully (deep) copied by columns"); System.Console.Out.WriteLine(e.Message); } bpacked = B.RowPackedCopy; try { check(bpacked, rowwise); try_success("getRowPackedCopy... ", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "getRowPackedCopy... ", "data not successfully (deep) copied by rows"); System.Console.Out.WriteLine(e.Message); } try { tmp = B.GetElement(B.RowDimension, B.ColumnDimension - 1); errorCount = try_failure(errorCount, "get(int,int)... ", "OutOfBoundsException expected but not thrown"); } catch (System.IndexOutOfRangeException e) { System.Console.Out.WriteLine(e.Message); try { tmp = B.GetElement(B.RowDimension - 1, B.ColumnDimension); errorCount = try_failure(errorCount, "get(int,int)... ", "OutOfBoundsException expected but not thrown"); } catch (System.IndexOutOfRangeException e1) { try_success("get(int,int)... OutofBoundsException... ", ""); System.Console.Out.WriteLine(e1.Message); } } catch (System.ArgumentException e1) { errorCount = try_failure(errorCount, "get(int,int)... ", "OutOfBoundsException expected but not thrown"); System.Console.Out.WriteLine(e1.Message); } try { if (B.GetElement(B.RowDimension - 1, B.ColumnDimension - 1) != avals[B.RowDimension - 1][B.ColumnDimension - 1]) { errorCount = try_failure(errorCount, "get(int,int)... ", "GeneralMatrix entry (i,j) not successfully retreived"); } else { try_success("get(int,int)... ", ""); } } catch (System.IndexOutOfRangeException e) { errorCount = try_failure(errorCount, "get(int,int)... ", "Unexpected ArrayIndexOutOfBoundsException"); System.Console.Out.WriteLine(e.Message); } SUB = new GeneralMatrix(subavals); try { M = B.GetMatrix(ib, ie + B.RowDimension + 1, jb, je); errorCount = try_failure(errorCount, "GetMatrix(int,int,int,int)... ", "ArrayIndexOutOfBoundsException expected but not thrown"); } catch (System.IndexOutOfRangeException e) { System.Console.Out.WriteLine(e.Message); try { M = B.GetMatrix(ib, ie, jb, je + B.ColumnDimension + 1); errorCount = try_failure(errorCount, "GetMatrix(int,int,int,int)... ", "ArrayIndexOutOfBoundsException expected but not thrown"); } catch (System.IndexOutOfRangeException e1) { try_success("GetMatrix(int,int,int,int)... ArrayIndexOutOfBoundsException... ", ""); System.Console.Out.WriteLine(e1.Message); } } catch (System.ArgumentException e1) { errorCount = try_failure(errorCount, "GetMatrix(int,int,int,int)... ", "ArrayIndexOutOfBoundsException expected but not thrown"); System.Console.Out.WriteLine(e1.Message); } try { M = B.GetMatrix(ib, ie, jb, je); try { check(SUB, M); try_success("GetMatrix(int,int,int,int)... ", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "GetMatrix(int,int,int,int)... ", "submatrix not successfully retreived"); System.Console.Out.WriteLine(e.Message); } } catch (System.IndexOutOfRangeException e) { errorCount = try_failure(errorCount, "GetMatrix(int,int,int,int)... ", "Unexpected ArrayIndexOutOfBoundsException"); System.Console.Out.WriteLine(e.Message); } try { M = B.GetMatrix(ib, ie, badcolumnindexset); errorCount = try_failure(errorCount, "GetMatrix(int,int,int[])... ", "ArrayIndexOutOfBoundsException expected but not thrown"); } catch (System.IndexOutOfRangeException e) { System.Console.Out.WriteLine(e.Message); try { M = B.GetMatrix(ib, ie + B.RowDimension + 1, columnindexset); errorCount = try_failure(errorCount, "GetMatrix(int,int,int[])... ", "ArrayIndexOutOfBoundsException expected but not thrown"); } catch (System.IndexOutOfRangeException e1) { try_success("GetMatrix(int,int,int[])... ArrayIndexOutOfBoundsException... ", ""); System.Console.Out.WriteLine(e1.Message); } } catch (System.ArgumentException e1) { errorCount = try_failure(errorCount, "GetMatrix(int,int,int[])... ", "ArrayIndexOutOfBoundsException expected but not thrown"); System.Console.Out.WriteLine(e1.Message); } try { M = B.GetMatrix(ib, ie, columnindexset); try { check(SUB, M); try_success("GetMatrix(int,int,int[])... ", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "GetMatrix(int,int,int[])... ", "submatrix not successfully retreived"); System.Console.Out.WriteLine(e.Message); } } catch (System.IndexOutOfRangeException e) { errorCount = try_failure(errorCount, "GetMatrix(int,int,int[])... ", "Unexpected ArrayIndexOutOfBoundsException"); System.Console.Out.WriteLine(e.Message); } try { M = B.GetMatrix(badrowindexset, jb, je); errorCount = try_failure(errorCount, "GetMatrix(int[],int,int)... ", "ArrayIndexOutOfBoundsException expected but not thrown"); } catch (System.IndexOutOfRangeException e) { System.Console.Out.WriteLine(e.Message); try { M = B.GetMatrix(rowindexset, jb, je + B.ColumnDimension + 1); errorCount = try_failure(errorCount, "GetMatrix(int[],int,int)... ", "ArrayIndexOutOfBoundsException expected but not thrown"); } catch (System.IndexOutOfRangeException e1) { try_success("GetMatrix(int[],int,int)... ArrayIndexOutOfBoundsException... ", ""); System.Console.Out.WriteLine(e1.Message); } } catch (System.ArgumentException e1) { errorCount = try_failure(errorCount, "GetMatrix(int[],int,int)... ", "ArrayIndexOutOfBoundsException expected but not thrown"); System.Console.Out.WriteLine(e1.Message); } try { M = B.GetMatrix(rowindexset, jb, je); try { check(SUB, M); try_success("GetMatrix(int[],int,int)... ", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "GetMatrix(int[],int,int)... ", "submatrix not successfully retreived"); System.Console.Out.WriteLine(e.Message); } } catch (System.IndexOutOfRangeException e) { errorCount = try_failure(errorCount, "GetMatrix(int[],int,int)... ", "Unexpected ArrayIndexOutOfBoundsException"); System.Console.Out.WriteLine(e.Message); } try { M = B.GetMatrix(badrowindexset, columnindexset); errorCount = try_failure(errorCount, "GetMatrix(int[],int[])... ", "ArrayIndexOutOfBoundsException expected but not thrown"); } catch (System.IndexOutOfRangeException e) { System.Console.Out.WriteLine(e.Message); try { M = B.GetMatrix(rowindexset, badcolumnindexset); errorCount = try_failure(errorCount, "GetMatrix(int[],int[])... ", "ArrayIndexOutOfBoundsException expected but not thrown"); } catch (System.IndexOutOfRangeException e1) { try_success("GetMatrix(int[],int[])... ArrayIndexOutOfBoundsException... ", ""); System.Console.Out.WriteLine(e1.Message); } } catch (System.ArgumentException e1) { errorCount = try_failure(errorCount, "GetMatrix(int[],int[])... ", "ArrayIndexOutOfBoundsException expected but not thrown"); System.Console.Out.WriteLine(e1.Message); } try { M = B.GetMatrix(rowindexset, columnindexset); try { check(SUB, M); try_success("GetMatrix(int[],int[])... ", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "GetMatrix(int[],int[])... ", "submatrix not successfully retreived"); System.Console.Out.WriteLine(e.Message); } } catch (System.IndexOutOfRangeException e) { errorCount = try_failure(errorCount, "GetMatrix(int[],int[])... ", "Unexpected ArrayIndexOutOfBoundsException"); System.Console.Out.WriteLine(e.Message); } /// <summary>Various set methods: /// /// </summary> try { B.SetElement(B.RowDimension, B.ColumnDimension - 1, 0.0); errorCount = try_failure(errorCount, "set(int,int,double)... ", "OutOfBoundsException expected but not thrown"); } catch (System.IndexOutOfRangeException e) { System.Console.Out.WriteLine(e.Message); try { B.SetElement(B.RowDimension - 1, B.ColumnDimension, 0.0); errorCount = try_failure(errorCount, "set(int,int,double)... ", "OutOfBoundsException expected but not thrown"); } catch (System.IndexOutOfRangeException e1) { try_success("set(int,int,double)... OutofBoundsException... ", ""); System.Console.Out.WriteLine(e1.Message); } } catch (System.ArgumentException e1) { errorCount = try_failure(errorCount, "set(int,int,double)... ", "OutOfBoundsException expected but not thrown"); System.Console.Out.WriteLine(e1.Message); } try { B.SetElement(ib, jb, 0.0); tmp = B.GetElement(ib, jb); try { check(tmp, 0.0); try_success("set(int,int,double)... ", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "set(int,int,double)... ", "GeneralMatrix element not successfully set"); System.Console.Out.WriteLine(e.Message); } } catch (System.IndexOutOfRangeException e1) { errorCount = try_failure(errorCount, "set(int,int,double)... ", "Unexpected ArrayIndexOutOfBoundsException"); System.Console.Out.WriteLine(e1.Message); } M = new GeneralMatrix(2, 3, 0.0); try { B.SetMatrix(ib, ie + B.RowDimension + 1, jb, je, M); errorCount = try_failure(errorCount, "SetMatrix(int,int,int,int,GeneralMatrix)... ", "ArrayIndexOutOfBoundsException expected but not thrown"); } catch (System.IndexOutOfRangeException e) { System.Console.Out.WriteLine(e.Message); try { B.SetMatrix(ib, ie, jb, je + B.ColumnDimension + 1, M); errorCount = try_failure(errorCount, "SetMatrix(int,int,int,int,GeneralMatrix)... ", "ArrayIndexOutOfBoundsException expected but not thrown"); } catch (System.IndexOutOfRangeException e1) { try_success("SetMatrix(int,int,int,int,GeneralMatrix)... ArrayIndexOutOfBoundsException... ", ""); System.Console.Out.WriteLine(e1.Message); } } catch (System.ArgumentException e1) { errorCount = try_failure(errorCount, "SetMatrix(int,int,int,int,GeneralMatrix)... ", "ArrayIndexOutOfBoundsException expected but not thrown"); System.Console.Out.WriteLine(e1.Message); } try { B.SetMatrix(ib, ie, jb, je, M); try { check(M.Subtract(B.GetMatrix(ib, ie, jb, je)), M); try_success("SetMatrix(int,int,int,int,GeneralMatrix)... ", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "SetMatrix(int,int,int,int,GeneralMatrix)... ", "submatrix not successfully set"); System.Console.Out.WriteLine(e.Message); } B.SetMatrix(ib, ie, jb, je, SUB); } catch (System.IndexOutOfRangeException e1) { errorCount = try_failure(errorCount, "SetMatrix(int,int,int,int,GeneralMatrix)... ", "Unexpected ArrayIndexOutOfBoundsException"); System.Console.Out.WriteLine(e1.Message); } try { B.SetMatrix(ib, ie + B.RowDimension + 1, columnindexset, M); errorCount = try_failure(errorCount, "SetMatrix(int,int,int[],GeneralMatrix)... ", "ArrayIndexOutOfBoundsException expected but not thrown"); } catch (System.IndexOutOfRangeException e) { System.Console.Out.WriteLine(e.Message); try { B.SetMatrix(ib, ie, badcolumnindexset, M); errorCount = try_failure(errorCount, "SetMatrix(int,int,int[],GeneralMatrix)... ", "ArrayIndexOutOfBoundsException expected but not thrown"); } catch (System.IndexOutOfRangeException e1) { try_success("SetMatrix(int,int,int[],GeneralMatrix)... ArrayIndexOutOfBoundsException... ", ""); System.Console.Out.WriteLine(e1.Message); } } catch (System.ArgumentException e1) { errorCount = try_failure(errorCount, "SetMatrix(int,int,int[],GeneralMatrix)... ", "ArrayIndexOutOfBoundsException expected but not thrown"); System.Console.Out.WriteLine(e1.Message); } try { B.SetMatrix(ib, ie, columnindexset, M); try { check(M.Subtract(B.GetMatrix(ib, ie, columnindexset)), M); try_success("SetMatrix(int,int,int[],GeneralMatrix)... ", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "SetMatrix(int,int,int[],GeneralMatrix)... ", "submatrix not successfully set"); System.Console.Out.WriteLine(e.Message); } B.SetMatrix(ib, ie, jb, je, SUB); } catch (System.IndexOutOfRangeException e1) { errorCount = try_failure(errorCount, "SetMatrix(int,int,int[],GeneralMatrix)... ", "Unexpected ArrayIndexOutOfBoundsException"); System.Console.Out.WriteLine(e1.Message); } try { B.SetMatrix(rowindexset, jb, je + B.ColumnDimension + 1, M); errorCount = try_failure(errorCount, "SetMatrix(int[],int,int,GeneralMatrix)... ", "ArrayIndexOutOfBoundsException expected but not thrown"); } catch (System.IndexOutOfRangeException e) { System.Console.Out.WriteLine(e.Message); try { B.SetMatrix(badrowindexset, jb, je, M); errorCount = try_failure(errorCount, "SetMatrix(int[],int,int,GeneralMatrix)... ", "ArrayIndexOutOfBoundsException expected but not thrown"); } catch (System.IndexOutOfRangeException e1) { try_success("SetMatrix(int[],int,int,GeneralMatrix)... ArrayIndexOutOfBoundsException... ", ""); System.Console.Out.WriteLine(e1.Message); } } catch (System.ArgumentException e1) { errorCount = try_failure(errorCount, "SetMatrix(int[],int,int,GeneralMatrix)... ", "ArrayIndexOutOfBoundsException expected but not thrown"); System.Console.Out.WriteLine(e1.Message); } try { B.SetMatrix(rowindexset, jb, je, M); try { check(M.Subtract(B.GetMatrix(rowindexset, jb, je)), M); try_success("SetMatrix(int[],int,int,GeneralMatrix)... ", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "SetMatrix(int[],int,int,GeneralMatrix)... ", "submatrix not successfully set"); System.Console.Out.WriteLine(e.Message); } B.SetMatrix(ib, ie, jb, je, SUB); } catch (System.IndexOutOfRangeException e1) { errorCount = try_failure(errorCount, "SetMatrix(int[],int,int,GeneralMatrix)... ", "Unexpected ArrayIndexOutOfBoundsException"); System.Console.Out.WriteLine(e1.Message); } try { B.SetMatrix(rowindexset, badcolumnindexset, M); errorCount = try_failure(errorCount, "SetMatrix(int[],int[],GeneralMatrix)... ", "ArrayIndexOutOfBoundsException expected but not thrown"); } catch (System.IndexOutOfRangeException e) { System.Console.Out.WriteLine(e.Message); try { B.SetMatrix(badrowindexset, columnindexset, M); errorCount = try_failure(errorCount, "SetMatrix(int[],int[],GeneralMatrix)... ", "ArrayIndexOutOfBoundsException expected but not thrown"); } catch (System.IndexOutOfRangeException e1) { try_success("SetMatrix(int[],int[],GeneralMatrix)... ArrayIndexOutOfBoundsException... ", ""); System.Console.Out.WriteLine(e1.Message); } } catch (System.ArgumentException e1) { errorCount = try_failure(errorCount, "SetMatrix(int[],int[],GeneralMatrix)... ", "ArrayIndexOutOfBoundsException expected but not thrown"); System.Console.Out.WriteLine(e1.Message); } try { B.SetMatrix(rowindexset, columnindexset, M); try { check(M.Subtract(B.GetMatrix(rowindexset, columnindexset)), M); try_success("SetMatrix(int[],int[],GeneralMatrix)... ", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "SetMatrix(int[],int[],GeneralMatrix)... ", "submatrix not successfully set"); System.Console.Out.WriteLine(e.Message); } } catch (System.IndexOutOfRangeException e1) { errorCount = try_failure(errorCount, "SetMatrix(int[],int[],GeneralMatrix)... ", "Unexpected ArrayIndexOutOfBoundsException"); System.Console.Out.WriteLine(e1.Message); } /// <summary>Array-like methods: /// Subtract /// SubtractEquals /// Add /// AddEquals /// ArrayLeftDivide /// ArrayLeftDivideEquals /// ArrayRightDivide /// ArrayRightDivideEquals /// arrayTimes /// ArrayMultiplyEquals /// uminus /// /// </summary> print("\nTesting array-like methods...\n"); S = new GeneralMatrix(columnwise, nonconformld); R = GeneralMatrix.Random(A.RowDimension, A.ColumnDimension); A = R; try { S = A.Subtract(S); errorCount = try_failure(errorCount, "Subtract conformance check... ", "nonconformance not raised"); } catch (System.ArgumentException e) { try_success("Subtract conformance check... ", ""); System.Console.Out.WriteLine(e.Message); } if (A.Subtract(R).Norm1() != 0.0) { errorCount = try_failure(errorCount, "Subtract... ", "(difference of identical Matrices is nonzero,\nSubsequent use of Subtract should be suspect)"); } else { try_success("Subtract... ", ""); } A = R.Copy(); A.SubtractEquals(R); Z = new GeneralMatrix(A.RowDimension, A.ColumnDimension); try { A.SubtractEquals(S); errorCount = try_failure(errorCount, "SubtractEquals conformance check... ", "nonconformance not raised"); } catch (System.ArgumentException e) { try_success("SubtractEquals conformance check... ", ""); System.Console.Out.WriteLine(e.Message); } if (A.Subtract(Z).Norm1() != 0.0) { errorCount = try_failure(errorCount, "SubtractEquals... ", "(difference of identical Matrices is nonzero,\nSubsequent use of Subtract should be suspect)"); } else { try_success("SubtractEquals... ", ""); } A = R.Copy(); B = GeneralMatrix.Random(A.RowDimension, A.ColumnDimension); C = A.Subtract(B); try { S = A.Add(S); errorCount = try_failure(errorCount, "Add conformance check... ", "nonconformance not raised"); } catch (System.ArgumentException e) { try_success("Add conformance check... ", ""); System.Console.Out.WriteLine(e.Message); } try { check(C.Add(B), A); try_success("Add... ", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "Add... ", "(C = A - B, but C + B != A)"); System.Console.Out.WriteLine(e.Message); } C = A.Subtract(B); C.AddEquals(B); try { A.AddEquals(S); errorCount = try_failure(errorCount, "AddEquals conformance check... ", "nonconformance not raised"); } catch (System.ArgumentException e) { try_success("AddEquals conformance check... ", ""); System.Console.Out.WriteLine(e.Message); } try { check(C, A); try_success("AddEquals... ", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "AddEquals... ", "(C = A - B, but C = C + B != A)"); System.Console.Out.WriteLine(e.Message); } A = R.UnaryMinus(); try { check(A.Add(R), Z); try_success("UnaryMinus... ", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "uminus... ", "(-A + A != zeros)"); System.Console.Out.WriteLine(e.Message); } A = R.Copy(); O = new GeneralMatrix(A.RowDimension, A.ColumnDimension, 1.0); C = A.ArrayLeftDivide(R); try { S = A.ArrayLeftDivide(S); errorCount = try_failure(errorCount, "ArrayLeftDivide conformance check... ", "nonconformance not raised"); } catch (System.ArgumentException e) { try_success("ArrayLeftDivide conformance check... ", ""); System.Console.Out.WriteLine(e.Message); } try { check(C, O); try_success("ArrayLeftDivide... ", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "ArrayLeftDivide... ", "(M.\\M != ones)"); System.Console.Out.WriteLine(e.Message); } try { A.ArrayLeftDivideEquals(S); errorCount = try_failure(errorCount, "ArrayLeftDivideEquals conformance check... ", "nonconformance not raised"); } catch (System.ArgumentException e) { try_success("ArrayLeftDivideEquals conformance check... ", ""); System.Console.Out.WriteLine(e.Message); } A.ArrayLeftDivideEquals(R); try { check(A, O); try_success("ArrayLeftDivideEquals... ", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "ArrayLeftDivideEquals... ", "(M.\\M != ones)"); System.Console.Out.WriteLine(e.Message); } A = R.Copy(); try { A.ArrayRightDivide(S); errorCount = try_failure(errorCount, "ArrayRightDivide conformance check... ", "nonconformance not raised"); } catch (System.ArgumentException e) { try_success("ArrayRightDivide conformance check... ", ""); System.Console.Out.WriteLine(e.Message); } C = A.ArrayRightDivide(R); try { check(C, O); try_success("ArrayRightDivide... ", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "ArrayRightDivide... ", "(M./M != ones)"); System.Console.Out.WriteLine(e.Message); } try { A.ArrayRightDivideEquals(S); errorCount = try_failure(errorCount, "ArrayRightDivideEquals conformance check... ", "nonconformance not raised"); } catch (System.ArgumentException e) { try_success("ArrayRightDivideEquals conformance check... ", ""); System.Console.Out.WriteLine(e.Message); } A.ArrayRightDivideEquals(R); try { check(A, O); try_success("ArrayRightDivideEquals... ", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "ArrayRightDivideEquals... ", "(M./M != ones)"); System.Console.Out.WriteLine(e.Message); } A = R.Copy(); B = GeneralMatrix.Random(A.RowDimension, A.ColumnDimension); try { S = A.ArrayMultiply(S); errorCount = try_failure(errorCount, "arrayTimes conformance check... ", "nonconformance not raised"); } catch (System.ArgumentException e) { try_success("arrayTimes conformance check... ", ""); System.Console.Out.WriteLine(e.Message); } C = A.ArrayMultiply(B); try { check(C.ArrayRightDivideEquals(B), A); try_success("arrayTimes... ", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "arrayTimes... ", "(A = R, C = A.*B, but C./B != A)"); System.Console.Out.WriteLine(e.Message); } try { A.ArrayMultiplyEquals(S); errorCount = try_failure(errorCount, "ArrayMultiplyEquals conformance check... ", "nonconformance not raised"); } catch (System.ArgumentException e) { try_success("ArrayMultiplyEquals conformance check... ", ""); System.Console.Out.WriteLine(e.Message); } A.ArrayMultiplyEquals(B); try { check(A.ArrayRightDivideEquals(B), R); try_success("ArrayMultiplyEquals... ", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "ArrayMultiplyEquals... ", "(A = R, A = A.*B, but A./B != R)"); System.Console.Out.WriteLine(e.Message); } /// <summary>LA methods: /// Transpose /// Multiply /// Condition /// Rank /// Determinant /// trace /// Norm1 /// norm2 /// normF /// normInf /// Solve /// solveTranspose /// Inverse /// chol /// Eigen /// lu /// qr /// svd /// /// </summary> print("\nTesting linear algebra methods...\n"); A = new GeneralMatrix(columnwise, 3); T = new GeneralMatrix(tvals); T = A.Transpose(); try { check(A.Transpose(), T); try_success("Transpose...", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "Transpose()...", "Transpose unsuccessful"); System.Console.Out.WriteLine(e.Message); } A.Transpose(); try { check(A.Norm1(), columnsummax); try_success("Norm1...", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "Norm1()...", "incorrect norm calculation"); System.Console.Out.WriteLine(e.Message); } try { check(A.NormInf(), rowsummax); try_success("normInf()...", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "normInf()...", "incorrect norm calculation"); System.Console.Out.WriteLine(e.Message); } try { check(A.NormF(), System.Math.Sqrt(sumofsquares)); try_success("normF...", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "normF()...", "incorrect norm calculation"); System.Console.Out.WriteLine(e.Message); } try { check(A.Trace(), sumofdiagonals); try_success("trace()...", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "trace()...", "incorrect trace calculation"); System.Console.Out.WriteLine(e.Message); } try { check(A.GetMatrix(0, A.RowDimension - 1, 0, A.RowDimension - 1).Determinant(), 0.0); try_success("Determinant()...", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "Determinant()...", "incorrect determinant calculation"); System.Console.Out.WriteLine(e.Message); } SQ = new GeneralMatrix(square); try { check(A.Multiply(A.Transpose()), SQ); try_success("Multiply(GeneralMatrix)...", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "Multiply(GeneralMatrix)...", "incorrect GeneralMatrix-GeneralMatrix product calculation"); System.Console.Out.WriteLine(e.Message); } try { check(A.Multiply(0.0), Z); try_success("Multiply(double)...", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "Multiply(double)...", "incorrect GeneralMatrix-scalar product calculation"); System.Console.Out.WriteLine(e.Message); } A = new GeneralMatrix(columnwise, 4); QRDecomposition QR = A.QRD(); R = QR.R; try { check(A, QR.Q.Multiply(R)); try_success("QRDecomposition...", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "QRDecomposition...", "incorrect QR decomposition calculation"); System.Console.Out.WriteLine(e.Message); } SingularValueDecomposition SVD = A.SVD(); try { check(A, SVD.GetU().Multiply(SVD.S.Multiply(SVD.GetV().Transpose()))); try_success("SingularValueDecomposition...", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "SingularValueDecomposition...", "incorrect singular value decomposition calculation"); System.Console.Out.WriteLine(e.Message); } DEF = new GeneralMatrix(rankdef); try { check(DEF.Rank(), System.Math.Min(DEF.RowDimension, DEF.ColumnDimension) - 1); try_success("Rank()...", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "Rank()...", "incorrect Rank calculation"); System.Console.Out.WriteLine(e.Message); } B = new GeneralMatrix(condmat); SVD = B.SVD(); double[] singularvalues = SVD.SingularValues; try { check(B.Condition(), singularvalues[0] / singularvalues[System.Math.Min(B.RowDimension, B.ColumnDimension) - 1]); try_success("Condition()...", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "Condition()...", "incorrect condition number calculation"); System.Console.Out.WriteLine(e.Message); } int n = A.ColumnDimension; A = A.GetMatrix(0, n - 1, 0, n - 1); A.SetElement(0, 0, 0.0); LUDecomposition LU = A.LUD(); try { check(A.GetMatrix(LU.Pivot, 0, n - 1), LU.L.Multiply(LU.U)); try_success("LUDecomposition...", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "LUDecomposition...", "incorrect LU decomposition calculation"); System.Console.Out.WriteLine(e.Message); } X = A.Inverse(); try { check(A.Multiply(X), GeneralMatrix.Identity(3, 3)); try_success("Inverse()...", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "Inverse()...", "incorrect Inverse calculation"); System.Console.Out.WriteLine(e.Message); } O = new GeneralMatrix(SUB.RowDimension, 1, 1.0); SOL = new GeneralMatrix(sqSolution); SQ = SUB.GetMatrix(0, SUB.RowDimension - 1, 0, SUB.RowDimension - 1); try { check(SQ.Solve(SOL), O); try_success("Solve()...", ""); } catch (System.ArgumentException e1) { errorCount = try_failure(errorCount, "Solve()...", e1.Message); System.Console.Out.WriteLine(e1.Message); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "Solve()...", e.Message); System.Console.Out.WriteLine(e.Message); } A = new GeneralMatrix(pvals); CholeskyDecomposition Chol = A.chol(); GeneralMatrix L = Chol.GetL(); try { check(A, L.Multiply(L.Transpose())); try_success("CholeskyDecomposition...", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "CholeskyDecomposition...", "incorrect Cholesky decomposition calculation"); System.Console.Out.WriteLine(e.Message); } X = Chol.Solve(GeneralMatrix.Identity(3, 3)); try { check(A.Multiply(X), GeneralMatrix.Identity(3, 3)); try_success("CholeskyDecomposition Solve()...", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "CholeskyDecomposition Solve()...", "incorrect Choleskydecomposition Solve calculation"); System.Console.Out.WriteLine(e.Message); } EigenvalueDecomposition Eig = A.Eigen(); GeneralMatrix D = Eig.D; GeneralMatrix V = Eig.GetV(); try { check(A.Multiply(V), V.Multiply(D)); try_success("EigenvalueDecomposition (symmetric)...", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "EigenvalueDecomposition (symmetric)...", "incorrect symmetric Eigenvalue decomposition calculation"); System.Console.Out.WriteLine(e.Message); } A = new GeneralMatrix(evals); Eig = A.Eigen(); D = Eig.D; V = Eig.GetV(); try { check(A.Multiply(V), V.Multiply(D)); try_success("EigenvalueDecomposition (nonsymmetric)...", ""); } catch (System.SystemException e) { errorCount = try_failure(errorCount, "EigenvalueDecomposition (nonsymmetric)...", "incorrect nonsymmetric Eigenvalue decomposition calculation"); System.Console.Out.WriteLine(e.Message); } print("\nTestMatrix completed.\n"); print("Total errors reported: " + System.Convert.ToString(errorCount) + "\n"); print("Total warnings reported: " + System.Convert.ToString(warningCount) + "\n"); }
public void ZeroBasedGetElementMethodReturnsCorrectElement() { var matrix = new GeneralMatrix(avals); var actual = matrix.GetElement(matrix.RowDimension-1, matrix.ColumnDimension-1); Assert.AreEqual(12.0,actual); }
/// <summary>Set a submatrix.</summary> /// <param name="i0"> Initial row index /// </param> /// <param name="i1"> Final row index /// </param> /// <param name="j0"> Initial column index /// </param> /// <param name="j1"> Final column index /// </param> /// <param name="X"> A(i0:i1,j0:j1) /// </param> /// <exception cref="System.IndexOutOfRangeException"> Submatrix indices /// </exception> public virtual void SetMatrix(int i0, int i1, int j0, int j1, GeneralMatrix X) { try { for (int i = i0; i <= i1; i++) { for (int j = j0; j <= j1; j++) { A[i][j] = X.GetElement(i - i0, j - j0); } } } catch (System.IndexOutOfRangeException e) { throw new System.IndexOutOfRangeException("Submatrix indices", e); } }
private double[][] calculateR() { sourceCount = sourceASTL.abstractLatticeGraph.VertexCount; targetCount = targetASTL.abstractLatticeGraph.VertexCount; GeneralMatrix newR = new GeneralMatrix(sourceCount * targetCount, 1); GeneralMatrix myA = new GeneralMatrix(A); //initial R for (int i = 0; i < newR.RowDimension; i++) { newR.SetElement(i, 0, 1.0 / newR.RowDimension); } //printMatrix(vector2Matrix(newR.Array, 0)); //move similarity matrix to a double index vector GeneralMatrix newSim = new GeneralMatrix(sourceCount * targetCount, 1); for (int i = 0; i < sourceCount; i++) for (int j = 0; j < targetCount; j++) { newSim.SetElement(targetCount * i + j, 0, simMatrix[i, j]); } //move structure similarity matrix to a double index vector GeneralMatrix newStructSim = new GeneralMatrix(sourceCount * targetCount, 1); //for (int i = 0; i < sourceCount; i++) // for (int j = 0; j < targetCount; j++) // { // newStructSim.SetElement(targetCount * i + j, 0, structSimMatrix[i, j]); // } //calculate R using power method (eigen vector) //========================================== int count = 0; while (count < 50) { //R = aAR + (1-2a)E1 + (1-2a)E2 where a = 0.333 //newR = (((myA.Multiply(newR)).Multiply(0.333)).Add(newStructSim.Multiply(0.333))).Add(newSim.Multiply(0.333));//ommited to have name similarity only 17/4/2012 newR = (((myA.Multiply(newR)).Multiply(0.5)).Add(newSim.Multiply(0.5))); double sum = 0; for (int i = 0; i < newR.RowDimension; i++) for (int j = 0; j < newR.ColumnDimension; j++) sum += newR.GetElement(i, j); newR = newR.Multiply(1.0 / sum); count++; } return newR.Array; }