Searched refs:JTJ (Results 1 – 1 of 1) sorted by relevance
| /petsc/src/dm/dt/interface/ |
| H A D | dt.c | 432 PetscScalar *JTJ; in PetscDTJacobianInverse_Internal() local
436 PetscCall(PetscMalloc1(n * n, &JTJ)); in PetscDTJacobianInverse_Internal()
443 JTJ[i * n + j] = val; in PetscDTJacobianInverse_Internal()
447 PetscCallBLAS("LAPACKgetrf", LAPACKgetrf_(&bn, &bn, JTJ, &bn, pivots, &info)); in PetscDTJacobianInverse_Internal()
449 PetscCallBLAS("LAPACKgetri", LAPACKgetri_(&bn, JTJ, &bn, pivots, W, &bn, &info)); in PetscDTJacobianInverse_Internal()
455 for (k = 0; k < n; k++) val += JTJ[i * n + k] * Js[j * n + k]; in PetscDTJacobianInverse_Internal()
460 PetscCall(PetscFree(JTJ)); in PetscDTJacobianInverse_Internal()
|