From 506acf04e336d33e860200cd7fe0325a13e5fd6a Mon Sep 17 00:00:00 2001 From: Antonio De Lucreziis Date: Sat, 4 May 2024 12:18:54 +0200 Subject: [PATCH] chore: added some comments --- main.c | 3 +++ 1 file changed, 3 insertions(+) diff --git a/main.c b/main.c index 0921672..9fca88c 100644 --- a/main.c +++ b/main.c @@ -138,6 +138,7 @@ PetscErrorCode ArnoldiIteration(Mat A, Vec b, PetscInt n, Vec *Q, Mat H) { VecDuplicate(b, &v); MatMult(A, Q[k - 1], v); + // Reorthogonalization using modified Gram-Schmidt for (PetscInt j = 0; j < k; j++) { PetscScalar h; VecDot(Q[j], v, &h); @@ -145,10 +146,12 @@ PetscErrorCode ArnoldiIteration(Mat A, Vec b, PetscInt n, Vec *Q, Mat H) { VecAXPY(v, -h, Q[j]); } + // Normalize PetscScalar h; VecNorm(v, NORM_2, &h); MatSetValue(H, k, k - 1, h, INSERT_VALUES); + // Check for convergence if (h > eps) { VecNormalize(v, NULL); Q[k] = v;