Question

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的新手。谢谢 AADI

Answer 1

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

 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时离开循环，那么只需在if语句中添加break;命令。所以你有两个解决方案，第一个：将if语句更改为：

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

其次，检查错误后再次检查k是否大于100：

if(error<(10^(-8)))
    loop=false;
else if(k < 100)
    loop=true;

数组超出索引错误 - jacobi方法实现

1 个答案: