The HilbertTwistReduction(P, A, B) calling sequence returns the image of P under the Hilbert twist reduction.  If the (optional) third argument B is present, it is assigned to the Ore ring whose automorphism is the same as the A's and whose (pseudo) derivation sends everything to zero.

The InverseOfHilbertTwistReduction(P, A) calling sequence returns the pre-image of P under the Hilbert twist reduction. Note that A is the source ring of the Hilbert twist reduction.

Let A be the shift, q-shift, or differential algebra. The AccurateIntegration(L, A) calling sequence performs accurate integration, which solves the following problem: Let y satisfy L(y)=0 and g satisfy lambda(g)=y, where lambda means the usual derivative in the differential case, the difference operator in the shift case, and the q-difference operator in the q-shift case. The function builds an annihilator S (represented as an OrePoly structure) for g of the same degree as that of L, and an operator K such that g=K(y) if both exist. Otherwise, it returns .

A := SetOreRing(n, 'difference',
  'sigma' = proc(p, x) eval(p, x=x+1) end,
  'sigma_inverse' = proc(p, x) eval(p, x=x-1) end,
  'delta' = proc(p, x) eval(p, x=x+1) - p end,
  'theta1' = 0);

Abramov, S.A., and van Hoeij, M. "Integration of Solutions of Linear Functional Equations." Integral Transformations and Special Functions. Vol. 8 No. 1-2. (1999): 3-12.