Computing Solution Operators of Boundary-value Problems for Some Linear Hyperbolic Systems of PDEs

We discuss possibilities of application of Numerical Analysis methods to proving computability, in the sense of the TTE approach, of solution operators of boundary-value problems for systems of PDEs. We prove computability of the solution operator for a symmetric hyperbolic system with computable real coefficients and dissipative boundary conditions, and of the Cauchy problem for the same system (we also prove computable dependence on the coefficients) in a cube $Q\subseteq\mathbb R^m$. Such systems describe a wide variety of physical processes (e.g. elasticity, acoustics, Maxwell equations). Moreover, many boundary-value problems for the wave equation also can be reduced to this case, thus we partially answer a question raised in Weihrauch and Zhong (2002). Compared with most of other existing methods of proving computability for PDEs, this method does not require existence of explicit solution formulas and is thus applicable to a broader class of (systems of) equations.Comment: 31 page

Computability and analysis: the legacy of Alan Turing

We discuss the legacy of Alan Turing and his impact on computability and analysis.Comment: 49 page

Computability of differential equations

In this chapter, we provide a survey of results concerning the computability and computational complexity of differential equations. In particular, we study the conditions which ensure computability of the solution to an initial value problem for an ordinary differential equation (ODE) and analyze the computational complexity of a computable solution. We also present computability results concerning the asymptotic behaviors of ODEs as well as several classically important partial differential equations.info:eu-repo/semantics/acceptedVersio

An Algorithm for Computing Fundamental Solutions of Difference Operators

We propose an FFT-based algorithm for computing fundamental solutions of difference operators with constant coefficients. Our main contribution is to handle cases where the symbol has zeros

Numeriska beräkningar med fundamentallösningar

Two solution strategies for large, sparse, and structured algebraic systems of equations are considered. The first strategy is to construct efficient preconditioners for iterative solvers. The second is to reduce the sparse algebraic system to a smaller, dense system of equations, which are called the boundary summation equations. The proposed preconditioners perform well when applied to equations that are discretizations of linear first order partial differential equations. Analysis shows that also very simple iterative methods converge in a number of iterations that is independent of the number of unknowns, if our preconditioners are applied to certain scalar model problems. Numerical experiments indicate that this property holds also for more complicated cases, and a flow problem modeled by the nonlinear Euler equations is treated successfully. The reduction process is applicable to a large class of difference equations. There is no approximation involved in the reduction, so the solution of the original algebraic equations is determined exactly if the reduced system is solved exactly. The reduced system is well suited for iterative solution, especially if the original system of equations is a discretization of a first order differential equation. The technique is used for several problems, ranging from scalar model problems to a semi-implicit discretization of the compressible Navier-Stokes equations. Both strategies use the concept of fundamental solutions, either of differential or difference operators. An algorithm for computing fundamental solutions of difference operators is also presented

Numeriska beräkningar med fundamentallösningar

Two solution strategies for large, sparse, and structured algebraic systems of equations are considered. The first strategy is to construct efficient preconditioners for iterative solvers. The second is to reduce the sparse algebraic system to a smaller, dense system of equations, which are called the boundary summation equations. The proposed preconditioners perform well when applied to equations that are discretizations of linear first order partial differential equations. Analysis shows that also very simple iterative methods converge in a number of iterations that is independent of the number of unknowns, if our preconditioners are applied to certain scalar model problems. Numerical experiments indicate that this property holds also for more complicated cases, and a flow problem modeled by the nonlinear Euler equations is treated successfully. The reduction process is applicable to a large class of difference equations. There is no approximation involved in the reduction, so the solution of the original algebraic equations is determined exactly if the reduced system is solved exactly. The reduced system is well suited for iterative solution, especially if the original system of equations is a discretization of a first order differential equation. The technique is used for several problems, ranging from scalar model problems to a semi-implicit discretization of the compressible Navier-Stokes equations. Both strategies use the concept of fundamental solutions, either of differential or difference operators. An algorithm for computing fundamental solutions of difference operators is also presented