Question

下面的代码使用名为 interp 的递归函数，但我找不到一种方法来避免使用 iter 和 fxInterpolated 的全局变量。完整的代码清单（执行N维线性插值）直接编译：

gcc NDimensionalInterpolation.c -o NDimensionalInterpolation -Wall -lm

给出的示例的输出是2.05。代码工作正常，但我想找到全局变量的替代品。任何有关这方面的帮助将不胜感激。感谢。

#include <stdio.h>
#include <stdlib.h>
#include <math.h>

int linearInterpolation(double *, double **, double *, int);
double ** allocateDoubleMatrix(int, int);
double * allocateDoubleVector(int);
void interp(int, int, double *, double *, double *);
double mult(int, double, double *, double *);

/* The objectionable global 
variables that I want to get rid of! */

int iter=0;
double fxInterpolated=0;

int main(int argc, char *argv[]){

 double *fx, **a, *x;
 int dims=2;

 x=allocateDoubleVector(dims);
 a=allocateDoubleMatrix(dims,2);
 fx=allocateDoubleVector(dims*2);

 x[0]=0.25;
 x[1]=0.4;

 a[0][0]=0;
 a[0][1]=1;
 a[1][0]=0;
 a[1][1]=1;

 fx[0]=1;
 fx[1]=3;
 fx[2]=2;
 fx[3]=4;
 linearInterpolation(fx, a, x, dims);
 printf("%f\n",fxInterpolated);

 return (EXIT_SUCCESS);

 } 

 int linearInterpolation(double *fx, double **a, double *x, int dims){

  double *b, *pos;
  int i;

  b=allocateDoubleVector(dims);
  pos=allocateDoubleVector(dims);

  for (i=0; i<dims;i++)
   b[i] = (x[i] - a[i][0]) / (a[i][1] -  a[i][0]);

  interp(0,dims,pos,fx,b); 

  return (EXIT_SUCCESS);

}  

void interp(int j, int dims, double *pos, double *fx, double *b) {

int i;

if (j == dims){
  fxInterpolated+=mult(dims,fx[iter],pos,b);
  iter++;
  return;
}

 for (i = 0; i < 2; i++){
   pos[j]=(double)i; 
   interp(j+1,dims,pos,fx,b);
 }

}

double mult(int dims, double fx, double *pos, double *b){

 int i;
 double val=1.0; 

 for (i = 0; i < dims; i++){
  val *= fabs(1.0-pos[i]-b[i]); 
 } 
 val *= fx;

 printf("mult val= %f fx=%f\n",val, fx);
 return val;

}

double ** allocateDoubleMatrix(int i, int j){

 int k;
 double ** matrix;
 matrix = (double **) calloc(i, sizeof(double *));
 for (k=0; k< i; k++)matrix[k] = allocateDoubleVector(j);
 return matrix;
}

 double * allocateDoubleVector(int i){
  double *vector;
  vector = (double *) calloc(i,sizeof(double));
  return vector;
}

感谢您的评论到目前为止。我想避免使用静态。我删除了全局变量，并且建议尝试使用 iter 变量进行解析。但没有快乐。另外我收到一个编译警告：＆＃34;未使用计算值＆＃34;参考 * iter ++; 我做错了什么？

void interp(int j, int dims, double *pos, double *fx, double *b, int *iter) {

int i;

if (j == dims){
fxInterpolated+=mult(dims,fx[*iter],pos,b);
*iter++;
return; 
}

 for (i = 0; i < 2; i++){
  pos[j]=(double)i; 
  interp(j+1,dims,pos,fx,b,iter);
 }

}

Answer 1

在查看此问题时，我会考虑两种方法：

将状态保留在参数

中

您可以使用一个或多个传递给函数的变量（如果需要，作为指针）以保持函数调用之间的状态。

例如，

int global = 0;
int recursive(int argument) {
  // ... recursive stuff
  return recursive(new_argument);
}

可能会成为

int recursive(int argument, int *global) {
  // ... recursive stuff
  return recursive(new_argument, global);
}

或有时甚至

int recursive(int argument, int global) {
  // ... recursive stuff
  return recursive(new_argument, global);
}

使用静态变量

您还可以使用static关键字在函数中声明要保留的函数中的变量：

int recursive(int argument) {
  static int global = 0;
  // ... recursive stuff
  return recursive(argument);
}

请注意，由于static关键字，global = 0仅在程序启动时设置，而不是每次调用函数时都设置，因为它没有关键字。这意味着如果您更改global的值，它将在下次调用函数时保留此值。

如果您只在程序中使用一次递归函数，则可以使用此方法;如果你需要多次使用它，我建议你使用上面的替代方法。

Answer 2

解决方案是使用静态，然后通过我称之为初始化的标志在第一次调用时重置变量。这样你可以选择重置或不重置变量。

double interp(int j, int dims, double *pos, double *fx, double *b, int initialise) {

 static double fxInterpolated = 0.0;
 static int iter = 0; 
 int i;


 if (initialise){
  fxInterpolated = 0.0;
  iter = 0;
 }

 .....
......
}

在C中使用递归函数时，避免使用全局变量

2 个答案: