型 "double **" の引数は型 "double *" のパラメーターと互換性がありません

Question

###実現したいこと
関数onestep内で計算したe0の値をmain内で出力したい。
###問題点
関数の関数の部分でのやり取りがうまくいかない。
multiIterからonestepに値を渡すところでうまくいっていないのだと思います。

gyouretu.cc: In function 'void ex::multiIter(int, double (*)[200], int, double*)':
gyouretu.cc:74:25: error: cannot convert 'double**' to 'double*'
   74 |             onestep(n,A,&e0);
      |                         ^~~
      |                         |
      |                         double**
gyouretu.cc:31:47: note:   initializing argument 3 of 'void ex::onestep(int, double (*)[200], double*)'
   31 |     void onestep(int n,double A[N][N],double *e0){
      |                                       ~~~~~~~~^~

###試したこと
そこで、onestep内のe0をdouble **e0にしてみたらプログラムが強制終了してしまいました。
最大固有値の近似値を出すところです。
どこを直せばmain内までe0を持ってくることができるのか教えてください。
###コード

#include<cstdio>
#include<cstdlib>
#include<cmath>

#define N 200
double A[N][N];
double x0[N];

namespace ex{
    double y0[N],y1[N];
    
    //乱数発生
    double uran(){
        double v=rand();
        return v/RAND_MAX;
    }

    //ノルム計算
    double norm(int n,double x[]){
        double s=0;
        int i;
        for(i=0;i<n;i++){
            s=s+x[i]*x[i];
        }
        return sqrt(s);
    }

    //1回のべき乗法の計算
    void onestep(int n,double A[N][N],double *e0){
        double u0;
        int i,j;
        printf("べき乗法:\n");
        for(i=0;i<n;i++){
            printf("x[%d]=%f ",i,x0[i]);
        }
        printf("\n");
        for(i=0;i<n;i++){
            for(j=0;j<n;j++){
                printf("%f ",A[i][j]);
            }
            printf("\n");
        }
        //y=Ax
        for(i=0;i<n;i++){
            //printf("AXi=%d\n",i);
            y0[i]=0;
            for(j=0;j<n;j++){
                y0[i]+=A[i][j]*x0[j];
                //printf("y0=%f\n",y0[i]);
            }
        }
        for(i=0;i<n;i++){
            printf("y[%d]=%f ",i,y0[i]);
        }
        e0=0;
        for(i=0;i<n;i++){
            *e0+=x0[i]*y0[i]; //最大固有値の近似値
        }
        printf("e0=%.10g\n",e0);
        u0=norm(n,y0); //||y0||
        for(i=0;i<n;i++){
            x0[i]=y0[i]/u0; //x=y/||y||
        }
    }

    //べき乗法の指示
    void multiIter(int n,double A[N][N],int kmax,double *e0){
        int k;
        //反復してべき乗法
        for(k=0;k<kmax;k++){
            //printf("k=%d\n",k);
            onestep(n,A,&e0);
        }
        printf("e0=%f\n",e0);
    }
}

int main(int ac,char**av)
{
    int i,j,n,kmax=3000;
    double e0;
    if(ac>1){
        kmax=atoi(av[1]); //文字列をint型に変換
    }
    printf("n=");
    scanf("%d",&n);
    printf("行列式Aを入力:\n");
    for(i=0;i<n;i++){
        for(j=0;j<n;j++){
            scanf("%lf",&A[i][j]);
        }
    }
    printf("ベクトルbを入力:\n");
    for(i=0;i<n;i++){
        scanf("%lf",&x0[i]);
    }

    ex::multiIter(n,A,kmax,&e0);
    printf("最大固有値: %f\n",e0);
}

Answer 1

関数onestepの引数のdouble *e0はe0は「double型の変数のポインタ」を受け取るという意味になります。

関数multiIterの引数も同じくdouble *e0となっていますが、この中でonestepを呼び出す際に
onestep(n,A,&e0);としていますが、こうすると、「double型の変数のポインタ（e0）のポインタ」を渡していることになってしまいます。

おそらく、onestep(n,A,e0);とすれば動くと思います。

Comments

[langhtorn]

できました。無事にmain内で表示できました。ありがとうございます。

Answer 2

変数名を引きずってしまって混乱していたりしませんか? 変数の「名前」はプログラムの動作に影響を与えませんから、ちょっと書き換えてみましょう。

    void onestep(int n,double A[N][N],double *e0){
        //略
        e0=0;//これ、ダメでしょう
        //略
    }

    //べき乗法の指示
    //この関数のローカル変数e0の名前をeeに変える
    void multiIter(int n,double A[N][N],int kmax,double *ee){
        //略
            onestep(n,A,ee);
        }
        printf("e0=%f\n",*ee);
    }

multiIterのローカルのee変数はdouble型です。
なら、onestepの第3バラメータがdouble型を要求しているのですから、(値がどうなっているかは知りませんが)そのまま渡すしかないでしょう。逆に、最後のprintfでは浮動小数点型を要求しているのですから、doubleの変数にはをつけないと、ですね。

余計なことを言えば、そんなところにポインタを使って混乱していないで、素直に関数が返す型をdoubleにすればよかったりしませんか?

Comments

[langhtorn]

そうですね。ポインタを使わないやり方も考えてみます。アドバイスありがとうございます