| /petsc/src/tao/tutorials/ |
| H A D | ex4.c | 551 PetscReal *Js, *hs; in TaylorTest() local
570 PetscCall(PetscCalloc2(numValues, &Js, numValues, &hs)); in TaylorTest()
579 Js[i] = J; in TaylorTest()
583 temp = PetscLogReal(Js[j] / Js[j - 1]) / PetscLogReal(hs[j] / hs[j - 1]); in TaylorTest()
588 PetscCall(PetscFree2(Js, hs)); in TaylorTest()
|
| /petsc/src/tao/pde_constrained/tutorials/ |
| H A D | hyperbolic.c | 16 Mat Js, Jd, JsBlockPrec, JsInv, JsBlock; member
179 …PetscCall(TaoSetJacobianStateRoutine(tao, user.Js, user.Js, user.JsInv, FormJacobianState, &user)); in main()
1012 …reateShell(PETSC_COMM_WORLD, PETSC_DETERMINE, PETSC_DETERMINE, user->m, user->m, user, &user->Js)); in HyperbolicInitialize()
1013 PetscCall(MatShellSetOperation(user->Js, MATOP_MULT, (PetscErrorCodeFn *)StateMatMult)); in HyperbolicInitialize()
1014 PetscCall(MatShellSetOperation(user->Js, MATOP_DUPLICATE, (PetscErrorCodeFn *)StateMatDuplicate)); in HyperbolicInitialize()
1015 …PetscCall(MatShellSetOperation(user->Js, MATOP_MULT_TRANSPOSE, (PetscErrorCodeFn *)StateMatMultTra… in HyperbolicInitialize()
1016 …PetscCall(MatShellSetOperation(user->Js, MATOP_GET_DIAGONAL, (PetscErrorCodeFn *)StateMatGetDiagon… in HyperbolicInitialize()
1075 PetscCall(StateMatInvMult(user->Js, user->q, user->ytrue)); in HyperbolicInitialize()
1079 PetscCall(StateMatInvMult(user->Js, user->q, user->y)); in HyperbolicInitialize()
1127 PetscCall(MatDestroy(&user->Js)); in HyperbolicDestroy()
|
| H A D | parabolic.c | 19 Mat Js, Jd, JsBlockPrec, JsInv, JsBlock; member
185 …PetscCall(TaoSetJacobianStateRoutine(tao, user.Js, user.JsBlockPrec, user.JsInv, FormJacobianState… in main()
918 …WORLD, PETSC_DETERMINE, PETSC_DETERMINE, user->m * user->nt, user->m * user->nt, user, &user->Js)); in ParabolicInitialize()
919 PetscCall(MatShellSetOperation(user->Js, MATOP_MULT, (PetscErrorCodeFn *)StateMatMult)); in ParabolicInitialize()
920 PetscCall(MatShellSetOperation(user->Js, MATOP_DUPLICATE, (PetscErrorCodeFn *)StateMatDuplicate)); in ParabolicInitialize()
921 …PetscCall(MatShellSetOperation(user->Js, MATOP_MULT_TRANSPOSE, (PetscErrorCodeFn *)StateMatMultTra… in ParabolicInitialize()
922 …PetscCall(MatShellSetOperation(user->Js, MATOP_GET_DIAGONAL, (PetscErrorCodeFn *)StateMatGetDiagon… in ParabolicInitialize()
1045 PetscCall(StateMatInvMult(user->Js, user->q, user->ytrue)); in ParabolicInitialize()
1070 PetscCall(StateMatInvMult(user->Js, user->q, user->y)); in ParabolicInitialize()
1203 PetscCall(MatDestroy(&user->Js)); in ParabolicDestroy()
|
| H A D | elliptic.c | 19 Mat Js, Jd, JsPrec, JsInv, JsBlock; member
152 PetscCall(TaoSetJacobianStateRoutine(tao, user.Js, NULL, user.JsInv, FormJacobianState, &user)); in main()
1009 …WORLD, ysubnlocal * user->ns, ysubnlocal * user->ns, user->nstate, user->nstate, user, &user->Js)); in EllipticInitialize()
1010 PetscCall(MatShellSetOperation(user->Js, MATOP_MULT, (PetscErrorCodeFn *)StateMatMult)); in EllipticInitialize()
1011 PetscCall(MatShellSetOperation(user->Js, MATOP_MULT_TRANSPOSE, (PetscErrorCodeFn *)StateMatMult)); in EllipticInitialize()
1012 PetscCall(MatSetOption(user->Js, MAT_SYMMETRIC, PETSC_TRUE)); in EllipticInitialize()
1013 PetscCall(MatSetOption(user->Js, MAT_SYMMETRY_ETERNAL, PETSC_TRUE)); in EllipticInitialize()
1186 PetscCall(MatDestroy(&user->Js)); in EllipticDestroy()
|
| /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()
384 for (i = 0; i < m * n; i++) Js[i] = J[i]; in PetscDTJacobianInverse_Internal()
386 Js = (PetscReal *)J; in PetscDTJacobianInverse_Internal()
395 PetscCall(PetscArraycpy(Jinvs, Js, m * m)); in PetscDTJacobianInverse_Internal()
412 for (k = 0; k < n; k++) val += Js[i * n + k] * Js[j * n + k]; in PetscDTJacobianInverse_Internal()
425 for (k = 0; k < m; k++) val += Js[k * n + i] * JJT[k * m + j]; in PetscDTJacobianInverse_Internal()
442 for (k = 0; k < m; k++) val += Js[k * n + i] * Js[k * n + j]; in PetscDTJacobianInverse_Internal()
455 for (k = 0; k < n; k++) val += JTJ[i * n + k] * Js[j * n + k]; in PetscDTJacobianInverse_Internal()
464 PetscCall(PetscFree2(Js, Jinvs)); in PetscDTJacobianInverse_Internal()
|