griejipger

.vscode/settings.json

@@ -10,9 +10,11 @@
         "LAPACKE",
         "MATDENSE",
         "MATMPIAIJ",
+        "matwrite",
         "PETSC",
         "PETSCVIEWERHDF",
         "Reorthogonalization",
+        "spdiagm",
         "VECMPI"
     ]
 }

arnoldi.c

@@ -194,9 +194,10 @@ int main(int argc, char **argv) {
 PetscErrorCode ArnoldiIteration(Mat A, Vec b, PetscInt n, PetscInt m, Vec *Q, double *h) {
     PetscFunctionBeginUser;
 
-    int rank;
+    int rank, total;
     MPI_Comm_rank(PETSC_COMM_WORLD, &rank);
-    printf("Rank %d\n", rank);
+    MPI_Comm_size(PETSC_COMM_WORLD, &total);
+    printf("Process %d of %d\n", rank, total);
 
     PetscScalar eps = 1e-12;
 

matrices/laplacian/julia-laplaciano.jl

@@ -2,8 +2,8 @@ using SparseArrays
 using MAT
 
 # 11 x 16
-nx = 10
-ny = 15
+nx = 11 - 1
+ny = 16 - 1
 ex = fill(1, nx)
 ey = fill(1, ny)
 Dxx = spdiagm(-1 => ex, 0 => -2 * ex, +1 => ex)