Correction to “On a Class of Hermite-Obreshkov One-Step Methods with Continuous Spline Extension” [Axioms 7(3), 58, 2018]

: The authors of the above mentioned paper specify that the considered class of one-step symmetric Hermite-Obreshkov methods satisﬁes the property of conjugate-symplecticity up to order p + r , where r = 2 and p is the order of the method. This generalization of conjugate-symplecticity states that the methods conserve quadratic ﬁrst integrals and the Hamiltonian function over time intervals of length O ( h − r ) . Theorem 1 of the above mentioned paper is then replaced by a new one. All the other results in the paper do not change. Two new ﬁgures related to the already considered Kepler problem are also added.


Introduction
In paper [1] we analyzed the numerical solution of the first order Ordinary Differential Equation (ODE), y (t) = f(y(t)), t ∈ [t 0 , t 0 + T], associated with the initial condition: y(t 0 ) = y 0 , where f : IR m → IR m , m ≥ 1, is a C R−1 , R ≥ 1, function on its domain and y 0 ∈ IR m is assigned.
In particular, we considered the numerical solution of Hamiltonian problems which in canonical form can be written as follows: y = J∇H(y), with y = q p , q, p ∈ IR , where q and p are the generalized coordinates and momenta, H : IR 2 → IR is the Hamiltonian function and I stands for the identity matrix of dimension . Note that the flow ϕ t : y 0 → y(t) associated with the dynamical system (3) is symplectic; this means that its Jacobian satisfies: ∂ϕ t (y) ∂y J ∂ϕ t (y) ∂y = J, ∀ y ∈ IR 2 .
We recall that a one-step numerical method Φ h : IR 2 → IR 2 with stepsize h is symplectic if the discrete flow y n+1 = Φ h (y n ), n ≥ 0, satisfies: Two numerical methods Φ h , Ψ h are conjugate to each other if there exists a global change of coordinates χ h , such that: uniformly for y varying in a compact set and • denoting a composition operator [2]. A method which is conjugate to a symplectic method is said to be conjugate symplectic, this is a less strong requirement than symplecticity, which allows the numerical solution to have the same long-time behavior of a symplectic method. A more relaxed property, shared by a wider class of numerical schemes, is a generalization of the conjugate-symplecticity property, introduced in [3]. A method y 1 = Ψ h (y 0 ) of order p is conjugate-symplectic up to order p + r, with r ≥ 0, if a global change of coordinates A consequence of property (7) is that the method Ψ h (y) nearly conserves all quadratic first integrals and the Hamiltonian function over time intervals of length O(h −r ) (see [3]).
Recently, the class of Euler-Maclaurin Hermite-Obreshkov (EMHO) methods for the solution of Hamiltonian problems has been analyzed in [4] where the conjugate symplecticity up to order p + 2 of the p-th order methods was proven. In this paper, we fix Theorem 1 of [1] related to symmetric one-step BS Hermite-Obreshkov (BSHO) methods, proving that the conjugate-symplecticity property is satisfied by the R-th one-step symmetric Hermite-Obreshkov method up to order 2R + 2.
Let t i , i = 0, . . . , N, be an assigned partition of the integration interval [t 0 , t 0 + T], and let us denote by u i an approximation of y(t i ). We consider one-step symmetric BSHO method as follows, setting u 0 := y 0 , where h n := t n+1 − t n , β where D 0 = I is the identity operator and D k f(z) is defined as the k-th total derivative of f(y(t)) computed at y(t) = z, where for the computation of the total derivative it is assumed that y satisfies the differential equation in (1). Thus for example D 1 f(z) = ∂f ∂y (z) f(z), where ∂f ∂y is the m × m Jacobian matrix of f. Note that we use the subscript to define the Lie operator to avoid confusion with the same order classical derivative operator in the following denoted as D k . With this clarification on the definition of u (j) i , following the lines of the proof given in [4], we can actually prove that the R-th one-step symmetric BSHO method is conjugate symplectic up to order 2R + 2.
We show that the map Theorem 1. The map u 1 = Ψ h (u 0 ) associated with the one-step method (8) admits a B-series expansion and is conjugate to a symplectic B-series integrator up to order 2R + 2.
Proof. The existence of a B-series expansion for y 1 = Ψ h (y 0 ) is directly deduced from [5], where a B-series representation of a generic multi-derivative Runge-Kutta method has been obtained. By defining the two characteristic polynomials of the trapezoidal rule: and the shift operator E(u n ) := u n+1 , the R-th method described in (8) reads, We now consider a function v(t), a stepsize h and the shift operator E h (v(t)) := v(t + h), and we look for a continuous function v(t) that satisfies (10) in the sense of formal series (a series where the number of terms is allowed to be infinite), using the relation By multiplying both sides of the previous equation by Dρ(e hD ) −1 , we obtain: Now, since Bernoulli numbers define the Taylor expansion of the function z/(e z − 1) and b 0 = 1, b 1 = −1/2 and b j = 0 for the other odd j, we have: Thus, we can write (11) aṡ Adding and subtracting terms involving the classical derivative operator D 2k , D 2k−1 , we geṫ that we recast asv Since v(t) = y(t) + O(h 2R ), due to the regularity conditions on the function f, we see that . . , R − 1 and hence the solution v(t) of (12) is O(h 2R+2 )-close to the solution of the following initial value probleṁ with: that has been derived from (12) by neglecting the sums containing the derivatives D 2k , D 2k−1 . Observe that δ j = 0 for j = 1, . . . , R − 1, since the method is of order 2R (see [2], Theorem 3.1, page 340). We may interpret (13) as the modified equation of a one-step method y 1 = Φ h (y 0 ), where Φ h is evidently the time-h flow associated with (13). Expanding the solution of (13) in Taylor series, we get the modified initial value differential equation associated with the numerical scheme by coupling (13) with the initial condition w(t 0 ) = y 0 . Thus, Φ h is a B-series integrators. The proof of the conjugated symplecticity of Φ h follows exactly the same steps of the analogous proof in Theorem 1 of [4]. Since Ψ h (y) = Φ h (y) + O(h 2R+3 ) and Φ h is conjugate-symplectic, the result follows using the same global change of coordinates χ h (y) associated to Φ h .
We report in Figure 1 the bottom-rigth picture of Figure 2 of [1], related to the Kepler problem, where we noticed that the error in the second component of the Lenz vector was not correctly computed, for completness Figure 1 also reports the error in the first component of the Lenz vector. To stress that the methods show a good long time behavior for Hamiltonian problems, we report also, in Figure 2 the results using a longer integration interval of 10 5 periods and all the other parameters unchanged. In the pictures we report the maximum error in each period for the Hamiltonian function, the angular mument and the Lenz vector. The results remain consistent, showing a linear grows in the error and in the Lenz vector and a near conservation of the Hamiltonian and of the angular moment.
All the other results in the paper do not change.