65,315 research outputs found
Status of the differential transformation method
Further to a recent controversy on whether the differential transformation
method (DTM) for solving a differential equation is purely and solely the
traditional Taylor series method, it is emphasized that the DTM is currently
used, often only, as a technique for (analytically) calculating the power
series of the solution (in terms of the initial value parameters). Sometimes, a
piecewise analytic continuation process is implemented either in a numerical
routine (e.g., within a shooting method) or in a semi-analytical procedure
(e.g., to solve a boundary value problem). Emphasized also is the fact that, at
the time of its invention, the currently-used basic ingredients of the DTM
(that transform a differential equation into a difference equation of same
order that is iteratively solvable) were already known for a long time by the
"traditional"-Taylor-method users (notably in the elaboration of software
packages --numerical routines-- for automatically solving ordinary differential
equations). At now, the defenders of the DTM still ignore the, though much
better developed, studies of the "traditional"-Taylor-method users who, in
turn, seem to ignore similarly the existence of the DTM. The DTM has been given
an apparent strong formalization (set on the same footing as the Fourier,
Laplace or Mellin transformations). Though often used trivially, it is easily
attainable and easily adaptable to different kinds of differentiation
procedures. That has made it very attractive. Hence applications to various
problems of the Taylor method, and more generally of the power series method
(including noninteger powers) has been sketched. It seems that its potential
has not been exploited as it could be. After a discussion on the reasons of the
"misunderstandings" which have caused the controversy, the preceding topics are
concretely illustrated.Comment: To appear in Applied Mathematics and Computation, 29 pages,
references and further considerations adde
The automatic solution of partial differential equations using a global spectral method
A spectral method for solving linear partial differential equations (PDEs)
with variable coefficients and general boundary conditions defined on
rectangular domains is described, based on separable representations of partial
differential operators and the one-dimensional ultraspherical spectral method.
If a partial differential operator is of splitting rank , such as the
operator associated with Poisson or Helmholtz, the corresponding PDE is solved
via a generalized Sylvester matrix equation, and a bivariate polynomial
approximation of the solution of degree is computed in
operations. Partial differential operators of
splitting rank are solved via a linear system involving a block-banded
matrix in operations. Numerical
examples demonstrate the applicability of our 2D spectral method to a broad
class of PDEs, which includes elliptic and dispersive time-evolution equations.
The resulting PDE solver is written in MATLAB and is publicly available as part
of CHEBFUN. It can resolve solutions requiring over a million degrees of
freedom in under seconds. An experimental implementation in the Julia
language can currently perform the same solve in seconds.Comment: 22 page
- …