1

我试图理解用 C 编写的一段代码，但不确定我是否完全理解它。

这是用C编写的函数：

int
gsl_multimin_diff (const gsl_multimin_function * f,
                   const gsl_vector * x, gsl_vector * g)
{
  size_t i, n = f->n;

  double h = GSL_SQRT_DBL_EPSILON;


  gsl_vector * x1 = gsl_vector_alloc (n);  /* FIXME: pass as argument */

  gsl_vector_memcpy (x1, x);

  for (i = 0; i < n; i++)
    {
      double fl, fh;

      double xi = gsl_vector_get (x, i);
      double dx = fabs(xi) * h;

      if (dx == 0.0) dx = h;    

     (x1, i, xi + dx);
      fh = GSL_MULTIMIN_FN_EVAL(f, x1);

      gsl_vector_set (x1, i, xi - dx);
      fl = GSL_MULTIMIN_FN_EVAL(f, x1);

      gsl_vector_set (x1, i, xi);
      gsl_vector_set (g, i, (fh - fl) / (2.0 * dx));
    }

  gsl_vector_free (x1);

  return GSL_SUCCESS;
}

这段代码中有第 14 行，其中有：(x1, i, xi + dx) 它是做什么的？对于referenc：x1 是指向为新创建的向量分配内存的函数的指针。i - 循环迭代器 xi - 从位置 i dx 处的向量返回一个元素只是一个值。谢谢你的帮助！

4

1 回答 1

0

那里好像少了点什么。它应该是函数参数。否则它是一个无操作 - (x1, i, xi + dx) 是 C 中的一个有效表达式，但它什么都不做。它只提到 x1，然后提到 i，然后提到 xi 和 dx 的总和。

于 2011-09-27T00:41:15.377 回答