Home
last modified time | relevance | path

Searched refs:jacobian (Results 1 – 25 of 51) sorted by relevance

123

/petsc/src/dm/impls/moab/
H A Ddmmbfem.cxx110 … *phypts, PetscReal *jxw, PetscReal *phi, PetscReal *dphidx, PetscReal *jacobian, PetscReal *ijaco… in Compute_Lagrange_Basis_1D_Internal() argument
115 PetscAssertPointer(jacobian, 9); in Compute_Lagrange_Basis_1D_Internal()
133 jacobian[0] = ijacobian[0] = volume[0] = 0.0; in Compute_Lagrange_Basis_1D_Internal()
136 jacobian[0] += dNi_dxi[i] * vertices[0]; in Compute_Lagrange_Basis_1D_Internal()
141 *volume = jacobian[0]; in Compute_Lagrange_Basis_1D_Internal()
142 ijacobian[0] = 1.0 / jacobian[0]; in Compute_Lagrange_Basis_1D_Internal()
162 jacobian[0] = ijacobian[0] = volume[0] = 0.0; in Compute_Lagrange_Basis_1D_Internal()
165 jacobian[0] += dNi_dxi[i] * vertices[0]; in Compute_Lagrange_Basis_1D_Internal()
170 *volume = jacobian[0]; in Compute_Lagrange_Basis_1D_Internal()
171 ijacobian[0] = 1.0 / jacobian[0]; in Compute_Lagrange_Basis_1D_Internal()
[all …]
/petsc/src/tao/complementarity/impls/ssls/
H A Dssls.c47 PetscCall(TaoComputeJacobian(tao, tao->solution, tao->jacobian, tao->jacobian_pre)); in Tao_SSLS_FunctionGradient()
49 …PetscCall(MatDFischer(tao->jacobian, tao->solution, tao->constraints, tao->XL, tao->XU, ssls->t1, … in Tao_SSLS_FunctionGradient()
50 PetscCall(MatDiagonalScale(tao->jacobian, ssls->db, NULL)); in Tao_SSLS_FunctionGradient()
51 PetscCall(MatDiagonalSet(tao->jacobian, ssls->da, ADD_VALUES)); in Tao_SSLS_FunctionGradient()
52 PetscCall(MatMultTranspose(tao->jacobian, ssls->ff, G)); in Tao_SSLS_FunctionGradient()
/petsc/src/ts/tutorials/autodiff/adolc-utils/
H A Ddrivers.cxx40 else jacobian(tag, m, n, u_vec, J); in PetscAdolcComputeRHSJacobian()
78 else jacobian(tag, m, n, u_vec, J); in PetscAdolcComputeRHSJacobianLocal()
119 else jacobian(tag1, m, n, u_vec, J); in PetscAdolcComputeIJacobian()
135 else jacobian(tag2, m, n, u_vec, J); in PetscAdolcComputeIJacobian()
179 else jacobian(tag, m, n, u_vec, J); in PetscAdolcComputeIJacobianIDMass()
224 else jacobian(tag1, m, n, u_vec, J); in PetscAdolcComputeIJacobianLocal()
239 else jacobian(tag2, m, n, u_vec, J); in PetscAdolcComputeIJacobianLocal()
283 else jacobian(tag, m, n, u_vec, J); in PetscAdolcComputeIJacobianLocalIDMass()
422 else jacobian(tag1, m, n, u_vec, J); in PetscAdolcComputeIJacobianAndDiagonalLocal()
435 else jacobian(tag2, m, n, u_vec, J); in PetscAdolcComputeIJacobianAndDiagonalLocal()
/petsc/src/tao/complementarity/impls/asls/
H A Dasils.c77 PetscCall(TaoComputeJacobian(tao, tao->solution, tao->jacobian, tao->jacobian_pre)); in Tao_ASLS_FunctionGradient()
78 …PetscCall(MatDFischer(tao->jacobian, tao->solution, tao->constraints, tao->XL, tao->XU, asls->t1, … in Tao_ASLS_FunctionGradient()
80 PetscCall(MatMultTranspose(tao->jacobian, asls->t1, G)); in Tao_ASLS_FunctionGradient()
166 PetscCall(MatNorm(tao->jacobian, NORM_1, &asls->identifier)); in TaoSolve_ASILS()
207 PetscCall(MatMult(tao->jacobian, tao->stepdirection, asls->t1)); in TaoSolve_ASILS()
213 PetscCall(TaoMatGetSubMat(tao->jacobian, asls->free, asls->w, tao->subset_type, &asls->J_sub)); in TaoSolve_ASILS()
214 if (tao->jacobian != tao->jacobian_pre) { in TaoSolve_ASILS()
H A Dasfls.c76 PetscCall(TaoComputeJacobian(tao, tao->solution, tao->jacobian, tao->jacobian_pre)); in Tao_ASLS_FunctionGradient()
78 …PetscCall(MatDFischer(tao->jacobian, tao->solution, tao->constraints, tao->XL, tao->XU, asls->t1, … in Tao_ASLS_FunctionGradient()
80 PetscCall(MatMultTranspose(tao->jacobian, asls->t1, G)); in Tao_ASLS_FunctionGradient()
169 PetscCall(MatNorm(tao->jacobian, NORM_1, &asls->identifier)); in TaoSolve_ASFLS()
210 PetscCall(MatMult(tao->jacobian, tao->stepdirection, asls->t1)); in TaoSolve_ASFLS()
215 PetscCall(TaoMatGetSubMat(tao->jacobian, asls->free, asls->w, tao->subset_type, &asls->J_sub)); in TaoSolve_ASFLS()
216 if (tao->jacobian != tao->jacobian_pre) { in TaoSolve_ASFLS()
/petsc/src/binding/petsc4py/test/
H A Dtest_ts_py.py24 def jacobian(self, ts, t, u, du, a, J, P): member in MyODE
107 ts.setIJacobian(ode.jacobian, J, J)
149 ts.setIJacobian(ode.jacobian, J, J)
161 ode.jacobian(ts, 0.0, u, u, 1.0, J, J)
/petsc/src/ts/tutorials/network/
H A DpipeInterface.c107 pipe->jacobian = Jin; in PipeCreateJacobian()
156 pipe->jacobian = Jpipe; in PipeCreateJacobian()
162 Mat *Jpipe = pipe->jacobian; in PipeDestroyJacobian()
264 Mat *Jv = junc->jacobian; in JunctionDestroyJacobian()
H A Dwash.h37 Mat *jacobian; member
H A Dpipe.h31 Mat *jacobian; member
/petsc/src/tao/leastsquares/tutorials/matlab/more_wild_probs/
H A Djacobian.m1 function J = jacobian(m,n,x,nprob) function
15 % J = jacobian(m,n,x,nprob)
/petsc/src/snes/impls/vi/ss/
H A Dviss.c263 PetscCall(SNESComputeJacobian(snes, X, snes->jacobian, snes->jacobian_pre)); in SNESSolve_VINEWTONSSLS()
268 PetscCall(SNESVIComputeBsubdifferentialVectors(snes, X, F, snes->jacobian, vi->Da, vi->Db)); in SNESSolve_VINEWTONSSLS()
270 PetscCall(SNESVIComputeJacobian(snes->jacobian, snes->jacobian_pre, vi->Da, vi->Db)); in SNESSolve_VINEWTONSSLS()
272 PetscCall(SNESVIComputeMeritFunctionGradient(snes->jacobian, vi->phi, vi->dpsi)); in SNESSolve_VINEWTONSSLS()
273 PetscCall(KSPSetOperators(snes->ksp, snes->jacobian, snes->jacobian_pre)); in SNESSolve_VINEWTONSSLS()
317 PetscCall(SNESVICheckLocalMin_Private(snes, snes->jacobian, vi->phi, X, gnorm, &ismin)); in SNESSolve_VINEWTONSSLS()
/petsc/src/binding/petsc4py/src/petsc4py/PETSc/
H A Dpetsctao.pxi417 (jacobian, args, kargs) = context
418 jacobian(tao, x, J, P, *args, **kargs)
433 (jacobian, args, kargs) = context
434 jacobian(tao, x, J, P, *args, **kargs)
451 (jacobian, args, kargs) = context
452 jacobian(tao, x, J, Jp, Ji, *args, **kargs)
465 (jacobian, args, kargs) = context
466 jacobian(tao, x, J, *args, **kargs)
509 (jacobian, args, kargs) = context
510 jacobian(tao, x, J, P, *args, **kargs)
[all …]
H A DTS.pyx635 jacobian: TSRHSJacobian | None,
646 jacobian
649 The matrix into which the jacobian is computed.
653 Additional positional arguments for ``jacobian``.
655 Additional keyword arguments for ``jacobian``.
666 if jacobian is not None:
669 context = (jacobian, args, kargs)
800 cdef object jacobian = self.get_attr('__rhsjacobian__')
801 return (J, P, jacobian)
844 jacobian: TSIJacobian | None,
[all …]
H A Dpetscts.pxi399 (jacobian, args, kargs) = context
400 jacobian(Ts, toReal(t), Xvec, Jmat, Pmat, *args, **kargs)
459 (jacobian, args, kargs) = context
460 jacobian(Ts, toReal(t), Xvec, XDvec, toReal(a), Jmat, Pmat, *args, **kargs)
479 (jacobian, args, kargs) = context
480 jacobian(Ts, toReal(t), Xvec, XDvec, toReal(a), Jmat, *args, **kargs)
525 (jacobian, args, kargs) = context
526 jacobian(Ts, toReal(t), Xvec, XDvec, XDDvec, toReal(v), toReal(a), Jmat, Pmat, *args, **kargs)
/petsc/src/snes/impls/ls/
H A Dls.c217 PetscCall(SNESComputeJacobian(snes, X, snes->jacobian, snes->jacobian_pre)); in SNESSolve_NEWTONLS()
219 PetscCall(KSPSetOperators(snes->ksp, snes->jacobian, snes->jacobian_pre)); in SNESSolve_NEWTONLS()
225 if (PetscLogPrintInfo) PetscCall(SNESNEWTONLSCheckResidual_Private(snes, snes->jacobian, F, Y)); in SNESSolve_NEWTONLS()
260 PetscCall(SNESNEWTONLSCheckLocalMin_Private(snes, snes->jacobian, F, fnorm, &ismin)); in SNESSolve_NEWTONLS()
/petsc/src/ts/tutorials/
H A Dex8.py29 def jacobian(self,ts,t,x,xdot,shift,J,P): member in MyODE
57 ts.setIJacobian(ode.jacobian, J)
H A Dex8.c18 TSIJacobianFn *jacobian; member
99 p->jacobian = &RoberJacobian; in RoberCreate()
187 p->jacobian = &CEJacobian; in CECreate()
273 p->jacobian = &OregoJacobian; in OregoCreate()
381 PetscCall(TSSetIJacobian(ts, A, A, problem->jacobian, problem->data)); in main()
/petsc/src/snes/impls/ksponly/
H A Dksponly.c42 PetscCall(SNESComputeJacobian(snes, X, snes->jacobian, snes->jacobian_pre)); in SNESSolve_KSPONLY()
45 PetscCall(KSPSetOperators(snes->ksp, snes->jacobian, snes->jacobian_pre)); in SNESSolve_KSPONLY()
/petsc/src/snes/interface/ftn-custom/
H A Dzsnesf.c71 PetscFortranCallbackId jacobian; member
178 …PetscObjectUseFortranCallback(snes, _cb.jacobian, (SNES *, Vec *, Mat *, Mat *, void *, PetscError… in oursnesjacobian()
217 …tFortranCallback((PetscObject)*snes, PETSC_FORTRAN_CALLBACK_CLASS, &_cb.jacobian, (PetscFortranCal… in snessetjacobian_()
241 …bjectGetFortranCallback((PetscObject)*snes, PETSC_FORTRAN_CALLBACK_CLASS, _cb.jacobian, NULL, ctx); in snesgetjacobian_()
256 …PetscObjectUseFortranCallback(snes, _cb.jacobian, (SNES *, Vec *, Mat *, Mat *, void *, PetscError… in oursnespicardjacobian()
266 …tFortranCallback((PetscObject)*snes, PETSC_FORTRAN_CALLBACK_CLASS, &_cb.jacobian, (PetscFortranCal… in snessetpicard_()
/petsc/src/snes/interface/
H A Dsnes.c591 if (!snes->vec_func && (snes->jacobian || snes->jacobian_pre)) { in SNESSetUpMatrixFree_Private()
592 Mat A = snes->jacobian, B = snes->jacobian_pre; in SNESSetUpMatrixFree_Private()
612 if (snes->jacobian) { in SNESSetUpMatrixFree_Private()
613 PetscCall(MatGetNullSpace(snes->jacobian, &nullsp)); in SNESSetUpMatrixFree_Private()
626 if (!snes->jacobian) PetscCall(SNESSetJacobian(snes, J, NULL, NULL, NULL)); in SNESSetUpMatrixFree_Private()
767 if (!snes->jacobian && snes->mf && !snes->mf_operator && !snes->jacobian_pre) { in SNESSetUpMatrices()
776 } else if (snes->mf_operator && !snes->jacobian_pre && !snes->jacobian) { in SNESSetUpMatrices()
791 J = snes->jacobian; in SNESSetUpMatrices()
1098 …PetscCall(SNESSetJacobian(snes, snes->jacobian, snes->jacobian_pre, SNESComputeJacobianDefault, fu… in SNESSetFromOptions()
1112 …PetscCall(SNESSetJacobian(snes, snes->jacobian, snes->jacobian_pre, SNESComputeJacobianDefaultColo… in SNESSetFromOptions()
[all …]
/petsc/doc/changes/
H A D216.md71 - -dmmg_jacobian_period -1 => jacobian is frozen after initial
77 - Allow user preallocation of jacobian using DASetGetMatrix()
/petsc/src/snes/impls/gs/
H A Dgssecant.c41 if (!snes->jacobian) PetscCall(SNESSetUpMatrices(snes)); in SNESComputeNGSDefaultSecant()
42 PetscCall(MatColoringCreate(snes->jacobian, &mc)); in SNESComputeNGSDefaultSecant()
/petsc/src/snes/impls/vi/rs/
H A Dvirs.c357 PetscCall(SNESComputeJacobian(snes, X, snes->jacobian, snes->jacobian_pre)); in SNESSolve_VINEWTONRSLS()
381 …PetscCall(MatCreateSubMatrix(snes->jacobian, vi->IS_inact, vi->IS_inact, MAT_INITIAL_MATRIX, &jac_… in SNESSolve_VINEWTONRSLS()
407 …PetscCall(MatCreateSubMatrix(snes->jacobian, vi->IS_inact, vi->IS_inact, MAT_INITIAL_MATRIX, &jac_… in SNESSolve_VINEWTONRSLS()
446 …if (snes->jacobian != snes->jacobian_pre) PetscCall(MatCreateSubMatrix(snes->jacobian_pre, vi->IS_… in SNESSolve_VINEWTONRSLS()
508 if (snes->jacobian != snes->jacobian_pre) PetscCall(MatDestroy(&prejac_inact_inact)); in SNESSolve_VINEWTONRSLS()
561 PetscCall(SNESVICheckLocalMin_Private(snes, snes->jacobian, F, X, gnorm, &ismin)); in SNESSolve_VINEWTONRSLS()
/petsc/src/tao/leastsquares/tutorials/matlab/
H A DProblemInitialize.m17 jac = @(x)jacobian(m,n,x,nprob);
/petsc/src/snes/impls/ncg/
H A Dsnesncg.c38 if (!snes->jacobian) PetscCall(SNESSetUpMatrices(snes)); in SNESLineSearchApply_NCGLinear()
44 PetscCall(SNESComputeJacobian(snes, X, snes->jacobian, snes->jacobian_pre)); in SNESLineSearchApply_NCGLinear()
48 PetscCall(MatMult(snes->jacobian, Y, W)); in SNESLineSearchApply_NCGLinear()

123