BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Asif Mushtaq appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
We announce some Python classes for numerical solution of partial differential equations, or boundary value problems of ordinary differential equations. These classes are built on routines in \texttt{numpy} and \texttt{scipy.sparse.linalg} (or \texttt{scipy.linalg} for smaller problems).
Some explicit algorithms for higher order symplectic integration of a large class of Hamilton's equations have recently been discussed by Mushtaq \emph{et. al}. Here we present a Python program for automatic numerical implementation of these algorithms for a given Hamiltonian, both for double precision and multiprecision computations. We provide examples of how to use this program, and illustrate behaviour of both the code generator and the generated solver module(s).
We discuss systematic extensions of the standard (St{ö}rmer-Verlet) splitting method for differential equations of Hamiltonian mechanics, with relative accuracy of order $τ^2$ for a timestep of length $τ$, to higher orders in $τ$. We present some splitting schemes, with all intermediate timesteps real and positive, which increase the relative accuracy to order $τ^{N}$ (for N=4, 6, and 8) for a large class of Hamiltonian systems.
We show how the standard (St{ö}rmer-Verlet) splitting method for differential equations of Hamiltonian mechanics (with accuracy of order $τ^2$ for a timestep of length $τ$) can be improved in a systematic manner without using the composition method. We give the explicit expressions which increase the accuracy to order $τ^8$, and demonstrate that the method work on a simple anharmonic oscillator.
We investigate a method to solve a class of Schr{ö}dinger equation eigenvalue problems numerically to very high precision $P$ (from thousands to a million of decimals). The memory requirement, and the number of high precision algebraic operations, of the method scale essentially linearly with $P$ when only eigenvalues are computed. However, since the algorithms for multiplying high precision numbers scale at a rate between $P^{1.6}$ and $P\,\log P\,\log\log P$, the time requirement of our method increases somewhat faster than $P^2$.