 # Non-constant bounded holomorphic functions of hyperbolic numbers - Candidates for hyperbolic activation functions

The Liouville theorem states that bounded holomorphic complex functions are necessarily constant. Holomorphic functions fulfill the socalled Cauchy-Riemann (CR) conditions. The CR conditions mean that a complex z-derivative is independent of the direction. Holomorphic functions are ideal for activation functions of complex neural networks, but the Liouville theorem makes them useless. Yet recently the use of hyperbolic numbers, lead to the construction of hyperbolic number neural networks. We will describe the Cauchy-Riemann conditions for hyperbolic numbers and show that there exists a new interesting type of bounded holomorphic functions of hyperbolic numbers, which are not constant. We give examples of such functions. They therefore substantially expand the available candidates for holomorphic activation functions for hyperbolic number neural networks. Keywords: Hyperbolic numbers, Liouville theorem, Cauchy-Riemann conditions, bounded holomorphic functions

## Authors

##### This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

## 1 Introduction

For the sake of mathematical clarity, we first carefully review the notion of holomorphic functions in the two number systems of complex and hyperbolic numbers.

The Liouville theorem states that bounded holomorphic complex functions are necessarily constant . Holomorphic functions are functions that fulfill the socalled Cauchy-Riemann (CR) conditions. The CR conditions mean that a complex -derivative

 df(z)dz,z=x+iy∈C,x,y∈R,ii=−1, (1)

is independent of the direction with respect to which the incremental ratio, that defines the derivative, is taken . Holomorphic functions would be ideal for activation functions of complex neural networks, but the Liouville theorem means that careful measures need to be taken in order to avoid poles (where the function becomes infinite).

Yet recently the use of hyperbolic numbers

 z=x+hy,h2=1,x,y∈R,h∉R. (2)

lead to the construction of hyperbolic number neural networks. We will describe the generalized Cauchy-Riemann conditions for hyperbolic numbers and show that there exist bounded holomorphic functions of hyperbolic numbers, which are not constant. We give a new example of such a function. They are therefore excellent candidates for holomorphic activation functions for hyperbolic number neural networks [2, 3]. In 

it was shown, that hyperbolic number neural networks allow to control the angle of the decision boundaries (hyperplanes) of the real and the unipotent

-part of the output. But Buchholz argued in , p. 114, that

Contrary to the complex case, the hyperbolic logistic function is bounded. This is due to the absence of singularities. Thus, in general terms, this seems to be a suitable activation function. Concretely, the following facts, however, might be of disadvantage. The real and imaginary part have different squashing values. Both component functions do only significantly differ from zero around the lines111Note that we slightly correct the two formulas of Buchholz, because we think it necessary to delete in Buchholz’ original , etc. and .

Complex numbers are isomorphic to the Clifford geometric algebra

which is generated by a single vector

of negative square , with algebraic basis . The isomorphism is realized by mapping .

Hyperbolic numbers are isomorphic to the Clifford geometric algebra which is generated by a single vector of positive square , with algebraic basis . The isomorphism between hyperoblic numbers and is realized by mapping .

## 2 Complex variable functions

We follow the treatment given in . We assume a complex function given by an absolute convergent power series.

 w=f(z)=f(x+iy)=u(x,y)+iv(x,y), (3)

where are real functions of the real variables . Since are obtained in an algebraic way from the complex number , they cannot be arbitrary functions but must satisfy certain conditions. There are several equivalent ways to obtain these conditions. Following Riemann, we state that a function is a function of the complex variable

if its derivative is independent of the direction (in the complex plane) with respect to which the incremental ratio is taken. This requirement leads to two partial differential equations, named after Cauchy and Riemann (CR), which relate

and .

One method for obtaining these equations is the following. We consider the expression only as a function of , but not of , i.e. the derivative with respect to shall be zero. First we perform the bijective substitution

 x=12(z+¯z),y=−i12(z−¯z), (4)

based on . For computing the derivative

with the help of the chain rule we need the derivatives of

and of (4)

 x,¯z=12,y,¯z=12i. (5)

Using the chain rule we obtain

 w,¯z =u,xx,¯z+u,yy,¯z+i(v,xx,¯z+v,yy,¯z) =12u,x+12iu,y+i(12v,x+12iv,y) =12[u,x−v,y+i(v,x+u,y)]!=0. (6)

Requiring that both the real and the imaginary part of (6) vanish we obtain the Cauchy-Riemann conditions

 u,x=v,y,u,y=−v,x. (7)

Functions of a complex variable that fulfill the CR conditions are functions of and , but they are only functions of , not of .

It follows from (7), that both and fulfill the Laplace equation

 u,xx=v,yx=v,xy=−u,yy⇔u,xx+u,yy=0, (8)

and similarly

 v,xx+v,yy=0. (9)

The Laplace equation is a simple example of an elliptic partial differential equation. The general theory of solutions to the Laplace equation is known as potential theory. The solutions of the Laplace equation are called harmonic functions and are important in many fields of science, notably the fields of electromagnetism, astronomy, and fluid dynamics, because they can be used to accurately describe the behavior of electric, gravitational, and fluid potentials. In the study of heat conduction, the Laplace equation is the steady-state heat equation .

Liouville’s theorem  states, that any bounded holomorphic function , which fulfills the CR conditions is constant. Therefore for complex neural networks it is not very meaningful to use holomorphic functions as activation functions. If they are used, special measures need to be taken to avoid poles in the complex plane. Instead separate componentwise (split) real scalar functions for the real part , and for the imaginary part , are usually adopted. Therefore a standard split activation function in the complex domain is given by

 g(u(x,y)+iv(x,y))=gr(u(x,y))+igi(v(x,y)). (10)

## 3 Hyperbolic numbers

Hyperbolic numbers are also known as split-complex numbers. They form a two-dimensional commutative algebra. The canonical hyperbolic system of numbers is defined  by

 z=x+hy,h2=1,x,y∈R,h∉R. (11)

The hyperbolic conjugate is defined as

 ¯z=x−hy. (12)

Taking the hyperbolic conjugate corresponds in the isomorphic algebra to taking the main involution (grade involution), which maps .

The hyperbolic invariant (corresponding to the Lorentz invariant in physics for ), or modulus, is defined as

 z¯z=(x+hy)(x−hy)=x2−y2, (13)

which is not positive definite.

Hyperbolic numbers are fundamentally different from complex numbers. Complex numbers and quaternions are division algebras, every non-zero element has a unique inverse. Hyperbolic numbers do not always have an inverse, but instead there are idempotents and divisors of zero.

We can define the following idempotent basis

 n1=12(1+h),n2=12(1−h), (14)

which fulfills

 n21 =14(1+h)(1+h)=14(2+2h)=n1, n22 =n2,n1+n2=1, n1n2 =14(1+h)(1−h)=14(1−1)=0, ¯n1 =n2,¯n2=n1. (15)

The inverse basis transformation is simply

 1=n1+n2,h=n1−n2. (16)

Setting

 z=x+hy=ξn1+ηn2, (17)

we get the corresponding coordinate transformation

 x=12(ξ+η),y=12(ξ−η), (18)

as well as the inverse coordinate transformation

 ξ=x+y∈R,η=x−y∈R. (19)

The hyperbolic conjugate becomes, due to (15), in the idempotent basis

 ¯z=ξ¯n1+η¯n2=ηn1+ξn2. (20)

In the idempotent basis, using (20) and (15), the hyperbolic invariant becomes multiplicative

 z¯z =(ξn1+ηn2)(ηn1+ξn2) =ξη(n1+n2)=ξη=x2−y2. (21)

In the following we consider the product and quotient of two hyperbolic numbers both expressed in the idempotent basis

 zz′=(ξn1+ηn2)(ξ′n1+η′n2)=ξξ′n1+ηη′n2, (22)

and

 zz′ =ξn1+ηn2ξ′n1+η′n2=z¯z′z′¯z′ =(ξn1+ηn2)(η′n1+ξ′n2)(ξ′n1+η′n2)(η′n1+ξ′n2) =(ξη′n1+ηξ′n2)(η′n1+ξ′n2)ξ′η′ =ξξ′n1+ηη′n2. (23)

Because of (23) it is not possible to divide by if , or if . Moreover, the product of a hyperbolic number with (on the axis) times a hyperbolic number with (on the axis) is

 (ξn1+0n2)(0n1+ηn2)=ξηn1n2=0, (24)

due to (15). We repeat that in (24) the product is zero, even though the factors are non-zero. The numbers along the axis are therefore called divisors of zero. The divisors of zero have no inverse.

The hyperbolic plane with the diagonal lines of divisors of zero (b), and the pairs of hyperbolas with constant modulus (c), and (a) is shown in Fig. 1. Figure 1: The hyperbolic number plane  with horizontal x-axis and vertical yh-axis, showing: (a) Hyperbolas with modulus z¯z=−1 (green). (b) Straight lines with modulus z¯z=0⇔x2=y2 (red), i.e. divisors of zero. (c) Hyperbolas with modulus z¯z=1 (blue).

## 4 Hyperbolic number functions

We assume a hyperbolic number function given by an absolute convergent power series

 w =f(z)=f(x+hy)=u(x,y)+hv(x,y), h2 =1,h∉R. (25)

where are real functions of the real variables . An example of a hyperbolic number function is the exponential function

 ez =ex+hy=exehy=ex(coshy+hsinhy) =u(x,y)+hv(x,y), (26)

with

 u(x,y)=excoshy,v(x,y)=exsinhy. (27)

Since are obtained in an algebraic way from the hyperbolic number , they cannot be arbitrary functions but must satisfy certain conditions. There are several equivalent ways to obtain these conditions. A function is a function of the hyperbolic variable , if its derivative is independent of the direction (in the hyperbolic plane) with respect to which the incremental ratio is taken. This requirement leads to two partial differential equations, so called generalized Cauchy-Riemann (GCR) conditions, which relate and .

To obtain the GCR conditions we consider the expression only as a function of , but not of , i.e. the derivative with respect to shall be zero. First we perform the bijective substitution

 x=12(z+¯z),y=h12(z−¯z), (28)

based on . For computing the derivative with the help of the chain rule we need the derivatives of and of (28)

 x,¯z=12,y,¯z=−12h. (29)

Using the chain rule we obtain

 w,¯z =u,xx,¯z+u,yy,¯z+h(v,xx,¯z+v,yy,¯z) =12u,x−12hu,y+h(12v,x−12hv,y) =12[u,x−v,y+h(v,x−u,y)]!=0. (30)

Requiring that both the real and the -part of (30) vanish we obtain the GCR conditions

 u,x=v,y,u,y=v,x. (31)

Functions of a hyperbolic variable that fulfill the GCR conditions are functions of and , but they are only functions of , not of . Such functions are called (hyperbolic) holomorphic functions.

It follows from (31), that and fulfill the wave equation

 u,xx=v,yx=v,xy=u,yy⇔u,xx−u,yy=0, (32)

and similarly

 v,xx−v,yy=0. (33)

The wave equation is an important second-order linear partial differential equation for the description of waves – as they occur in physics – such as sound waves, light waves and water waves. It arises in fields like acoustics, electromagnetics, and fluid dynamics. The wave equation is the prototype of a hyperbolic partial differential equation .

Let us compute the partial derivatives , for the exponential function of (26):

 u,x =excoshy,u,y=exsinhy, v,x =exsinhy=u,y,v,y=excoshy=u,x. (34)

We clearly see that the partial derivatives (34) fulfill the GCR conditions (31) for the exponential function , as expected by its definition (26). The exponential function is therefore a manifestly holomorphic hyperpolic function, but it is not bounded.

In the case of holomorphic hyperbolic functions the GCR conditions do not imply a Liouville type theorem like for holomorphic complex functions. This can most easily be demonstrated with a counter example

 f(z) =u(x,y)+hv(x,y), u(x,y) =v(x,y)=11+e−xe−y. (35)

The function is pictured in Fig. 2.

Let us verify that the function of (35) fulfills the GCR conditions

 u,x =−1(1+e−xe−y)2(−e−xe−y) =e−xe−y(1+e−xe−y)2, (36)

where we repeatedly applied the chain rule for differentiation. Similarly we obtain

 u,y=v,x=v,y=e−xe−y(1+e−xe−y)2. (37)

The GCR conditions (31) are therefore clearly fulfilled, which means that the hyperbolic function of (35) is holomorphic. Since the exponential function has a range of , the product also has values in the range of . Therefore the function has values in , and the components of the function of (35) have values

 0<11+e−xe−y<1. (38)

We especially have

 limx,y→−∞11+e−xe−y=0, (39)

and

 limx,y→∞11+e−xe−y=1. (40) Figure 2: Function u(x,y)=1/(1+e−xe−y). Horizontal axis −3≤x≤3, from left corner into paper plane −3≤y≤3. Vertical axis 0≤u≤1. (Figure produced with .)

The function (35) is representative for how to turn any real neural node activation function into holomorphic hyperbolic activation function via

 f(x)=r(x+y)(1+h). (41)

We note that in [3, 4] another holomorphic hyperbolic activation function was studied, namely

 f′(z)=11+e−z, (42)

but compare the quote from , p. 114, given in the introduction. The split activation function used in 

 f′′(x,y)=11+e−x+h11+e−y, (43)

is clearly not holomorphic, because the real part depends only on and not on , and the -part depends only on and not on , thus the GCR conditions (31) can not be fulfilled.

## 5 Geometric interpretation of multiplication of hyperbolic numbers

In order to geometrically interpret the product of two complex numbers, it proves useful to introduce polar coordinates in the complex plane. Similarly, for the geometric interpretation of the product of two hyperbolic numbers, we first introduce hyperbolic polar coordinates for with radial coordinate

 ρ=√|z¯z|=√|x2−y2|. (44)

The hyperbolic polar coordinate transformation  is then given as

1. , :

 θ=artanh(y/x),z=ρehθ,

i.e. the quadrant in the hyperbolic plane of Fig. 1 limitted by the diagonal idempotent lines, and including the positive -axis (to the right).

2. , :

 θ=artanh(y/x),z=−ρehθ,

i.e. the quadrant in Fig. 1 including the negative -axis (to the left).

3. , :

 θ=artanh(x/y),z=hρehθ,

i.e. the quadrant in Fig. 1 including the positive -axis (top).

4. , :

 θ=artanh(x/y),z=−hρehθ,

i.e. the quadrant in Fig. 1 including the negative -axis (bottom).

The product of a constant hyperbolic number (assuming )

 a =ax+hay=ρaehθa, ρa =√a2x−a2y,θa=artanh(ay/ax), (45)

with a hyperbolic number (assuming , ) in hyperbolic polar coordinates is

 az=ρaehθaρehθ=ρaρeh(θ+θa). (46)

The geometric interpretation is a scaling of the modulus and a hyperbolic rotation (movement along a hyperbola) .

In the physics of Einstein’s special relativistic space-time [11, 12], the hyperbolic rotation corresponds to a Lorentz transformation from one inertial frame with constant velocity to another inertial frame with constant velocity . Neural networks based on hyperbolic numbers (dimensionally extended to four-dimensional space-time) should therefore be ideal to compute with electromagnetic signals, including satellite transmission.

## 6 Conclusion

We have compared complex numbers and hyperbolic numbers, as well as complex functions and hyperbolic functions. We saw that according to Liouville’s theorem bounded complex holomorphic functions are necessarily constant, but non-constant bounded hyperbolic holomorphic functions exist. One such function has already beeng studied in [3, 4]. We have studied a promising example of a hyperbolic holomorphic function

 f(z)=1+h1+e−x−y, (47)

in some detail. The distinct notions of idempotents and divisors of zero, special to hyperbolic numbers, were introduced. After further introducing hyperbolic polar coordinates, a geometric interpretation of the hyperbolic number multiplication was given.

Hyerbolic neural networks offer, compared to complex neural networks, therefore the advantage of suitable bounded non-constant hyperbolic holomorphic activation functions. It would certainly be of interest to study convergence, accuracy and decision boundaries of hyperbolic neural networks with the activation function (35), similar to [3, 4].

## Acknowledgment

I want to acknowledge God :

In the beginning was the Word222Greek term: logos. Note: A Greek philosopher named Heraclitus first used the term Logos around 600 B.C. to designate the divine reason or plan which coordinates a changing universe. , and the Word was with God, and the Word was God. He was with God in the beginning. Through him all things were made; without him nothing was made that has been made. In him was life, and that life was the light of all mankind.

I want to thank my dear family, as well as T. Nitta and Y. Kuroe.

## References

•  K. Guerlebeck et al, Holomorphic Functions in the Plane and n-dimensional Space, Birkhauser, 2008, chp. 7.3.3.
•  S. Buchholz, G. Sommer,

A hyperbolic multilayer perceptron

, Proceedings of the International Joint Conference on Neural Networks, Como, Italy, vol. 2, 129/133 (2000).
•  T. Nitta, S. Buchholz,

On the Decision Boundaries of Hyperbolic Neurons

, Proceedings of the International Joint Conference on Neural Networks, IJCNN’08-HongKong, June 1-6, 2973/2979(2008).
•  S. Buchholz, PhD Thesis, A Theory of Neural Computation with Clifford Algebras, University of Kiel, 2005.
•  F. Catoni et al, The Mathematics of Minkowski Space-Time, Birkhauser, 2008.
•  Laplace’s equation, Wikipedia, accessed 24 August 2011, http://en.wikipedia.org/wiki/Laplace’s_equation
•  Wave equation, Wikipedia, accessed 24 August 2011, http://en.wikipedia.org/wiki/Wave_equation
•  Online 3D function grapher, http://www.livephysics.com/ptools/online-3d-function-grapher.php?
•  Split-complex number, Wikipedia, accessed 29 August 2011, http://en.wikipedia.org/wiki/Split-complex_number
•  Notes of collaboration with H. Ishi, Feb. 2011, p. 15.
•  C. Doran and A. Lasenby, Geometric Algebra for Physicists, Cambridge University Press, Cambridge (UK), 2003.
•  E. Hitzer, Relativistic Physics as Application of Geometric Algebra, in K. Adhav (ed.), Proceedings of the International Conference on Relativity 2005 (ICR2005), University of Amravati, India, January 2005, 71/90(2005).
•  The Bible, New International Version (NIV), The Gospel according to John, chapter 1, verses 1-4, http://www.biblegateway.com/
• 

Strong’s Bible lexicon entry G3056 for

logos, available online at Blue Letter Bible. http://www.blueletterbible.org/lang/lexicon/lexicon.cfm?Strongs=G3056&t=KJV