-4
public class testrun {


    public static  double ellipticalequn(float x, float y){
         double exp;
         exp=(-2)*Math.exp(Math.pow((y-(.35*6.5)),2)-Math.pow((x-(.75*8)),2)/(2*.05*6.5));
         return exp;
     }
     public static void jacobi(){
         int ny=80;
        int nx=80;
         double[][][] p=new double[100][100][100];
         for(int i=0;i<100;i++){p[i][i][0]=0;} //initialise
         for(int i=0;i<nx;i++){
             p[0][i][0]=10;                 //Boundary conditions
             p[79][i][0]=10;
             p[i][79][0]=50;
             p[i][0][0]=50;

        }
         boolean loop=true;
         int k=0;
         while(loop){
             double error=0;
             for(int i=1;i<nx-1;i++){
                 for(int j=1;j<ny-1;j++){
                     p[i][j][k]=0.25*((p[i-1][j][k]+p[i+1][j][k]+p[i][j+1][k]+p[i][j-1][k])-ellipticalequn((float)i,(float)j));
                     error+=Math.pow(p[i][j][k],2);
                 }
             }
             error=Math.sqrt(error);
             k++;
             if(k>50){
                 loop=false;}

             if(error<(10^(-8))){loop=false;}
             else { loop=true;}
         }
     }

     public static void main(String[] args){
            jacobi();

        }
}

我已经花了 3 个小时,无法找出错误!请在这里帮助我!错误出现了,但我不知道为什么会这样。我是 java 新手。谢谢阿迪

4

1 回答 1

0

你得到数组索引超出范围异常,因为在你的循环中:

 boolean loop=true;
 int k=0;
 while(loop){
     double error=0;
     for(int i=1;i<nx-1;i++){
         for(int j=1;j<ny-1;j++){
             p[i][j][k]=0.25*((p[i-1][j][k]+p[i+1][j][k]+p[i][j+1][k]+p[i][j-1][k])-ellipticalequn((float)i,(float)j));
             error+=Math.pow(p[i][j][k],2);
         }
     }
     error=Math.sqrt(error);
     k++;
     if(k>50){
         loop=false;}

     if(error<(10^(-8))){loop=false;}
     else { loop=true;}
 }

您没有检查 k 永远不会大于 99(如果事实如此,并且当它达到 100 时,您会得到异常)。我不熟悉 Jacobi 算法,但是如果您打算在 k 达到 50 的值时离开循环,那么只需break;在 if 语句中添加一个命令。所以你有两个解决方案,首先:将 if 语句更改为:

if(k>50){
   loop=false;
   break;
}

其次,检查错误后再次检查 k 是否不大于 100:

if(error<(10^(-8)))
    loop=false;
else if(k < 100)
    loop=true;
于 2013-06-07T15:23:07.957 回答