| /petsc/src/ts/trajectory/utils/ |
| H A D | reconstruct.c | 64 …tf(tj->monitor, "Reconstructing at time %g, order %" PetscInt_FMT "\n", (double)t, tj->lag.order)); in TSTrajectoryReconstruct_Private()
65 if (!tj->lag.T) { in TSTrajectoryReconstruct_Private()
66 PetscInt o = tj->lag.order + 1; in TSTrajectoryReconstruct_Private()
67 …PetscCall(PetscMalloc5(o, &tj->lag.L, o, &tj->lag.T, o, &tj->lag.WW, 2 * o, &tj->lag.TT, o, &tj->l… in TSTrajectoryReconstruct_Private()
68 for (i = 0; i < o; i++) tj->lag.T[i] = PETSC_MAX_REAL; in TSTrajectoryReconstruct_Private()
69 PetscCall(VecDuplicateVecs(U ? U : Udot, o, &tj->lag.W)); in TSTrajectoryReconstruct_Private()
72 PetscCall(PetscArrayzero(tj->lag.TT, 2 * (tj->lag.order + 1))); in TSTrajectoryReconstruct_Private()
76 PetscInt up = PetscMin(nid + tj->lag.order / 2 + 1, tshn); in TSTrajectoryReconstruct_Private()
77 PetscInt low = PetscMax(up - tj->lag.order - 1, 0); in TSTrajectoryReconstruct_Private()
78 up = PetscMin(PetscMax(low + tj->lag.order + 1, up), tshn); in TSTrajectoryReconstruct_Private()
[all …]
|
| /petsc/src/ts/trajectory/interface/ |
| H A D | traj.c | 67 if (tj->lag.caching) tj->lag.Udotcached.time = PETSC_MIN_REAL; in TSTrajectorySet()
178 if (U && tj->lag.caching) { in TSTrajectoryGetVecs()
185 …if (id == tj->lag.Ucached.id && *time == tj->lag.Ucached.time && state == tj->lag.Ucached.state) U… in TSTrajectoryGetVecs()
187 …if (id == tj->lag.Ucached.id && stepnum == tj->lag.Ucached.step && state == tj->lag.Ucached.state)… in TSTrajectoryGetVecs()
196 if (Udot && tj->lag.caching) { in TSTrajectoryGetVecs()
203 …if (id == tj->lag.Udotcached.id && *time == tj->lag.Udotcached.time && state == tj->lag.Udotcached… in TSTrajectoryGetVecs()
205 …if (id == tj->lag.Udotcached.id && stepnum == tj->lag.Udotcached.step && state == tj->lag.Udotcach… in TSTrajectoryGetVecs()
246 PetscCall(PetscObjectStateGet((PetscObject)U, &tj->lag.Ucached.state)); in TSTrajectoryGetVecs()
247 PetscCall(PetscObjectGetId((PetscObject)U, &tj->lag.Ucached.id)); in TSTrajectoryGetVecs()
248 tj->lag.Ucached.time = *time; in TSTrajectoryGetVecs()
[all …]
|
| /petsc/src/dm/dt/dualspace/impls/lagrange/ |
| H A D | dspacelagrange.c | 858 PetscDualSpace_Lag *lag = (PetscDualSpace_Lag *)sp->data; in PetscDualSpaceDestroy_Lagrange() local
861 if (lag->symperms) { in PetscDualSpaceDestroy_Lagrange()
862 PetscInt **selfSyms = lag->symperms[0]; in PetscDualSpaceDestroy_Lagrange()
865 PetscInt i, **allocated = &selfSyms[-lag->selfSymOff]; in PetscDualSpaceDestroy_Lagrange()
867 for (i = 0; i < lag->numSelfSym; i++) PetscCall(PetscFree(allocated[i])); in PetscDualSpaceDestroy_Lagrange()
870 PetscCall(PetscFree(lag->symperms)); in PetscDualSpaceDestroy_Lagrange()
872 if (lag->symflips) { in PetscDualSpaceDestroy_Lagrange()
873 PetscScalar **selfSyms = lag->symflips[0]; in PetscDualSpaceDestroy_Lagrange()
877 PetscScalar **allocated = &selfSyms[-lag->selfSymOff]; in PetscDualSpaceDestroy_Lagrange()
879 for (i = 0; i < lag->numSelfSym; i++) PetscCall(PetscFree(allocated[i])); in PetscDualSpaceDestroy_Lagrange()
[all …]
|
| /petsc/src/dm/dt/dualspace/impls/lagrange/tests/ |
| H A D | ex1.c | 222 PetscDualSpace_Lag *lag = (PetscDualSpace_Lag *)sp->data; in testLagrange() local
224 …PetscCall(PetscLagNodeIndicesGetData_Internal(lag->allNodeIndices, &nodeIdxDim, &nodeVecDim, &nNod… in testLagrange()
260 …PetscCall(PetscLagNodeIndicesGetData_Internal(lag->intNodeIndices, &intNodeIdxDim, &intNodeVecDim,… in testLagrange()
270 …PetscCheck(lag->intNodeIndices == lag->allNodeIndices, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Interior … in testLagrange()
|
| /petsc/src/snes/tests/ |
| H A D | ex12f.F90 | 18 PetscInt :: its, lag
112 if (mod(snesm%its, snesm%lag) == 0) then
263 PetscCallA(SNESGetLagJacobian(snes, snesm%lag, ierr))
|
| /petsc/src/snes/interface/ |
| H A D | snes.c | 962 PetscInt i, indx, lag, grids, max_its, max_funcs; in SNESSetFromOptions() local
1008 …ften to rebuild preconditioner", "SNESSetLagPreconditioner", snes->lagpreconditioner, &lag, &flg)); in SNESSetFromOptions()
1010 …PetscCheck(lag != -1, PetscObjectComm((PetscObject)snes), PETSC_ERR_USER, "Cannot set the lag to -… in SNESSetFromOptions()
1011 PetscCall(SNESSetLagPreconditioner(snes, lag)); in SNESSetFromOptions()
1015 …_jacobian", "How often to rebuild Jacobian", "SNESSetLagJacobian", snes->lagjacobian, &lag, &flg)); in SNESSetFromOptions()
1017 …PetscCheck(lag != -1, PetscObjectComm((PetscObject)snes), PETSC_ERR_USER, "Cannot set the lag to -… in SNESSetFromOptions()
1018 PetscCall(SNESSetLagJacobian(snes, lag)); in SNESSetFromOptions()
3608 PetscErrorCode SNESSetLagPreconditioner(SNES snes, PetscInt lag) in SNESSetLagPreconditioner() argument
3612 …PetscCheck(lag >= -2, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Lag must be -2, -1, 1 or greater… in SNESSetLagPreconditioner()
3613 PetscCheck(lag, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Lag cannot be 0"); in SNESSetLagPreconditioner()
[all …]
|
| /petsc/src/snes/impls/nasm/ |
| H A D | nasm.c | 673 PetscInt i, lag = 1; in SNESNASMComputeFinalJacobian_Private() local
709 else if (subsnes->lagjacobian > 1) lag = subsnes->lagjacobian; in SNESNASMComputeFinalJacobian_Private()
712 if (lag > 1) subsnes->lagjacobian = lag; in SNESNASMComputeFinalJacobian_Private()
|
| /petsc/src/ts/impls/rosw/ |
| H A D | rosw.c | 1141 PetscInt lag; in TSStep_RosW() local
1174 PetscCall(SNESGetLagJacobian(snes, &lag)); in TSStep_RosW()
1175 …if (lag == 1) { /* use did not set a nontrivial lag, so lag over all st… in TSStep_RosW()
1180 …if (!ros->recompute_jacobian && i == s - 1 && lag == 1) PetscCall(SNESSetLagJacobian(snes, lag)); … in TSStep_RosW()
|
| /petsc/include/petsc/private/ |
| H A D | tsimpl.h | 119 } lag; member
|
| /petsc/share/petsc/datafiles/meshes/ |
| H A D | testcase3D.cas | 4672 (pdf-transport-lag? #t)
7291 (phaselag/auto-create/default-name-prefix "phase-lag-intf")
|