Exemple #1
0
        public static void MainCycle(object obj)
        {
            AllParameters p = (AllParameters)obj;

            for (int i = p.St; i < p.End; i++)
            {
                for (int j = 0; j < p.End2; j++)
                {
                    p.res[i][j] = -p.mulH[i] - p.mulV[j];
                    for (int k = 0; k < p.End3 / 2; k++)
                    {
                        p.res[i][j] = p.res[i][j] + (p.matr1[i][2 * k + 1] + p.matr2[2 * k][j]) * (p.matr1[i][2 * k] + p.matr2[2 * k + 1][j]);
                    }
                }
            }
        }
Exemple #2
0
        // This addresation works faster
        public static void MainCycleOptimize(object obj)
        {
            AllParameters p = (AllParameters)obj;

            int[][] res = p.res, matr1 = p.matr1, matr2 = p.matr2;
            int[]   mulH = p.mulH, mulV = p.mulV;
            int     st = p.St, end = p.End;
            int     end2 = p.End2, end3 = p.End3;

            for (int i = st; i < end; i++)
            {
                for (int j = 0; j < end2; j++)
                {
                    res[i][j] = -mulH[i] - mulV[j];
                    for (int k = 0; k < end3 / 2; k++)
                    {
                        res[i][j] = res[i][j] + (matr1[i][2 * k + 1] + matr2[2 * k][j]) * (matr1[i][2 * k] + matr2[2 * k + 1][j]);
                    }
                }
            }
        }
Exemple #3
0
        public static int[][] ParallelMultVin(int[][] matr1, int[][] matr2, int nThreads)
        {
            int n1 = matr1.Length;
            int n2 = matr2.Length;

            if (n1 == 0 || n2 == 0)
            {
                return(null);
            }

            int m1 = matr1[0].Length;
            int m2 = matr2[0].Length;

            if (m1 != n2)
            {
                return(null);
            }

            int[] mulH = new int[n1];
            int[] mulV = new int[m2];

            int[][] res = new int[n1][];
            for (int i = 0; i < n1; i++)
            {
                res[i] = new int[m2];
            }

            Thread[] t = new Thread[nThreads];

            /*
             * // Count MulV
             * int colForThread = m2 / nThreads;
             * int m2St = 0;
             * for (int i = 0; i < nThreads; i++)
             * {
             *  int m2End = m2St + colForThread;
             *  if (i == nThreads - 1)
             *      m2End = m2;
             *  Parameters p = new Parameters(matr2, mulV, m2St, m2End, n2);
             *
             *  t[i] = new Thread(CountMulV);
             *  t[i].Start(p);
             *
             *  m2St = m2End;
             * }
             *
             * foreach (Thread thread in t)
             * {
             *  thread.Join();
             * }
             *
             * // Count MulH
             * int rowsForThread = n1 / nThreads;
             * int n1St = 0;
             * for (int i = 0; i < nThreads; i++)
             * {
             *  int n1End = n1St + rowsForThread;
             *  if (i == nThreads - 1)
             *      n1End = n1;
             *  Parameters p = new Parameters(matr1, mulH, n1St, n1End, m1);
             *
             *  t[i] = new Thread(CountMulH);
             *  t[i].Start(p);
             *
             *  n1St = n1End;
             * }
             * foreach (Thread thread in t)
             * {
             *  thread.Join();
             * }
             */

            for (int i = 0; i < n1; i++)
            {
                for (int j = 0; j < m1 / 2; j++)
                {
                    mulH[i] = mulH[i] + matr1[i][j * 2] * matr1[i][j * 2 + 1];
                }
            }

            for (int i = 0; i < m2; i++)
            {
                for (int j = 0; j < n2 / 2; j++)
                {
                    mulV[i] = mulV[i] + matr2[j * 2][i] * matr2[j * 2 + 1][i];
                }
            }


            int rowsForThread = n1 / nThreads;
            int n1St          = 0;

            for (int i = 0; i < nThreads; i++)
            {
                int n1End = n1St + rowsForThread;
                if (i == nThreads - 1)
                {
                    n1End = n1;
                }
                AllParameters p = new AllParameters(res, matr1, matr2, mulV, mulH, n1St, n1End, m2, m1);

                t[i] = new Thread(MainCycleOptimize);
                t[i].Start(p);

                n1St = n1End;
            }
            foreach (Thread thread in t)
            {
                thread.Join();
            }


            if (m1 % 2 == 1)
            {
                for (int i = 0; i < n1; i++)
                {
                    for (int j = 0; j < m2; j++)
                    {
                        res[i][j] = res[i][j] + matr1[i][m1 - 1] * matr2[m1 - 1][j];
                    }
                }
            }

            return(res);
        }