Home
last modified time | relevance | path

Searched refs:Jac (Results 1 – 18 of 18) sorted by relevance

/petsc/src/ts/tests/
H A Dex24.c187 Mat Jac, Jac2; in main() local
209 PetscCall(DMCreateMatrix(da, &Jac)); in main()
226 PetscCall(TSSetRHSJacobian(ts, Jac, Jac, RHSJacobian, &appctx)); in main()
239 PetscCall(TSComputeIJacobian(ts, t, U, Udot, shift, Jac, Jac, PETSC_FALSE)); in main()
241 PetscCall(TSComputeIJacobian(ts, t, U, Udot, shift, Jac, Jac, PETSC_FALSE)); in main()
242 PetscCall(MatAXPY(Jac, -1, Jac2, SAME_NONZERO_PATTERN)); in main()
243 PetscCall(MatNorm(Jac, NORM_INFINITY, &norm)); in main()
245 PetscCall(MatDestroy(&Jac)); in main()
/petsc/src/ts/tutorials/
H A Dex23fwdadj.c20 Mat Jac; /* Jacobian matrix */ member
114 PetscCall(MatCreate(PETSC_COMM_WORLD, &user.Jac)); in main()
115 PetscCall(MatSetSizes(user.Jac, PETSC_DECIDE, PETSC_DECIDE, 1, 1)); in main()
116 PetscCall(MatSetFromOptions(user.Jac)); in main()
117 PetscCall(MatSetUp(user.Jac)); in main()
118 PetscCall(MatCreateVecs(user.Jac, &user.x, NULL)); in main()
123 PetscCall(TSSetIJacobian(ts, user.Jac, user.Jac, IJacobian, &user)); in main()
158 PetscCall(MatCreateVecs(user.Jac, &user.lambda[0], NULL)); in main()
174 PetscCall(MatDestroy(&user.Jac)); in main()
H A Dex20fwd.c26 Mat Jac; /* Jacobian matrix */ member
157 PetscCall(MatCreate(PETSC_COMM_WORLD, &user.Jac)); in main()
158 PetscCall(MatSetSizes(user.Jac, PETSC_DECIDE, PETSC_DECIDE, 2, 2)); in main()
159 PetscCall(MatSetFromOptions(user.Jac)); in main()
160 PetscCall(MatSetUp(user.Jac)); in main()
161 PetscCall(MatCreateVecs(user.Jac, &user.x, NULL)); in main()
169 PetscCall(TSSetIJacobian(ts, user.Jac, user.Jac, IJacobian, &user)); in main()
221 PetscCall(MatDestroy(&user.Jac)); in main()
H A Dex44.c94 PetscScalar Jac; in I2Jacobian() local
100 Jac = shiftA + shiftV * app->Cd * v[0]; in I2Jacobian()
106 PetscCall(MatSetValue(P, i, i, Jac, INSERT_VALUES)); in I2Jacobian()
H A Dex31.c1178 Mat Jac = NULL; /* Jacobian matrix */ in SolveODE() local
1212 PetscCall(MatCreate(PETSC_COMM_WORLD, &Jac)); in SolveODE()
1213 PetscCall(MatSetSizes(Jac, PETSC_DECIDE, PETSC_DECIDE, N, N)); in SolveODE()
1214 PetscCall(MatSetFromOptions(Jac)); in SolveODE()
1215 PetscCall(MatSetUp(Jac)); in SolveODE()
1216 PetscCall(TSSetRHSJacobian(ts, Jac, Jac, RHSJacobian, &ptype[0])); in SolveODE()
1221 PetscCall(MatCreate(PETSC_COMM_WORLD, &Jac)); in SolveODE()
1222 PetscCall(MatSetSizes(Jac, PETSC_DECIDE, PETSC_DECIDE, N, N)); in SolveODE()
1223 PetscCall(MatSetFromOptions(Jac)); in SolveODE()
1224 PetscCall(MatSetUp(Jac)); in SolveODE()
[all …]
H A Dex30.c1456 …ivate(DM dm, PetscReal time, Vec locX, Vec locX_t, PetscReal X_tShift, Mat Jac, Mat JacP, void *us… in DMPlexTSComputeIJacobianFEM_Private() argument
1464 PetscCall(DMPlexTSComputeIJacobianFEM(dm, time, locX, locX_t, 0.0, Jac, JacP, user)); in DMPlexTSComputeIJacobianFEM_Private()
1471 PetscCall(DMPlexTSComputeIJacobianFEM(dm, time, locX, locX_t, X_tShift, Jac, JacP, user)); in DMPlexTSComputeIJacobianFEM_Private()
/petsc/src/ts/tutorials/power_grid/
H A Dex3sa.c55 PetscCall(MatCreate(PETSC_COMM_WORLD, &ctx.Jac)); in main()
56 PetscCall(MatSetSizes(ctx.Jac, n, n, PETSC_DETERMINE, PETSC_DETERMINE)); in main()
57 PetscCall(MatSetType(ctx.Jac, MATDENSE)); in main()
58 PetscCall(MatSetFromOptions(ctx.Jac)); in main()
59 PetscCall(MatSetUp(ctx.Jac)); in main()
60 PetscCall(MatCreateVecs(ctx.Jac, &U, NULL)); in main()
127 PetscCall(TSSetRHSJacobian(ts, ctx.Jac, ctx.Jac, (TSRHSJacobianFn *)RHSJacobian, &ctx)); in main()
146 PetscCall(MatCreateVecs(ctx.Jac, &lambda[0], NULL)); in main()
246 PetscCall(MatDestroy(&ctx.Jac)); in main()
H A Dex3opt.c100 PetscCall(MatCreate(PETSC_COMM_WORLD, &ctx.Jac)); in main()
101 PetscCall(MatSetSizes(ctx.Jac, 2, 2, PETSC_DETERMINE, PETSC_DETERMINE)); in main()
102 PetscCall(MatSetType(ctx.Jac, MATDENSE)); in main()
103 PetscCall(MatSetFromOptions(ctx.Jac)); in main()
104 PetscCall(MatSetUp(ctx.Jac)); in main()
109 PetscCall(MatCreateVecs(ctx.Jac, &ctx.U, NULL)); in main()
122 PetscCall(TSSetRHSJacobian(ctx.ts, ctx.Jac, ctx.Jac, (TSRHSJacobianFn *)RHSJacobian, &ctx)); in main()
126 PetscCall(MatCreateVecs(ctx.Jac, &lambda[0], NULL)); in main()
202 PetscCall(MatDestroy(&ctx.Jac)); in main()
H A Dex3.h16 Mat Jac; /* Jacobian matrix */ member
/petsc/src/snes/tutorials/ex10d/
H A Dex10.c76 Mat Jac; /* Jacobian matrix */ in main() local
369 PetscCall(MatCreate(MPI_COMM_WORLD, &Jac)); in main()
370 PetscCall(MatSetSizes(Jac, PETSC_DECIDE, PETSC_DECIDE, N, N)); in main()
371 PetscCall(MatSetFromOptions(Jac)); in main()
372 PetscCall(MatSetUp(Jac)); in main()
378 PetscCall(MatSetLocalToGlobalMapping(Jac, isl2g, isl2g)); in main()
394 PetscCall(SNESSetJacobian(snes, Jac, Jac, FormJacobian, (void *)&user)); in main()
400 PetscCall(FormJacobian(snes, x, Jac, Jac, &user)); in main()
402 PetscCall(MatColoringCreate(Jac, &mc)); in main()
407 PetscCall(MatFDColoringCreate(Jac, iscoloring, &matfdcoloring)); in main()
[all …]
/petsc/src/dm/interface/
H A Ddmceed.c93 const CeedScalar *x = in[0], *Jac = in[1], *w = in[2]; in Geometry2D() local
99 {Jac[i + Q * 0], Jac[i + Q * 2]}, in Geometry2D()
100 {Jac[i + Q * 1], Jac[i + Q * 3]} in Geometry2D()
117 const CeedScalar *Jac = in[1], *w = in[2]; in Geometry3D() local
123 {Jac[i + Q * 0], Jac[i + Q * 3], Jac[i + Q * 6]}, in Geometry3D()
124 {Jac[i + Q * 1], Jac[i + Q * 4], Jac[i + Q * 7]}, in Geometry3D()
125 {Jac[i + Q * 2], Jac[i + Q * 5], Jac[i + Q * 8]} in Geometry3D()
/petsc/src/binding/petsc4py/demo/legacy/wrap-f2py/
H A DBratu2D.F90125 subroutine JacobianLocal(grd, x, Jac, lambda, ierr) argument
129 Mat Jac
155 PetscCall(MatSetValuesLocal(Jac, ione, row, ione, col, v, INSERT_VALUES, ierr))
168 PetscCall(MatSetValuesLocal(Jac, ione, row, ifive, col, v, INSERT_VALUES, ierr))
/petsc/src/snes/tutorials/network/
H A Dex1.c465 Mat Jac; in main() local
713 PetscCall(SNESGetJacobian(snes, &Jac, NULL, NULL, NULL)); in main()
714 PetscCall(MatView(Jac, PETSC_VIEWER_DRAW_WORLD)); in main()
724 PetscCall(MatView(Jac, PETSC_VIEWER_DRAW_WORLD)); in main()
738 PetscCall(DMCreateMatrix(networkdm, &Jac)); in main()
739 PetscCall(SNESSetJacobian(snes_power, Jac, Jac, FormJacobian_subPower, &user)); in main()
748 PetscCall(SNESGetJacobian(snes_power, &Jac, NULL, NULL, NULL)); in main()
749 PetscCall(MatView(Jac, PETSC_VIEWER_DRAW_WORLD)); in main()
823 PetscCall(MatDestroy(&Jac)); in main()
/petsc/src/ts/utils/
H A Ddmplexts.c191 …anFEM(DM dm, PetscReal time, Vec locX, Vec locX_t, PetscReal X_tShift, Mat Jac, Mat JacP, PetscCtx… in DMPlexTSComputeIJacobianFEM() argument
225 if (hasJac && hasPrec) PetscCall(MatZeroEntries(Jac)); in DMPlexTSComputeIJacobianFEM()
228 …PetscCall(DMPlexComputeJacobianByKey(plex, key, cellIS, time, X_tShift, locX, locX_t, Jac, JacP, c… in DMPlexTSComputeIJacobianFEM()
/petsc/src/ts/tutorials/network/
H A Dpipes.c798 Mat Jac; in main() local
800 PetscCall(SNESGetJacobian(snes, &Jac, NULL, NULL, NULL)); in main()
801 PetscCall(MatView(Jac, PETSC_VIEWER_DRAW_WORLD)); in main()
/petsc/src/snes/utils/
H A Ddmplexsnes.c541 PetscErrorCode DMPlexSNESComputeJacobianFEM(DM dm, Vec X, Mat Jac, Mat JacP, PetscCtx ctx) in DMPlexSNESComputeJacobianFEM() argument
575 if (hasJac && hasPrec) PetscCall(MatZeroEntries(Jac)); in DMPlexSNESComputeJacobianFEM()
578 PetscCall(DMPlexComputeJacobianByKey(plex, key, cellIS, 0.0, 0.0, X, NULL, Jac, JacP, ctx)); in DMPlexSNESComputeJacobianFEM()
/petsc/src/dm/impls/plex/
H A Dplexfem.c4618 …alSection, IS cellIS, PetscReal t, PetscReal X_tShift, Vec X, Vec X_t, Mat Jac, Mat JacP, PetscCtx… in DMPlexComputeJacobian_Patch_Internal() argument
4629 Mat J = Jac, JP = JacP; in DMPlexComputeJacobian_Patch_Internal()
4665 assembleJac = hasJac && hasPrec && (Jac != JacP) ? PETSC_TRUE : PETSC_FALSE; in DMPlexComputeJacobian_Patch_Internal()
4783 …PlexMatSetClosure_Internal(dm, section, globalSection, mesh->useMatClPerm, Jac, cell, &elemMat[(c … in DMPlexComputeJacobian_Patch_Internal()
4797 PetscCall(MatAssemblyBegin(Jac, MAT_FINAL_ASSEMBLY)); in DMPlexComputeJacobian_Patch_Internal()
4798 PetscCall(MatAssemblyEnd(Jac, MAT_FINAL_ASSEMBLY)); in DMPlexComputeJacobian_Patch_Internal()
5906 …etIS, Vec locX, Vec locX_t, PetscReal t, DMField coordField, PetscReal X_tShift, Mat Jac, Mat JacP) in DMPlexComputeBdJacobianSingleByLabel() argument
6050 …exMatSetClosure_Internal(plex, section, globalSection, mesh->useMatClPerm, Jac, support[0], &elemM… in DMPlexComputeBdJacobianSingleByLabel()
6057 …exMatSetClosure_Internal(plex, section, globalSection, mesh->useMatClPerm, Jac, support[0], &elemM… in DMPlexComputeBdJacobianSingleByLabel()
6097 …alues[], PetscInt fieldI, Vec locX, Vec locX_t, PetscReal t, PetscReal X_tShift, Mat Jac, Mat JacP) in DMPlexComputeBdJacobianSingle() argument
[all …]
/petsc/include/
H A Dpetscts.h754 …ts, PetscReal t, Vec U, Vec U_t, Vec U_tt, PetscReal v, PetscReal a, Mat J, Mat Jac, PetscCtx ctx);