Searched refs:Jinvs (Results 1 – 1 of 1) sorted by relevance
| /petsc/src/dm/dt/interface/ |
| H A D | dt.c | 374 PetscScalar *Js, *Jinvs; in PetscDTJacobianInverse_Internal() local
383 PetscCall(PetscMalloc2(m * n, &Js, m * n, &Jinvs)); in PetscDTJacobianInverse_Internal()
387 Jinvs = Jinv; in PetscDTJacobianInverse_Internal()
395 PetscCall(PetscArraycpy(Jinvs, Js, m * m)); in PetscDTJacobianInverse_Internal()
396 PetscCallBLAS("LAPACKgetrf", LAPACKgetrf_(&bm, &bm, Jinvs, &bm, pivots, &info)); in PetscDTJacobianInverse_Internal()
398 PetscCallBLAS("LAPACKgetri", LAPACKgetri_(&bm, Jinvs, &bm, pivots, W, &bm, &info)); in PetscDTJacobianInverse_Internal()
426 Jinvs[i * m + j] = val; in PetscDTJacobianInverse_Internal()
456 Jinvs[i * m + j] = val; in PetscDTJacobianInverse_Internal()
463 for (i = 0; i < m * n; i++) Jinv[i] = PetscRealPart(Jinvs[i]); in PetscDTJacobianInverse_Internal()
464 PetscCall(PetscFree2(Js, Jinvs)); in PetscDTJacobianInverse_Internal()
|