Beispiel #1
0
  @Override
  public Matrix run(Matrix coef, Matrix free) {
    int n = coef.getHeight();
    double[][] alpha = new double[n][n];
    double[][] beta = new double[n][1];
    for (int i = 0; i < n; i++) {
      for (int j = 0; j < n; j++) {
        alpha[i][j] = (i == j ? 0. : -coef.get(i, j) / coef.get(i, i));
      }
      beta[i][0] = free.get(i, 0) / coef.get(i, i);
    }
    Matrix alphaMatrix = new Matrix(alpha);
    Matrix betaMatrix = new Matrix(beta);
    alphaNorm = alphaMatrix.getNorm();

    Matrix cur = new Matrix(approximate);
    Matrix prev;
    iteration = 0;
    do {
      iteration++;
      prev = new Matrix(cur);
      cur = alphaMatrix.multiply(prev).add(betaMatrix);
    } while (check(cur, prev) && iteration < allowedIterations);

    return new Matrix(cur);
  }
Beispiel #2
0
 private boolean check(Matrix cur, Matrix prev) {
   return cur.subtract(prev).getVectorNorm() > (1 - alphaNorm) * eps / alphaNorm;
 }