jacobi方法覆盖旧的x

出于某种原因，我的雅可比方法覆盖了x_old变量。jacobi方法覆盖旧的x

public static double[] Jacobi(double[][] A, double[] b, double tol) { 
    int m = b.length; 
    double[] x = b; 
    double err = tol*100; 
    while(err > tol) { 
     double[] x_old = x; 
     for(int i = 0; i < m; i++) { 
      double sum = 0.0; 
      for(int j = 0; j < m; j++) { 
       if(i != j){ 
        sum += A[i][j]*x_old[j]; 
       } 
      } 
      System.out.println(sum); 
      x[i] = (b[i]-sum)/A[i][i]; 
     } 
     printVector(x); 
     printVector(x_old); 
     err = norm(subtract(x,x_old)); 

    } 
    return x; 
}

测试代码

double[][] A = {{1.48 , 5.244, 0, -2},{4, 2, 4, 7}, {9, 2, 11, -3}, {-1, 0.2, 3, 12}}; 
    double[][] B = {{92 , 1.3, 0.5, 0.5},{2, 23.3, 1, 0.3}, {0, -2, 28, 3.3}, {-1, 0.2, 3, 12}}; 
    double[][] v = {{2 , 1, 1, 0},{4, 3, 3, 1}, {8, 7, 9, 5}, {6, 7, 9, 8}}; 
    double[] x = {-2, 4, 13.2, 0.22}; 
    double[] y = {1.5, -3, 8.87, 0.6}; 

    double[] c = MyMath.Jacobi(B,y,1e-10);

导致

0.8349999999999993 9.06445652173913 3.015575667102071 0.516473922373577

[0.007228260869565225 -0.5177878335510356 0.20908658331778313 0.0069605064688685785] [0.00722826 0869565225 -0.5177878335510356 0.20908658331778313 0.0069605064688685785]

显然程序退出是因为错误为零。我不明白这是如何工作的。有任何想法吗？

另外一个侧面说明，我有一个覆盖A.

public static double[][] cholesky(double[][] A) { 
    int m = A.length; // rows 
    double akk, akjkk; 
    for(int k = 0 ; k < m ; k++){ 
     akk = A[k][k]; 
     for(int j = k+1 ; j < m ; j++){ 
      akjkk = A[k][j]/akk; 
      for(int i = j ; i < m ; i++){ 
       A[j][i] -= A[k][i]*akjkk; 
      } 
     } 
     for(int i = k ; i < m ; i++){ 
      A[k][i] /= Math.sqrt(akk); 
     } 
    } 
    return A; 
}

我敢肯定有这种明显的解决方案，但是当我通过讲矩阵U到函数，然后尝试另一个功能稍后访问U，我会收到覆盖的版本。有没有办法改变这一点，而不是在函数中复制矩阵？

来源

2016-02-26 Marius Popescu

您每次通过'while'循环重新声明'old_x'。那是你的意图吗？ – bradimus

哈哈不，我被告知要在之前发布的问题中这样做。我认为这至少在解决这个问题方面不会造成伤害。 –

，如果任何人的好奇此代码的工作：

public static double[] Jacobi(double[][] A, double[] b, double tol) { 
    int m = b.length; 
    double[] x_old = new double[m]; 
    double[] x = new double[m]; 
    double err = tol*10; 
    while(err > tol) { 
     System.arraycopy(x,0,x_old,0,m); 
     for(int i = 0; i < m; i++) { 
      double sum = 0.0; 
      for(int j = 0; j < m; j++) { 
       if(i != j){ 
        sum += A[i][j]*x_old[j]; 
       } 
      } 
      x[i] = (b[i]-sum)/A[i][i]; 
     } 
     err = norm(subtract(x,x_old)); 
    } 
    return x; 
}

设置X = x_old时，这可能是一个内存分配的问题。 IIRC在java中使用等号实际上引用了旧变量（我也使用了这个旧变量）。用于Matlab。胆怯的功能，从我的理解，我相信唯一真正的解决方案是在函数范围内调用一个新的变量，我希望这可以避免内存问题，但如果我想存储原始A并创建一个相同大小的新矩阵，我相信整体内存消耗不会有太大的改变。我唯一的问题是我是否可以设置double[][] MatrixInsideFunction = A;或者如果我必须再次执行复制数组。

来源

2016-02-26 16:35:25

jacobi方法覆盖旧的x

回答

相关问题