diff --git a/arnoldi.c b/arnoldi.c index 7c4a52a..eab1b31 100644 --- a/arnoldi.c +++ b/arnoldi.c @@ -240,7 +240,6 @@ int main(int argc, char **argv) { PetscCall(MatDestroy(&A)); PetscCall(VecDestroy(&b)); // PetscCall(MatDestroy(&H)); - PetscFinalize(); return 0; } @@ -281,8 +280,7 @@ PetscErrorCode ArnoldiIteration(Mat A, Vec b, PetscInt n, PetscInt m, Vec *Q, do // Reorthogonalization using modified Gram-Schmidt for (PetscInt j = 0; j < k - 1; j++) { // anche solo 3 PetscScalar h_ij; - - // h_ij = q_j . v + // h_(j, k - 1) = q_j . v PetscCall(VecDot(Q[j], v, &h_ij)); // v -= h_ij * q_j