6 research outputs found

    The roots to peace in the Democratic Republic of Congo:conservation as a platform for green development

    Get PDF
    Several existing algorithms for computing the matrix cosine employ polynomial or rational approximations combined with scaling and use of a double angle formula. Their derivations are based on forward error bounds. We derive new algorithms for computing the matrix cosine, the matrix sine, and both simultaneously, that are backward stable in exact arithmetic and behave in a forward stable manner in floating point arithmetic. Our new algorithms employ both Pad\'e approximants of sinx\sin x and new rational approximants to cosx\cos x and sinx\sin x obtained from Pad\'e approximants to exe^x. The amount of scaling and the degree of the approximants are chosen to minimize the computational cost subject to backward stability in exact arithmetic. Numerical experiments show that the new algorithms have backward and forward errors that rival or surpass those of existing algorithms and are particularly favorable for triangular matrices

    A Meshfree Splitting Method for Soliton Dynamics in Nonlinear Schrödinger Equations

    No full text
    A new method for the numerical simulation of the so-called solitondynamics arising in a nonlinear Schroedinger equation in the semi-classical regime is proposed. For the time discretization a classical fourth-order splitting method is used. For the spatial discretization, however, a meshfree method is employed in contrast to the usual choice of (pseudo) spectral methods. This approach allows to keep the degrees of freedom almost constant as the semi-classical parameter becomes small. This behavior is confirmed by numerical experiments
    corecore