conatus.expert.hill_from_plan¶
- conatus.expert.hill_from_plan(build, base, *, Nh, omega, label, group=None, linearity_tol=1e-12)[source]¶
Hill operator of an IR eigenvalue problem whose base flow is T-periodic.
build(ncc)must return the plan (discretize(problem, ncc=ncc)) of an EVP in the convention (L + σ M) x = 0, σ being the growth rate; the base-flow fields are NCC.base[field][m]gives the harmonic U_m of each base field (callable on the grid, like any NCC); a missing harmonic is zero. Harmonic m of L is L(U := U_m) - L(U := 0); linearity in U is checked (L(2U_m) - L(0) = 2 (L(U_m) - L(0))) and a non-linear dependence is refused (harmonic_matrices()). The diagonal block of harmonic n is L_0 + (μ + i n ω) M. An IVP plan works too: M X’ + L X = 0 is the same pencil.- Parameters:
build (Callable[[Mapping[Any, Any]], Any])
base (Mapping[Any, Mapping[int, Any]])
Nh (int)
omega (float)
label (str)
group (Any)
linearity_tol (float)
- Return type: