Searched refs:minMN (Results 1 – 1 of 1) sorted by relevance
| /petsc/src/mat/tests/ |
| H A D | ex116.c | 129 PetscBLASInt minMN, maxMN, im, in; in main() local
135 PetscCall(PetscBLASIntCast(PetscMin(m, n), &minMN)); in main()
137 PetscCall(PetscBLASIntCast(5 * minMN + maxMN, &lwork)); in main()
138 …PetscCall(PetscMalloc4(m * minMN, &arrayU, m * minMN, &arrayVT, m * minMN, &arrayErr, lwork, &work… in main()
142 PetscCall(MatSetSizes(Err, PETSC_DECIDE, PETSC_DECIDE, m, minMN)); in main()
148 PetscCall(PetscArraycpy(arrayErr, arrayA, m * minMN)); in main()
153 …LAPACKgesvd_("S", "S", &im, &in, arrayA, &im, evals, arrayU, &minMN, arrayVT, &minMN, work, &lwork… in main()
156 …cCall(PetscPrintf(PETSC_COMM_SELF, " 1st 10 of %" PetscBLASInt_FMT " singular values: \n", minMN)); in main()
164 for (j = 0; j < minMN; j++) { /* U[:,j] = sigma[j]*U[:,j] */ in main()
168 …BLASgemm_("N", "N", &im, &minMN, &minMN, &alpha, arrayU, &im, arrayVT, &minMN, &beta, arrayErr, &i… in main()
|