运行代码时出现奇怪的错误/结果... C fdtd方法

时间:2017-09-15 15:54:35

标签: c simulation

每当我启动for循环为数组赋值时,它会在另一个循环中修改结果,尽管没有在该循环中调用...问题来自epsR ...某种程度上注释/取消注释它会改变结果时间步进循环(计算电极/磁场)。为什么?当espR输入的for循环和if语句被注释掉时,结果很好。如果没有注释掉,我会得到无意义的结果!

是否会像内存损坏一样?我不知道可能导致问题的原因。

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

#define SIZE 200

int main()
{
    double ez[SIZE]={0.}, hy[SIZE]={0.}, imp0=377.0;
    double epsR[SIZE];
    int qTime, maxTime=650, mm;
    char basename[80]="sim", filename[100];
    int frame=0;
    FILE *snapshot;
    /* Init */
    for (mm=0;mm<SIZE;mm++){
        ez[mm]=0.0;
    }
    for (mm=0;mm<SIZE-1;mm++){
        hy[mm]=0.0;
    }
    /* Leaving this loop uncommented modifies the results
     of the loop time stepping*/
    for (mm=0;mm<SIZE;mm++){
        if (mm<100){
            epsR[mm]=1.0;
        }
        else {
            epsR[mm]=9.0;
        }
    }

    /* time stepping*/
    for (qTime=0; qTime<maxTime;qTime++){

        /*Mag field*/
        hy[SIZE-1]=hy[SIZE-2];

        for (mm=0;mm<SIZE-1;mm++){
            hy[mm]=hy[mm]+(ez[mm+1]-ez[mm])/imp0;
        }
        hy[49]-=exp(-(qTime-30.)*(qTime-30.)/100.)/imp0;
        /*Elec field*/
        ez[0]=ez[1];
        ez[SIZE-1]=ez[SIZE-2];

        for (mm=0;mm<SIZE;mm++){
            ez[mm]=ez[mm]+(hy[mm]-hy[mm-1])*imp0/9.0;
        }
        ez[50]+=exp(-(qTime+0.5-(-0.5)-30.)*(qTime+0.5-(-0.5)-30.0)/100.);

        /*Write*/
        if (qTime % 10==0){
            sprintf(filename,"%s.%d", basename, frame++);
            snapshot=fopen(filename,"w");
            for (mm=0;mm<SIZE;mm++){
                fprintf(snapshot,"%g\n",ez[mm]);
            }
        fclose(snapshot);
        }
    }
    return 0;
}

1 个答案:

答案 0 :(得分:1)

在这个循环中:

for (mm=0;mm<SIZE;mm++){
    ez[mm]=ez[mm]+(hy[mm]-hy[mm-1])*imp0/9.0;
}

在第一次迭代中,mm == 0

所以hy[mm-1]hy[-1],这是超出界限的未定义行为。意味着任何事情都可能发生,因为它正在访问随机数据。

如果您希望它访问最后一个元素,请按以下方式检查:

for (mm=0;mm<SIZE;mm++){
    if(mm==0) {
        ez[0]=ez[0]+(hy[0]-hy[SIZE - 1])*imp0/9.0;
    } else {    
        ez[mm]=ez[mm]+(hy[mm]-hy[mm-1])*imp0/9.0;
    }
}