protected override GeneralMatrix CalculateNextHessianApproximation(GeneralMatrix previousH, double[]prevX, double[]curX, double[]prevGrad, double[]curGrad) { GeneralMatrix currentH = new GeneralMatrix(_nDim,_nDim); GeneralMatrix cX = new GeneralMatrix(curX,_nDim); GeneralMatrix pX = new GeneralMatrix(prevX,_nDim); GeneralMatrix cG = new GeneralMatrix(curGrad,_nDim); GeneralMatrix pG = new GeneralMatrix(prevGrad,_nDim); GeneralMatrix dX = cX.Subtract(pX); GeneralMatrix dG = cG.Subtract(pG); double aK1 = 1/(dX.Transpose().Multiply(dG).GetElement(0,0)); GeneralMatrix aK2 = dX.Multiply(dX.Transpose()); GeneralMatrix aK = aK2.Multiply(aK1); double bK1 = -1/(dG.Transpose().Multiply(previousH).Multiply(dG).GetElement(0,0)); GeneralMatrix bK2 = previousH.Multiply(dG).Multiply(dG.Transpose()).Multiply(previousH.Transpose()); GeneralMatrix bK =bK2.Multiply(bK1); currentH = previousH.Add(aK).Add(bK); return currentH; }
public void CholeskyDecomposition2() { 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}}; GeneralMatrix A = new GeneralMatrix(pvals); CholeskyDecomposition chol = A.chol(); GeneralMatrix X = chol.Solve(GeneralMatrix.Identity(3, 3)); Assert.IsTrue(GeneralTests.Check(A.Multiply(X), GeneralMatrix.Identity(3, 3))); }
protected override GeneralMatrix CalculateNextHessianApproximation(GeneralMatrix pH, double[]prevX, double[]curX, double[]prevGrad, double[]curGrad) { GeneralMatrix cH = new GeneralMatrix(_nDim,_nDim); GeneralMatrix cX = new GeneralMatrix(curX,_nDim); GeneralMatrix pX = new GeneralMatrix(prevX,_nDim); GeneralMatrix cG = new GeneralMatrix(curGrad,_nDim); GeneralMatrix pG = new GeneralMatrix(prevGrad,_nDim); GeneralMatrix sigma = cX.Subtract(pX); GeneralMatrix gamma = cG.Subtract(pG); double sigmaTGamma = sigma.Transpose().Multiply(gamma).GetElement(0,0); GeneralMatrix hGammaSigmaT = pH.Multiply(gamma.Multiply(sigma.Transpose())); GeneralMatrix sigmaGammaTH = sigma.Multiply(gamma.Transpose().Multiply(pH)); double gammaTHGamma = (gamma.Transpose().Multiply(pH.Multiply(gamma))).GetElement(0,0); GeneralMatrix sigmaSigmaT = sigma.Multiply(sigma.Transpose()); GeneralMatrix term1 = (hGammaSigmaT.Add(sigmaGammaTH)).Multiply(1/sigmaTGamma); GeneralMatrix term2 = (sigmaSigmaT.Multiply(1/sigmaTGamma)).Multiply(1+gammaTHGamma/sigmaTGamma); return pH.Subtract(term1).Add(term2); }
public void EigenValueDecomposition2() { 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}}; GeneralMatrix A = new GeneralMatrix(evals); EigenvalueDecomposition Eig = A.Eigen(); GeneralMatrix D = Eig.D; GeneralMatrix V = Eig.GetV(); Assert.IsTrue(GeneralTests.Check(A.Multiply(V), V.Multiply(D))); }
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"); }
//encrypt line of input static void Encrypt(string input, GeneralMatrix key, Dictionary<char, int> hillDict, StringBuilder encText) { string ct = ""; string[] pt = input.Trim().Split(' '); for (int i = 0; i < pt.Length; i += 3) { //encrypt 3 letters at a time double[] temp = new double[3]; temp[0] = Convert.ToInt32(pt[i]); temp[1] = Convert.ToInt32(pt[i + 1]); temp[2] = Convert.ToInt32(pt[i + 2]); //create plain text matrix, transpose and encrypt it GeneralMatrix ptMat = new GeneralMatrix(new double[] { temp[0], temp[1], temp[2] }, 3); GeneralMatrix trnasPTMat = ptMat.Transpose(); GeneralMatrix ctMat = key.Multiply(ptMat); for (int x = 0; x < ctMat.RowDimension; x++) { for (int y = 0; y < ctMat.ColumnDimension; y++) { var tempElement = Convert.ToInt32(ctMat.GetElement(x, y)) % 31; ctMat.SetElement(x, y, tempElement); } } for (int x = 0; x < ctMat.RowDimension; x++) { for (int y = 0; y < ctMat.ColumnDimension; y++) { ct += hillDict.FirstOrDefault(z => z.Value == ctMat.GetElement(x, y)).Key; } } } //append to string builder encText.AppendLine(ct); }
public void Rotate(double sineAngle, double cosineAngle) { double[][] temp = new double[2][]; temp[0] = new double[2]; temp[1] = new double[2]; temp[0][0] = cosineAngle; temp[0][1] = -sineAngle; temp[1][0] = sineAngle; temp[1][1] = cosineAngle; GeneralMatrix rotationMatrix = new GeneralMatrix(temp); GeneralMatrix newVector = rotationMatrix.Multiply(this.ToMatrix()); this.X = newVector.GetElement(0, 0); this.Y = newVector.GetElement(1, 0); }
public void Rotate(double angleInRadians) { /* http://mathworld.wolfram.com/RotationMatrix.html */ /* матрица на ротација * |- -| * | cos(a) -sin(a) | * R(a) = | | * | sin(a) cos(a) | * |- -| */ //V = R(a) * V0 double[][] temp = new double[2][]; temp[0] = new double[2]; temp[1] = new double[2]; temp[0][0] = Math.Cos(angleInRadians); temp[0][1] = -Math.Sin(angleInRadians); temp[1][0] = Math.Sin(angleInRadians); temp[1][1] = Math.Cos(angleInRadians); GeneralMatrix rotationMatrix = new GeneralMatrix(temp); GeneralMatrix newVector = rotationMatrix.Multiply(this.ToMatrix()); this.X = newVector.GetElement(0, 0); this.Y = newVector.GetElement(1, 0); }
/// <summary> /// multiply normalized priority matrix by sum of average rows /// </summary> /// <param name="argMatrix"></param> /// <param name="selection"></param> /// <returns></returns> private GeneralMatrix FCalc(GeneralMatrix argMatrix, GeneralMatrix selection) { GeneralMatrix matrix = argMatrix.Multiply(selection); return (matrix.ArrayRightDivide(selection)); }
private void GetFiles_Decompose_WriteToArff() { try { // sampleFactor is the amount to divide by the total size of the // data set // when determining the subsample that will be used in svd Console.WriteLine("Creating arff file..."); DisplayMessage("Creating arff file..."); // Create folder System.IO.Directory.CreateDirectory(OutputDir + "Results/"); StreamWriter output = new StreamWriter(OutputDir + "Results/" + "resultsGalaxy.arff"); output.Write("@relation 'galaxy'\n"); output.Write("@attribute class real\n"); output.Write("@attribute colorF real\n"); output.Write("@attribute bulgeF real\n"); output.Write("@attribute constF real\n"); for (int i = 0; i < sv * 3; i++) { output.Write("@attribute " + i + " real\n"); } output.Write("@data\n"); Console.WriteLine("Begin galaxy sampling"); DisplayMessage("Begin galaxy sampling"); // Initialize a three matrices that will hold all of the images // (r,g,b of each image where each row is an image) dataRed = new GeneralMatrix(galaxyData.Count() / sampleFactor, imageScaleSize * imageScaleSize); dataGreen = new GeneralMatrix(galaxyData.Count() / sampleFactor, imageScaleSize * imageScaleSize); dataBlue = new GeneralMatrix(galaxyData.Count() / sampleFactor, imageScaleSize * imageScaleSize); // subsample from galaxydata System.Threading.Tasks.Parallel.For(0, galaxyData.Count() / sampleFactor, (int index) => { Bitmap tempImage = getImage(OutputDir + "galaxies/" + galaxyData[sampleFactor * index][0] + ".jpg", imageScaleSize); for (int i = 0; i < imageScaleSize; i++) { for (int j = 0; j < imageScaleSize; j++) { int pixelColor = tempImage.GetPixel(i, j).ToArgb(); int[] rgb = new int[3]; rgb[0] += ((pixelColor & 0x00ff0000) >> 16); rgb[1] += ((pixelColor & 0x0000ff00) >> 8); rgb[2] += (pixelColor & 0x000000ff); dataRed.SetElement(index, i * imageScaleSize + j, rgb[0]); dataGreen.SetElement(index, i * imageScaleSize + j, rgb[1]); dataBlue.SetElement(index, i * imageScaleSize + j, rgb[2]); } } //if (index % 10 == 0) //DisplayImage(index); //Console.WriteLine("Galaxy " + (sampleFactor * index) + " finished"); }); Console.WriteLine("Galaxy sampling finished\nBegin R, G and B channel SVD"); DisplayMessage("Galaxy sampling finished, Begin R, G and B channel SVD"); // Perform svd on subsample: var redWorker = System.Threading.Tasks.Task.Factory.StartNew(() => svdR = new SingularValueDecomposition(dataRed)); var greenWorker = System.Threading.Tasks.Task.Factory.StartNew(() => svdG = new SingularValueDecomposition(dataGreen)); var blueWorker = System.Threading.Tasks.Task.Factory.StartNew(() => svdB = new SingularValueDecomposition(dataBlue)); System.Threading.Tasks.Task.WaitAll(redWorker, greenWorker, blueWorker); // Create the basis for each component GeneralMatrix rV = svdR.GetV(); Console.Write("dim rU: " + rV.RowDimension + ", " + rV.ColumnDimension); GeneralMatrix gV = svdG.GetV(); GeneralMatrix bV = svdB.GetV(); rV = GetSVs(rV, sv); Console.Write("Svs: " + sv); Console.Write("Dim SsV: " + rV.RowDimension + ", " + rV.ColumnDimension); gV = GetSVs(gV, sv); bV = GetSVs(bV, sv); // Perform the pseudoinverses frV = pinv(rV.Transpose()); fgV = pinv(gV.Transpose()); fbV = pinv(bV.Transpose()); // Stores frV, fgV and fbV to file WriteToFile(); Console.WriteLine("SVD finished"); DisplayMessage("SVD finished, load full dataset for testing"); Console.WriteLine("Begin filling dataset for testing"); // fill the 'full' datasets dataRed = new GeneralMatrix(galaxyData.Count(), imageScaleSize * imageScaleSize); dataGreen = new GeneralMatrix(galaxyData.Count(), imageScaleSize * imageScaleSize); dataBlue = new GeneralMatrix(galaxyData.Count(), imageScaleSize * imageScaleSize); System.Threading.Tasks.Parallel.For(0, galaxyData.Count(), (int index) => { Bitmap tempImage = getImage(OutputDir + "galaxies/" + galaxyData[index][0] + ".jpg", imageScaleSize); for (int i = 0; i < imageScaleSize; i++) { for (int j = 0; j < imageScaleSize; j++) { int pixelColor = tempImage.GetPixel(i, j).ToArgb(); int[] rgb = new int[3]; rgb[0] += ((pixelColor & 0x00ff0000) >> 16); rgb[1] += ((pixelColor & 0x0000ff00) >> 8); rgb[2] += (pixelColor & 0x000000ff); dataRed.SetElement(index, i * imageScaleSize + j, rgb[0]); dataGreen.SetElement(index, i * imageScaleSize + j, rgb[1]); dataBlue.SetElement(index, i * imageScaleSize + j, rgb[2]); } } }); Console.WriteLine("Finished filling dataset for testing"); DisplayMessage("Finished filling dataset for testing, begin projecting galaxies to U coordinate system"); Console.WriteLine("Begin projecting galaxies to U coordinate system, writing to ARFF file"); // Do the coordinate conversion rV = dataRed.Multiply(frV); gV = dataGreen.Multiply(fgV); bV = dataBlue.Multiply(fbV); Console.Write("Dim Final rU: " + rV.ColumnDimension + ", " + rV.RowDimension); Console.WriteLine("galaxyData.Count(): " + galaxyData.Count()); // write to the output file here: for (int imageIndex = 0; imageIndex < galaxyData.Count(); imageIndex++) { Bitmap tempImage = getImage(OutputDir + "galaxies/" + galaxyData[imageIndex][0] + ".jpg", imageScaleSize); float colorFactor = (GetColor(tempImage)[0] / GetColor(tempImage)[2]); float centralBulgeFactor = getCentralBulge(tempImage); float consistencyFactor = GetConsistency(tempImage); output.Write(galaxyData[imageIndex][1] + ", "); output.Write(colorFactor + ", "); output.Write(centralBulgeFactor + ", "); output.Write(consistencyFactor + ", "); // output data (r,g,b) for (int i = 0; i < rV.ColumnDimension; i++) { output.Write(rV.GetElement(imageIndex, i) + ", "); output.Write(gV.GetElement(imageIndex, i) + ", "); if (i == rV.ColumnDimension - 1) { output.Write(bV.GetElement(imageIndex, i) + "\n"); } else { output.Write(bV.GetElement(imageIndex, i) + ", "); } } //if (imageIndex % (galaxyData.Count() / 100) == 0) // DisplayImage(imageIndex); DisplayMessage("Finished galaxy " + imageIndex.ToString() + " - " + (100 * imageIndex / galaxyData.Count()).ToString() + "%"); } output.Flush(); output.Close(); output.Dispose(); Console.Write("Finished creating arff file..."); DisplayMessage("Finished creating arff file..."); } catch (Exception ex) { Console.Write(ex.ToString()); } }
private void ClassifyGroup(string[] images, string filePath, StreamWriter output, int groupNum) { // fill the 'full' datasets dataRed = new GeneralMatrix(images.Count(), imageScaleSize * imageScaleSize); dataGreen = new GeneralMatrix(images.Count(), imageScaleSize * imageScaleSize); dataBlue = new GeneralMatrix(images.Count(), imageScaleSize * imageScaleSize); //create the entire size of the dataset System.Threading.Tasks.Parallel.For(0, images.Count(), (int imageIndex) => { Bitmap tempImage = getImage(images[imageIndex], imageScaleSize); for (int i = 0; i < imageScaleSize; i++) { for (int j = 0; j < imageScaleSize; j++) { int pixelColor = tempImage.GetPixel(i, j).ToArgb(); int[] rgb = new int[3]; rgb[0] += ((pixelColor & 0x00ff0000) >> 16); rgb[1] += ((pixelColor & 0x0000ff00) >> 8); rgb[2] += (pixelColor & 0x000000ff); dataRed.SetElement(imageIndex, i * imageScaleSize + j, rgb[0]); dataGreen.SetElement(imageIndex, i * imageScaleSize + j, rgb[1]); dataBlue.SetElement(imageIndex, i * imageScaleSize + j, rgb[2]); } } Console.WriteLine("SDSS Galaxy " + imageIndex + " put in main dataset"); }); //then convert the whole dataset to the same coordinate // Do the coordinate conversion GeneralMatrix rV = dataRed.Multiply(frV); GeneralMatrix gV = dataGreen.Multiply(fgV); GeneralMatrix bV = dataBlue.Multiply(fbV); Console.WriteLine("Dim Final rU: " + rV.ColumnDimension + ", " + rV.RowDimension); Console.WriteLine("images.length: " + images.Count()); // write to the output file here: for (int imageIndex = 0; imageIndex < images.Count(); imageIndex++) { Bitmap tempImage = getImage(images[imageIndex], imageScaleSize); float colorFactor = (GetColor(tempImage)[0] / GetColor(tempImage)[2]); float centralBulgeFactor = getCentralBulge(tempImage); float consistencyFactor = GetConsistency(tempImage); output.Write(-999 + ", "); output.Write(colorFactor + ", "); output.Write(centralBulgeFactor + ", "); output.Write(consistencyFactor + ", "); // output data (r,g,b) for (int i = 0; i < rV.ColumnDimension; i++) { output.Write(rV.GetElement(imageIndex, i) + ", "); output.Write(gV.GetElement(imageIndex, i) + ", "); if (i == rV.ColumnDimension - 1) { output.Write(bV.GetElement(imageIndex, i) + "\n"); } else { output.Write(bV.GetElement(imageIndex, i) + ", "); } } Console.WriteLine("Creating ARFF classification file - " + (100 * imageIndex / images.Count()).ToString() + "%"); } }
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; }
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; }