Cauchy Formulae and Hardy Spaces in Discrete Octonionic Analysis

In this paper, we continue the development of a fundament of discrete octonionic analysis that is associated to the discrete first order Cauchy–Riemann operator acting on octonions. In particular, we establish a discrete octonionic version of the Borel–Pompeiu formula and of Cauchy’s integral formula. The latter then is exploited to introduce a discrete monogenic octonionic Cauchy transform. This tool in hand allows us to introduce discrete octonionic Hardy spaces for upper and lower half-space together with Plemelj projection formulae.


Introduction
Complex analytic tools together with their related discretisation methods provide a very powerful toolkit to compute numerically and some cases even analytically the solutions to many partial differential equations with given boundary data that arise in the context of harmonic analysis.In the focus there are null-solutions to a discrete version of the Cauchy-Riemann operator that satisfy a discretised version of the Cauchy integral and the Borel-Pompeiu formula.These integral formulas provide us with the key ingredients in the representation of the solutions.Now, there are several possibilities to generalise complex function theory together with discretisations to higher dimensional settings.Many branches of engineering particularly focus on three-dimensional generalisations.To address the three-dimensional setting the function theory related to associative Clifford algebras has proven to be a very effective tool.In this context one considers a three-dimensional Cauchy-Riemann operator defined by D := ∂ xi e i where e 1 and e 2 are two different imaginary units generating the Clifford algebra Cℓ 0,2 which is spanned as an R-vector space by the basis elements 1, e 1 , e 2 , e 1 e 2 .Actually the particular case Cℓ 0,2 represents the Hamiltonian quaternions which for a non-commutative skew-field.
Using the multiplicative structure of the Clifford algebra, the generalised first order Cauchy-Riemann operator factorises the second-order Laplacian in R 3 .In fact, this approach could be generalised very easily to all associative Clifford algebras.Its function theory now is broadly known under the term Clifford analysis while functions in the kernel of the generalised Cauchy-Riemann operator are often called monogenic, hyperholomorphic or Clifford holomorphic functions.Classical references are for instance the books [2,19] as well as also J. Ryan's edited volume containing the milestone contribution [25] among others which all provided a boost in the development of this function theory including applications to Calderon-Zygmund type operators.Particularly, [19] exhibits how related integral operators, in particular Plemelj-Sokhotzkij type formulas related to the Cauchy transform, can successfully be applied to treat boundary value problems.Additionally, J. Ryan and his co-authors have also addressed unbounded domains with this function theoretical toolkit in [18].
Over the last two decades there has been a growing interest in the discretisation of the continuous Clifford analysis for developing numerical algorithms for the higher dimensional associative framework, cf. for instance [3,5,6,8,9,13,14,17] among others.However, instead of embedding the vector space R n+1 into associative Clifford algebras Cℓ 0,n , which have zero-divisors for n > 2, there are many more ways to generalise complex analysis to higher dimensions.If the Cayley-Dickson duplication process to the complex numbers is applied, then we first arrive at the four-dimensional Hamiltonian quaternions, which, as mentioned is still a particular case of a Clifford algebra.However, if we apply, as the next step, the Cayley-Dickson doubling to the quaternions, then we obtain a different sort of algebra, namely the octonions or Cayley numbers, denoted by O in all that follows, cf. for instance [1] among others.The octonions turn out to be non-associative.Consequently they do not belong to the Clifford algebras and they are not representable with matrices.But, in contrast to the eight-dimensional Clifford algebra Cℓ 0,3 octonions still form a normed division algebra (in the wider sense of being non-associative without zero-divisors).So, every non-zero element has a multiplicative inverse.Accord-ing to the famous theorem of Hurwitz, they form the largest real normed zero divisor-free algebra over R.
Surprisingly, although the octonions are not associative, according to recently published research results in particle physics, see for instance [4,16,26], they seem to offer a much more adequate model for a unified description of particle physics including gravity, which was also already proposed in the book by F. Gürsey and H. Tze in 1996, see [20].This provides one further motivation to develop analytic and discretised function theoretic tools in the non-associative context of octonions.On the continuous level at least the fundaments of an octonionic function theory are already well-developed, see for example [12,27,32,33].In [34] generalisations of the Cauchy transform together with Plemelj projection formulas and with some basic applications to Calderon-Zygmund type operators have been presented particularly.Recently one also managed to introduce meaningful octonionic generalisations of Bergman and Hardy spaces on the continuous level, see for example [30,31,11,10], as well as the papers [15,28] addressing octonionic Hilbert spaces on a more general level.
Nevertheless, if we want to solve practical octonionic boundary value problems numerically, then we are in need of discretised versions of these octonionic operators.As far as we know, the development of a discrete octonionic function theory is still a rather open research field.In our recent paper [23], we developed a fundament for research in this direction, namely, we introduced a discretised version of octonionic Cauchy-Riemann operator in terms of appropriate forward and backward operators in the Hermitian sense, and established a discrete Stokes' formula for that operator.
In this paper, we depart from the discretised Stokes' formula and exploit that one further in order to obtain a discrete octonionic version of Borel-Pompeiu's formula.The latter in turn produces a discrete Cauchy formula in a special case.As a next step, we establish a discretised version of the octonionic Cauchy transform that has been introduced on the continuous level in [34].The related Plemelj-Sokhotzkij projection formulae then give rise to the definition of a discrete octonionic monogenic generalisation of Hardy spaces.In particular, the focus of this paper is put on the consideration of the upper octonionic half-space (resp.its related upper half-lattice) where the last component satisfies x 7 > 0. In this case, we have a simple geometric description of the inner product that has to be used here.An alternative way to introduce discrete octonionic monogenic Hardy spaces consists in defining appropriate extension operators using the Fourier transform on the complexified octonions and its related Fourier symbols, which in our framework are also complexified octonions.Also this approach will be carefully developed in this paper as it has the advantage of being fully explicit.
As a consequence of the lack of associativity, the signs in the discrete octonionic Borel-Pompeiu and Cauchy differ from the results in discrete Clifford analysis.Surprisingly, the discretisation of the continuous octonionic analysis that we develop here has the interesting effect that the associator expressions that we intrinsically have in the Stokes formula of the continuous case (see [34]), disappear in our discretised constructions.This effect also shows up when considering instead of the forward and backward operator its centralised version to which we also pay a particular attention in this paper to get a more complete view on this particular phenomenon.
The discrete octonionic setting thus turns out to be really different as well from its associative version in Clifford analysis but also from continuous octonionic monogenic function theory.It exposes new intrinsic peculiarities that have to be taken really carefully into account.

Continuous octonionic analysis
In this section, we briefly recall basic notions and the most fundamental results on continuous octonionic analysis.The basic context is the 8-dimensional Euclidean vector space R 8 , where the standard vectors are denoted by e k , k = 0, 1, . . ., 7. A vector from R 8 can be expressed as usual in terms of its real coordinates in the way x = (x 0 , x 1 , . . ., x 7 ).Next, R 8 -vectors can also be described as octonions x = x 0 e 0 + x 1 e 1 + x 2 e 2 + x 3 e 3 + x 4 e 4 + x 5 e 5 + x 6 e 6 + x 7 e 7 , where we now additionally identify e 4 = e 1 e 2 , e 5 = e 1 e 3 , e 6 = e 2 e 3 and e 7 = e 4 e 3 = (e 1 e 2 )e 3 .Moreover, we have e 2 i = −1 and e 0 e i = e i e 0 for all i = 1, . . ., 7, and e i e j = −e j e i for all mutual distinct i, j ∈ {1, . . ., 7}, as well as e 0 is the neutral element and, therefore, often will be omitted.This definition endows R 8 additionally with a multiplicative closed structure.
Table 1 fully describes the multiplication rules for real octonions.This table clearly indicates in particular that the octonionic multiplication actually is closed but not associative, for instance we have (e i e j )e k = −e i (e j e k ).Here, we use the same labelling of the basis elements as used in [1].Additionally, we also need to consider the octonionic conjugate which is given by x = x 0 e 0 − x 1 e 1 − x 2 e 2 − x 3 e 3 − x 4 e 4 − x 5 e 5 − x 6 e 6 − x 7 e 7 .
The Euclidean norm from R 8 then is reproduced by x • x = 7 j=0 x 2 i = |x| 2 , so every non-zero octonion x is invertible via x −1 = x |x| 2 .While there are several possibilities to extend the classical function theory to octonions, we recall here the definition in the sense of the Riemannapproach, following the classical development of P. Dentoni and M. Sce [12], K. Nono [27], the school of Xingmin-Li and Zhong Peng, see for instance [32] and others.In this sense we recall: is the octonionic first order Cauchy-Riemann operator.If f satisfies Df = 0 (resp.f D = 0), then we call f left (right) octonionic anti-monogenic.
Octonionic analysis has two crucial differences to classical Clifford analysis: (i) Octonionic analysis considers functions from R 8 back into R 8 , while Clifford analysis addresses null-solutions to the Cauchy-Riemann operator defined on the paravector space R ⊕ R 7 with values in the Clifford algebra Cℓ 7 .Cℓ 7 however is a real vector space being isomorphic to R 128 .(ii) Left(right) octonionic monogenic functions do neither form a right nor a left O-module.See for example J. Kauhanen and H. Orelma in [22] for concrete examples.This fact significantly complicates the development of a consistent function theory and the theory of generalised Hilbert function spaces in octonionic settings, see also [10,11,28].
Additionally, the lack of associativity leads to modifications of the classical integral formulae, such as for example the Stokes' formula [34]: Here, the expression [a, b, c] := (ab)c − a(bc) called the associator appears.This is an intrinsic feature, which is cancelled out in the cases of associativity.It is important to remark, that although the associator appears in most of octonionic constructions, it is nevertheless possible to introduce specific structures, where the associator would vanish.For example, it has been pointed out in [32], that considering the two functions being octonionic monogenic and Stein-Weiss conjugate harmonics, i.e. ∂g j ∂x i = ∂g i ∂x j for all 0 ≤ i < j ≤ 7, the associator will vanish.Further, a generalisation of the Cauchy's integral formula to octonionic setting has been presented, see for instance [27,33]: open and G ⊆ U be an 8-D compact oriented manifold with a strongly Lipschitz boundary ∂G.
However, if the parenthesis were set differently, then one would obtain the different formula which involves the associator again.

Discretisation of octonionic analysis
Let us consider the unbounded uniform lattice hZ 8 with the lattice constant h > 0, which is defined in the classical way as follows Next, we define the classical forward and backward differences ∂ ±j h as for discrete functions f (mh) with mh ∈ hZ 8 .In the sequel, we consider functions defined on Ω h ⊂ hZ 8 and taking values in octonions O.As usual, all important properties such as the l p -summability (1 ≤ p < ∞) are defined component-wisely.A next step is to introduce discretisations of the Cauchy-Riemann operators in octonions.In contrast to the classical discrete Clifford analysis presented for instance in [3,14], where fundamental ideas of the Weyl calculus are used, we follow the alternative approach approach presented in [13] which uses a direct discretisation of the continuous Cauchy-Riemann (or Dirac) operators in terms of forward and backward finite difference operators.Following this approach, the non-associativity of octonionic multiplication can be respected, see [23] for details.Hence, by using the finite difference operators (2), we introduce a discrete forward Cauchy-Riemann operator Direct computations show that the star-Laplacian ∆ h is represented in this setting as follows: where D + h and D − h are the discrete conjugated forward and backward Cauchy-Riemann operators, respectively: Since we consider discrete forward and backward Cauchy-Riemann operators, it is also necessary to distinguish between discrete forward and backward monogenic functions:

Discrete fundamental solution
To construct discrete versions of the Borel-Pompeiu and Cauchy formulae, it is necessary to work with discrete fundamental solutions of the discrete Cauchy-Riemann operator and of the discrete Laplace operator.Thus, the following definition introduces a discrete fundamental solution to the discrete Cauchy-Riemann operators (3): , for mh = 0, 0, for mh = 0, for all grid points mh of hZ 8 .Analogously, the function for all grid points mh of hZ 8 .
As usual, a discrete fundamental solution can be constructed by means of the discrete Fourier transform of u ∈ l p hZ 8 , O , 1 ≤ p < +∞, where mh, ξ = h 8 j=1 m j ξ j .It is worth to underline, that the last expression implicitly introduces a complexified octonionic structure in the sense of O C := O ⊗ R C, implying that the complex imaginary unit i commutes with all real octonions.The inverse transform is given by and R h denotes its restriction to the lattice hZ 8 .Note again that f is O-valued, so the expression Next, we recall the known symbols for the forward and backward differences ∂ ±j h , namely ξ ±j h = ∓h −1 1 − e ∓ihξj , as well as the symbol for the star-Laplacian: Thus, applying the discrete Fourier transform to forward and backward Cauchy-Riemann operators leads to implying that the operators D ± h have complexified octonionic symbols ξ ± = 7 j=0 e j ξ ±j h , respectively.Hence, the fundamental solutions E ± h can be expressed by Evidently, the discrete fundamental solutions E ± h constructed above are different to the discrete fundamental solutions typically considered in the framework of discrete Clifford analysis, see for example [6].The difference comes from the fact, that we work with a direct discretisation of the Cauchy-Riemann operators without involving the splitting of basis unit vectors, as it is done in discrete Clifford analysis.
Next, we will provide two basic properties of the discrete fundamental solutions E ± h .For proving one of these properties, we need to recall before the following theorem [29]: Theorem 1 Let n be the dimension of the Euclidean space and assume that p 1 , p 2 are two positive integers with p 2 < p 1 + n.For a positive integer N > 0 let κ N consider the set of functions of the form where T j (Θ) are trigonometric polynomials which satisfy the following conditions: (i) there are ordinary homogeneous polynomials P j (Θ) of degree p j , j = 1, 2, such that For any N > 0 satisfying (ii)-(iv) there is a constant C such that for all µ (with integer components) and These tools in hand now allow us to establish the following theorem: Theorem 2 The discrete fundamental solutions E ± h to the discrete forward and backward Cauchy-Riemann operators satisfy: Proof The proof of property (i) can be done by a straightforward calculations.
To prove (ii), we are going to use the integral representation of the discrete fundamental solution E + h : Taking into account the definition of ξ+ and that and after using known trigonometric identity, we get sin(hξ j ) hd 2 e −i mh,ξ dξ .
The conditons of Theorem 1 are satisfied for the two integrals above if p 1 = p 2 = 2 for the first integral, and if p 1 = 1, p 2 = 2 for the second integral, see also [24] for a detailed discussion.Thus, the following estimate is obtained: Hence, it can be concluded that the fundamental solution E + h belongs to l p (hZ 8 , O) for p > 8  7 .The proof for the formula of E − h can be performed analogously.

Discrete Stokes' and Borel-Pompeiu formulae
We start this section by recalling the discrete octonionic Stokes' formulae from [23].Again, for shortening the notations, the long list of indices m 0 , m 1 , . . ., m 7 , will be omitted from the argument, i.e. we will simply write f (mh) instead of f (m 0 h, m 1 h, . . ., m 7 h).
We are particularly interested in studying upper and lower half-spaces (or half-lattices), which are defined as follows: The following two theorems express the discrete octonionic Stokes' formulae for the upper half-lattice and for the lower half-lattice, respectively [23]: Theorem 3 The discrete Stokes' formula for the upper half-lattice hZ 8  + is given by for all discrete functions f and g such that the series converge.

Theorem 4
The discrete Stokes' formula for the lower half-lattice hZ 8 − is given by for all discrete functions f and g such that the series converge.
Remark 1 We would like to remark that the discrete Stokes' formulae introduced above do not contain the associator, which appears in the continuous case as an intrinsic term.The non-appearance of the associator therefore represents a surprising effect of the discrete setting.Moreover, the cancellation of the associator is not related to the consideration of the forward or backward discrete Cauchy-Riemann operators, either.Note that if we consider instead for example the central discrete Cauchy-Riemann operator: then the discrete octonionic formula for the whole space would have the form For the sake of completeness, let us briefly outline the proof of this formula.We start with the first summand (h is omitted for the sake of abbreviation) By using the relation (e i e j )e k = −e i (e j e k ) and the definitions of finite differences, we get We just have proved the discrete octonionic Stokes' formula for the whole space in the case of the central discrete Cauchy-Riemann operator.This highlights that the constructions in the discrete settings exhibit an essentially different nature than the constructions in the continuous case, where the associstor appears as an intrinsic ingredient.
By means of the discrete discrete octonionic Stokes' formule ( 6)-( 7), the discrete octononionic Borel-Pompeiu formulae can be introduced: h be the discrete fundamental solution to the discrete Cauchy-Riemann operator D + h .Then the discrete octonionic Borel-Pompeiu formula for the upper half-lattice hZ 8  + is given by for any discrete function f such that the series converge.
Proof To prove the discrete octonionic Borel-Pompeiu formula for the upper half-lattice, we use at first the discrete octonionic Stokes's formula for the upper half-lattice (6), and replace therein the function g by the shifted discrete fundamental solution E + h (• − mh) with m ∈ Z 8 + .Next, by taking into account the definition of the discrete fundamental solution, we note that E + h (nh − mh)D + h = 0 for n = m and E + h (nh − mh)D + h = h −8 for n = m.Thus, the formula is established.
By using the discrete octonionic Borel-Pompeiu formula (8) and requiring that the function f is discrete left backward monogenic in hZ 8 + , we immediately arrive at the discrete octonionic Cauchy formula: Theorem 6 Let f be a discrete left backward monogenic function with respect to the operator D − h , and let E + h be the discrete fundamental solution to the operator D + h .Then the discrete octonionic Cauchy formula for the upper halflattice hZ 8  + is given by which holds for any discrete function f such that the series converge.
Analogously we may introduce discrete octonionic Borel-Pompeiu and Cauchy formulae for the lower half-lattice hZ + − : Corollary 1 Let f be a discrete left backward monogenic function with respect to the operator D − h , and let E + h be the discrete fundamental solution to operator D + h .Then the discrete octonionic Borel-Pompeiu formula for the lower halflattice hZ 8  − is given by for any discrete function f such that the series converge.In the case when f is a discrete left backward monogenic function with respect to the operator D − h , we obtain the discrete octonionic Cauchy formula for the lower half-lattice hZ 8  − in the form Theorem 6 and Corollary 1 immediately lead to the definition of discrete octonionic Cauchy transforms for the upper and lower half-lattices: Definition 4 For a discrete l p -function f , 1 ≤ p < +∞, defined on the boundary layers (m, 0) and (m, 1) the discrete octonionic Cauchy transform for the upper half-lattice hZ 8  + is defined by Analogously, for a discrete l p -function f , 1 ≤ p < +∞, defined on the boundary layers (m, 0) and (m, −1) the discrete octonionic Cauchy transform for the lower half-lattice hZ 8 − is defined by It is worth to mention, that alike in the continuous case also in discrete Clifford analysis, the discrete octonionic Cauchy formulae and Cauchy transform indicate the dependence of (discrete) monogenic functions on their boundary values.However, as a consequence of the lack of associativity, the signs in the discrete octonionic Borel-Pompeiu and Cauchy formulae ( 8)-( 11) differ from the results in discrete Clifford analysis, compare for example with [6].Additionally, in line with the discussion in [23], the discretisation of the continuous octonionic analysis has the interesting effect that the associator disappears from the constructions.
Let us now present some properties of the discrete octonionic Cauchy transforms: Theorem 7 Let us consider the discrete upper half-lattice hZ 8  + and the lower half-lattice hZ 8 − .Then the discrete Cauchy transforms ( 12)-( 13) satisfy the following properties: (i) The interior and exterior Cauchy transforms have the following mapping properties: Proof The proof of this theorem can be performed along the same ideas as in [6,9], and, therefore, we will only mention that the first statement follows from a direct application of Hölder's inequality and the properties of the discrete fundamental solution E + h .The proof of (ii) and (iii) is done by straightforward application of the discrete forward Cauchy-Riemann operator D + h and applying its properties.

Discrete octonionic Hardy spaces
The aim of this section is to construct discrete octonionic Hardy spaces.Looking at the results in discrete Clifford analysis related to Hardy spaces, see for example [6,8,9], it becomes evident that a typical approach to introduce discrete Hardy spaces will be to work with Fourier transforms on boundary layers of the discrete fundamental solution.In this way, discrete Riesz kernels can be defined, and, hence, discrete Plemelj (or Hardy) projections are then introduced.Alternatively, one could follow the classical continuous approach for defining Hardy spaces and adapt it to the discrete setting.For the purpose of a better understanding the effect of non-associativity, both approaches will be presented in this section.
At first, we recall from [10] the following definition of an octonionic Hilbert space, which straightforwardly extends to the discrete setting: Definition 5 An octonionic Hilbert space H is a left O-module with an octonion-valued inner product •, • : H × H → O such that (H, •, • 0 ) is a real Hilbert space, where •, • 0 := Re •, • .The octonion-valued inner product is supposed to satisfy for all f, g, h ∈ H and for all α ∈ O the following rules: Before introducing the notion of a discrete octonionic Hardy space, let us fix the notation γ h for the discrete boundary of Ω h (0-layer in the case of half-lattices).Now we introduce the following notion of a discrete octonionic Hardy space, which is adapted from the definition presented in [10]: Definition 6 The discrete octonionic Hardy space H 2 h (Ω h , O) is the closure of the set of l 2 (γ h , O) functions that are discrete left monogenic inside of Ω h and extendable to the boundary γ h .
It is important to remark, that we have not distinguished in this definition between discrete forward and backward monogenic functions as it has been introduced in Definition 2. The reason for this is to avoid too many specific definitions, which do not generally change the notion of a discrete octonionic Hardy space, since both types of discrete monogenic functions can be used.Further, Definition 6 requires a discrete monogenic function being extendable to the boundary γ h , implying that an extension operator must be defined.Although there are several possibilities to define such an operator, arguably the most straightforward way is to work with the Fourier transform of the discrete fundamental solution E + h , as it is done in the discrete Clifford analysis, see again [6,8,9].Hence, we will introduce an extension operator based on this idea, and thus, establish a clear connection between the discrete Clifford analysis approach and the "continuous theory-based" approach.
Next, we need to introduce an octonion-valued inner product on the discrete boundary γ h : Definition 7 For discrete octonionic function f, g ∈ l 2 (γ h ), we introduce the following O-valued inner products: if the upper half-lattice is considered, and f, g if the lower half-lattice is considered.
It is then easy to verify the following statement: The next step is to establish the theory of discrete Hardy spaces in the sense of discrete Clifford analysis.This requires a definition of discrete Riesz kernels (convolution kernels) implying that the behaviour of the discrete fundamental solution E + h on boundary layers needs to be studied.Let us recall the integral representation of the discrete fundamental solution e j ξ +j h .To study Fourier symbols on the boundary layers, i.e. for m 7 ∈ {−1, 0, 1}, we apply the 7-dimensional discrete Fourier transform to the discrete fundamental solution: dξ.
Let us study now the integral (I): The integrals above need to be calculated for m 7 ∈ {−1, 0, 1}, and these calculations have been presented in [6] for the case of splitting of basis unit elements.Hence, the Fourier symbols of the discrete fundamental solutions on the layers m 7 = −1, m 7 = 0, and m 7 = 1 are given by F Next, following ideas in [6], we introduce the pair of operators: which fulfil the condition (H + ) 2 = (H − ) 2 = I.By help of these operators, we can formulate conditions for a function to be a discrete boundary value of a discrete octonionic monogenic function in hZ 8 + or hZ 8 − : Finally, we define two extension operators, which extend a discrete function from layers m 7 = −1 and m 7 = 1 to the boundary layer m 7 = 0, see [7] for the details: Definition 9 The upper extension operator, denoted as A + , is an operator extending a function given on the boundary layer m 7 = 1 to the boundary layer m 7 = 0, i.e. it is a mapping A + : l p (hZ 7 ) → l p (hZ 7 ) given by A + := F Similarly, the lower extension operator, denoted by A − , is an operator extending a function given on the boundary layer m 7 = −1 to the boundary layer m 7 = 0, i.e. it is a mapping A − : l p (hZ 7 ) → l p (hZ 7 ) given by This definition rounds off the discussion around discrete octonionic Hardy spaces, which we have initiated in Definition 6.Now it is clear how an extension of a discrete function to the boundary layer m 7 = 0 can be explained.

Summary
In this paper, we have continued the development of discrete octonionic analysis by introducing discrete Borel-Pompeiu and Cauchy formulae, and defined discrete octonionic Hardy spaces for half-spaces.Moreover, we have discussed two approaches for constructing discrete Hardy spaces: (i) by a direct discretisation of the continuous case; and (ii) by using an approach used in discrete Clifford analysis by studying Fourier symbols of the discrete fundamental solution of the discrete Cauchy-Riemann operator on the boundary layers.Both approaches complement each other and contribute to a better understanding of the discrete octonionic setting.
A very surprising result is that the associator, which appears in the continuous case, does not appear in the discrete setting.For a better understanding of this effect, the discrete Stokes' formula for the whole space has been proved also for the central discrete Cauchy-Riemann operator.Nevertheless, the associator does not appear even in that case, which underlines the particularity of the constructions in the discrete setting.
The results presented in this paper provide us with a powerful basic toolkit for a further development of discrete octonionic analysis.In particular, the discrete versions of octonionic Hardy spaces that we introduced allow us to study concrete boundary value problems for monogenic functions in the discrete octonionic setting.This also opens the door to study Calderon-Zygmund type operators in this context.Furthermore, the discrete toolkit allow us us develop numerical algorithms for solving some particular physical problems arising in the unification of particle physics and gravity as illustrated in a series of recent works in this direction, see again [4,26], such as already mentioned roughly in the introductory text.
Further, after understanding the difficulties arising on the way discretising octonionic analysis for the half-space settings, the results of this paper can subsequentially be extended to the case of considering bounded domains in R 8 .

8 .
g i (m + e j ) − g i (m)) f k (m)e i (e j e k )− − (g i (m) − g i (m − e j )) f k (m)e i (e j e k )] h 8 i (m + e j )e i f k (m) + g i (m − e j )e i f k (m)] (e j e k )hPerforming a change of variables in the latter expression, we get i (m)e i f k (m − e j ) + g i (m)e i f k (m + e j )] (e j e k )h 8 i (m)e i (f k (m + e j ) + f k (m − e j ))] (e j e k )h8 spaces in the sense of Definition 5.

Definition 8 8 + 8 −
This enables us now to introduce another definition of discrete octonionic Hardy spaces: The space of discrete functions f ∈ l p (hZ 8 + , O) whose discrete 2D-Fourier transform fulfills f = H + f for m 7 = 1 is called the upper discrete octonionic Hardy space and it is denoted by h + p,hZ 8 + .Analogously, the space of discrete functions f ∈ l p (hZ 8 − , O) whose discrete 2D-Fourier transform fulfills f = H − f for m 7 = −1 is called the lower discrete octonionic Hardy space and it is denoted by h − p,hZ 8 − .By means of the operators H + and H − , the discrete Plemelj or Hardy projections can be now introduced as follows P + = 1 2 (I + H + ) and P − = 1 2 (I + H − ) .Additionally, combining these projections with the previous definition we get f ∈ h + p,hZ ⇐⇒ P + f = f, and f ∈ h − p,hZ ⇐⇒ P − f = f.

Table 1
Multiplication table for real octonions O