0

あるクラスでは、微分方程式を実装して、両端が氷浴に浸された 100 度の棒の温度を求めるプログラムを作成することになっています。それを解くために微分方程式が与えられます。ロッドのセグメントとして x 方向に分割され、時間の反復として y 方向に分割された配列があります。r 値は、ロッドが冷却する速度を決定します。実際のプログラムがどのように実装されるかはそれほど重要ではないことに注意してください。このプログラムはコンピューター クラスター上で実行することを意図しており、割り当てのポイントは、クラスター上でジョブを発行する方法を紹介することです。ただし、プログラムは (当然のことながら) 正しくなければなりません。

私の配列は時間の繰り返しで（一見）ランダムなポイントでガベージデータを取得し始めているように見えるため、プログラムに問題があります。私は Netbeans を持っていて、これをデバッグしようとしましたが、最後に実行したときは、配列変数のさまざまな要素を覗くことができず (デバッグ プロセスが実質的に役に立たなくなりました)、今ではブレーク ポイントで停止することさえできます。私はしばらくの間、この問題に頭を悩ませてきましたが、私よりもずっと賢い誰かが単に問題を「見て」助けてくれることを望んでいます。

ありがとう。

これが私のコードです：

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

int main(void)
{
    int segmentsl;
    int segmentst;
    float rvalue;
    int i,j;
    float k;

    //Pull in initial data
    printf("Enter as integers without spaces: Number of Segments - Length, Number of Segments - time ,value\n");
    scanf("%d,%d,%f", &segmentsl,&segmentst,&rvalue);

    float tempvstime[segmentsl][segmentst];


    //at t0 -> temp = 100sin(pi*x)
    for(i = 0; i < segmentsl; i++){
          k = (float)(i/(segmentsl-1));
          tempvstime[i][0] = 100*sin(M_PI*k); 
          printf("%f,",tempvstime[i][0]); 
    }

    printf("\n\n\nEND OF INITIALIZATION \n\n\n"); 

    for(j = 0; j < (segmentst-1); j++){
      for(i = 0; i < segmentsl; i++){
        if(i == 0 || i == (segmentsl - 1)){tempvstime[i][j] = 0;}
        else{
        tempvstime[i][j+1] = (rvalue*tempvstime[i-1][j]) + ((1-(2*rvalue))*tempvstime[i][j]) + (rvalue*tempvstime[i+1][j]); 
        }
      }
    }


    for(j = 0; j < segmentst; j++){
            for(i = 0; i < segmentsl; i++){
            printf("%f,",tempvstime[i][j]);
                if(i == segmentsl - 1){printf("\n");}
            }
    }

}

出力例を次に示します。

Enter as integers without spaces: Number of Segments - Length, Number of Segments - time ,value
6,10,0.5
0.000000,58.778526,95.105652,95.105652,58.778522,0.000000,


END OF INITIALIZATION 


0.000000,58.778526,95.105652,95.105652,58.778522,0.000000,
0.000000,47.552826,76.942093,76.942085,47.552826,0.000000,
0.000000,38.471046,62.247456,62.247459,38.471043,0.000000,
0.000000,31.123728,50.359253,50.359249,31.123730,0.000000,
0.000000,25.179626,40.741489,40.741493,25.179625,0.000000,
0.000000,20.370745,32.960560,32.960556,-nan,0.000000,
0.000000,16.480280,26.665649,-nan,-nan,0.000000,
0.000000,13.332825,-nan,-nan,-nan,0.000000,
0.000000,-nan,-nan,-nan,-nan,0.000000,
0.000000,-nan,-nan,-nan,-nan,0.000000,
4

2 に答える 2

4

の整数除算とは別に

k = (float)(i/(segmentsl-1));

間違ったエンドポイントを設定しました:

for(j = 0; j < (segmentst-1); j++){
  for(i = 0; i < segmentsl; i++){
    if(i == 0 || i == (segmentsl - 1)){tempvstime[i][j] = 0;}
    else{
    tempvstime[i][j+1] = (rvalue*tempvstime[i-1][j]) + ((1-(2*rvalue))*tempvstime[i][j]) + (rvalue*tempvstime[i+1][j]); 
    }
  }
}

iが 0 または最後のインデックスの場合は を設定tempvstime[i][j]し、それ以外のi場合は を設定しtempvstime[i][j+1]ます。

tempvstime[i][j+1]エンドポイントにも設定する必要があります。そうしないと、次の繰り返しで、初期化されていない (ガベージ) 値が使用され、未定義の動作が呼び出されます。

その変更と、設定時の浮動小数点除算により、k妥当な値が得られます。

6,10,0.5
0.000000,58.778526,95.105652,95.105652,58.778522,0.000000,


END OF INITIALIZATION 


0.000000,58.778526,95.105652,95.105652,58.778522,0.000000,
0.000000,47.552826,76.942093,76.942085,47.552826,0.000000,
0.000000,38.471046,62.247456,62.247459,38.471043,0.000000,
0.000000,31.123728,50.359253,50.359249,31.123730,0.000000,
0.000000,25.179626,40.741489,40.741493,25.179625,0.000000,
0.000000,20.370745,32.960560,32.960556,20.370747,0.000000,
0.000000,16.480280,26.665649,26.665653,16.480278,0.000000,
0.000000,13.332825,21.572968,21.572964,13.332827,0.000000,
0.000000,10.786484,17.452894,17.452896,10.786482,0.000000,
0.000000,8.726447,14.119690,14.119688,8.726448,0.000000,
于 2013-06-15T22:28:06.150 に答える
0

この行に気付きました：

k = (float)(i/(segmentsl-1));

整数除算を行っています。i<segmentslこれは 0の場合になるためですi/segmentsl。これはおそらく意図したものではありません。

また、プログラムの後半で and を使用tempvstime[i][j+1]しtempvstime[i][j]、両方とも である必要がある場合もありますj+1。

正しく実行され、変数名が変更され、数学の括弧が少なくなったプログラムの修正版が続きます。

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

int main(void){
    int width, maxt;
    float rvalue;
    int i,t;
    float k;

    //Pull in initial data
    printf("Enter as integers without spaces: Number of Segments - Length, Number of Segments - time ,value\n");
    scanf("%d,%d,%f", &width,&maxt,&rvalue);

    float tempvstime[width][maxt];

    fprintf(stderr,"Initializing with %d segments...",width);
    //at t0 -> temp = 100sin(pi*x)
    for(i = 0; i < width; ++i)
      tempvstime[i][0] = 100*sin((M_PI*i)/(width-1));
    fprintf(stderr,"done.\n");

    for(t = 0; t <(maxt-1); ++t){
      for(i = 0; i < width; ++i){
        if( i == 0 || i == (width-1) ){
          tempvstime[i][t+1] = 0;
        } else {
          tempvstime[i][t+1] = rvalue*tempvstime[i-1][t] + (1-2*rvalue)*tempvstime[i][t] + rvalue*tempvstime[i+1][t];
        }
      }
    }

    for(t = 0; t < maxt; t++){
      for(i = 0; i < width; i++)
        printf("%f,",tempvstime[i][t]);
      printf("\n");
    }

}
于 2013-06-15T22:24:15.327 に答える