Searched refs:JJT (Results 1 – 1 of 1) sorted by relevance
| /petsc/src/dm/dt/interface/ |
| H A D | dt.c | 402 PetscScalar *JJT; in PetscDTJacobianInverse_Internal() local
406 PetscCall(PetscMalloc1(m * m, &JJT)); in PetscDTJacobianInverse_Internal()
413 JJT[i * m + j] = val; in PetscDTJacobianInverse_Internal()
417 PetscCallBLAS("LAPACKgetrf", LAPACKgetrf_(&bm, &bm, JJT, &bm, pivots, &info)); in PetscDTJacobianInverse_Internal()
419 PetscCallBLAS("LAPACKgetri", LAPACKgetri_(&bm, JJT, &bm, pivots, W, &bm, &info)); in PetscDTJacobianInverse_Internal()
425 for (k = 0; k < m; k++) val += Js[k * n + i] * JJT[k * m + j]; in PetscDTJacobianInverse_Internal()
430 PetscCall(PetscFree(JJT)); in PetscDTJacobianInverse_Internal()
|