Lines Matching refs:nmat
18 PetscInt nmat = 3; /* number of matrices */ in main() local
27 PetscCall(PetscOptionsGetInt(NULL, NULL, "-nmat", &nmat, NULL)); in main()
32 PetscCall(PetscMalloc1(nmat + 3, &A)); in main()
35 …for (i = 1; i < nmat + 1; i++) PetscCall(MatCreateAIJ(PETSC_COMM_WORLD, PETSC_DECIDE, PETSC_DECIDE… in main()
36 …CreateAIJ(PETSC_COMM_WORLD, PETSC_DECIDE, PETSC_DECIDE, n / 2, n, 3, NULL, 3, NULL, &A[nmat + 1])); in main()
37 for (i = 0; i < nmat + 2; i++) PetscCall(MatSetRandom(A[i], rctx)); in main()
51 for (i = 2; i < nmat + 1; i++) { in main()
58 PetscCall(MatCreateComposite(PETSC_COMM_WORLD, nmat, A + 1, &B)); in main()
85 for (i = 1; i < nmat; i++) { in main()
92 PetscCall(MatCreateComposite(PETSC_COMM_WORLD, nmat, A, &B)); in main()
112 for (i = 3; i < nmat + 1; i++) { in main()
116 PetscCall(MatMult(A[nmat + 1], z, v)); in main()
118 PetscCall(MatCreateComposite(PETSC_COMM_WORLD, nmat, A + 2, &B)); in main()
133 PetscCall(MatCreateComposite(PETSC_COMM_WORLD, nmat, A, &B)); in main()
135 …if (nmat != n) PetscCall(PetscPrintf(PETSC_COMM_WORLD, "Error with GetNumberMat %" PetscInt_FMT " … in main()
136 PetscCall(MatCompositeGetMat(B, 0, &A[nmat + 2])); in main()
137 if (A[0] != A[nmat + 2]) PetscCall(PetscPrintf(PETSC_COMM_WORLD, "Error with GetMat\n")); in main()
153 for (i = 0; i < nmat + 2; i++) PetscCall(MatDestroy(&A[i])); in main()