4 research outputs found
Detection of permutation symmetries in numerical constraint satisfaction problems
Technical reportThis technical report summarizes the work done by Ieva Dauzickaite during her traineeship at IRI from 2016-10-03 to 2017-03-31. The main objectives of the stay at IRI were to study detection of variable symmetries in numerical constraint satisfaction problems and to implement a method that detects such symmetries.Peer ReviewedPostprint (published version
Recommended from our members
On the preconditioning for weak constraint four-dimensional variational data assimilation
Data assimilation is used to obtain an improved estimate (analysis) of the state of a dynamical system by combining a previous estimate with observations of the system. A
weak constraint four-dimensional variational assimilation (4D-Var) method accounts for
the dynamical model error and is of large interest in numerical weather prediction. The
analysis can be approximated by solving a series of large sparse symmetric positive definite (SPD) or saddle point linear systems of equations. The iterative solvers used for
these systems require preconditioning for a satisfactory performance. In this thesis, we
use randomised numerical methods to construct effective preconditioners that are cheap
to construct and apply. We employ a randomised eigenvalue decomposition to construct
limited memory preconditioners (LMPs) for a forcing formulation of 4D-Var independently of the previously solved systems. This preconditioning remains effective even if the
subsequent systems change significantly. We propose a randomised approximation of a
control variable transform technique (CVT) to precondition the SPD system of the state
formulation, which preserves potential for a time-parallel model integration. A new way
to include the observation information in the approximation of the inverse Schur complement in the block diagonal preconditioner for the saddle point formulations is introduced,
namely applying the randomised LMPs. Numerical experiments with idealised systems
show that the proposed preconditioners improve the performance of the iterative solvers.
We provide theoretical results describing the change of the extreme eigenvalues of the
unpreconditioned and preconditioned coefficient matrices when new observations of the
dynamical system are added. These show that small positive eigenvalues can cause convergence issues. New eigenvalue bounds for the SPD and saddle point coefficient matrices
in the state formulation emphasize their sensitivities to the observations. These results
can guide the design of other preconditioners
On time-parallel preconditioning for the state formulation of incremental weak constraint 4D-Var
Using a high degree of parallelism is essential to perform data assimilation efficiently. The state formulation of the incremental weak constraint four-dimensional variational data assimilation method allows parallel calculations in the time dimension. In this approach, the solution is approximated by minimising a series of quadratic cost functions using the conjugate gradient method. To use this method in practice, effective preconditioning strategies that maintain the potential for parallel calculations are needed. We examine approximations to the control variable transform (CVT) technique when the latter is beneficial. The new strategy employs a randomised singular value decomposition and retains the potential for parallelism in the time domain. Numerical results for the Lorenz 96 model show that this approach accelerates the minimisation in the first few iterations, with better results when CVT performs well
Detection of permutation symmetries in numerical constraint satisfaction problems
Technical reportThis technical report summarizes the work done by Ieva Dauzickaite during her traineeship at IRI from 2016-10-03 to 2017-03-31. The main objectives of the stay at IRI were to study detection of variable symmetries in numerical constraint satisfaction problems and to implement a method that detects such symmetries.Peer Reviewe