In this paper, we report the application of two parameter coupled alternating group explicit (CAGE) iteration and Newton-CAGE iteration methods for the cubic spline solution of non-linear differential equation u" = f(r,u,u') subject to given natural boundary conditions. The error analysis for CAGE iteration method is discussed in details. We compared the results of proposed CAGE iteration metho...

A mathematical model of Calcium Induced Calcium Release in a heart cell has been developed that consists of three coupled non-linear advection-diffusion-reaction equations. A program in C with MPI based on matrixfree Newton-Krylov method gives very good scalability, but still requires large run times for fine meshes. A programming model with CUDA and MPI that utilizes GPUs on a hybrid node can ...

This article develops fast numerical methods for the practical solution of the famous EIT and DC-resistivity problems in the presence of discontinuities and potentially many experiments or data. Based on a Gauss-Newton (GN) approach coupled with preconditioned conjugate gradient (PCG) iterations, we propose two algorithms. One determines adaptively the number of inner PCG iterations required to...

We present a method for the quantum mechanical inverse scattering problem at xed energy for coupled channels in reactions with particles having internal degrees of freedom. The scattered particles can be excited by a local interaction between the relative motion and the internal dynamics which can be expanded in multipoles. The inverse scattering problem is solved by an extension of the modiied...

In this paper, the numerical algorithm for solving the state and output feedback H∞-constrained LQG control problem is investigated. Although the equations that have to be solved to design the controller consist of the nonlinear crosscoupled algebraic Riccati equations (CAREs), it is newly proven that both the uniqueness and the positive semidefiniteness of the iterative solutions can be guaran...

We present a Monte Carlo synthetic-acceleration method for solving the time-dependent, nonlinear, equilibrium diffusion equation in three-dimensions. The new scheme uses the adjoint Monte Carlo method as a relaxation step that accelerates standard Jacobi iteration. Results show that this method is 40% faster than regular Jacobi iteration and is competitive with Jacobi-preconditioned Conjugate G...

In this paper we introduce a new framework for image registration. Our formulation is based on consistent discretization of the optimization problem coupled with a multigrid solution of the linear system which evolve in a Gauss-Newton iteration. We show that our discretization is h-elliptic independent of parameter choice and therefore a simple multigrid implementation can be used. To overcome ...