compute_nu_jets#
- t3toolbox.backend.sampling_derivatives.compute_nu_jets(tt_cores, xi_jets, trs)#
def compute_nu_jets( tt_cores: typ.Sequence[NDArray], # len=d, elm_shape=C+(rLi,nUi,rR(i+1)) xi_jets: typ.Sequence[NDArray], # input jets, len=d, elm_shape=(2,)+W+C+(nUi,) trs: NDArray, # binomial tensor -- only its shape (order) is read here ) -> typ.Tuple[NDArray, ...]: # nu_jets. len=d, elm_shape=(order+1,)+W+C+(rR(i+1),). nu_jets[i][t]=nu_i^(t)
Right derivative-pushthrough jets (standard recurrence form).
The mirror image of
compute_mu_jets(): reverse the tensor train (tt_reverseswaps bonds and core order), run the left banded-recurrence sweep, reverse the result.nu_jets[i]is the right edge variable entering corei(nu_i), stacked over derivative orders. Equal to the densecompute_nu_jets_trs()to tolerance; see it for the binomial-tensor reference form.