For more information see the official workshop home page.

Registration deadline: **August 1, 2017**

Time: 14:00 - 15:30

Address: Sokolovská 83, Praha

Room: K3

Speaker: Esmond G. Ng (Berkeley Lab)

Abstract: pdf file

A lecture organized in cooperation with the Mathematical Institute of the Czech Academy of Sciences.

]]>the Institute There is no age limit for applicants; however eligible

candidates in Mathematics must have earned their PhD not earlier than

January 1, 2012.

Candidates who are preparing their doctoral thesis are eligible to apply;

however, they must have obtained their Ph.D degree before taking up their

appointment with GSSI.

______________________________

Research Grants The annual gross salary is € 36.000,00. Each post doctoral

research grant is intended for a duration of two years.

______________________________

Topics Stochastic Differential Equations. Particles systems and

macroscopic limits, statistical mechanics and phase transitions.

Mathematical theory in classical and nonclassical fluids, Nonlinear

dispersive PDEs. Dynamics and Control on Networks. Mathematical models of

collective behaviour, in biology and social sciences. Mathematical Models in

continuum mechanics of solids, fluids and biological matter. Numerical

methods for nonlinear PDEs.

http://www.gssi.it/postdoc/

]]>**Abstract:**

In this talk we study the numerical behavior of the generalized minimal residual (GMRES) method for solving singular linear systems. GMRES determines a solution without breakdown in theory in the two cases: the coefficient matrix is symmetric in its range space (EP); its range space and null spaces are disjoint (GP). We show how the inconsistency of a linear system and the principal angles between the range o A and the range o A^T affect the conditioning of the extended Hessenberg matrix in the Arnoldi decomposition and the accuracy of computed iterates. We compare GMRES with the range restricted GMRES (RR-GMRES) method and the simpler GMRES method. Numerical experiments show typical behavior of GMRES in the EP and GP cases.

]]>The core topics include a posteriori error estimation and adaptive methods for nonlinear partial differential equations and evolutionary problems, for a broad range of numerical methods including high-order methods.

For the full details on the workshop, we invite you to look at our dedicated webpage https://project.inria.fr/gatipor/events/workshop/ where you may find a list of speakers and some further information.

Registration is free but compulsory, and it gives access to the talks and coffee breaks. You may register directly by emailing your name and institution to the organisers at organisers_aposterioriworkshop2017@inria.fr

]]>**Abstract:**

The efficient solution of large linear least squares problems in which the system matrix A contains rows with very different densities is challenging. There have been many classical contributions to solving this problem that focus on direct methods; they can be found in the monograph by Ake Bjorck. Such solvers typically perform a splitting of the rows of A into two row blocks, As and Ad. The block As is such that the sparse factorization of the normal matrix A^T_s As is feasible, while the rows in the block Ad have a relatively large number of nonzero entries. These dense rows are initially ignored, a factorization of the sparse part is computed using a sparse direct solver and then the solution updated to take account of the omitted dense rows.

There are two potential weaknesses of this approach. First, in practical applications the number of rows that contain a significant number of entries may not be small. Processing some of the denser rows separately may improve performance. Furthermore, large-scale problems require the use of preconditioned iterative solvers. A straightforward proposal to precondition the iterative solver using only an incomplete factorization of the sparse block while discarding the dense block may not lead to any success. In this presentation, we propose processing As separately within a conjugate gradient method using an incomplete factorization preconditioner combined with the factorization of a dense matrix of size equal to the number of rows in Ad. Problems arising from practical applications are used to demonstrate the potential of the new approach.

**Abstract:**

We consider a convection-diffusion boundary value problem (1D, convection dominated) with Dirichlet boundary conditions. Because of the occurrence of boundary layers in the solution, such problems are difficult to solve numerically. Standard discretization techniques typically cannot resolve the layers and have to be stabilized in order to yield an acceptable numerical solution. Here we consider discretizations using a Shishkin mesh, which clusters mesh points in the layer instead of putting them equidistantly over the whole region, leading to linear algebraic systems with highly nonnormal matrices. In this talk we are interested in solving such systems using the GMRES method and the multiplicative Schwarz method.

]]>