Home
last modified time | relevance | path

Searched refs:yk (Results 1 – 5 of 5) sorted by relevance

/petsc/src/tao/bound/impls/bncg/
H A Dbncg.c131 if (!cg->yk) PetscCall(VecDuplicate(tao->gradient, &cg->yk)); in TaoSetUp_BNCG()
141 PetscCall(MatLMVMAllocate(cg->B, cg->sk, cg->yk)); in TaoSetUp_BNCG()
162 PetscCall(VecDestroy(&cg->yk)); in TaoDestroy_BNCG()
448 PetscCall(VecWAXPY(cg->yk, -1.0, cg->G_old, tao->gradient)); in TaoBNCGStepDirectionUpdate()
449 PetscCall(VecNorm(cg->yk, NORM_2, &ynorm)); in TaoBNCGStepDirectionUpdate()
451 PetscCall(VecDot(cg->yk, tao->stepdirection, &dk_yk)); in TaoBNCGStepDirectionUpdate()
519 PetscCall(VecDot(cg->yk, tao->gradient, &gkp1_yk)); in TaoBNCGStepDirectionUpdate()
525 PetscCall(VecDot(cg->yk, cg->g_work, &gkp1_yk)); in TaoBNCGStepDirectionUpdate()
533 PetscCall(VecWAXPY(cg->yk, -1.0, cg->G_old, tao->gradient)); in TaoBNCGStepDirectionUpdate()
534 PetscCall(VecNorm(cg->yk, NORM_2, &ynorm)); in TaoBNCGStepDirectionUpdate()
[all …]
H A Dbncg.h14 Vec sk, yk; member
/petsc/src/tao/unconstrained/impls/bmrm/
H A Dbmrm.h32 PetscReal *g, *y, *tempv, *d, *Qd, *t, *xplus, *tplus, *sk, *yk; member
H A Dbmrm.c17 PetscReal *xplus = df->xplus, *tplus = df->tplus, *sk = df->sk, *yk = df->yk; in solve() local
157 yk[i] = tplus[i] - t[i]; in solve()
181 bk += sk[i] * yk[i]; in solve()
541 PetscCall(PetscMalloc1(n, &df->yk)); in init_df_solver()
627 PetscCall(PetscFree(df->yk)); in ensure_df_space()
628 PetscCall(PetscMalloc1(n, &df->yk)); in ensure_df_space()
666 PetscCall(PetscFree(df->yk)); in destroy_df_solver()
/petsc/src/dm/dt/interface/
H A Ddt.c2417 PetscReal yk; /* Quadrature point 1 - x_k on reference domain [-1, 1] */ in PetscDTTanhSinhIntegrate() local
2443 …yk = 1.0 / (PetscExpReal(0.5 * PETSC_PI * PetscSinhReal(k * h)) * PetscCoshReal(0.5 * PETSC_PI * P… in PetscDTTanhSinhIntegrate()
2444 lx = -alpha * (1.0 - yk) + beta; in PetscDTTanhSinhIntegrate()
2445 rx = alpha * (1.0 - yk) + beta; in PetscDTTanhSinhIntegrate()
2482 mpfr_t yk; /* Quadrature point 1 - x_k on reference domain [-1, 1] */ in PetscDTTanhSinhIntegrateMPFR() local
2492 … * PetscLogReal(10.) / PetscLogReal(2.)), alpha, beta, h, sum, osum, psum, yk, wk, lx, rx, tmp, ma… in PetscDTTanhSinhIntegrateMPFR()
2533 mpfr_set_d(yk, 1.0, MPFR_RNDZ); in PetscDTTanhSinhIntegrateMPFR()
2535 mpfr_div(yk, yk, tmp, MPFR_RNDZ); in PetscDTTanhSinhIntegrateMPFR()
2537 mpfr_div(yk, yk, tmp, MPFR_RNDZ); in PetscDTTanhSinhIntegrateMPFR()
2539 mpfr_sub_d(lx, yk, 1.0, MPFR_RNDZ); in PetscDTTanhSinhIntegrateMPFR()
[all …]