conatus.expert.find_relative_periodic_orbit

conatus.expert.find_relative_periodic_orbit(solver, x0, T0, ell0, newton, *, steps, reflect=False, log=None)

S τ_ℓ Φ_T(u) = u with unknowns (u, T, ℓ) and phase conditions ⟨v(u*), δu⟩ = 0 (time), ⟨∂x u*, δu⟩ = 0 (translation). reflect=True gives pre-periodic orbits R τ_ℓ Φ_T u = u (the ℓ freedom there is the choice of reflection point, pinned by the translation phase). solver: a Simulation or an IVPSolver (its scheme defines the discrete map); scope as find_equilibrium() (SCOPE).

Parameters:

• x0 (ndarray)

• T0 (float)

• ell0 (float)

• newton (NewtonKrylovHookstep)

• steps (int)

• reflect (bool)

Return type:

ECSSolution