Past

For large--scale saddle point problems, the application of exact iterative schemes and preconditioners may be computationally expensive. In practical situations, only approximations to the inverses of the diagonal block or the related cross-product matrices are considered, giving rise to inexact versions of various solvers. Therefore, the approximation effects must be carefully studied. In this talk we study numerical behavior of several iterative Krylov subspace solvers applied to the solution of large-scale saddle point problems. Two main representatives of the segregated solution approach are analyzed: the Schur complement reduction method, based on an (iterative) elimination of primary variables and the null-space projection method which relies on a basis for the null-space for the constraints. We concentrate on the question what is the best accuracy we can get from inexact schemes solving either Schur complement system or the null-space projected system when implemented in finite precision arithmetic. The fact that the inner solution tolerance strongly influences the accuracy of computed iterates is known and was studied in several contexts.

In particular, for several mathematically equivalent implementations we study the influence of inexact solving the inner systems and estimate their maximum attainable accuracy. When considering the outer iteration process our rounding error analysis leads to results similar to ones which can be obtained assuming exact arithmetic. The situation is different when we look at the residuals in the original saddle point system. We can show that some implementations lead ultimately to residuals on the the roundoff unit level independently of the fact that the inner systems were solved inexactly on a much higher level than their level of limiting accuracy. Indeed, our results confirm that the generic and actually the cheapest implementations deliver the approximate solutions which satisfy either the second or the first block equation to the working accuracy. In addition, the schemes with a corrected direct substitution are also very attractive. We give a theoretical explanation for the behavior which was probably observed or it is already tacitly known. The implementations that we pointed out as optimal are actually those which are widely used and suggested in applications.

I will introduce the moving frame approach to the analysis of invariant variational problems and the evolution of differential invariants under invariant submanifold flows. Applications will include differential geometric flows, integrable systems, and image processing.

Regions of a link diagram can be colored in black and white in a checkerboard manner. Putting a vertex in each black region and connecting two vertices by an edge if the corresponding regions share a crossing yields a planar graph. In 1987 Thistlethwaite proved that the Jones polynomial of the link can be obtained by a specialization of the Tutte polynomial of this planar graph. The goal of my talk will be an explanation of a generalization of Thistlethwaite's theorem to virtual links. In this case graphs will be embedded into a (higher genus, possibly non-oriented) surface. For such graphs we used a generalization of the Tutte polynomial called the Bollobas-Riordan polynomial. For graphs on
surfaces the natural duality can be generalized to a duality with respect to a subset of edges. The generalized dual graph might be embedded into a different surface. I will explain a relation between the Bollobas-Riordan polynomials of dual graphs. This relation unifies various Thistlethwaite type theorems.

The theory of risk measurement has been extensively developed over the past ten years or so, but there has been comparatively little effort devoted to using this theory to inform portfolio choice. One theme of this paper is to study how an investor in a conventional log-Brownian market would invest to optimize expected utility of terminal wealth, when subjected to a bound on his risk, as measured by a coherent law-invariant risk measure. Results of Kusuoka lead to remarkably complete expressions for the solution to this problem.

The second theme of the paper is to discuss how one would actually manage (not just measure) risk. We study a principal/agent problem, where the principal is required to satisfy some risk constraint. The principal proposes a compensation package to the agent, who then optimises selfishly ignoring the risk constraint. The principal can pick a compensation package that induces the agent to select the principal's optimal choice.

We report on work joint with Scott Parsell in which estimates are obtained for the set of real numbers not closely approximated by a given form with real coefficients. "Slim"

technology plays a role in obtaining the sharpest estimates.

This is joint work with Diaz, Glesser and Park.

In Proc. Instructional Conf, Oxford 1969, G. Glauberman shows that

several global properties of a finite group are determined by the properties

of its p-local subgroups for some prime p. With Diaz, Glesser and Park, we

reviewed these results by replacing the group by a saturated fusion system

and proved that the ad hoc statements hold. In this talk, we will present

the adapted versions of some of Glauberman and Thompson theorems.

Multicore chips are nearly ubiquitous in modern machines, and to fully exploit this continuation of Moore's Law, numerical algorithms need to be able to exploit parallelism. We describe recent approaches to both dense and sparse parallel Cholesky factorization on shared memory multicore systems and present results from our new codes for problems arising from large real-world applications. In particular we describe our experiences using directed acyclic graph based scheduling in the dense case and retrofitting parallelism to a

sparse serial solver.

We introduce the notion of $G$-Higgs bundle from studying the representations of the fundamental group of a closed connected oriented surface $X$ in a Lie group $G$. If $G$ turns to be the isometry group of a Hermitian symmetric space, much more can be said about the moduli space of $G$-Higgs bundles, but this also implies dealing with exceptional cases. We will try to face all these subjects intuitively and historically, when possible!

In this talk, I shall endeavour to explain to the uneducated and uninitiated the joys and pleasures one can have studying automata.