Next Article in Journal / Special Issue
Cohomology Theory of Nonassociative Algebras with Metagroup Relations
Previous Article in Journal
Generalized Hyers–Ulam Stability of the Additive Functional Equation
Previous Article in Special Issue
Unification Theories: New Results and Examples
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dual Numbers and Operational Umbral Methods

1
Institut de Recherche en Informatique Fondamentale (IRIF), Université de Paris, Bâtiment Sophie Germain, Case Courier 7014, 8 Place Aurélie Nemours, CEDEX 13, 75205 Paris, France
2
ENEA—Frascati Research Center, Via Enrico Fermi 45, 00044 Rome, Italy
3
H. Niewodniczański Institute of Nuclear Physics, Polish Academy of Science, ul. Radzikowskiego 152, 31-342 Kraków, Poland
*
Author to whom correspondence should be addressed.
Submission received: 22 May 2019 / Revised: 25 June 2019 / Accepted: 26 June 2019 / Published: 2 July 2019
(This article belongs to the Special Issue Non-associative Structures and Other Related Structures)

Abstract

:
Dual numbers and their higher-order version are important tools for numerical computations, and in particular for finite difference calculus. Based on the relevant algebraic rules and matrix realizations of dual numbers, we present a novel point of view, embedding dual numbers within a formalism reminiscent of operational umbral calculus.

1. Introduction

Dual numbers (DNs), introduced during the second half of the 19th century [1,2,3,4,5], can be viewed as abstract entities, similar to ordinary complex numbers, and are defined as
z = x + ϵ y , ( x , y ) R
where the corresponding “imaginary” unit or dual number unit (DNU) ϵ is a nilpotent number,
ϵ 2 = 0 a n d ϵ 0 .
Dual numbers were originally introduced within the context of geometrical studies, and later exploited to deal with problems in pure and applied mechanics [6,7]. For instance, it has been demonstrated in [8,9,10] how to formulate the equations of rigid body motion in terms of just three “dual” equations instead of their six “real” counterparts (thereby realizing an equivalence between spherical and spatial kinematics). This approach has been extended in [11,12,13] to a treatment of rigid body motion in terms of a certain variant of “hyper-dual” numbers, implemented in contrast to our approach via sets of “ordinary” dual numbers together with certain algebraic relations that are motivated from the specific requirements within the relevant field of robotics and of mechanics. More recently, as further discussed in the present paper, the importance of a different kind of higher-order dual numbers has been recognized in numerical analysis to reduce round-off errors [14]. We believe that the use of dual numbers in the applied sciences is not as widespread as it could be, and that many new fields of research would benefit from their relevant introduction. An important domain in which they may bring significant novelties is that of the perturbative techniques in classical and quantum mechanics.
The main contribution of this paper consists in fixing the underlying algebraic rules of the dual numbers in the wider context of umbral and operational calculus. The paper is organized as follows: Section 2 delivers a basic mathematical introduction to dual numbers. Section 3 is devoted to the description of the computational procedure based upon dual numbers and umbral calculus. In Section 4, we provide insight into how this powerful method can be applied to deal with problems arising in different contexts. For illustration, we consider the Schrödinger and the heat equation, cornerstones in their respective fields of physics. Section 5 provides a conclusion with further considerations for future works.

2. Higher-Order Dual Numbers

The DN algebraic rules [15,16], summarized below, are a straightforward consequence of the previous identity in Equation (2) (with z = x + ϵ y and w = u + ϵ v ):
Component - wise algebraic addition z + w = x + u + ϵ ( y + ν ) Product z · w = x u + ϵ ( x ν + y u ) Inverse z 1 = 1 x 1 ϵ y x ( x 0 ) Power z n = x n 1 + n ϵ y x ( n Z 0 , x 0 )
While the addition operation is entirely analogous to the component-wise addition operation on two-dimensional vectors, the last three operations (product, inverse and power) characterize the distinguishing special algebraic properties of dual numbers (DNs). The multiplication is commutative, associative and distributive, thus DNs form a two-dimensional associative and commutative algebra over the real numbers.
We now extend this traditional dual number formalism as motivated by the following type of problem. Consider the Taylor expansion up to some order k (denoted k ) of an at least k-fold continuously differentiable function f around a point x,
f ( x + y ) k m = 0 k y m m ! f ( m ) ( x ) .
Following the automatic differentiation paradigm [17,18,19], since in practice the function f will be implemented in some algorithmic from, it may be advantageous to formulate truncations such as Equation (4) in terms of generalized (or higher-order) dual numbers. To this end, let us introduce the families of square matrices k ϵ ^ ± , k 1 ^ and k 0 ^ with entries (for i , j = 1 , , k )
k ϵ ^ ± i , j : = δ j , i ± 1 , k 1 ^ i , j : = δ i , j , k 0 ^ i , j : = 0 ,
where δ i , j denotes the Kronecker symbol. It is straightforward to verify that for all k 2 and 0
k ϵ ^ ± i , j = δ j , i ± k ϵ ^ ± k = 0 .
Then, under the assumptions in Equation (4), endowing the function f ( x ) suitably with a component-wise action on square matrices, we find (for k 2 )
f x k 1 ^ + y k ϵ ^ ± = m = 0 k 1 1 m ! y m f ( m ) ( x ) k ϵ ^ ± m .
For example, setting k = 2 , reproduces the well-known dual number identity [19]
f x 2 1 ^ + y 2 ϵ ^ + = f ( x ) 2 1 ^ + y f ( x ) 2 ϵ ^ + = f ( x ) y f ( x ) 0 f ( x ) .
For illustration, setting k = 3 , we obtain
f x * [ 3 ] 1 ^ + y * [ 3 ] ϵ ^ + = f ( x ) 3 1 ^ + y f ( x ) 3 ϵ ^ + + 1 2 y 2 f ( x ) 3 ϵ ^ + 2 = f ( x ) y f ( x ) 1 2 y 2 f ( x ) 0 f ( x ) y f ( x ) 0 0 f ( x ) .
It may be verified that, e.g., for the choice “+” in Equation (5), the first row of the resulting matrices in Equation (7) contain the terms of the Taylor expansion up to order k 1 . More explicitly, introducing the auxiliary notations for the row vector k e 1 and the column vector k 1 of length k 2 ,
k e 1 : = ( 1 , 0 , , 0 ) , k 1 : = 1 1 ,
allows us to define the order k evaluation operation acting on some function F ( k ϵ ^ + ) (which is itself assumed to act entry-wise) depending on a generalized dual number k ϵ ^ + as
F ( k ϵ ^ + ) k : = k e 1 F ( k ϵ ^ + ) k 1 .
We thus find that
f x k 1 ^ + y k ϵ ^ + k = m = 0 k 1 y m m ! f ( m ) ( x ) .
Recently, expansions such as Equation (7) have received considerable interest in the field of numerical analysis [20]. Referring to Fike [19] for an overview, various alternative types of “numbers” have been studied for the purpose of finding optimized numerical schemes for computing kth-order derivatives of functions. For example, it has been demonstrated that the use of so-called hyper-dual numbers results in first- and second-derivative calculations that are exact, regardless of the step size [14].
For later convenience, motivated by the identity (for k 2 )
exp k ϵ ^ + x = r = 0 k 1 x r r ! k ϵ ^ + r ,
we may introduce the so-called truncated exponential polynomials [21] e n ( x ) defined through the series
e n ( x ) : = r = 0 n x r r ! ,
which may be expressed in terms of generalized dual numbers as
e n ( x ) e n 1 ( x ) e 1 ( x ) 1 : = exp * [ n + 1 ] ϵ ^ + x n + 1 1 .
One may thus easily verify (via Equations (13) and (14)) the property
e n ( x ) = e n 1 ( x ) .
Having provided a matrix-based extension of ordinary to kth-order dual numbers of arbitrary order k 2 , we now proceed to develop a computational procedure embedding dual numbers with other techniques inspired by the operational umbral formalism.

3. Umbral-Type Methods and Dual Numbers

Starting from this section, we employ the notational simplification of writing ϵ for the dual number unit (DNU) k ϵ ^ ± of generalized dual numbers (cf. Equation (5)), making the order k 2 of the DN explicit only via the analog of the notation in Equation (11), and masking the matrix nature of k ϵ ^ ± . Thus, for some function F F ( ϵ ) , we write
F k G : G = F | ϵ k + 1 0 | ϵ 1
for the truncation of F via setting ϵ k + 1 = 0 and afterwards ϵ = 1 . It is straightforward to verify that this formal definition may be implemented in terms of the matrix representations introduced in Section 2 via use of Equation (11) as
G = F k + 1 ϵ ^ + ) k + 1 .
Consider then the dual complex parameter
z ^ z ^ ( a , b ) : = a + ϵ b .
Following the principles of umbral calculus, we treat the dual complex parameter z ^ as an ordinary algebraic quantity in calculations of integrals, derivatives and other operations, delaying the evaluation of z ^ via performing the operation k to the very end of the computations. Note that albeit the term umbral calculus has been introduced in the seminal papers by Roman and Rota [22], in the following we make reference to the formalism developed in [23] which enriches the original formalism with the wealth of techniques derived from the operational calculus [22,23]. We now illustrate the computational benefits of this approach via a number of examples.

3.1. Dual Shifted Gaussians

We first consider a Gaussian-type function explicitly containing in its argument the dual complex parameter in Equation (19), whence the dual-shifted Gaussian function
f ( x ) = e α x 2 + z ^ ( a , b ) x .
Assuming for instance third-order dual numbers (i.e., ϵ 3 = 0 ), we may write the above function in more conventional terms as
f ( x ) 2 e α x 2 + a x 1 + b x + 1 2 ( b x ) 2 ,
which is easily recognized as the product of a shifted Gaussian with a second-degree polynomial.
In full analogy to the umbral operational methods of Licciardi [23], it is then straightforward to calculate the following integral of the function f of Equation (20) via the standard Gaussian integral formula
+ f ( x ) d x = π α e z ^ ( a , b ) 2 4 α = π α e a 2 4 α + a b 2 α ϵ + b 2 4 α ϵ 2 .
The term on the right has in fact a definite meaning, since the use of the generating function of the two-variable Hermite polynomials [24]
n = 0 t n n ! H n ( x , y ) = e x t + y t 2
H n ( x , y ) = e y x 2 x n = n ! r = 0 n 2 x n 2 r y r ( n 2 r ) ! r !
permits to cast the right-hand side of Equation (22) into the form
π α e a 2 4 α + a b 2 α ϵ + b 2 4 α ϵ 2 = π α e a 2 4 α m 0 ϵ m m ! H m a b 2 α , b 2 4 α k π α e a 2 4 α H e k a b 2 α , b 2 4 α .
Here, H e k ( x , y ) denotes the Hermite-based truncated exponential polynomial [25,26,27] defined as
H e k ( x , y ) : = r = 0 k 1 r ! H r x , y .

3.2. Another Form of Dual Gaussian

Let us consider as a further example
g ( x ) : = e z ^ ( a , b ) x 2
and the following infinite integral (for Re ( a ) > 0 )
+ g ( x ) d x = π z ^ = π a + ϵ b k π a r = 0 k 1 2 r b a r .
Here, by invoking the operation k , we obtain a finite series, thus obviating the need to impose any condition on the relevant convergence range.

3.3. Examples From Symbolic Calculus

The calculus of higher-order dual numbers may be further refined via combining it with the wealth of techniques available from the theory of special functions and symbolic calculus as put forward in [23,28,29,30,31,32]. Consider for illustration the following identity, known from the theory of two-variable Hermite polynomials [33],
x n e α x 2 = H n ( 2 α x , α ) e α x 2
which allows to simplify the task of calculating successive derivatives of the dual Gaussian introduced in Equation (26), such as in the computation
x n e z ^ x 2 = H n ( 2 z ^ x , z ^ ) e z ^ x 2 = n ! r = 0 n 2 s 0 ( 1 ) n r + s 2 n 2 r x n 2 ( r s ) ( n 2 r ) ! r ! s ! z ^ n r + s .
Here, the first step follows from (28) and the second by invoking (23b).
Another interesting type of calculus concerns infinite integrals involving rational functions such as
Φ ( x ; a , b ) : = 1 1 + z ^ x 2 k 1 1 + a x 2 r = 0 k b x 2 1 + a x 2 r .
For example, the infinite integral
+ 1 1 + z ^ x 2 d x = π z ^
may be easily transformed into truncated form in full analogy to the calculation summarized in Equation (27).

3.4. Umbral Image Type Techniques

Referring to Behr et al. [34] for the precise technical details (compare also [33]), suffice it here to provide the following definition for the action of the formal integration operator I ^ on the formal variable v (for α C ):
I ^ ( v α ) : = 1 Γ ( α ) .
Then, an interesting variant of the example presented in Equation (29) may be obtained as
I ^ + v Φ ( x ; a , v b ) d x = I ^ v π z ( a , v β ) k π a r = 0 k 1 Γ ( 1 2 r ) ( r ! ) 2 b a r .
In summary, the combination of the concept of higher order dual numbers with techniques from symbolic and umbral-image type calculus appears to offer a large potential in view of novel tools of computation. To corroborate this claim, we now present some first high-level results in this direction.

4. Dual Numbers and Solution of Heat- and Schrödinger-Type Equations

Before entering the main topic of this section, let us recall a few useful “operational rules”, starting with the Glaisher identity [27,35]
e τ d 2 d x 2 e α x 2 = 1 1 + 4 τ α e α x 2 1 + 4 τ α ,
which can also be understood as the solution of the heat equation with a Gaussian as initial function. It proves particularly useful in the following to note that, according to the definition of the Hermite polynomials H n ( x , y ) as given in Equation (23b), an alternative interpretation of Equation (33) is provided in terms of the double-lacunary exponential generating function H 2 , 0 ( λ ; x , y ) of the polynomials H n ( x , y ) , where we employ notations as in [36]
e τ d 2 d x 2 e α x 2 = n 0 ( α ) n n ! H 2 n ( x , τ ) = H 2 , 0 ( α ; x , τ ) .
By specializing Equation (33) to α = z ^ (with z ^ = a + ϵ b the dual complex parameter in Equation (19)), we obtain the operational identity
e τ d 2 d x 2 e z ^ x 2 = 1 1 + 4 τ z ^ e z ^ x 2 1 + 4 τ z ^ .
Via the simple factorizations
1 + 4 z ^ τ = γ ( a , τ ) γ b ϵ γ ( a , τ ) , τ , z ^ 1 + 4 z ^ τ = a γ ( a , τ ) + b ϵ [ γ ( a , τ ) ] 2 γ b ϵ γ ( a , τ ) , τ , γ ( c , τ ) = 1 + 4 c τ ,
we may transform the identity in Equation (35) as
e τ d 2 d x 2 e z ^ x 2 = H 2 , 0 b ϵ [ γ ( a , τ ) ] 2 ; x , τ γ ( a , τ ) H 2 , 0 ( a ; x , τ ) .
By re-inserting the definition of the first double-lacunary exponential generating function, using the Glaisher-identity in Equation (33) for the second one and finally truncating to order k, we eventually arrive at the compact result
e τ d 2 d x 2 e z ^ x 2 k e α x 2 γ ( a , τ ) γ ( a , τ ) n = 0 k 1 n ! b [ γ ( a , τ ) ] 2 n H 2 n x , τ γ ( a , τ ) , γ ( a , τ ) = 1 + 4 a τ .
For example, by evaluating the above expression for second order dual numbers, one finds
e τ d 2 d x 2 e z ^ x 2 2 e a x 2 γ ( a , τ ) γ ( a , τ ) 1 b γ ( a , τ ) 2 H 2 ( x , τ γ ( a , τ ) ) + b 2 2 γ ( a , τ ) 4 H 4 ( x , τ γ ( a , τ ) ) .
The above result may be interpreted as the solution of the heat-type equation
τ F ( x , τ ) = x 2 F ( x , τ ) , F ( x , 0 ) = e z ^ x 2 .
An analogous problem has been addressed in [33] within the framework of a different method. The techniques we have envisaged may be further exploited to treat the paraxial propagation of the so-called flattened distributions, introduced in [37] to study the laser field evolution in optical cavities employing super-Gaussian mirrors [38]. These cavities shape beams whose transverse distribution is not reproduced by a simple Gaussian, but by a function exhibiting a quasi-constant flat-top, expressible through a function of the type
E ( x ; p ) : = e | x | p , p Z > 0 .
The paraxial propagation of these beams has less obvious properties than, e.g., Laguerre or Hermite Gauss modes [38]. To overcome this drawback, Gori introduced the so-called flattened beams [37], which permit a fairly natural expansion in terms of Gauss Laguerre/Hermite modes, thus providing a straightforward solution to the corresponding paraxial wave equation.
Invoking our formalism as developed thus far, we may approximate the aforementioned Gori beams in the form
E ( x ; p ) Y ( x ; α | m ) : = e α x 2 e m ( x 2 ) .
Here, e m ( x ) denotes the truncated exponential polynomials introduced in Equation (14), and both parameters α and m depend on p (see [33] for further details). Recalling from Equation (19) the definition z ^ ( a , b ) : = a + b ϵ of the dual complex parameter, the r.h.s. of Equation (42) may be equivalently expressed as
e z ^ ( α , 1 ) x 2 m Y ( x ; α | m ) ,
whence as an instance of a dual Gaussian as described in Section 3.2. The problem of the relevant propagation can accordingly be reduced to that of an ordinary Gaussian mode, namely to the solution of the Schrödinger type equation
i τ Ψ ( x , τ ) = x 2 Ψ ( x , τ ) , Ψ ( x , 0 ) = Y ( x ; α | m ) .
Consequently, by invoking the operational identity in Equation (35), the paraxial evolution of a flattened beam may be expressed in the form
Ψ ( x , τ ) = e i τ x 2 e z ^ ( α , 1 ) x 2 = 1 1 + 4 i τ z ^ ( α , 1 ) e z ^ ( α , 1 ) x 2 1 + 4 i τ z ^ ( α , 1 ) ,
which reproduces indeed the known solution of our problem (compare [33]).
In a forthcoming paper, we will discuss this specific application in further detail by applying the method to the problem of designing super-Gaussian optical systems.

5. Weyl Formula and Modified Hermite Polynomials

The wide flexibility of the method we propose is corroborated by the following further example, relevant to the use of operational ordering tools. Let us consider an evolution equation of the form
τ F ( x , τ ) = [ γ x z ^ x ] F ( x , τ ) , F ( x , 0 ) = f ( x ) .
The relevant procedure for combining differential calculus with the umbral formalism is described in [39]. Following this approach, the solution of Equation (46) can be expressed as
F ( x , τ ) = e τ ( γ x z ^ x ) f ( x ) .
To evaluate the solution of Equation (47) explicitly, we need to suitably “factorize” the exponential operator. This so-called disentanglement operation may be implemented via the Weyl formula [40]
e X ^ + Y ^ = e 1 2 [ X ^ , Y ^ ] e X ^ e Y ^ ,
which is applicable whenever the identities [ X ^ , [ X ^ , Y ^ ] ] = [ [ X ^ , Y ^ ] , Y ^ ] ] = 0 hold. Applying the Weyl formula in Equation (48) to our solution in Equation (47), if we let X ^ = τ z ^ x and Y ^ = τ γ x (resulting in [ X ^ , Y ^ ] = τ 2 γ z ^ , and with higher-order commutators vanishing), we obtain
F ( x , τ ) = e 1 2 τ 2 γ z ^ e z ^ x τ f ( x + γ τ ) .
Thus, the solution at any desired truncation order k may be obtained by invoking the dual number evaluation operation k of Equation (17).
As already mentioned above, the Weyl formula applies in the example presented because the algebraic structure of the argument of the exponential in Equation (47) satisfies a special property: the commutators of the associated generators reduce to a constant after the first commutation bracket. A more interesting extension is given by the case in which the generators are embedded into a solvable Lie algebra. In this case, the combined use of the dual number formalism and of the Wei-Norman ordering method [41] leads to new and interesting results. They deserve a separate treatment that will be reported in a forthcoming paper.
As a final example, we define modified Hermite polynomials H n ( x , z ^ ) , whence ordinary two-variable Hermite polynomials H n ( x , y ) as introduced in Equation (23b) evaluated at y = z ^ , with z ^ z ^ ( a , b ) the dual complex parameter of Equation (19),
H n ( x , z ^ ) = e z ^ x 2 x n .
It is straightforward to verify that these modified polynomials inherit all the relevant properties from the polynomials H n ( x , y ) , such as the recurrences
x H n ( x , z ^ ) = n H n 1 ( x , z ^ ) , H n + 1 ( x , z ^ ) = x H n ( x , z ^ ) + 2 z ^ x H n ( x , z ^ ) ,
and we find that they satisfy the second order differential equation
2 z ^ x 2 H n ( x , z ^ ) + x x H n ( x , z ^ ) = n H n ( x , z ^ ) .
The explicit form of these truncated polynomials is easily obtained. For example, by using third order dual numbers, which implies
e z ^ x 2 2 e a x 2 1 + b x 2 + 1 2 b 2 x 4 ,
we find the explicit formula
H n ( x , z ^ ) 2 H n ( x , a ) + b a H n ( x , a ) + 1 2 b 2 a 2 H n ( x , a ) ,
where we have invoked the well-known identity
x 2 H n ( x , y ) = y H n ( x , y ) .

6. Final Comments

The method we have outlined in this paper offers many computational advantages to treat problems where truncated expansions (not necessarily of Taylor type) of functions are involved. At its core, the umbral formalism and the notion of higher order dual numbers allow delaying the explicit expansions to later stages in a given calculation, thus opening the possibility to exploit numerous efficient computation strategies from the theory of operational calculus and special functions.
The technique we have introduced in this paper is amenable for new applications in various different fields. We have presented herein the solution of parabolic equations in transport problems, and within such a context a fairly important example has been provided by treating the propagation of flattened beams [33,37] in optics. For brevity, we have just outlined the procedure in terms of a one-dimensional computation. The relevant extension to the three-dimensional case does not require any particular conceptual effort, but only a consistent numerical implementation. In a forthcoming investigation, we will further extend the method and study its potential for treating perturbative problems in classical and quantum mechanics.

Author Contributions

Conceptualization, G.D.; methodology, N.B. and G.D.; validation, N.B., G.D. and S.L.; formal analysis, N.B., G.D., A.L., and S.L.; writing—original draft preparation, N.B., G.D., and A.L.; and writing—review and editing, N.B., A.L., and S.L.

Funding

The work of N.B. is supported by funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No. 753750. A.L. was supported by the NCN research project OPUS 12 No. UMO-2016/23/B/ST3/01714 and by the NAWA project: Program im. Iwanowskiej PPN/IWA/2018/1/00098. S.L. was supported by an Enea-Research Center Individual Fellowship.

Acknowledgments

N.B. would like to thank the LPTMC (Paris 06) and ENEA Frascati for warm hospitality.

Conflicts of Interest

The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

References

  1. Clifford, W.K. Preliminary Sketch of Biquaternions. Proc. Lond. Math. Soc. 1871, s1–s4, 381–395. [Google Scholar] [CrossRef]
  2. Grünwald, J. Über duale Zahlen und ihre Anwendung in der Geometrie. Monatshefte für Mathematik Physik 1906, 17, 81–136. [Google Scholar] [CrossRef]
  3. Segre, C. Le Geometrie Projettive nei Campi di Numeri Duali; Vincenzo Bona: Torino, Italy, 1912. [Google Scholar]
  4. Yaglom, I.M. A Simple Non-Euclidean Geometry and Its Physical Basis: An Elementary Account of Galilean Geometry and the Galilean Principle of Relativity; Springer: Berlin, Germany, 2012. [Google Scholar]
  5. Klawitter, D. Clifford Algebras: Geometric Modelling and Chain Geometries with Application in Kinematics; Springer: Berlin, Germany, 2014. [Google Scholar]
  6. Kotelnikov, A.P. Screw Calculus and Some Applications to Geometry and Mechanics. Ann. Imp. Univ. Kazan 1895, 24. (In Russian) [Google Scholar]
  7. Study, E. Die Geometrie der Dynamen. Jahresbericht der Deutschen Mathematiker-Vereinigung 1900, 8, 204–216. [Google Scholar]
  8. Rooney, J. On the principle of transference. In Proceedings of the Fourth World Congress on the Theory of Machines and Mechanisms; Institution of Mechanical Engineers: London, UK, 1975; pp. 1088–1092. [Google Scholar]
  9. Hsia, L.; Yang, A. On the principle of transference in three-dimensional kinematics. J. Mech. Des. 1981, 103, 652–656. [Google Scholar] [CrossRef]
  10. Martínez, J.M.R.; Duffy, J. The principle of transference: History, statement and proof. Mech. Mach. Theory 1993, 28, 165–177. [Google Scholar] [CrossRef]
  11. Cohen, A.; Shoham, M. Application of Hyper-Dual Numbers to Multibody Kinematics. J. Mech. Rob. 2015, 8, 011015. [Google Scholar] [CrossRef]
  12. Cohen, A.; Shoham, M. Application of hyper-dual numbers to rigid bodies equations of motion. Mech. Mach. Theory 2017, 111, 76–84. [Google Scholar] [CrossRef]
  13. Cohen, A.; Shoham, M. Principle of transference—An extension to hyper-dual numbers. Mech. Mach. Theory 2018, 125, 101–110. [Google Scholar] [CrossRef]
  14. Fike, J.; Alonso, J. The Development of Hyper-Dual Numbers for Exact Second-Derivative Calculations. In Proceedings of the 49th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition; American Institute of Aeronautics and Astronautics: Reston, VA, USA, 2011. [Google Scholar] [CrossRef] [Green Version]
  15. Harkin, A.A.; Harkin, J.B. Geometry of Generalized Complex Numbers. Math. Mag. 2004, 77, 118–129. [Google Scholar] [CrossRef]
  16. Özdemir, M. Introduction to Hybrid Numbers. Adv. Appl. Clifford Algebr. 2018, 28. [Google Scholar] [CrossRef]
  17. Rall, L.B. (Ed.) Automatic Differentiation: Techniques and Applications; Springer: Berlin/Heidelberg, Germany, 1981. [Google Scholar] [CrossRef]
  18. Rall, L.B.; Corliss, G.F. An introduction to automatic differentiation. In Computational Differentiation: Techniques, Applications, and Tools; SIAM: Philadelphia, PA, USA, 1996; Volume 89. [Google Scholar]
  19. Fike, J.A. Derivative Calculations Using Hyper-Dual Numbers; Technical Report; Sandia National Lab. (SNL-NM): Albuquerque, NM, USA, 2016. [Google Scholar]
  20. Berland, H.; (Department of Mathematical Sciences, NTNU, Trondheim, Norway). Personal communication, 2006.
  21. Dattoli, G.; Cesarano, C.; Sacchetti, D. A note on truncated polynomials. Appl. Math. Comput. 2003, 134, 595–605. [Google Scholar] [CrossRef]
  22. Roman, S.M.; Rota, G.C. The umbral calculus. Adv. Math. 1978, 27, 95–188. [Google Scholar] [CrossRef] [Green Version]
  23. Licciardi, S. Umbral Calculus, a Different Mathematical Language. arXiv 2018, arXiv:1803.03108. [Google Scholar]
  24. Appel, P.; De Feriet, J.K. Fonctions hypergéométriques et hypersphériques. In Polynômes d’Hermite; Gauthier-Villars: Paris, France, 1926. [Google Scholar]
  25. Dattoli, G. Hermite-Bessel and Laguerre-Bessel functions: A by-product of the monomiality principle. In Proceedings of the Melfi School on Advanced Topics in Mathematics and Physics, Advanced Special Functions and Applications; Aracne Editrice: Melfi, Italy, 2000; pp. 147–164. [Google Scholar]
  26. Dattoli, G. Generalized polynomials, operational identities and their applications. J. Comput. Appl. Math. 2000, 118, 111–123. [Google Scholar] [CrossRef] [Green Version]
  27. Dattoli, G.; Khan, S.; Ricci, P. On Crofton—Glaisher type relations and derivation of generating functions for Hermite polynomials including the multi-index case. Integral Transf. Spec. Funct. 2008, 19, 1–9. [Google Scholar] [CrossRef]
  28. Babusci, D.; Dattoli, G.; Górska, K.; Penson, K. Repeated derivatives of composite functions and generalizations of the Leibniz rule. Appl. Math. Comput. 2014, 241, 193–199. [Google Scholar] [CrossRef] [Green Version]
  29. Babusci, D.; Dattoli, G. On Ramanujan Master Theorem. arXiv 2011, arXiv:1103.3947. [Google Scholar]
  30. Górska, K.; Babusci, D.; Dattoli, G.; Duchamp, G.; Penson, K. The Ramanujan master theorem and its implications for special functions. Appl. Math. Comput. 2012, 218, 11466–11471. [Google Scholar] [CrossRef] [Green Version]
  31. Dattoli, G.; Di Palma, E.; Sabia, E.; Gorska, K.; Horzela, A.; Penson, K. Operational versus umbral methods and the Borel transform. Int. J. Appl. Comput. Math. 2017, 3, 3489–3510. [Google Scholar] [CrossRef]
  32. Babusci, D.; Dattoli, G.; Górska, K.; Penson, K. Lacunary generating functions for the Laguerre polynomials. Séminaire Lotharingien Combinatoire 2017, 76, B76b. [Google Scholar]
  33. Dattoli, G.; Migliorati, M. The truncated exponential polynomials, the associated Hermite forms and applications. Int. J. Math. Math. Sci. 2006, 2006, 98175-1. [Google Scholar] [CrossRef]
  34. Behr, N.; Dattoli, G.; Duchamp, G.; Licciardi, S.; Penson, K. Operational Methods in the Study of Sobolev-Jacobi Polynomials. Mathematics 2019, 7, 124. [Google Scholar] [CrossRef]
  35. Crofton, M. On operative symbols in the differential calculus. Proc. Lond. Math. Soc. 1880, 1, 122–134. [Google Scholar] [CrossRef]
  36. Behr, N.; Duchamp, G.H.; Penson, K.A. Explicit formulae for all higher order exponential lacunary generating functions of hermite polynomials. arXiv 2018, arXiv:1806.08417. [Google Scholar]
  37. Gori, F. Flattened gaussian beams. Opt. Commun. 1994, 107, 335–341. [Google Scholar] [CrossRef]
  38. Siegman, A.E. Lasers university science books. Mill Valley CA 1986, 37, 208. [Google Scholar]
  39. Babusci, D.; Dattoli, G. Umbral methods and operator ordering. arXiv 2011, arXiv:1112.1570. [Google Scholar]
  40. Dattoli, G.; Ottaviani, P.L.; Torre, A.; Vázquez, L. Evolution operator equations: Integration with algebraic and finite difference methods. Applications to physical problems in classical and quantum mechanics and quantum field theory. La Rivista del Nuovo Cimento 1997, 20, 3–133. [Google Scholar] [CrossRef]
  41. Wei, J.; Norman, E. Lie algebraic solution of linear differential equations. J. Math. Phys. 1963, 4, 575–581. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Behr, N.; Dattoli, G.; Lattanzi, A.; Licciardi, S. Dual Numbers and Operational Umbral Methods. Axioms 2019, 8, 77. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms8030077

AMA Style

Behr N, Dattoli G, Lattanzi A, Licciardi S. Dual Numbers and Operational Umbral Methods. Axioms. 2019; 8(3):77. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms8030077

Chicago/Turabian Style

Behr, Nicolas, Giuseppe Dattoli, Ambra Lattanzi, and Silvia Licciardi. 2019. "Dual Numbers and Operational Umbral Methods" Axioms 8, no. 3: 77. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms8030077

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop