并行程序的求和算法

我正在尝试编写一个并行算法，使其比执行基本相同操作的顺序算法快三倍。请看糊状物。

http://pastebin.com/3DDyxfPP

粘贴：

大家好。我正在做一个课堂作业，并且大部分都完成了，但是我在数学方面遇到了一些问题。我正在尝试计算表达式：

100000000
∑ (9999999/10000000)^i * i^2
i = 1

我从1000万到1000万。给出了一个快速顺序算法：

  double sum = 0.0;
  double fact1 = 0.9999999;
  for (int i = 1; i <= 10000000; i++)
  {
     sum += (fact1 * i * i);
     fact1 *= 0.9999999;
  }

我们应该实现它并验证它是否有效，并在发布模式下计时。我已经完成了这项工作并正常工作。然后，时间将显示在控制台上。

 DateTime t = DateTime.Now;
 long saveticks = t.Ticks;
 double sum = 0.0;
 double fact1 = 0.9999999;
 for (int i = 1; i <= 100000000; i++)
 {
    sum += (fact1 * i * i);
    fact1 *= 0.9999999;
 }
 t = DateTime.Now;

然后，我们必须编写一个定时并行算法来击败时间，并且应该在示例并行程序之后对其进行建模。它必须至少比顺序算法快 3 倍。我们将为并行程序使用 4 个处理元素。

有一个提示，"在你弄清楚每个处理元素将要做的工作后，你可能需要从耗时的 Pow 函数开始处理元素"。

例如：Math.Pow(x，y(

"不要在每次迭代时对并行代码使用 pow 函数，因为它不会打败时间。">

这是我的并行程序代码。这既执行顺序算法，也执行并行算法，并将它们都计时。

 const int numPEs = 4;
 const int size = 100000000;
 static double pSum;
 static int numThreadsDone;
 static int nextid;
 static object locker1 = new object();
 static object locker2 = new object();
 static long psaveticks;
 static DateTime pt;
 static void Main(string[] args)
 {
     DateTime t = DateTime.Now;
     long saveticks = t.Ticks;
     double sum = 0.0;
     double fact1 = 0.9999999;
     for (int i = 1; i <= 100000000; i++)
     {
         sum += (fact1 * (i * i));
         fact1 *= 0.9999999;
     }
     t = DateTime.Now;
     Console.WriteLine("sequential: " + ((t.Ticks - saveticks) / 100000000.0) + " seconds");
     Console.WriteLine("sum is " + sum);
     // time it
     pt = DateTime.Now;
     psaveticks = pt.Ticks;
     for (int i = 0; i < numPEs; i++)
     new Thread(countThreads).Start();
     Console.ReadKey();
 }

 static void countThreads()
 {
     int id;
     double localcount = 0;
     lock (locker1)
     {
         id = nextid;
         nextid++;
     }
     // assumes array is evenly divisible by the number of threads
     int granularity = size / numPEs;
     int start = granularity * id;
     for (int i = start; i < start + granularity; i++)
         localcount += (Math.Pow(0.9999999, i) * (i * i));
     lock (locker2)
     {
         pSum += localcount;
         numThreadsDone++;
         if (numThreadsDone == numPEs)
         {
             pt = DateTime.Now;
             Console.WriteLine("parallel: " + ((pt.Ticks - psaveticks) / 10000000.0) + " seconds");
             Console.WriteLine("parallel count is " + pSum);
         }
     }
 }

我的

问题是我的顺序程序比并行程序快得多。我使用的算法一定有问题。

谁能帮忙？

Console.WriteLine("sequential: " + ((t.Ticks - saveticks) / 100000000.0) + " seconds");

// Inner loop of sequential code:
sum += (fact1 * (i * i));
fact1 *= 0.9999999;
// Inner loop of parallel code:
localcount += (Math.Pow(0.9999999, i) * (i * i));

double fact1 = Math.Pow(0.9999999, start + 1);
for (int i = start + 1; i <= start + granularity; i++)
{
    localcount += (fact1 * (i * i));
    fact1 *= 0.9999999;
}