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 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");
        }
예제 #3
0
        public static double[,] Accord_MMult(double[,] m1, double[,] m2)
        {
            DiagnosticsMenu.StartStopwatch();
            var retVal = Matrix.Dot(m1, m2);

            DiagnosticsMenu.StopStopwatch();
            return(retVal);
        }
        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)}");
        }
 private void ThisAddIn_Startup(object sender, System.EventArgs e)
 {
     DiagnosticsMenu.PrimeDiagnostics();          //initializes any diagnostic objects early so that they don't affect the results
     MyApp       = Globals.ThisAddIn.Application; //Grab Excel at startup.
     TempAppList = new List <Excel.Application>();
 }