Distance function wavelets – Part III: “Exotic” transforms and series

W. Chen
Simula Research Laboratory, P. O. Box. 134, 1325 Lysaker, Norway
E-mail: wenc@simula.no
(8 June 2002)

Summary
Part III of the reports consists of various unconventional distance function wavelets
(DFW). The dimension and the order of partial differential equation (PDE) are first used
as a substitute of the scale parameter in the DFW transform and series, especially with the
space and time-space potential problems. It is noted that the recursive multiple
reciprocity formulation is the DFW series. The Green second identity is used to avoid the
singularity of the zero-order fundamental solution in creating the DFW series. The
fundamental solutions of various composite PDEs are found very flexible and efficient to
handle a broad range of problems. We also discuss the underlying connections between
the crucial concepts of dimension, scale and the order of PDE through the analysis of
dissipative acoustic wave propagation. The shape parameter of the potential problems is
also employed as the “scale parameter” to create the non-orthogonal DFW. This paper
also briefly discusses and conjectures the DFW correspondences of a variety of
coordinate variable transforms and series. Practically important, the anisotropic and
inhomogeneous DFW’s are developed by using the geodesic distance variable. The DFW
and the related basis functions are also used in making the kernel distance sigmoidal
functions, which are potentially useful in the artificial neural network and machine
learning. As or even worse than the preceding two reports, this study scarifies
mathematical rigor and in turn unfetter imagination. Most results are intuitively obtained
without rigorous analysis. Follow-up research is still under way. The paper is intended to
inspire more research into this promising area.

1

Keywords: distance function wavelets, potential DFW, multiple reciprocity DFW,
recursive multiple reciprocity, composite multiple reciprocity, differentiation smoothing,
Green second identity, dimension, scale, the order of PDE, multifractal, fractional
derivative, fractional Laplacian, fractional distance, fractional algebra, Possion kernel,
shift parameter DFW, kernel geodesic function, geodesic DFW, inhomogeneous,
anisotropic, distance variable, kernel distance sigmoidal function, machine learning.

1. Introduction
The report is the third in series [1,2] about my latest advances on the distance function
wavelets (DFW). One of the main motivations behind this research is to extend the DFW
to various scale-invariant potential problems where unlike the Helmholtz and convectiondiffusion problems, the standard scale parameter does not exist. In section 2, we instead
use the dimension as a substitute of the scale parameter to develop the space and timespace potential DFW’s, whose transforms involve continuous dimension variation (fractal
and multifractal). The links between these DFW’s and the Kontorovich-Lebedev and
Mellin transforms are pointed out. In section 3, it is noted that the boundary particle
method (BPM) [3,4] and recursive multiple reciprocity BEM have the DFW series
formulation, where the “scale parameter” is interpreted as the order of the high-order
general solutions and fundamental solutions of partial differential equations (PDE). Its
corresponding DFW transforms involve the fractional derivative and complex-order
derivative. The singularity of the zero-order fundamental solution is a troublesome
specter in creating the DFW series in terms of the order of PDE. Section 4 applies the
Green second identity to circumvent this thorny problem. It is also noted that the
recursive multiple reciprocity method can become very efficient via the scale
orthogonality. The fundamental solutions of various composite PDEs are also presented
to handle a wide range of problems. The technique is called the composite multiple
reciprocity. However, in some cases like the convection-diffusion equation, the
reciprocity principle does not hold. Thus, the differentiation smoothing is a better term

2

for this approach. The time-space DFW series is also proposed. It is discussed that the
composite multiple reciprocity DFW could lead to a boundary-only DFW technique not
only for numerical PDE but also for general data processing. In section 5, the underlying
connections between the essential concepts of dimension, scale and the order of PDE are
critically discussed through the analysis of the frequency-dependent attenuation of
acoustic wave propagation. Furthermore, section 6 employs the shape parameter of
various rotational and translate invariant MQ-type (the MQ is the abbreviation of the
multiquadratic) distance functions as the “scale parameter” to construct the pre-wavelets.
Section 7 mirrors the DFW on a variety of existing standard integral transforms and
series. In section 8, we introduce practically important geodesic DFW, which uses the
geodesic distance variable and will be potentially very useful in handling the anisotropic
and inhomogeneous problems. Section 9 employs the DFW and the related basis
functions to make the kernel distance sigmoidal functions for the artificial neural network
and machine learning. Finally, section 10 gives a brief remark on the DFW promise.
“Exotic” in the title means that the DFW’s introduced in this report use the
unconventional dimension, the order of PDE, shape parameter as a substitute of the scale
parameter. As or even worse than the preceding reports [1,2], this research lacks
mathematical rigor and is contemplated to be future-orientated in the best hope leading to
an enormous territory rich in open problems.

2. Potential DFW transform using dimension parameter
Unlike the Helmholtz and convection-diffusion problems, the kernel solutions of various
potential PDEs, e.g. the Laplacian, are scale invariant and have not explicit scale
provision within them. Notwithstanding, in terms of dimensional parameter as a
substitute of the scale parameter, this section aims to create the potential DFW via the
PDE kernel solutions of space and time-space potential problems, followed by section 3
to present the multiple reciprocity DFW using the order of PDE as the scale argument,
and then, section 4 derives the multiple reciprocity (MR) DFW series. Section 5 provides
a unified mathematical and physical view on dimensionality, PDE order and scale

3

parameters and their effects. The DFW using the dimension argument may find some
applications in high-dimensional problems such as the neural network.
2.1. Space potential DFW transform
Without a loss of generality, let us consider the fundamental solution of the Laplace
equation [5]

1
 − rk ,
S1

1

u L*n (rk ) =  − ln rk ,
 S 2 2−n
rk

− (n − 2 )S ,

n

n =1
n=2

(1)

n≥3

where rk = x − x k , Sn is the surface area of the unit sphere given by

S1=2,

Sn =

S2=2π,

2π n 2
.
Γ(n 2 )

(2)

Γ is the gamma function. We can verify (1) via the radial Laplacian
d 2 u L*n
dr 2

*

+

n − 2 du Ln
= −∆ i .
r dr

(3)

In terms of the fractal geometry and fractal Hausdorff dimension, n is not necessarily an
integer. We note that
1

 − S ln rk ,
2
u L*n (rk ) = 
rk2−n
−
,
 (n − 2 )S n

n=2

(4)
n≠2

4

could satisfy equation (3) with arbitrary real-valued n. It is the author’s guess that two
dimensions are special since the Laplacian happens to be the second order of partial
differential equation model. A fraction n may suggest the fractal field and potential or the
fractional Laplace operator [6,7].
The Laplacian is scale invariance, which happens to agree with the definition of the selfsimilarity and fractal dimension. In order to create the continuous DFW transform, it is
necessary to use the concept of continuous dimension variation instead of conventional
integer quantum jumps of dimensions. In terms of Laplacian fundamental solution (4), we
have the continuous DFW:

P (n,ξ ) = ∫

f (x )

IR n

u L*n ( ξ − x

ξ−x

)

−1

dx

(5a)

and

f (x ) =

1
CP

+∞

∫ ∫
0

IR n

P (n,ξ )u L*n ( x − ξ )ndξdn .

(5b)

where Cp can be determined by the continuous wavelets theory. It is known that negative
dimension occurs in the exploration of some hidden physics secrets. It is therefore
reasonable to extend dimension n to negative real number space. Furthermore, the
complex space has been presented in literature. That means n could even be a complex
number. It is also noted that the derivative of the Laplacian fundamental solutions of any
dimensions with respect to the distance has the same form
du L*n (rk )
drk

=

1
.
S n rkn −1

(6)

(6) can be used to replace the Laplacian fundamental solution in (5) to construct the
double layer potential DFW.

5

In the Helmholtz-Laplace transforms [1,2], we use the fundamental and general solutions
of the modified Helmholtz equation

w1 (µrk ) =

wn (µrk ) =

2

e − µrk ,

µ n −1 2
(2πµrk )−(n 2 )+1 K (n 2 )−1 (µrk ) ,
2π

wˆ 1 (µrk ) =

wˆ n (µrk ) =

µ1 2

µ1 2
2

(7a)

n≠1;

e µrk ,

µ n −1 2
(2πµrk )−(n 2 )+1 I (n 2 )−1 (µrk ) ,
2π

(7b)

(8a)

n≠1,

(8b)

where I and K denote the modified Bessel function of the first and second kinds. Instead
of using the parameter µ, we use the dimension parameter n as the “scale parameter” and
have the DFW transform for a suitable function f(x)
H (n,ξ ) = ∫

IR

f (x ) =

1
Cw

γ + i∞

∫γ ∫

−i∞ IR n

n

f (x )

wn (µ ξ − x

ξ −x

)

dx ,

H (n,ξ )wˆ n (µ x − ξ )ndξdn .

(9a)

(9b)

Note that the dimensionality here is the complex number. In contrast, the classic forward
and inverse Kontorovich-Lebedev transforms [8] are defined by

∞

L(µ ,α ) = ∫ K α (µx ) f ( x )
0

6

dx
,
x

(10a)

f (x ) =

1
iπ

∫

i∞

−i∞

αI α (µx )L(µ ,α )dα ,

(10b)

where the imaginary index α of Kα matches the dimensionality in the DFW transform (9).
As µ→0, the modified Helmholtz equation degenerates into the Laplace equation and the
Kontorovich-Lebedev transforms become the Mellin transforms, i.e.
∞

M (α ) = ∫ x α −1 f ( x )dx ,
0

f (x ) =

1 c +i∞ −α
x M (α )dα .
2πi ∫c −i∞

(11a)

(11b)

It is now clear that the Mellin transform has the lurking background of the Laplace
equation. In terms of the complex dimensionality, the inverse transform (5b) is
alternatively replaced by
f (x ) =

1
ˆ
CP

c + i∞

∫ ∫

c + i∞ IR n

P (n,ξ )u L*n ( x − ξ )ndξdn .

(12)

The Laplacian potential DFW transform (5a) and (12) may be considered the DFW
counterpart of the Mellin transform which finds applications in scale-invariant image and
speech recognization. By analogy with the Mellin-Fourier series [9], we can construct the
expansion series of the Laplacian DFW. It is the author’s guess that the existence
condition for the Laplacian DFW is similar to that of the Mellin transform, i.e.

∫

IR n

f ( x )u L*n ( ξ − x )dx p ∞ .

(13)

We may consider the Kontorovich-Lebedev and Mellin transforms are the degenerate
DFW transforms respectively using the solutions of the modified Helmholtz and Laplace
equations as the Hankle transform is to the Helmholtz-Fourier transform. In the later
section 3, we will show that the solutions of high-order Laplacian can also be used to
create the DFW transforms of this kind.
7

2.2. Time-space potential DFW transform

The time-space potential can be seen as an extension of the spatial potential under the
general Euclidean-like space. The origin of the universe, as stated by the big bang theory,
starts from a singularity (zero dimension) to time dimension to one spatial dimension to
two and three spatial dimensions. If we speculate this cosmos dimension evolution is a
continuous (or finite quantum jumps but integrable) process, the continuous time-space
DFW could come into play. Unfortunately, as far as the author knows, the temporal
fractal fundamental solution of potential PDE is not available albeit in principle existing
and solvable. So, the time-space potential DFW transforms developed below are in fact
the space potential DFW augmented with one time dimension. It is stressed that besides
the finite four time-space dimensions, the dimensionality could reach much higher and
even infinite in terms of other physical parameters. The dimensional DFW also has close
ties with the multifractal concept.
For time-dependent problems, we typically meet the so-called “heat potential” and “wave
potential” [5]. Consider the diffusion equation
∂u
− κ∇ 2 u = f (x, t ) ,
∂t

(14)

where κ is the conductivity parameter of media in terms of heat problems, the heat
potential is defined by
u ( x, t ) = ∫ ∫ n uh*n ( x − ξ , t − τ ) f (ξ , τ )dξdτ ,
t

0 IR

(15)

where f is a distribution that vanish for τ<0, and
uh*n ( x − ξ , t − τ ) =

2
H (t − τ )
− x −ξ
e
n2
[4πκ (t − τ )]

8

4κ (t −τ )

,

(16)

where H is the Heaviside step function, which is added to ensure the causality. In terms
of the heat potential (15), we can make the time-space DFW transforms

H (n,ξ ,τ ) = ∫

+∞

0

∫

IR n

f ( x, t )

u h*n (ξ − x,τ − t )

−1

d xdt

ξ−x

(17a)

and
f ( x, t ) =

1
CH

c + i∞ t

∫ ∫∫

c −i∞ 0 IR n

H (n,ξ ,τ )u h*n ( x − ξ , t − τ )ndξdτdn .

(17b)

It is known that the Black-Scholes PDE model for financial option is also a diffusion-type
equation usually involving multiple variables (up to thousands dimensions). The Monte
Carlo method is now the dominant simulation technique in this field. The dimensional
DFW may play a role for analyzing this high dimension problem. Through the
dimensional DFW transforms, we may find which dimensions in very high-dimensional
problems play the key role in particular objects so that a “dimension compression” may
be achieved.

If κ in (14) is purely imaginary such that κ =

ih
, where m is the mass of the quantum
2m

particle and h the Plank’s constant, equation (14) becomes the Schrodinger equation [10]
∂u h 2 2
ih
+
∇ u = f ( x, t )
∂t 2m

(18)

to define a particle density. u ( x, t ) ∆u denotes the probability of the particle being in
2

the neighborhood u(x,t) [10]. The fundamental solution of the Schrodinger equation is [5]

1+ i
u (x − ξ , t − τ ) = − H (t − τ )
h 2
*
Sn

9



m
 2πh(t − τ )



n2

e

im x −ξ

2

2 h (t −τ )

.

(19)

We have the corresponding DFW transforms

S (n,ξ ,τ ) = ∫

+∞

∫

IR n

0

V ( x, t )

u S*n (ξ − x,τ − t )

−1

ξ −x

d xdt

(20a)

and
V ( x, t ) =

1
CS

∞

t

0

0 IR n

∫ ∫∫

S (n,ξ ,τ )u S*n ( x − ξ , t − τ )ndξdτdn .

(20b)

Dimension n in (20) also can be the complex number. Now we go to the hyperbolic wave
equation
∂ 2u
− c 2 ∇ 2 u = f ( x, t ) ,
2
∂t

(21)

where c is the media wave velocity. The wave potential with density f is defined by [5]
Wn = u w* n * f ,

(22)

where * denotes the convolution operation and u w* n is the time-harmonic fundamental
solution of wave equation (21). In terms of (22), we can create the corresponding wave
potential DFW transforms in terms of dimensionality in the same way as did with the
heat potential.
We can also make the DFW series in terms of discrete dimension parameters. However,
the basis functions should not be singular, e.g. nonsingular general solution of the
Helmholtz equation. It is noted that the higher-dimension basis functions tend to be more
compactly supported. Thus, the corresponding interpolation matrices tend to be sparser in
the higher dimensions which offsets the so-called the curse of dimensionality.

10

3. Multiple reciprocity DFW transforms

This section is concerned with the DFW transforms using the order of PDE as a substitute
of the scale parameter. Since this kind of the DFW is essentially related to the boundary
particle method [3,4], a boundary-only numerical technique based on the multiple
reciprocity principle [11] to solve a PDE problem, we call it the multiple reciprocity
(MR) DFW. To clearly illustrate our idea, consider the following example without a loss
of generality

where ℜ

ℜ{u} = f (x ) ,

x∈Ω ,

(23)

u (x ) = R(x ) ,

x ∈ Su ,

(24)

∂u ( x )
= N (x ) ,
∂n

x ∈ ST ,

(25)

is a linear differential operator, x means multi-dimensional independent

variable, and n is the unit outward normal. The solution can be expressed as
u = u h0 + u 0p ,

(26)

where u h0 and u 0p are the zero-order homogeneous and particular solutions, respectively.
The multiple reciprocity method evaluates the particular solution in Eq. (26) by a sum of
higher-order homogeneous solution, namely,

∞

u 0p = ∑ uhm ,

(27)

m =1

where superscript m is the order index of the homogeneous solution. Thus, the BPM
solution of inhomogeneous equation (23) can be expressed as

11

∞

u = uh0 + u 0p = ∑ uhm .

(28)

m =0

We do not want to plunge the readers into details of the BPM and recursive multiple
reciprocity solution procedure which can be found in [3,4]. The final solution is given by
M

L

u ( xi ) = ∑∑ β km um# ( x − xk ) ,

(29)

m =0 k =1

where u m# is the general or fundamental solution of operator ℜm{} which denotes the mth order of ℜ{}, say ℜ0{}=ℜ{}, ℜ1{}=ℜ{ℜ{}}. For instance, the general and
fundamental solutions of the m-order Helmholtz equation ℜ m { } = (∇ 2 + γ ) are given by
m

u H# m (r ) = Am (γr )

− n 2 +1+ m

J n 2−1+ m (γr ) ,

(30a)

u H* m (r ) = Am (γr )

− n 2 +1+ m

H n(12) −1+ m (γr ) ,

(30b)

n

n

where Am=Am-1/(2*m*γ2), A0=1, and n is the dimension of the problem; J represents the
Bessel function of the first kind and H(1) is the first kind of the Hankle function. It is
noted that the recursive multiple reciprocity method described in [3,4] can be directly
used with the boundary integral formulation in the same way as in the collocation BPM
and will greatly improve computing efficiency of the so-called multiple reciprocity BEM.
It is obvious that (29) is a DFW series in terms of integer “scale parameter” m. To
introduce the continuous DFW transforms, we need to use the concept of the fractional
order of derivative [12], i.e., m could be a real number rather than an integer. Thus, we
could create the multiple reciprocity DFW transforms with the high-order fundamental
solutions of the operator ℜ{}
12

M n (m,ξ ) = ∫

f ( x )u H* m (ξ − x )d x

(31a)

M n (m,ξ )u H* m ( x − ξ )mdξdm ,

(31b)

IR n

n

and

f (x ) =

1
CM

+∞

∫ ∫
0

IR n

n

The multiple reciprocity DFW transforms (31) reveal the distribution of function f(x) with
respect to the order of the Helmholtz operator. Observing the MR DFW series (29), the
MR DFW can also alternatively be given by
Π n (m,ξ ) = ∫

IR

n

{(∇

}

+ λ2 ) f (x ) u H* 0 (ξ − x )d x

(32a)

Π n (m,ξ )u H* 0 ( x − ξ )mdξdm ,

(32b)

2

m

n

and

f (x ) =

1
CΠ

+∞

∫ ∫
0

IR n

n

where u H* 0 denotes the fundamental solution of the zero-order Helmholtz ∇ 2 + γ 2 . (32)
n

may be better in some practical uses than (31).
[11,13,14] give the high-order general and fundamental solutions of the Laplace,
modified Helmholtz, convection-diffusion, Berger and Winkler equations. Consider the
scale-invariant Laplacian, we have its 2D and 3D high-order fundamental solutions [11]:

u L*m (r ) =
2

1 2m
r ( Am ln r − Bm ) ,
2π

(33)

where Am +1 = Am 4(m + 1) and Bm+1 = (Bm (m + 1) + Bm ) 4(m + 1) ;
2

2

u L*m (r ) =
3

1
r 2 m −1 ,
4π (2m )!

13

(34)

where m is the Laplacian order. In section 2, we have constructed the potential Laplacian
DFW in terms of dimensionality parameter. Here we instead have the Laplacian MR
DFW transforms
Gn (m,ξ ) = ∫

IR n

f (x )

u L*m ( ξ − x
n

)

−1

ξ−x

dx ,

(35a)

and
f (x ) =

1
CG

Gn (m,ξ )u L*m ( x − ξ )mdξdm ,

c + i∞

∫ ∫

c −i∞ IR n

n

(35b)

where complex order m corresponds to the complex-order Laplacian. The alternatives of
the MR DFW (35) are given by

Gn (m,ξ ) = ∫ n ∇

2m

IR

f (x )

u L*0 ( ξ − x
n

)

−1

ξ−x

dx ,

(36a)

and
f (x ) =

1
CL

c + i∞

∫ ∫

c −i∞ IR n

Gn (m,ξ )u L*0 ( x − ξ )mdξdm .
n

(36b)

Comparing the MR Laplacian DFW transforms (35) and (36) with the preceding potential
Laplacian DFW transforms (5) and (12), one sees that both look very much similar. Zahle
[6] pointed out that the spectral dimension of the Laplacian agrees with the Hausdorff
dimension underlying fractal through the fractional derivative. In other words, the
dimension and the order of the Laplace equation can be reciprocally connected under
certain conditions through the Green function of the Reize potential on a fractal. For a
detailed discussion see section 5. We also can develop the high-order potential Laplacian
DFW transform in terms of dimension parameter, i.e.
P (n,ξ ) = ∫
m

IR n

14

f (x )

u L*m ( ξ − x
n

ξ−x

)

−1

dx ,

(37a)

f (x ) =

1
CP

c + i∞

∫ ∫

c + i∞ IR n

P m (n,ξ )u L*m ( x − ξ )ndξdn .

(37b)

n

In particular, the Mellin transform can be seen as the degenerate form of the onedimension high-order Laplacian DFW transforms (37).

4. Composite multiple reciprocity DFW series

The singularity of the fundamental solutions of multidimensional PDEs does not cause
troubles in the creating of various continuous DFW transforms, but makes it seemingly
impossible to construct a complete multiple reciprocity DFW series. We have good
reason to allot a separate section to deal with this thorny issue. The basic strategy is to
use the Green second identity to evaluate the zero-order singular fundamental solution
term and leave the other nonsingular high-order fundamental solution terms to be handled
in the normal fashion. The title “composite multiple reciprocity” suggests involving more
than one PDE operators in some cases, e.g. composite Laplace and Helmholtz operators,
and thus is different from the traditional multiple reciprocity [11] where the mere one
type of PDEs of different orders is concerned.
4.1. Laplacian MR DFW series

To eliminate the edge effect of the Helmholtz-Fourier (HF) series, section 3.3 of [1]
developed the HF series of the form:

∞

∞

f ( x ) = f 0 ( x ) + ∑∑ β jk φ n (λ j x − x k ) ,

n≥2,

(38)

j =1 k =1

where f0(x) and the sum term respectively represent the solutions corresponding to
nonzero boundary condition and zero boundary condition, and n is the dimensionality and

15

λ the eigenvalues; φn denote the HF basis functions based on the nonsingular kernel
general solution of the Helmholtz equation. f0(x) is evaluated by a Laplace equation
(degenerate Helmholtz equation with a zero eigenvalue)

∇ 2 f 0 (x ) = 0 .

(39)

In terms of the Green second identity, we have

(

∂u L*0 x − x j
 ∂f (x j ) *
u L0 x − x j − f (x j ) n
f 0 ( x ) = ∫ n −1 
n
S
∂
∂n
n


(

)

)

dx j ,


(40)

where u L*0 is the fundamental solution of the zero-order Laplacian ∇2u, and Sn-1 is the
n

surface of finite domains. (40) can be easily evaluated by the boundary element method
(BEM) with boundary f(x) data, and then f0(x) at any inner locations can be calculated via
(40) again. It is well known that the boundary is the dominant features in various data
sets. The Green second identity should be the right tool to capture it. The boundary knot
method (BKM) [3,13] is an alternative to the BEM for this task.
[56] discusses the essential concept of the complete fundamental solution. For instance,
the 2D Laplacian has the essential fundamental solution
fundamental solution

−1
ln(r ) and the complete
2π

−1
(ln(r ) + C ) , where C is a constant. The standard BEM only uses
2π

the former. For the DFW series, this may lead to the missing of the constant term and
incomplete basis functions. Thus, in the following DFW series, we may need to add a
constant term whenever the completeness concern arises, or instead use the complete
fundamental solution in (40). On the other hand, the complete fundamental solution also
provides an alternative explanation of the multiquadratic (MQ) type radial basis

16

functions. For instance, the shifted TPS ln (rk2 + ck2 ) has something to do with the 2D
Laplacian complete fundamental solution
inverse MQ 1

−1
(ln(r ) + C ) , while the MQ
2π

rk2 + c k2 and

rk2 + ck2 may respectively underlie the 1D and 3D Laplacian complete

fundamental solutions

1 1
1

r + C and
 +C .
4π  r
2


The above HF series provides an insight to construct the multiple reciprocity Laplacian
DFW series

∞

∞

f ( x ) = f 0 ( x ) + ∑∑α mk u L*m ( x − x k ) , n≥2,
m =1 k =1

n

(41)

where f0(x) is evaluated in the same way as (40), and u L*m is the fundamental solution of
n

the m-th order Laplacian ∇2(m+1)u. Observing (33) and (34) of the 2D and 3D Laplacian
high-order fundamental solutions, it is found that when the order m≥1, these kernel
solutions are no longer singular at the origin. Thus, (41) kills the singularity culprit and is
called the multiple reciprocity Laplacian DFW series. Note that except of evaluating f0(x)
with the boundary elements or BKM, the nodes could be anywhere inside domain and on
the boundary. (41) can be understood splitting f(x) into two parts of f0(x) and the
remainder, where the latter is approximated by the high-order Laplacians.
The MR Laplacian DFW series could be an efficient tool to deal with the aperiodic
signals which do not exhibit the periodic behavior. It is also very interesting to note that
the order m of the Laplacian plays an analogous role of the polynomial order in the
standard polynomial interpolation. In fact, expansion series (41) is a multidimensional
analog of the univariate polynomial approximation, where the order of monomials is
interpreted as the order of the fundamental solution of the one-dimension high-order
Laplacian. It is worth pointing out that u L*m of different orders are orthogonal, which
n

17

greatly reduces computing effort. The Gram-Schmidt orthogonality method may be
helpful to enforce the orthogonalization of the Laplacian basis functions of the same
order over translates. The translation invariant polynomial DFW series presented in [2]
can also be efficiently calculated in the same recursive way as the MR DFW series for
data processing and numerical PDEs.
4.2. Composite multiple reciprocity DFW series with the solutions of composite PDE

[1] noted that the nonsingular general solution of the convection-diffusion problem was
not suitable to handle exterior unbounded domain problems, while its fundamental
solution capable for this task could not be used to create the DFW series due to its
singularity at the origin. A solution of this perplexity is to construct the composite MR
DFW series with the high-order solutions of a composite PDE.
To illustrate the basic strategy clearly, let us start with the equation
v
∇ 2 (D∇ 2 u + v • ∇u − ku ) = − ∆ i ,

(42)

v
where v denotes velocity vector, D is the diffusivity coefficient, and k represents the
reaction coefficient. (42) combines the convection-diffusion and Laplace equations. Its
fundamental solutions are
v

v ⋅( x − xk )
−

−1 
 ln (rk ) + e 2 D K 0 (ρrk ) ,
u (rk ) =

2π 

*
C20

−
rk2−n
u (rk ) =
+e
(n − 2 )S n (1)
*
Cn0

v
v ⋅( x − xk )
2D

(2πρrk )− (n 2 )+1 K

where dot denotes the inner product of two vectors;

18

2π

( n 2 )−1

(43a)

(ρrk ) ,

3≤ n, (43b)

1

 vv  2 k  2
 +  .
ρ = 

D
 2 D 


(44)

The fundamental solution (43) is singular at the origin. The high-order fundamental
solutions of (42) are a sum of those of the convection-diffusion equation and the Laplace
equation. It is found [14] that the high-order fundamental solutions of the convectiondiffusion equations except in 2D case are also singular in contrast to the nonsingular
high-order fundamental solutions of the Laplacian. Due to this observation, we construct
a composite multiple reciprocity series for a continuously differential function Q(x)

M

N

Q (x ) ≅ Q0 ( x ) + ∑∑ c mk u L*m ( x − x k ) ,
m =1 k =1

n

n≥2,

(45)

where Q0 satisfies (42) and is evaluated by the Green second identity as in (40). In data
processing, the Dirichlet boundary data are often only accessible. We can use the Green
second identity twice, namely, first with the Laplacian to get the Neumann data, and then
with the fundamental solution (43) for all necessary high-order boundary data. The
expansion (45) is in fact equal to the iterative differentiation of function Q(x) by

v
∇ 2 M (D∇ 2 + v • ∇ − k )Q = ∇ 2 M F ( x ) = q( x ) ,

(46)

where F(x) is the inhomogeneous term of the convection-diffusion equation. If the
residues q(x) continuously tends to zero with increasing m, (45) will converge. It is noted
v
that expansion (45) includes the direction vector v and thus is also well suitable to
handle the track data. If Q(x) is due to a quasi-periodic source P(x) consisting of a finite
number of periodic components [15], a composite operator of the convection-diffusion,
Laplace, and Helmholtz equations of the form

v
∇ 2 M (∇ 2 + λ12 )K(∇ 2 + λ2S )(D∇ 2 + v • ∇ − k )Q = ∇ 2 M (∇ 2 + λ12 )K(∇ 2 + λ2S )P ( x ) = g ( x ) (47)

19

will be preferred, where g(x) should be smooth and close to a constant or zero. The
corresponding DFW series is

M

N

S

N

Q (x ) ≅ Q0 ( x ) + ∑∑α mk u L*m ( x − x k ) + ∑∑ β ik u H# 0 (λi x − x k ) ,
m =1 k =1

n

i =1 k =1

n

n≥2,

(48)

where u H# 0 is the general solution of the zero-order Helmholtz equation. The strategy can
n

be extended to handle other problems via a varied combining of different PDEs. We may
call it the differentiation smoothing. Following are some composite PDE examples of
maybe highly interest:
∇ 2 M (∇ 2 ± µ 2 )u ,

(49a)

∇ 2 M (∇ 4 ± κ 2 )u ,

(49b)

∂N
∂t N

 1 ∂

 2 − ∇ 2 u ,
 η ∂t


(49c)

∂N
∂t N

 1 ∂2

 2 2 − ∇ 2 u .
 c ∂t


(49d)

(49a) under m=1 is the Burger equation governing large deflections of plate without
planar displacements which combines the Laplacian and the modified Helmholtz
equation. [13] gives some composite PDEs and their high-order fundamental and general
solutions. (49c,d) are the temporal multiple reciprocity method which aims to eliminate
the inhomogeneous terms due to time-dependent sources. We have a great deal of
freedom to combine any PDEs such as the convection-diffusion, Helmholtz, modified
Helmholtz, Laplace, elastostatics, time-dependent diffusion and wave and transport
equations to tackle a very broad variety of type data processing and PDE problems. The
20

corresponding fundamental solutions will be simply a sum of all individual fundamental
solutions. The composite multiple reciprocity approach is also very useful in the kernel
distance function [3] to handle domain integral and particular solution of PDE subjected
to different types of inhomogeneous function terms.
4.3. Time-space MR DFW series

The preceding DFW transform and series are only concerned with spatial variables and
equilibrium data. It is feasible to develop the time-space MR DFW series and transforms
due to the time-harmonic solution of the diffusion and wave equations, which behaves
similarly as the spatial Laplacian harmonic functions. For brevity, we simply display an
analogous time-space MR DFW series for wave problems below

∞

∞

(

)

f ( x ) = f 0 ( x, t ) + ∑∑α mk u L*m  c 2 (t − t k ) − rk2  H c 2 (t − t k ) − rk2 , n≥2,
n 

m =1 k =1
2

2

(50)

where c is the wave velocity, H is the Heaviside step function augmented for the
causality, and rk = x − x k . f0(x,t) can be evaluated by the time-space boundary element
as in (40).
4.4. Boundary-only MR DFW for data processing and numerical PDE

Inverse problems are much often common in engineering computation than direct
problems. The boundary element method is found advantageous over the domain-type
numerical techniques such as the finite elements and finite differences in handling inverse
problems since the boundary data in most cases are dominant in determining the
systematic behavior and much more easily accessible than the inside-domain data.
However, we unfortunately find that the BEM is seldom used in practice to handle
inverse problem and data processing problems. This may be due to a few factors: 1)
singular integration, 2) full and asymmetric interpolation matrix, 3) costly domain
integral, 4) mathematical difficulty. The preceding composite MR DFW series seems a
21

promise to remedy all these defects of the BEM and could be very efficient if the proper
composite PDE is found. As well, the boundary particle method may be more promising
to handle data processing by only using the boundary data. [3,4,14] shows that the BPM
could be symmetric and nonsingular and truly boundary-only for inhomogeneous PDE
problems. The wavelet threshold shooting of the composite MR DFW series will lead to a
sparse interpolation matrix. On the other hand, the advantages of the MR-BEM over the
BPM are twofold: 1) the integral formulation of the MR-BEM is generally more stable
than the collocation BPM; 2) the MR-BEM can use the singular fundamental solution
without requiring the fictitious boundary outside physical domain. Thus, the BPM and
BEM have their respective relative strengths and weaknesses. In conclusion, the BEM
and BPM with the help of the recursive multiple reciprocity DFW series are two
boundary-only techniques not only for PDE problems but also for data processing.

The basic procedure in the multiple reciprocity solution of inhomogeneous PDE is a
reversely recursive iteration [3,4]. Different from the solution of PDE problems, the data
processing, however, is a directly recursive iteration, namely, we start to interpolate the
accessible boundary data with the zero-order fundamental solution or general solution,
and attain the derivative boundary data of one order higher, and then interpolate these
obtained data with high-order fundamental or general solution and evaluate the higherorder derivative boundary data, and repeat this processing upward until the residues norm
is under acceptable magnitude. For details on the BPM and recursive multiple reciprocity
method see [3,4,14]. In this way, we can effectively interpolate the whole domain data
only by the boundary data. If the composite multiple reciprocity method presented in
section 4.2 is used, we need to evaluate a finite number of interpolation matrices
corresponding to each different element PDE. In addition, [57] presents a generalized
boundary element method to transfer the domain integral caused by the inhomogeneous
media to the boundary integral. The basic strategy behind this method can be applied
with the DFW in handling the inhomogeneous media problems.

22

One may argue that the MR DFW series using the solutions of high-order PDE may not
be suitable for discontinuous problems, e.g. functions F(x) in (46) and Q(x) in (47)
exhibit discontinuous or low continuous property. The author thinks that this
discontinuity issue could be solved by the subregion approach. In cases of concentrated
loading, the distance function can be used in the same very simple fashion as in the BEM
[16]. It is also stressed that the MR DFW as well as other DFW’s are not necessarily
global approximation in correspondence to local and global boundary integral equations.

5. Connections among dimension, scale and PDE order

In many scientific disciplines, the dimension, scale (frequency) and PDE order are treated
as independent parameters. Observing the dimension, PDE order and scale DFW’s
developed in the present report series, one may wonder if these three basic parameters
exhibit inherent connections just like what relationships between mass and energy, time
and space are. Recent years witness intense research on applications of fractional
derivatives [12], fractional Fourier transform [17] and fractal to many real world
problems. Now it is known that there is a definite connection between fractional Fourier
and fractional derivative, namely,
 ∂ s p (t ) 
s
 = (− iω ) P (ω ) ,
FT+ 
s
 ∂t 

(51)

where FT+ indicates the Fourier transform, ω represents frequency, and s is the derivative
order, not necessarily an integer, and P(ω) the Fourier transform of p(t). In the following
we will establish the links between fractional derivative and fractal through the analysis
of the acoustic frequency dependent attenuation. The power law of attenuation [18] is
described by
E = E0 e −α (ω )Λx ,

23

(52)

α (ω ) = α 0ω y ,

y∈[0,2],

(53)

where E denotes signal energy, α0 and y are media-dependent parameters obtained
through a fitting of measured data. Note that y is real. (53) can be restated as

y=

ln α (ω ) α 0
.
ln ω

(54)

(54) can be seen as a self-similar object, i.e. scale invariance [19]. Self-similar fractals
with parameters N and n are normally described by a power law
N = qn ,

(55)

ln N
ln q

(56)

where
n=

is known as the Hausdorff dimension. Note that q here could be any physical parameters
besides time and space. Comparing (54) with (56), it is clear that y represents fractal
dimensionality in terms of frequency. In other words, the power law attenuation (54) on
different frequencies underlies the invariant parameter y. The author is very curious why
fractal is between 0 and 2.
[20] presented the fractional temporal derivative PDE model to describe the power law
(53) of frequency dependent attenuation:

∇2 p =

1 ∂ 2 p 2α 0 ∂ ∂ y p
,
+
c 2 ∂t 2 c1+2 y ∂t ∂t y

24

(57)

where p denotes acoustic pressure and c is acoustic speed. When y=0, (57) is reduced to
the standard damped wave equation. In 1D case, y=2 leads (57) to an acoustic equation
very similar to equation (13) in [21]. Here we find fractal y is mirrored by the y-order
temporal derivative. Namely, (57) suggests that y represents temporal fractal. However, it
is observed from (54) that y is independent of frequency (time scale) ω. It is therefore
much reasonable to think that y may in fact mean the spatial fractal. For example, y varies
with different human body tissues, which have different spatial microstructures. [20] also
proposed a spatial fractional partial derivative model:
1 ∂ 2 p 2α 0 ∂ y
∇ p = 2 2 + 1+ y
∇ •p.
c ∂t
c ∂t
2

(58)

When y=2, (58) is the known augmented wave equation. It is not difficult to derive the
fundamental solution of the above fractional wave equation via the Fourier transform.
The absolute value of the fractional derivative in (58) can be calculated via the positive
definite Laplacian, i.e.
∇ y • p = (− ∇ 2 p )

y 2

.

(59)

Both attenuation temporal model (57) and spatial model (58) involving fractional order
derivatives show that the damping behavior has much to do with the fractal structure of
media, while the differential orders of inertia and diffusion terms are independent of
fractional dimensional effect, no matter either integer or fractal dimensions. It is clearly
recognized from the above analysis that the dimensional effect on some physical
mechanism, e.g. frequency-dependent attenuation, can be exactly reflected by the
corresponding derivative order in the PDE model. Now we can conclude that the
dimension, scale (frequency) and derivative order are inherently related to describe the
attenuation mechanism.
To more clearly clarify the connection between the dimension and derivative order, let us
consider the Riesz potential [6,7]
25

I s f (x ) =

f (ξ )
Γ((n − s ) 2 )
dξ ,
n
s2 s
n−s
∫
IR
π 2 Γ(s 2 )
x −ξ

(60)

where s is the order of Riesz potential, and n is the dimension. By analogy with the
Riemann-Liouville definitions of the fractional derivative [22], we define the fractional
Laplacian

Lns − y { f ( x )} = Lns + y { f ( x )} =


f (ξ )
Γ((n − s ) 2 ) 2 
d
ξ
∇

,
π s 2 2 s Γ(s 2 ) ∫IR n x − ξ n − s + y 

0p y p 2.

(61)

(61) clearly shows that the dimension and the order of the Laplacian are inherently
connected. By using the multiple reciprocity expansion series, we can have the expansion
definition of (62) just like Caputo’s expansion definition of the fraction derivative. Two
intuitive definitions of the fractional distance variable are given by
1s

s
 n
r =  ∑ (x i − x ki )  ,
 i =1

s
k

(62a)

and
s− y 
 n
rks − y =  ∑ (x i − x ki ) 
 i =1


1 (s − y )

,

(62b)

where n is the topological dimensionality, and s could even be a complex number. It is
noted that s=2 leads to the Newtonian potential and Euclidean distance.
Varying y means the multifractal, i.e. continuous or quantum dimension variation, which
is visible in some physics problems such as acoustic energy absorption variations over
different human body tissues, while a broadband excitation will cause multiscale
frequency components. Therefore, a complete description of these complex physical
phenomena should involve multifractal, multiscale and multivariate (3M) systems. It is

26

also interesting to see from the attenuation power law (53) that the time scale (frequency)
and space dimension are reflected on media acoustic attenuation as what the general
relativity reads dependence of time and space on mass.
[20] also noted that the algebra counterpart of the fractional derivative and fractal should
be the fractional order of a matrix. For instance, the finite element (FE) numerical
discretization (59) is given by

(− ∇ p )
2

y 2

v
= Ay 2 p ,

(63)

v
where A is the positive definite FE interpolation matrix of the Laplacian, and p is the

discrete value vector of p. From this point of view, the fractional Haar transforms should
be the fractional order of the original transform matrices, where the fraction may indicate
the fractal dimension or the fractional order derivative. Similarly, we can get the
fractional Welsh and Hadamard transforms corresponding to a fractal. Now we have a
complete theory of fractional calculus, fractal geometry, and fractional algebra in place to
model and analyze complex mathematical physics problems.
Now we try to use the above insights to analyze the electrical current in a transmission
line governed by the telegraph equation [10]

∇2u = C

∂ 2u
∂u
+ (CR + GL ) + GRu ,
2
∂t
∂t

(64)

where C, G, R, and L are the capacitance, leakage resistance, resistance, and selfinductance per unit length, respectively. These different resistances in the electric current
play the same role as does the viscoelastic damping in acoustic wave propagation and
mechanical vibration. By analogy with the previous acoustic attenuation PDE model (58),
(64) may be reformulated as

27

∇2u = C

∂ 2u
∂
+ γ ∇τ u + GRu ,
2
∂t
∂t

(65)

where τ is a positive real number to describe the resistance behavior of media and γ the
alternative collective parameter of C, G, R, and L. (65) is different from the fractional
telegraph equation given in [23] in that the spatial fractional derivative is used instead of
the temporal fractional derivative. It is interesting to imagine that the electric current
resistance may obey a power law similar to (52) and (53), where a temperature-related
parameter may act like the frequency in damped acoustics.
Furthermore, [24] presented the new complex-order derivative model for frequency
dependent storage and damping loss mechanisms more accurately. For example, the
following force-displacement model was given for the viscoelastic vibration [24]

q(t ) + λ

d r q(t )
dq(t )
=ν0
,
r
dt
dt

(66)

where parameters λ, r and υ0 are all complex-valued. The Fourier transform of the
complex-order derivative is defined by [24]
 d r q(t )
r
= (iω ) FT+ [q(t )] .
FT+ 
r 
 dt 

(67)

(67) is valid for Re(r)>0 and for n times differentiable function q(t) where n-Re(r)>0.
According the preceding analysis, the complex-order r may mean the complex
dimension. This naturally raises a question about the definition of the distance variable
under complex dimension. In research of wave equation
1
utt = u xx ,
c2

28

(68)

[25] introduced the pseudo-Euclidean distance in terms of Minkowisk 4-dimension space
as opposed to the Euclidean distance. Let
s = ict ,

(69)

we have the Laplacian-alike wave potential
u xx + u ss = 0 .

(70)

Its Euclidean distance variable is

rˆk =

(x − xk )2 + (s − sk )2

=

(x − xk )2 − c 2 (t − tk )2

.

(71)

The author wonders whether the fundamental solution of the complex-order PDE can be
expressed in terms of the pseudo-Euclidean distance.

6. Shape parameter DFW’s

The MQ

2

rk2 + ck2 and inverse MQ [26], Gaussian e rk

ck2

, the thin plate spline (TPS)

rk2 m+1 and rk2 m ln rk [27] as well as the shifted TPS rk2 m ln (rk2 + ck2 ) [28], a MQ-type TPS,
are among the most popular distance functions (radial basis function, RBF), where ck is
called the shape parameter, a pseudo-scale parameter. In general, replacing the Euclidean
distance variable rk by the MQ in the rotational invariant fundamental solution and
general solution of PDEs will yield a variety of the MQ-type kernel distance functions
[3]. Enormous numerical experiments [29,30] show that the MQ-type distance functions
can achieve spectral accuracy if the shape parameter ck are optimized. Despite intense
research has been devoted to analyzing and determining of the shape parameter,
unfortunately, ck is found to be problem dependent and there are not general approaches
optimally determining it a priori. But nevertheless [31] proves that the MQ and TPS are
29

conditional positive definite function and [32] shows the MQ could be seen as the prewavelets.
On the other hand, Hon and Wu [33] applied the translate-invariant 2D Laplacian
harmonic function e −αx

2

+αy 2

cos 2αxy as the basis function to construct very simple and

efficient numerical scheme for solving the 2D Laplace problems. Just like the MQ, their
Laplacian harmonic function also requires determining a problem-dependent parameter α
to get the optimal accuracy and convergence speed. In this study, α is also interpreted as
the shape parameter. Now we have rotational and translation invariants Laplacian
distance functions, both of which include a shape parameter. The purpose of this section
is to create these two-type shape parameter distance function wavelets.
6.1. Rotational invariant shape parameter DFW

Beside the known MQ, the Possion kernel [34] is also the MQ-type distance function:
Psn (x − x k , s ) =

where n is dimensionality, and x∈Rn.

(x−x

cn s
2
k

+s

2

)

(n +1) 2

, s>0,

(72)

∫ P (x )dx = 1 for all s>0, which mimics the scale
n
s

function of the multiresolution analysis (MRA). It is, however, stressed that the DFW is
not a MRA. Given a function f ∈Lp(Rn) (p≥1), the solution of the Dirichlet problem with
boundary value f is [35]
u ( x, s ) = (Psn * f )( x ) =

2s

ω n +1

∫

IR

n

(ξ − x

f (ξ )
2

+ s2

( n +1) 2 dξ

)

.

(73)

It is worth mentioning that since s and ck play the analogous role under the same potential
theory backdrops, the Possion kernel gives an explicit physical explanation of the shape
parameter of the MQ. It is also obvious
30

u ( x,0) = f ( x ) .

(74)

A basic result of Littlewood-Paley theory [34] is that
f

p

≈ g1 ( f ) p ,

1<p<∞,

(75)

where
12

 ∞ ∂u ( x, s ) 2 ds 
 .
g1 ( f )( x ) =  ∫ s
 0

s
s
∂



(76)

Based on the above observations, now we construct the Possion DFW transform
H (ξ , s ) = ∫ n ∇ 2 Psn (ξ − x, s ) f ( x )dx .
IR

(77)

Parameter s is understood the scale parameter in the sense of the pre-wavelets. Unlike the
MQ, the Possion kernel function (72) has a rapid decay at infinity and is unconditional
positive definite function. As such, the MQ itself can also be employed to construct the
DFW. Very interestingly, [36] uses the complex shape parameter to stabilize the illconditioning MQ interpolation matrix. This reminds us that we could have the complex
shape parameter DFW, i.e. s in (77) is the complex number.
To get the compactly supported DFW, by analogy with the RBF wavelets given in [37]
we introduce the following basis function
 Psn ( x − x k , s ),
 0,

ψ ( x − xk , s ) = 

31

x − x k p s,
x − x k ≥ s.

(78)

∇ 2ψ could serve the compactly-supported DFW basis function. In addition, ∂Psn ∂s
may be also a good DFW basis function. [38] researches the MQ behavior with an
infinite shape parameter.
The classic Gaussian RBF can be stated in the following fashion

φ (rk ) =

1
2 πβ

2

e −rk

4β

,

(79)

where shape parameter β is related to the time translate in the heat potential (16). φ(rk,β)
approaches the Dirac function δ(rk) as β→0 [5]. Thus, the shape parameter in (79) is in
fact a special form of the preceding time-space heat potential DFW (17). As in the MQ
case, β can be the complex number.
The diffusion radial basis function presented in [3]

ϕ (rk ) = e −αr ,
k

(80)

is related to the fundamental solution of the 1D diffusion equation. Due to its simplicity,
it is very convenient to use (80) creating the shape parameter prewavelet. It is also
straightforward to have the corresponding shape parameter DFW series. [39] gave the
MQ DFW series. For the Possion kernel DFW, we have the corresponding series
f ( x ) ≅ ∑∑ A jk Psn ( x − x k , s j ) .
N

M

j =1 k =1

6.2. Translation invariant shape parameter DFW

Before carrying on, we list some lemmas [33,35] below:

32

(81)

Lemma 7.1 The shift of harmonic solution (Laplacian solution) is a harmonic function.

Any linear combination of harmonic functions is also a harmonic function.
Lemma 7.2 The derivatives of any harmonic function are also harmonic functions.
Corollary 7.1. A sequence of harmonic functions on Ω which converges locally

uniformly on Ω.. Then, the limit function is again harmonic.
[33] uses the 2D translate invariant harmonic function

ϕ ( x − x k , y − y k ) = e −α ( x − x

k

)2 +α ( y − yk )2 −2 iα ( x − xk )( y − yk )

(82)

with a shape parameter α for numerical solution of the Laplace equation. In fact, there are
numerous such translate invariant harmonic functions, e.g.

ψ ( x − x k , y − y k ) = eα ( x − x

k

)+iα ( y − yk )

.

(83)

For higher-dimension Laplacian, we also have similar harmonic functions. These
translate invariant harmonic functions can be used as the basis function to construct the
DFW transform and series for the function bounded within a finite domain. The author
has not a clear idea of physical background of (82) and (83) and wonders if there are the
underlying connections between the shape parameters of the rotational and translation
invariant Laplacian basis functions.
[2] introduced the harmonic polynomial DFW series using the translate invariant
monomial solutions of the high-order Laplacian. It is worth mentioning that such
polynomial DFW is not compactly-supported if we do not enforce some additional
condition as did (78) to the Possion kernel.

33

7. DFW transforms versus common integral transforms

Besides the foremost outstanding Fourier and Laplace transforms and series, there are
many other theoretically and practically important transforms such as the Hartley,
Hilbert, Mellin, Abel, Radon, Riemann-Lioville, Weyl, inverse scattering transforms, and
their corresponding expansion series, some of which use the eigenfunctions of PDEs. [1]
has developed the Helmholtz-Fourier and Helmholtz-Laplace DFW’s as the counterparts
of the Fourier and Laplace transforms. Section 2 introduced the DFW correspondences of
the Mellin and Kontorovich-Lebedev transforms. This section will try to build the DFW
correspondences of a few more normal transforms. Furthermore, the DFW methodology
might be extendable to group transforms of irregular high-dimensional domains.

It is noted that both the integral equation and integral transform have the similar
expression in terms of the kernel function w(t,a), i.e.
g (a ) = ∫ f (t )w(a − t )dt .
b

(84)

a

Similarly, the distance function wavelet transforms and Green integral are expressed as
w(λ ,ξ ) = ∫

IR n

f ( x )g (ξ − x, λ )dx

(85)

Thus, it is not difficult to make the DFW transforms in connecting to the traditional
integral transforms if we know the kernel solution of partial differential equation behind
the latter. As examples, we first intuitively conjecture a DFW transform mimicking Weyl
transform [40] without any mathematical analysis:

Wn (ξ , γ ) =

[

2Γ(n 2 )
2
f (x ) ξ − x − γ 2
n
∫
 n − 1  IR
π Γ

 2 

34

]

(n −3 ) 2

(

2

)

ξ − x H ξ − x − γ 2 dx ,

(86)

where H is the Heaviside step function. (86) looks like the wave potential DFW
transforms discussed in section 2. In terms of the Hilbert transform, a DFW transform is
proposed by

H n (ξ ,η ) =

1
Cn

1
Cn

∞

f (x ) =

∫

IR

∫ ∫
0

IR n

f (x )
dx ,
ξ − x −η

(87a)

H n (ξ ,η )
dξdη ,
η − x −ξ

(87b)

n

By analogy with the cylinder Abel transform [41], the corresponding DFW transform is
given by

An (ξ ,κ ) =

1
CA

∫

IR

f (x ) ξ − x
n

2

ξ − x −κ

2

(

)

2

H ξ − x − κ 2 dx ,

(88)

The basis functions of the DFW correspondence of the Hartley transform are

g n (λrk ) =

λn −1 2
4

(2πλrk )−(n 2 )+1 [J (n 2 )−1 (λrk ) + Y(n 2 )−1 (λrk )] ,

n≥2,

(89)

where Y is the Bessel function of the second kind. It is straightforward to construct the
DFW transforms based on (89).
In terms of the Stieltjes Transforms, we have the DFW transform
S (ξ , t ) = ∫ n Γ( p )( ξ − x + t ) dx .
−p

IR

(90)

We need to mention again that the present study is more intuitive conjecture than the
rigorous mathematical inference.

35

In the radial function wavelets (RFW) [40,42,43], the generalized translation operator Tx
(x≥0) for smooth function φ (x) on (0,+∞] is defined by

Txφ ( y ) =

∞
Γ(n 2 )
φ
∫
π Γ(n 2 − 1 2 ) 0

(x

2

)

+ y 2 + 2 xy cosθ (sin θ )

n−2

dθ ,

(91)

where x and y represents two different nodes coordinates. If the above translation
operator is redefined by
T xφ ( y ) = φ ( x − y ) ,

(92)

and then, all the related results of the RFW [40,42,43] can be used with the distance
function wavelets.
Further, in the sense of translation invariant (φ (x) on (-∞,+∞)) [44]
T xφ ( y ) = φ ( x − y ) ,

(93)

we may get more general results. In addition, θ in the translation operator (91) may be
expressed via both the inner product distance variable and the Euclidean distance
variable, namely,

θ = χ ( x − y , x , y , x ⋅ y),
2

2

2

2x ⋅ y = x + y − x − y ,

(94)

(95)

where dot denotes the inner product of two vectors for the ridge distance variable. The
other radial transforms and series [45] may be feasible to be transferred into the DFW
transforms and series.
36

8. Geodesic DFW

The processing of geodesic data is one of the important applications of the radial basis
function [46]. It is, however, noted that there is no particular kernel distance function
available now in literature to deal with the geodesic problems. This section will give a
few geodesic kernel distance functions based on the fundamental solutions of the related
PDEs.
For anisotropic and inhomogeneous objects, [5] lists the fundamental solutions of the ndimensional time-dependent diffusion equation
2
H (t − τ )κ
u (x − ξ , t − τ ) =
e − R 4 (t −τ ) ,
n 2
[4π (t − τ )]

−1 2

*
hn

(96)

and of the n-dimensional Laplace equation
 κ −1 2
 − S ln R,
*
2
u Ln ( x − ξ ) = 
−1 2 2 − n
κ
R
−
,
 (n − 2 )S n

n = 2,
(97)
n ≠ 2,

where the coefficient matrix κ=[κij] represents the parameters in different directions (j)
and locations (i) of anisotropic and inhomogeneous media, and the geodesic distance R is
defined by

R2 =

n

∑κ (x
−1

ij

i , j =1

i

− ξ i )(x j − ξ j ) .

(98)

By analogy with the normal Laplacian high-order solutions (33) and (34), it is straight
forward to have the high-order fundamental solutions of the anisotropic and

37

inhomogeneous Laplacians. It is also not difficult to get the geodesic fundamental
solutions of other PDEs such as the convection-diffusion, Berger plate, Winkler plate,
and wave equations. For instance, we can easily derive the fundamental solution of the
Helmholtz equation

iκ −1 2
ϕ n (x − ξ ) =
4

 λ 


 2πR 

(n 2 )−1

H ((n1)2 )−1 (λR ) , n≥2.

(99)

In addition, [58] presents the general approach to obtain the Laplacian fundamental
solutions in the isotropic heterogeneous media. It is straightforward to construct various
geodesic (inhomogeneous) DFW transforms and series by the PDE kernel solutions of
(96), (97), (99) and those given in [58].

9. Kernel distance sigmoidal functions

In artificial neural network and machine learning, each node has a transfer function which
consists of an activation function (a summation of inputs) and an output function. The
type of activation function depends strongly on its task [47]. Among various possibilities,
a variety of distance variables construct the simplest and efficient activation function.
Note that we here use the widest definition of the distance function [1], which includes
the inner product in the ridgelet function, various distance functions listed in [47], and a
combination of both. The readers have seen from these reports I, II and III that the
translation invariant distance vector, inner product, Euclidean distance function and their
combinations appeared in the foregoing elliptical, hyperbolic, and parabolic distance
function basis functions. [47] argues that the symbolic (nominal) variables can not be in
any liner order. A color variable example is discussed, where each different color is
assigned an identity value and a linear distance thus makes little sense in this type of data.
However, if we use a different digital definition of colors, situations will alter. For
instance, we can assign 0, 0.5, 1 for blue, green, red. Then, any other colors can be
something between 0 and 1. Under this definition of variables, the distance variable does

38

make sense. The basic requirement is that if the distance measure between two values is
small, both should belong to similar classifications [48].
As of the output function, the S-shaped sigmoidal functions may be the most popular
[49,50]. The sigmoidal function yields a value between 0 and 1 in terms of the activation
function value A and a slope (scale) parameter s. The logistic function [49]

σ (s ) =

1
.
1 + e − As

(100)

illustrates the typical form of the sigmoidal functions.
This section will focus on developing of the kernel sigmoidal functions with the help of
the kernel distance functions. It is stressed that although the kernel distance function has
the origin of PDEs, any activation function can be used with them to create a transfer
function. The singularity at the origin of some kernel distance function is not an issue for
the sigmoidal functions.
As an extension of the hyperbolic functions, we first define the multidimensional
hyperbolic functions
sin g ( x ) = ϕ n# (r ) − ϕ n* (r ) ,

(101a)

cos g ( x ) = ϕ n# (r ) + ϕ n* (r ) ,

(101b)

tan g ( x ) =

ϕ n# (r ) − ϕ n* (r )
,
ϕ n# (r ) + ϕ n* (r )

(101c)

csc g ( x ) =

1
,
ϕ (r ) − ϕ n* (r )

(101d)

sec g ( x ) =

1
,
ϕ (r ) + ϕ n* (rk )

(101e)

ϕ n# (r ) + ϕ n* (r )
,
ϕ n# (r ) − ϕ n* (r )

(101f)

#
n

coth g ( x ) =

#
n

39

where ϕn with superscripts # and * are respectively the general solution and the
fundamental solution of the modified Helmholtz equation

ϕ 1# (srk ) =
ϕ n# (srk ) =

(102a)

1
(2πsrk )−(n 2 )+1 I (n 2 )−1 (srk ) ,
2π

ϕ 1* (srk ) =
ϕ n* (srk ) =

1 srk
e ,
2s
n≥2,

1 − srk
e ,
2s

(102b)

(103a)

1
(2πsrk )−(n 2 )+1 K (n 2 )−1 (srk ) ,
2π

n≥2;

(103b)

where rk = x − x k , n is dimensionality, and I denotes the modified Bessel function of
the first kind which grows exponentially as rk→∞. In contrast, the modified Bessel
function of the second kind K decays exponentially and has singularity at the origin.
(101) degenerate into the hyperbolic function in 1D case.
The modified Helmholtz sigmoidal function using the fundamental solution is given by

σ MF (s ) =

1

1 + ϕ n* (sA)

,

(104)

where A is the activation function for a summation of inputs. The natural choice of the
distance variable is the Euclidean one for the modified Helmholtz equation. This,
however, is not necessary. σMF can be simplified as

σˆ MF (s ) =

1+ e

− sA

1
.
(sA)n 2−1 ln( A)

[

]

The modified Helmholtz sigmoidal function using the general solution is given by

40

(105)

σ MG (s ) =

ϕ n# (sA) − 1
.
ϕ n# (sA) + 1

(106)

The modified Helmholtz equation is also called the diffusion equation and exhibits
closely relationship with the stochastic dynamics. Thus, the above modified Helmholtz
sigmoidal functions are also underpinned by the statistics.
As seen before, the DFW basis functions using the general and fundamental solutions of
the convection-diffusion equation naturally combine the inner product for direction
vector and Euclidean distance variables, i.e.

1 −
φ n (v ⋅ ( x − x k ), ρ x − x k ) =
e
2π

v
v ⋅( x − xk )
2D

1 −
φ (v ⋅ ( x − x k ), ρ x − x k ) =
e
2π
#
n

(2πρrk )−(n 2 )+1 K (n 2 )−1 (ρrk ) ,

v
v ⋅( x − xk )
2D

(2πρrk )−(n 2 )+1 I (n 2 )−1 (ρrk ) .

(107)

(108)

It is straightforward to have the convection-diffusion kernel sigmoidal function of the
form

σ CF (w, s ) =

1
,
1 + φ (w ⋅ x, s x − x k )

(109)

σ CG (w, s ) =

φ n# (wT ⋅ x, s x − x k ) − 1
,
φ n# (w T ⋅ x, s x − x k ) + 1

(110)

*
n

T

where w is the direction vector. The kernel solution u h*n (16) of the time-dependent heat
equation can be used to create the time-space sigmoidal function

41

σ H (s,α ) =

1
.
1 + u (s x − ξ ,α (t − τ ))

(111)

*
h

In the same way, the fundamental solutions (96) of heterogeneous transient diffusion
problem using the geodesic distance variable (98) can be employed in making the timespace sigmoidal functions with anisotropic and inhomogeneous parameters. In addition,
the kernel solution of the time-dependent convection-diffusion equation is

φ n (v ⋅ ( x − x k ), ρ x − x k ) = e

−

v
v ⋅( x − xk )
2D

2
H (t − τ )
− x − xk
e
[4πκ (t − τ )]n 2

4κ (t −τ )

.

(112)

In terms of (112), it is straightforward to construct the convection-diffusion temporal
sigmoidal function.
The fundamental solutions of the potential Laplace equations are entitled in creating the
cheap sigmoidal functions

σ (s ) =

1 + (sA)

1

2−n

(n − 2 )S n

,

n≥3,

(113)

which are much more computationally efficient since no costly Bessel functions and
exponents are involved. In addition, the fundamental solutions of wave equations are also
eligible to make the sigmoidal function. We also can use the fundamental solutions of the
Helmholtz equations to construct the periodic sigmoidal functions

σ H (s ) =

1
,
1 + H n* (sA)

(114)

σ H (s ) =

1 − H n# (sA)
,
1 + H n# (sA)

(115)

and

42

where Hn with superscript * and # respectively represent the fundamental and general
solutions of the n-dimensional Helmholtz equation. We can also create variants of the
above distance sigmoidal functions by analogy with the sigmoidal functions (46-51) in
[49]. Among them, the geodesic kernel distance functions (96) and (97) are especially
attractive. In addition, most of DFW’s and their basis functions in this report series such
as the Possion kernel can be directly used as the transfer function in their own right.
It is observed that the S-shaped sigmoidal function behaves very much like what the
shock wave does. Therefore, various distance sigmoidal functions presented above can
also be used in the numerical solution of shock wave and boundary layer problems. It is
very interesting to compare the shock wave, boundary layer and neuron transfer
behaviors. For instance, the Peclet number in the convection-diffusion equation, which is
the ratio of the bulk heat transfer and conductive heat transfer, plays an analogous role as
dose the slope parameter s in the sigmoidal function. The Peclet number is decided by the
heat capacity Cp, thermal conductivity k, density µ, characteristic length D and density ρ
of media [51], i.e.

Pe =

D ⋅Cp ⋅ µ ⋅ ρ
k

.

(116)

The author guesses that the electric signals traveling within neuron may also exhibit
similar mechanism and obey the similar mathematical equations such as the NavierStokes equation or simpler Burgers equation. Thus, the slope parameter s could be
decided via a few physical parameters parallel to those in (116) for the Peclet number. On
the other hand, the Reynolds number is proportional to the ratio of inertial force and
viscous force and is used in momentum, heat, and mass transfer to describe dynamic
similarity [52]. It is conceivable that the solutions of the mathematical physics PDE
models of brain neuron may depend crucially on some parameters like the Peclet or
Reynolds numbers.

43

10. Concluding remarks

A distributed parameter system is represented by a PDE, its counterpart is lumped
parameter systems or simply data set. This suggests any discrete signal, no matter
whether linear or nonlinear, memoryless and memory, random or deterministic, having or
having not PDE model, must have their respective PDE solution structures. So, we should
be aware of the lurking PDE backdrop in processing of specific data such as those in
social sciences, linguistics, and literature. It is well known that all PDEs can be
categorized into elliptical, parabolic, hyperbolic or their mixture types. As such, we must
recognize which kind of PDE type the data we are handling belong to. Majority of data
approximation techniques available now, e.g. spline, however, have their roots on
elliptical Laplace and biharmonic equations and are not suitable and efficient in handling
the hyperbolic and parabolic data.
One of the major drawbacks in the standard wavelets is less incapable of including
inherently the features of the particular type objects. In contrast, the DFW can be easily
designed with the characteristics of certain type problems. For instance, the fundamental
solution of Laplacian under an axisymmetric region has the form [53,54]

w=

4 E (s )
,
( p + q )1 2

(117)

where E is the complete elliptic integral of the first kind, and
p = xi2 + x k2 + ( y i − y k ) ,

(118a)

q = 2 xi x k ,

(118b)

2q
,
p+q

(118c)

2

s=

where x and y are the Cartesian coordinates. With the Laplacian fundamental solution
(117), we can easily make the axisymmetric Laplacian DFW transforms and series as in
sections 2 and 3.
44

In dealing with the infinite domain problems, the distance function presented in [54,55]

u (rk ) =

2C − rk
(rk + C )4

(119)

satisfies the regularity conditions of diffusion problems at infinity, where C is a problemdependent constant. It is obvious that (119) and its variants will be nice DFW basis
functions for diffusion type problems within infinite domains.
This report as well as the other two [1,2] explored many DFW-related issues intuitively
without enough rigorous mathematical justifications. These DFW’s share all of the major
merits of the standard wavelets. Unlike the latter, the DFW in principle, however, is
mathematically very simple and computationally efficient to perform meshfree highdimensional data processing and numerical PDEs. The essential differences between the
DFW and the traditional integral transform and series lie in that the DFW includes the
dimensional effect and keeps the translate and/or rotational invariance.
As complements, there are another two reports [59,60] on the kernel distance function
and its applications. The readers are advised to read them for a complete comprehension
of the distance function. The DFW is expected to find numerous applications in data
processing and numerical PDE. For instance, the composite multiple reciprocity and the
differentiation smoothing are very effective and efficient techniques in evaluating the
domain integral in terms of the BEM and particular solution in terms of the BKM and
BPM, while the Helmholtz-Fourier series is a competitive alternative to the Fourier series
method for the domain integral. Some important PDEs and their solutions were not
mentioned in the three reports. For instance, the kernel solutions of the delay PDEs and
the Black-Scholes equation can also be used to construct the important DFW’s.
A generalized distance function wavelet transform with multiple arguments is stated as

45

v
F (λ , w, n, m,ξ ) = ∫

IR

n

v
f ( x )ℜ n (λ , w, m,ξ − x )d x ,

(120)

v
where arguments λ, w , n, m and ξ represent scale, direction, dimension, order of high-

order PDE, and translate, respectively. The wavelet basis function ℜn could be the kernel
distance function solution of a partial differential equation, which best fits the system
behavior of interest. The shape parameter discussed in section 7 is considered an expedite
substitute of the scale, dimension, or order of PDE and thus is not considered an
independent parameter to be included in (120). Multiple parameters DFW transform
(120) may be useful to analyze very complex systems, where some or all of those
parameters manifest locally.
The ultimate goal of this research is to show some recent advances in building a complete
framework of the distance function. Given their immature nature, it is inevitable to have
more or less errors in these reports. Some cautions on the given results are necessary
based on readers’ own discretion to avoid any potential misleading, and the readers may
better regard the reports as a gathering of intuitive ideas rather than the established
theory.
References

1. W. Chen, Distance function wavelets – Part I: Helmholtz and convection-diffusion
transforms and series, CoRR preprint, http://xxx.lanl.gov/abs//0205019,
Research report of Simula Research Laboratory, May, 2002.
2. W. Chen, Distance function wavelets - Part II: extended results and conjectures,
CoRR preprint, http://xxx.lanl.gov/abs//0205063, Research report of Simula
Research Laboratory, May, 2002.
3. W. Chen, New RBF collocation schemes and kernel RBF with applications, in Int.
Workshop for Meshfree Methods for PDEs, Bonn, Germany, 2001 (also see CoRR
preprint, http://xxx.lanl.gov/abs/cs.NA/0111063).

46

4. W. Chen, Meshfree boundary particle method applied to Helmholtz problems, Engng.
Anal. Bound. Elem., (in press), 2002
5. P.K. Kythe, Fundamental Solutions for Differential Operators and Applications,
Birkhauser, Boston, 1996.
6. M. ZÄHLE: Harmonic calculus on fractals -- A measure--geometric approach II.
(Preprint), July, 2000.
7. M. ZÄHLE: Fractional differentiation in the self-affine case. V - The local degree of
differentiability. Math. Nachr. 185, 279-306, 1997.
8. B. Davies, Integral Transforms and Their Applications, Springer-Verlag, 1978.
9. P.L. Butzer and S. Janshe, Mellin-Fourier series and the classical Mellin transform,
Comput. Math. Appl., 40, 49-62, 2000.
10. W.A. Strauss, Partial Differential Equations, an Introduction, John Wiely & Sons,
1992.
11. Nowak AJ, Neves AC (ed.). The Multiple Reciprocity Boundary Element Method.
Southampton: Comput. Mech. Publ., 1994.
12. S.G. Samko, A.A. Kilbas and O.I. Marichev, Fractional integrals and derivatives:
theory and applications, Gordon and Breach Science Publishers, 1993.
13. W. Chen, Boundary knot method for Laplace and biharmonic problems. The Proc. Of
the 14th Nordic Seminar on Computational Mechanics, Lund, Sweden, Oct. 2001.
14. W. Chen, High-order fundamental and general solutions of convection-diffusion
equation and their applications with boundary particle method, Engng. Anal. Bound.
Elem., (in press), 2002.
15. O. Erosy, Fourier-Related Transforms, Fast Algorithms and Applications, PrenticeHall, 1997.
16. W. Chen and J. He, A study on radial basis function and quasi-Monte Carlo methods,
Int. J. Nonlinear Sci. & Numer. Simulation, 1(4), 343-342, 2000.
17. L.M. Ozaktas, Z. Zalevsky, M.A. Kutay, The Fractional Fourier Transform with
Applications in Optics and Signal Processing, John Wiley & Sons, Ltd, 2001.
18. P. He, Simulation of Ultrasound Pulse Propagation in Lossy Media Obeying a
Frequency Power Law, IEEE Trans. Ultra. Ferro. Freq. Contr., 45(1), 1998
19. Mathematics World Website, http://mathworld.wolfram.com/Self-Similarity.html.

47

20. W. Chen, A note on fractional derivative modeling of broadband frequencydependent absorption: Model III, Research report of Simula Research Laboratory,
April, 2002.
21. T. L. Szabo, Time domain wave equations for lossy media obeying a frequency
power law, J. Acoust. Soc. Amer. 96(1), 491-500, 1994.
22. K. M. Kolwankar and A. D. Gangal, Fractional differentiability of nowhere
differentiable functions and dimensions, CoRR preprint, http://arxiv.org/abs/chaodyn/9609016, Sept. 1996.
23. R. C. Cascaval, E. C. Eckstein, C. L. Frota and J. A. Goldstein, Fractional Telegraph
Equations, http://www.msci.memphis.edu/faculty/jmp_cefg.pdf.
24. N. Makris and M.C. Constantinous, Models of viscoelasticity with complex-order
derivatives, J. Engng. Mech., 119(7), 1993.
25. W. Chen, Definitions of distance function in radial basis function approach, preprint,
http://mpej.unige.ch/mp_arc/a/00-304, 2000.
26. R.L. Hardy, Multiquadratic equations for topography and other irregular surfaces, J.
Geophys. Res., 176 1905-1915 (1971).
27. J. Duchon, Splines minimizing rotation invariant semi-norms in Sobolov spaces, in:
Constructive Theory of Functions of Several Variables, Springer-Verlag, Berlin,
1976.
28. M.D. Buhmann, Radial basis functions. Acta Numerica. 1-38. 2000.
29. R. Franke, Scattered data interpolation: tests of some methods, Math. Comp., 38, 181199, 1982.
30. A. I. Fedoseyev, M. J. Friedman and E. J. Kansa, Continuation for nonlinear elliptic
partial differential equation discretized by the multiquadric method, Int. J. Bifurcat
Chaos, 10, 481-492, 2000.
31. C.A. Micchelli, Interpolation of scattered data: distance matrices and conditionally
positive definite functions, Constr. Approx. 1, 11-22, 1986.
32. M.D. Buhmann, Multiquadric prewavelets on nonequally spaced knots in one
dimension, Math. Comp., 64, 1611-1625, 1995.
33. Y.C. Hon and Z. Wu, A numerical computation for inverse boundary determination
problem, Engng. Anal. Bound. Elem., 24, 599-606, 2000.

48

34. M. Frazier, B. Jawerth and G. Weiss, Littlewood-Paley Theory and the Study of
Function Spaces, CBMS Regional Conference Series in Mathematics, 79, AMS.
35. http://texas.math.ttu.edu/~gilliam/ttu/ode_pde_pdf/Ch8.pdf.
36. B. Fornberg and G. Wright, Stable computation of multiquadric interpolants for all
values of the shape parameter. Submitted to BIT.
37. C. Li and N. Zheng, A theory of constructing locally supported radial wavelet frame
and its application, Science in China (Series E), 42(6), 584-594, 1999.
38. T. A. Driscoll and B. Fornberg, Interpolation in the limit of increasingly flat radial
basis functions, Comp. Math. Appl, 43(3-5), 413-422, 2002.
39. W. Chen, Orthonormal RBF Wavelet and Ridgelet-like Series and Transforms for
High-Dimensional Problem, Int. J. Nonlinear Sci. & Numer. Simulation, 2(2), 155160, 2001.
40. M.A. Mourou and K. Trimeche, Inversion of the Weyl integral transform and the
radon transform on Rn using generalized wavelets, C.R. Math. Rep. Acad. Sci.
Canada, XVIII(2-3), 80-84, 1996.
41. R. Bracewell, The Fourier Transform and Its Applications, 3rd ed., New York:
McGraw-Hill, pp. 262-266, 1999
42. K. Trimeche, Generalized Harmonic Analysis and Wavelet Packets, Gordon &
Breach Sci. Publ., 2001.
43. M. Rosler, Radial wavelets and Bessel-Kingman Hypergroups, preprint.
44. W. Light and H. Wayne, Spaces of distributions, interpolation by translates of a basis
function and error estimates, Numer. Math., 81, 415-450, 1999.
45. M. Ahmed and A.R. Leyman, Tissue characterization using radial transform and
higher order statistics, NORSIG 2000, pp.13-16, Kolmarden, Sweden, 2000.
46. T. A. Foley, The map and blend scattered data interpolant on a sphere, Comp. Math.
Appl. 24(2), 49-60, 1992.
47. D. R. Wilson and T. R. Martinez, Improved heterogeneous distance functions, J.
Artificial Intelligence Research, 6, 1-34, 1997.
48. F. Girosi, Some extensions of radial basis functions and their applications in artificial
intelligence, Comp. Math. Appl. 24(12), 61-80, 1992.

49

49. W. Duch and N. Jankowski, Survey of neural transfer functions, Neural Computing
Surveys, 2, 163-212, 1999.
50. W. Light, Ridge functions, sigmoidal functions and neural networks. In
Approximation Theory VII, E.W. Cheney, C.K. Chui and L.L. Schumaker (eds.),
163-206, 1992.
51. http://www.processassociates.com/process/dimen/dn_pec.htm
52. http://www.processassociates.com/process/dimen/dn_rey.htm
53. C.A. Brebbia, J.C.F. Telles and L.C. Wrobel, Boundary Element Techniques,
Springer-Verlag, Berlin, 1984.
54. P.W. Partridge, C.A. Brebbia, L.W. Wrobel, The Dual Reciprocity Boundary Element
Method, Southampton: Comput. Mech. Publ., 1992.
55. C.F. Loeffler and W.J. Mansur, Dual reciprocity boundary element formulation for
potential problems in infinite domains. In Boundary Elements X, Vol. 2, Comput.
Mech. Publ., Southampton and Springer-Verlag, Berlin, 1988.
56. T. Westphal Jr. and C.S. de Barcellos, On general fundamental solutions of some
linear elliptic differential operators, Engng. Anal. Bound. Elem. 17, 279-285, 1996.
57. L. Chen, A.J. Kassab, D.W. Nicholson and M.B. Chopra, Generalized boundary
element method for solids exhibiting nonhomogeneities, Engng. Anal. Bound. Elem.,
25, 407-422, 2001.
58. E. Divo and A. Kassab, A generalized BEM for steady and transient heat conduction
in media with spatially varying thermal conductivity, Chapter 2, in Advances in
Boundary Elements: Numerical and Mathematical Aspects, M.A. Golberg, (ed.),
Comput. Mech. Publ., Boston, 1998. pp. 37-76.
59. W. Chen, Kernel distance function – Part I: Green second identity and harmonic
function, CoRR preprint, Research report of Simula Research Laboratory, (in
preparation), 2002.
60. W. Chen, Kernel distance function – Part II: Data processing and numerical PDEs
CoRR preprint, Research report of Simula Research Laboratory, (in preparation),
2002.

50

