Lines Matching refs:R
22 Mat block[4], H, R, C, M; in main() local
37 PetscCall(MatCreate(PETSC_COMM_WORLD, &R)); in main()
38 PetscCall(MatSetSizes(R, PETSC_DECIDE, PETSC_DECIDE, n, n)); in main()
39 PetscCall(MatSetFromOptions(R)); in main()
45 PetscCall(MatGetOwnershipRange(R, &Istart, &Iend)); in main()
47 if (i > 1) PetscCall(MatSetValue(R, i, i - 2, a, INSERT_VALUES)); in main()
48 if (i > 0) PetscCall(MatSetValue(R, i, i - 1, b, INSERT_VALUES)); in main()
49 PetscCall(MatSetValue(R, i, i, c, INSERT_VALUES)); in main()
50 if (i < n - 1) PetscCall(MatSetValue(R, i, i + 1, PetscConj(b), INSERT_VALUES)); in main()
51 if (i < n - 2) PetscCall(MatSetValue(R, i, i + 2, PetscConj(a), INSERT_VALUES)); in main()
61 PetscCall(MatAssemblyBegin(R, MAT_FINAL_ASSEMBLY)); in main()
63 PetscCall(MatAssemblyEnd(R, MAT_FINAL_ASSEMBLY)); in main()
66 block[0] = R; in main()
72 PetscCall(MatCreateTranspose(R, &block[3])); in main()
75 PetscCall(MatCreateNest(PetscObjectComm((PetscObject)R), 2, NULL, 2, NULL, block, &H)); in main()
83 PetscCall(MatDestroy(&R)); in main()