45   TSAdapt_DSP *dsp = (TSAdapt_DSP *)adapt->data;  in TSAdaptRestart_DSP()  local
48   dsp->cerror[0] = dsp->hratio[0] = 1.0;  in TSAdaptRestart_DSP()
49   dsp->cerror[1] = dsp->hratio[1] = 1.0;  in TSAdaptRestart_DSP()
50   dsp->cerror[2] = dsp->hratio[2] = 1.0;  in TSAdaptRestart_DSP()
56   TSAdapt_DSP *dsp = (TSAdapt_DSP *)adapt->data;  in TSAdaptRollBack_DSP()  local
59   dsp->cerror[0] = dsp->cerror[1];  in TSAdaptRollBack_DSP()
60   dsp->cerror[1] = dsp->cerror[2];  in TSAdaptRollBack_DSP()
61   dsp->cerror[2] = 1.0;  in TSAdaptRollBack_DSP()
62   dsp->hratio[0] = dsp->hratio[1];  in TSAdaptRollBack_DSP()
63   dsp->hratio[1] = dsp->hratio[2];  in TSAdaptRollBack_DSP()
64   dsp->hratio[2] = 1.0;  in TSAdaptRollBack_DSP()
70   TSAdapt_DSP *dsp   = (TSAdapt_DSP *)adapt->data;  in TSAdaptChoose_DSP()  local
111   if (!ts->steprollback) dsp->rollback = PETSC_FALSE;  in TSAdaptChoose_DSP()
112   else if (!dsp->rollback) {  in TSAdaptChoose_DSP()
113     dsp->rollback = PETSC_TRUE;  in TSAdaptChoose_DSP()
121     PetscReal b1 = dsp->kBeta[0];  in TSAdaptChoose_DSP()
122     PetscReal b2 = dsp->kBeta[1];  in TSAdaptChoose_DSP()
123     PetscReal b3 = dsp->kBeta[2];  in TSAdaptChoose_DSP()
124     PetscReal a2 = dsp->Alpha[0];  in TSAdaptChoose_DSP()
125     PetscReal a3 = dsp->Alpha[1];  in TSAdaptChoose_DSP()
128     PetscReal ctr1 = dsp->cerror[0];  in TSAdaptChoose_DSP()
129     PetscReal ctr2 = dsp->cerror[1];  in TSAdaptChoose_DSP()
131     PetscReal rho1 = dsp->hratio[0];  in TSAdaptChoose_DSP()
132     PetscReal rho2 = dsp->hratio[1];  in TSAdaptChoose_DSP()
152       dsp->cerror[2] = dsp->cerror[1];  in TSAdaptChoose_DSP()
153       dsp->cerror[1] = dsp->cerror[0];  in TSAdaptChoose_DSP()
154       dsp->cerror[0] = ctr0;  in TSAdaptChoose_DSP()
155       dsp->hratio[2] = dsp->hratio[1];  in TSAdaptChoose_DSP()
156       dsp->hratio[1] = dsp->hratio[0];  in TSAdaptChoose_DSP()
157       dsp->hratio[0] = rho0;  in TSAdaptChoose_DSP()
158       dsp->rollback  = PETSC_FALSE;  in TSAdaptChoose_DSP()
181   TSAdapt_DSP *dsp = (TSAdapt_DSP *)adapt->data;  in TSAdaptView_DSP()  local
187     double a2 = (double)dsp->Alpha[0], a3 = (double)dsp->Alpha[1];  in TSAdaptView_DSP()
188     double b1 = (double)dsp->kBeta[0], b2 = (double)dsp->kBeta[1], b3 = (double)dsp->kBeta[2];  in TSAdaptView_DSP()
227   TSAdapt_DSP      *dsp = (TSAdapt_DSP *)adapt->data;  in TSAdaptDSPSetFilter_DSP()  local
241   dsp->kBeta[0] = tab->kBeta[0] / tab->scale;  in TSAdaptDSPSetFilter_DSP()
242   dsp->kBeta[1] = tab->kBeta[1] / tab->scale;  in TSAdaptDSPSetFilter_DSP()
243   dsp->kBeta[2] = tab->kBeta[2] / tab->scale;  in TSAdaptDSPSetFilter_DSP()
244   dsp->Alpha[0] = tab->Alpha[0] / tab->scale;  in TSAdaptDSPSetFilter_DSP()
245   dsp->Alpha[1] = tab->Alpha[1] / tab->scale;  in TSAdaptDSPSetFilter_DSP()
251   TSAdapt_DSP *dsp = (TSAdapt_DSP *)adapt->data;  in TSAdaptDSPSetPID_DSP()  local
254   dsp->kBeta[0] = kkI + kkP + kkD;  in TSAdaptDSPSetPID_DSP()
255   dsp->kBeta[1] = -(kkP + 2 * kkD);  in TSAdaptDSPSetPID_DSP()
256   dsp->kBeta[2] = kkD;  in TSAdaptDSPSetPID_DSP()
257   dsp->Alpha[0] = 0;  in TSAdaptDSPSetPID_DSP()
258   dsp->Alpha[1] = 0;  in TSAdaptDSPSetPID_DSP()
264   TSAdapt_DSP *dsp = (TSAdapt_DSP *)adapt->data;  in TSAdaptSetFromOptions_DSP()  local
283 …PetscCall(PetscOptionsRealArray("-ts_adapt_dsp_kbeta", "Filter parameters", "", dsp->kBeta, (n = 3…  in TSAdaptSetFromOptions_DSP()
286     for (i = n; i < 3; i++) dsp->kBeta[i] = 0;  in TSAdaptSetFromOptions_DSP()
288 …PetscCall(PetscOptionsRealArray("-ts_adapt_dsp_alpha", "Filter parameters", "", dsp->Alpha, (n = 2…  in TSAdaptSetFromOptions_DSP()
291     for (i = n; i < 2; i++) dsp->Alpha[i] = 0;  in TSAdaptSetFromOptions_DSP()
373   TSAdapt_DSP *dsp;  in TSAdaptCreate_DSP()  local
376   PetscCall(PetscNew(&dsp));  in TSAdaptCreate_DSP()
379   adapt->data                = (void *)dsp;  in TSAdaptCreate_DSP()