Question

我正在尝试转换一个函数，使其在C＃中并行执行而不是顺序执行，但我不确定我做错了什么：

// sequential
static void Romberg(double a, double b, int n, double[,] R)
{
    int i, j, k;
    double h, sum;

    h = b - a;
    R[0, 0] = (h / 2) * (f(a) + f(b));

    for (i = 1; i <= n; i++)
    {
        h = h / 2;
        sum = 0;

        for (k = 1; k <= (Math.Pow(2.0, i) - 1); k += 2)
        {
            sum += f(a + k * h);
        }

        R[i, 0] = R[i - 1, 0] / 2 + sum * h;

        for (j = 1; j <= i; j++)
        {
            R[i, j] = R[i, j - 1] + (R[i, j - 1] - R[i - 1, j - 1]) / (Math.Pow(4.0, j) - 1);
        }
    }
}


// parallel
static void RombergCP(double a, double b, int n, double[,] R)
{
    int i,j, k;
    double h, sum;

    h = b - a;
    R[0, 0] = (h / 2) * (f(a) + f(b));

    Parallel.For(0, n, options, i =>
    {
        h = h / 2;
        sum = 0;

         for (k = 1; k <= (Math.Pow(2.0, i) - 1); k += 2)
        {
            sum += f(a + k * h);
        };

        R[i, 0] = R[i - 1, 0] / 2 + sum * h;

        for (j = 1; j <= i; j++)
        {
            R[i, j] = R[i, j - 1] + (R[i, j - 1] - R[i - 1, j - 1]) / (Math.Pow(4.0, j) - 1);
        }
    });
}

我得到的错误是“i”无法声明，因为它会给“i”赋予不同的含义，“i”用于“父或当前”范围。我尝试在并行函数中将其重命名为i2但它给出了相同的错误。提前谢谢！

Answer 1

在最顶层删除int i的声明。它由下面的lambda声明。

Answer 2

有几个问题：

在尽可能小的范围内声明变量。
您的外部循环从for (i = 1; i <= n; i++)转到Parallel.For(0, n, options, ...)，这意味着R[i-1, ...]会抛出并行版本。
h = h / 2;不是线程安全的。

// parallel
static void RombergCP(double a, double b, int n, double[,] R)
{
  //int i,j, k;
  //double h, sum;

  double h0 = b - a;
  R[0, 0] = (h0 / 2) * (f(a) + f(b));

  Parallel.For(1, n, options, i =>   // start at 1
  {
     //h = h / 2;
     double h = (b - a) / Math.Pow(2, i);    // derive from i
     double sum = 0;

     for (int k = 1; k <= (Math.Pow(2.0, i) - 1); k += 2)   // keep k local
       ...

C＃将顺序函数转换为并行执行

2 个答案: