Как максимально ускорить программу в MVS?

for (int w1 = 0; w1 <= 51; w1++)
            {
                for (int w2 = w1 + 1; w2 <= 51; w2++)
                {
                    for (int w8 = w1; w8 <= 51; w8++)
                    {
                        for (int w9 = w8 + 1; w9 <= 51; w9++)
                        {
                            for (int w6 = 0; w6 <= 51; w6++)
                            {
                                for (int w7 = w6 + 1; w7 <= 51; w7++)
                                {
                                    if (Hstr1[w1,w2,w6,w7] == 0 | Hstr1[w8,w9,w6,w7] == 0)
                                    {
                                        continue;
                                    }
                                    if (Hstr1[w1, w2, w6, w7] > Hstr1[w8, w9, w6, w7])
                                    {
                                        mC[w1, w2, w8, w9]++;
                                    }
                                    else if (Hstr1[w1, w2, w6, w7] < Hstr1[w8, w9, w6, w7])
                                    {
                                        mC[w8, w9, w1, w2]++;
                                    }
                                    else
                                    {
                                        mC[w1, w2, w8, w9]+=0.5;
                                        mC[w8, w9, w1, w2]+=0.5;
                                    }
                                }
                            }
                        }
                    }
                }
            }

for (int w1 = 0; w1 <= 51; w1++)
{
    for (int w2 = w1 + 1; w2 <= 51; w2++)
    {
        for (int w8 = w1; w8 <= 51; w8++)
        {
            for (int w9 = w8 + 1; w9 <= 51; w9++)
            {
                for (int w6 = 0; w6 <= 51; w6++)
                {
                    for (int w7 = w6 + 1; w7 <= 51; w7++)
                    {
                        var x = Hstr1[w1, w2, w6, w7];
                        var y = Hstr1[w8, w9, w6, w7];
                        if (x == 0 || y == 0)
                        {
                            continue;
                        }
                        if (x > y)
                        {
                            mC[w1, w2, w8, w9]++;
                        }
                        else if (x < y)
                        {
                            mC[w8, w9, w1, w2]++;
                        }
                        else
                        {
                            mC[w1, w2, w8, w9]+=0.5;
                            mC[w8, w9, w1, w2]+=0.5;
                        }
                    }
                }
            }
        }
    }
}

for (int w1 = 0; w1 <= 51; w1++)
{
    for (int w2 = w1 + 1; w2 <= 51; w2++)
    {
        for (int w8 = w1; w8 <= 51; w8++)
        {
            for (int w9 = w8 + 1; w9 <= 51; w9++)
            {
                ref double mC1 = ref mC[w1, w2, w8, w9];
                ref double mC2 = ref mC[w8, w9, w1, w2];
                for (int w6 = 0; w6 <= 51; w6++)
                {
                    for (int w7 = w6 + 1; w7 <= 51; w7++)
                    {
                        var x = Hstr1[w1, w2, w6, w7];
                        var y = Hstr1[w8, w9, w6, w7];
                        if (x == 0 || y == 0)
                        {
                            continue;
                        }
                        if (x > y)
                        {
                            mC[w1, w2, w8, w9]++;
                        }
                        else if (x < y)
                        {
                            mC[w8, w9, w1, w2]++;
                        }
                        else
                        {
                            mC1+=0.5;
                            mC2+=0.5;
                        }
                    }
                }
            }
        }
    }
}

            Parallel.For(fromInclusive: 0, toExclusive: 52, body: (w1) =>
            {
                for (int w2 = w1 + 1; w2 <= 51; w2++)
                {
                    for (int w8 = w1; w8 <= 51; w8++)
                    {
                        for (int w9 = w8 + 1; w9 <= 51; w9++)
                        {
                            for (int w6 = 0; w6 <= 51; w6++)
                            {
                                for (int w7 = w6 + 1; w7 <= 51; w7++)
                                {
                                    if (Hstr1[w1,w2,w6,w7] == 0 | Hstr1[w8,w9,w6,w7] == 0)
                                    {
                                        continue;
                                    }
                                    if (Hstr1[w1, w2, w6, w7] > Hstr1[w8, w9, w6, w7])
                                    {
                                        //lock (locker)
                                        //{
                                        //    mC[w1, w2, w8, w9]++;
                                        //}
 
                                        InterlockedEx.Add(ref mC[w1, w2, w8, w9], 1);
                                    }
                                    else if (Hstr1[w1, w2, w6, w7] < Hstr1[w8, w9, w6, w7])
                                    {
                                        //lock (locker)
                                        //{
                                        //    mC[w8, w9, w1, w2]++;
                                        //}
 
                                        InterlockedEx.Add(ref mC[w8, w9, w1, w2], 1);
                                    }
                                    else
                                    {
                                        //lock (locker)
                                        //{
                                        //    mC[w1, w2, w8, w9] += 0.5;
                                        //    mC[w8, w9, w1, w2] += 0.5;
                                        //}
 
                                        InterlockedEx.Add(ref mC[w1, w2, w8, w9], 0.5);
                                        InterlockedEx.Add(ref mC[w8, w9, w1, w2], 0.5);
                                    }
                                }
                            }
                        }
                    }
                }
            });
 
 
 
        public static class InterlockedEx
        {
            //[MethodImpl(MethodImplOptions.AggressiveInlining)]
            public static double Add(ref double location1, double value)
            {
                Thread.MemoryBarrier();
                double snapshot = location1;
                double snapshot2;
                double result;
                do
                {
                    result = snapshot + value;
                    snapshot2 = Interlocked.CompareExchange(ref location1, result, snapshot);
                    if (snapshot2 == snapshot)
                    {
                        return result;
                    }
                    else
                    {
                        snapshot = snapshot2;
                        //Thread.SpinWait(300);
                    }
                } while (true);
            }
        }