public void btnTestEigenvalues_Click(IRibbonControl e) { Random rando = new Random(); double[,] TwoArray = new double[500, 500]; for (int i = 0; i < TwoArray.GetLength(0); i++) { for (int j = 0; j < TwoArray.GetLength(1); j++) { TwoArray[i, j] = rando.NextDouble(); } } List <double> ticks1 = new List <double>(); List <double> ticks2 = new List <double>(); for (int i = 0; i < 10; i++) { DiagnosticsMenu.StartStopwatch(); double[] result1 = MatrixOps.RealEigenvalues_Accord(TwoArray); ticks1.Add(DiagnosticsMenu.StopStopwatch(TimeUnit.milliseconds)); DiagnosticsMenu.StartStopwatch(); result1 = null; } MessageBox.Show($"{ticks1.Average().ToString()} Accord"); }
public void btnTestTranspose_Click(IRibbonControl e) { Random rando = new Random(); double[,] TwoArray = new double[1000, 1000]; for (int i = 0; i < TwoArray.GetLength(0); i++) { for (int j = 0; j < TwoArray.GetLength(1); j++) { TwoArray[i, j] = rando.NextDouble(); } } List <double> ticks1 = new List <double>(); List <double> ticks2 = new List <double>(); //int equalTo = 16; for (int i = 0; i < 10; i++) { DiagnosticsMenu.StartStopwatch(); var result1 = MatrixOps.Transpose(TwoArray); ticks1.Add(DiagnosticsMenu.StopStopwatch(TimeUnit.milliseconds)); DiagnosticsMenu.StartStopwatch(); var result2 = ThisAddIn.MyApp.WorksheetFunction.Transpose(TwoArray); ticks2.Add(DiagnosticsMenu.StopStopwatch(TimeUnit.milliseconds)); //int places = TestEquality(result1, result2); //if (places < equalTo) // equalTo = places; result1 = null; result2 = null; } MessageBox.Show($"{ticks1.Average().ToString()} Accord"); MessageBox.Show($"{ticks2.Average().ToString()} VBA"); }
public void btnTestInverse_Click(IRibbonControl e) { Random rando = new Random(); double[,] TwoArray = new double[100, 100]; double[,] result1 = new double[100, 100]; double[,] result2 = new double[100, 100]; List <double> ticks1 = new List <double>(); List <double> ticks2 = new List <double>(); //int equalTo = 16; for (int i = 0; i < 50; i++) { for (int i2 = 0; i2 < TwoArray.GetLength(0); i2++) { for (int j = 0; j < TwoArray.GetLength(1); j++) { TwoArray[i2, j] = rando.NextDouble(); } } DiagnosticsMenu.StartStopwatch(); result1 = MatrixOps.MatrixInverse(TwoArray); ticks1.Add(DiagnosticsMenu.StopStopwatch(TimeUnit.milliseconds)); DiagnosticsMenu.StartStopwatch(); result2 = Utilities.Conversions.ConvertObjectArrayToDouble(ThisAddIn.MyApp.WorksheetFunction.MInverse(TwoArray)); //result2 = ThisAddIn.MyApp.WorksheetFunction.MInverse(TwoArray); ticks2.Add(DiagnosticsMenu.StopStopwatch(TimeUnit.milliseconds)); } ThisAddIn.MyApp.ActiveSheet.Range("A1").Value = "Original Matrix"; ThisAddIn.MyApp.ActiveSheet.Range("A2:CV101").Value = TwoArray; ThisAddIn.MyApp.ActiveSheet.Range("A102").Value = "Accord Inverse"; ThisAddIn.MyApp.ActiveSheet.Range("A103:CV202").Value = result1; ThisAddIn.MyApp.ActiveSheet.Range("A203").Value = "VBA Matrix"; ThisAddIn.MyApp.ActiveSheet.Range("A204:CV303").Value = result2; MessageBox.Show($"ACC v VBA Precision: {TestEquality(result1, result2, 9)}"); MessageBox.Show($"{ticks1.Average().ToString()} Accord"); MessageBox.Show($"{ticks2.Average().ToString()} VBA"); }
public void btnTestDeterminant_Click(IRibbonControl e) { Random rando = new Random(); double[,] TwoArray = new double[5, 5]; List <double> ticks1 = new List <double>(); List <double> ticks2 = new List <double>(); int equalTo = 16; for (int i = 0; i < 10; i++) { for (int i2 = 0; i2 < TwoArray.GetLength(0); i2++) { for (int j = 0; j < TwoArray.GetLength(1); j++) { TwoArray[i2, j] = rando.NextDouble(); } } DiagnosticsMenu.StartStopwatch(); double result1 = MatrixOps.Determinant(TwoArray); ticks1.Add(DiagnosticsMenu.StopStopwatch(TimeUnit.ticks)); DiagnosticsMenu.StartStopwatch(); double result2 = ThisAddIn.MyApp.WorksheetFunction.MDeterm(TwoArray); ticks2.Add(DiagnosticsMenu.StopStopwatch(TimeUnit.ticks)); int places = TestEquality(result1, result2); if (places < equalTo) { equalTo = places; } result1 = 0; result2 = 0; } MessageBox.Show($"{ticks1.Average().ToString()} Accord"); MessageBox.Show($"{ticks2.Average().ToString()} VBA"); MessageBox.Show($"Equal to {equalTo} decimal places."); }
public void btnTestMMult_Click(IRibbonControl e) { Random rando = new Random(); double[,] m1 = new double[5, 5]; double[,] m2 = new double[5, 5]; List <long> ticks1 = new List <long>(); List <long> ticks2 = new List <long>(); List <long> ticks3 = new List <long>(); for (int i = 0; i < 100; i++) //measure timing { for (int i2 = 0; i2 < m2.GetLength(0); i2++) { for (int j = 0; j < m2.GetLength(1); j++) { m1[i2, j] = rando.NextDouble(); m2[i2, j] = rando.NextDouble(); } } DiagnosticsMenu.StartStopwatch(); double[,] result1 = MatrixOps.MMult(m1, m2); DiagnosticsMenu.StopStopwatch(TimeUnit.ticks); ticks1.Add(DiagnosticsMenu.stopwatch.ElapsedTicks); DiagnosticsMenu.StartStopwatch(); double[,] result2 = MatrixOps.Accord_MMult(m1, m2); DiagnosticsMenu.StopStopwatch(TimeUnit.ticks); ticks2.Add(DiagnosticsMenu.stopwatch.ElapsedTicks); DiagnosticsMenu.StartStopwatch(); dynamic result3 = ThisAddIn.MyApp.WorksheetFunction.MMult(m1, ThisAddIn.MyApp.WorksheetFunction.Transpose(m2)); DiagnosticsMenu.StopStopwatch(TimeUnit.ticks); ticks3.Add(DiagnosticsMenu.stopwatch.ElapsedTicks); //int places = TestEquality(result1, result2); //if (places < equalTo) // equalTo = places; result1 = null; result2 = null; result3 = null; } MessageBox.Show($"{ticks1.Average().ToString()} Manual"); MessageBox.Show($"{ticks2.Average().ToString()} Accord"); MessageBox.Show($"{ticks3.Average().ToString()} VBA"); double[,] vba_out = Utilities.Conversions.ConvertObjectArrayToDouble(ThisAddIn.MyApp.WorksheetFunction.MMult(m1, m2)); double[,] man_out = MatrixOps.MMult(m1, m2); double[,] acc_out = MatrixOps.Accord_MMult(m1, m2); ThisAddIn.MyApp.ActiveSheet.Range("A1").Value = "Manual"; ThisAddIn.MyApp.ActiveSheet.Range("A2:E6").Value = man_out; ThisAddIn.MyApp.ActiveSheet.Range("A7").Value = "Accord"; ThisAddIn.MyApp.ActiveSheet.Range("A8:E12").Value = acc_out; ThisAddIn.MyApp.ActiveSheet.Range("A13").Value = "VBA Function"; ThisAddIn.MyApp.ActiveSheet.Range("A14:E18").Value = vba_out; ThisAddIn.MyApp.ActiveSheet.Range("G1").Value = "M1"; ThisAddIn.MyApp.ActiveSheet.Range("G2:K6").Value = m1; ThisAddIn.MyApp.ActiveSheet.Range("G7").Value = "M2"; ThisAddIn.MyApp.ActiveSheet.Range("G8:K12").Value = m2; MessageBox.Show($"MAN v ACC Precision: {TestEquality(man_out, acc_out)}"); MessageBox.Show($"MAN v VBA Precision: {TestEquality(man_out, vba_out)}"); MessageBox.Show($"ACC v VBA Precision: {TestEquality(acc_out, vba_out)}"); }