Lines Matching refs:mats

108       Mat mats[9] = {NULL, NULL, A, NULL, A, NULL, A, NULL, NULL};  in main()  local
115 PetscCall(MatCreateNest(PETSC_COMM_WORLD, 3, NULL, 3, NULL, mats, &B)); in main()
118 Mat mats[4]; in main() local
127 PetscCall(CreateIdentity(m, &mats[0])); in main()
128 PetscCall(CreateRandom(m, n, &mats[1])); in main()
129 mats[3] = A; in main()
136 if (!flg) PetscCall(MatCreateTranspose(mats[1], &mats[2])); in main()
137 else PetscCall(MatTranspose(mats[1], MAT_INITIAL_MATRIX, &mats[2])); in main()
143 PetscCall(MatDuplicate(mats[1], MAT_COPY_VALUES, &B)); in main()
144 PetscCall(MatCreateHermitianTranspose(B, &mats[2])); in main()
147 PetscCall(MatScale(mats[2], PetscCMPLX(4.0, -2.0))); in main()
148 … PetscCall(MatShift(mats[2], PetscCMPLX(-2.0, 1.0))); // mats[2] = (4 - 2i) B* - (2 - i) I in main()
149 PetscCall(MatCreateHermitianTranspose(mats[2], &B)); in main()
150 PetscCall(MatDestroy(mats + 2)); in main()
153 PetscCall(MatCreateHermitianTranspose(B, &mats[2])); // mats[2] = B* = (2 - i) B* in main()
155 PetscCall(MatScale(mats[1], PetscCMPLX(2.0, 1.0))); // mats[1] = (2 + i) B = mats[2]* in main()
157 } else PetscCall(MatHermitianTranspose(mats[1], MAT_INITIAL_MATRIX, &mats[2])); in main()
163 PetscCall(MatDuplicate(mats[1], MAT_COPY_VALUES, &B)); in main()
164 PetscCall(MatCreateTranspose(B, &mats[2])); in main()
167 PetscCall(MatScale(mats[2], 4.0)); in main()
168 PetscCall(MatShift(mats[2], -2.0)); // mats[2] = 4 B' - 2 I in main()
169 PetscCall(MatCreateTranspose(mats[2], &B)); in main()
170 PetscCall(MatDestroy(mats + 2)); in main()
173 PetscCall(MatCreateTranspose(B, &mats[2])); // mats[2] = B' = 2 B' in main()
175 PetscCall(MatScale(mats[1], 2.0)); // mats[1] = 2 B = mats[2]' in main()
177 } else PetscCall(MatTranspose(mats[1], MAT_INITIAL_MATRIX, &mats[2])); in main()
179 PetscCall(MatCreateNest(PETSC_COMM_WORLD, 2, NULL, 2, NULL, mats, &B)); in main()
180 PetscCall(MatDestroy(&mats[0])); in main()
181 PetscCall(MatDestroy(&mats[1])); in main()
182 PetscCall(MatDestroy(&mats[2])); in main()