1,204 research outputs found
Numerical Verification of Affine Systems with up to a Billion Dimensions
Affine systems reachability is the basis of many verification methods. With
further computation, methods exist to reason about richer models with inputs,
nonlinear differential equations, and hybrid dynamics. As such, the scalability
of affine systems verification is a prerequisite to scalable analysis for more
complex systems. In this paper, we improve the scalability of affine systems
verification, in terms of the number of dimensions (variables) in the system.
The reachable states of affine systems can be written in terms of the matrix
exponential, and safety checking can be performed at specific time steps with
linear programming. Unfortunately, for large systems with many state variables,
this direct approach requires an intractable amount of memory while using an
intractable amount of computation time. We overcome these challenges by
combining several methods that leverage common problem structure. Memory is
reduced by exploiting initial states that are not full-dimensional and safety
properties (outputs) over a few linear projections of the state variables.
Computation time is saved by using numerical simulations to compute only
projections of the matrix exponential relevant for the verification problem.
Since large systems often have sparse dynamics, we use Krylov-subspace
simulation approaches based on the Arnoldi or Lanczos iterations. Our method
produces accurate counter-examples when properties are violated and, in the
extreme case with sufficient problem structure, can analyze a system with one
billion real-valued state variables
Multiphase modelling of vascular tumour growth in two spatial dimensions
In this paper we present a continuum mathematical model of vascular tumour growth which is based on a multiphase framework in which the tissue is decomposed into four distinct phases and the principles of conservation of mass and momentum are applied to the normal/healthy cells, tumour cells, blood vessels and extracellular material. The inclusion of a diffusible nutrient, supplied by the blood vessels, allows the vasculature to have a nonlocal influence on the other phases. Two-dimensional computational simulations are carried out on unstructured, triangular meshes to allow a natural treatment of irregular geometries, and the tumour boundary is captured as a diffuse interface on this mesh, thereby obviating the need to explicitly track the (potentially highly irregular and ill-defined) tumour boundary. A hybrid finite volume/finite element algorithm is used to discretise the continuum model: the application of a conservative, upwind, finite volume scheme to the hyperbolic mass balance equations and a finite element scheme with a stable element pair to the generalised Stokes equations derived from momentum balance, leads to a robust algorithm which does not use any form of artificial stabilisation. The use of a matrix-free Newton iteration with a finite element scheme for the nutrient reaction-diffusion equations allows full nonlinearity in the source terms of the mathematical model. Numerical simulations reveal that this four-phase model reproduces the characteristic pattern of tumour growth in which a necrotic core forms behind an expanding rim of well-vascularised proliferating tumour cells. The simulations consistently predict linear tumour growth rates. The dependence of both the speed with which the tumour grows and the irregularity of the invading tumour front on the model parameters are investigated
A bibliography on parallel and vector numerical algorithms
This is a bibliography of numerical methods. It also includes a number of other references on machine architecture, programming language, and other topics of interest to scientific computing. Certain conference proceedings and anthologies which have been published in book form are listed also
Implementing vertex dynamics models of cell populations in biology within a consistent computational framework
The dynamic behaviour of epithelial cell sheets plays a central role during development, growth, disease and wound healing. These processes occur as a result of cell adhesion, migration, division, differentiation and death, and involve multiple processes acting at the cellular and molecular level. Computational models offer a useful means by which to investigate and test hypotheses about these processes, and have played a key role in the study of cell–cell interactions. However, the necessarily complex nature of such models means that it is difficult to make accurate comparison between different models, since it is often impossible to distinguish between differences in behaviour that are due to the underlying model assumptions, and those due to differences in the in silico implementation of the model. In this work, an approach is described for the implementation of vertex dynamics models, a discrete approach that represents each cell by a polygon (or polyhedron) whose vertices may move in response to forces. The implementation is undertaken in a consistent manner within a single open source computational framework, Chaste, which comprises fully tested, industrial-grade software that has been developed using an agile approach. This framework allows one to easily change assumptions regarding force generation and cell rearrangement processes within these models. The versatility and generality of this framework is illustrated using a number of biological examples. In each case we provide full details of all technical aspects of our model implementations, and in some cases provide extensions to make the models more generally applicable
GPU Accelerated Approach to Numerical Linear Algebra and Matrix Analysis with CFD Applications
A GPU accelerated approach to numerical linear algebra and matrix analysis with CFD applications is presented. The works objectives are to (1) develop stable and efficient algorithms utilizing multiple NVIDIA GPUs with CUDA to accelerate common matrix computations, (2) optimize these algorithms through CPU/GPU memory allocation, GPU kernel development, CPU/GPU communication, data transfer and bandwidth control to (3) develop parallel CFD applications for Navier Stokes and Lattice Boltzmann analysis methods. Special consideration will be given to performing the linear algebra algorithms under certain matrix types (banded, dense, diagonal, sparse, symmetric and triangular). Benchmarks are performed for all analyses with baseline CPU times being determined to find speed-up factors and measure computational capability of the GPU accelerated algorithms. The GPU implemented algorithms used in this work along with the optimization techniques performed are measured against preexisting work and test matrices available in the NIST Matrix Market. CFD analysis looked to strengthen the assessment of this work by providing a direct engineering application to analysis that would benefit from matrix optimization techniques and accelerated algorithms. Overall, this work desired to develop optimization for selected linear algebra and matrix computations performed with modern GPU architectures and CUDA developer which were applied directly to mathematical and engineering applications through CFD analysis
Past, Present, and Future of Simultaneous Localization And Mapping: Towards the Robust-Perception Age
Simultaneous Localization and Mapping (SLAM)consists in the concurrent
construction of a model of the environment (the map), and the estimation of the
state of the robot moving within it. The SLAM community has made astonishing
progress over the last 30 years, enabling large-scale real-world applications,
and witnessing a steady transition of this technology to industry. We survey
the current state of SLAM. We start by presenting what is now the de-facto
standard formulation for SLAM. We then review related work, covering a broad
set of topics including robustness and scalability in long-term mapping, metric
and semantic representations for mapping, theoretical performance guarantees,
active SLAM and exploration, and other new frontiers. This paper simultaneously
serves as a position paper and tutorial to those who are users of SLAM. By
looking at the published research with a critical eye, we delineate open
challenges and new research issues, that still deserve careful scientific
investigation. The paper also contains the authors' take on two questions that
often animate discussions during robotics conferences: Do robots need SLAM? and
Is SLAM solved
Parallelisation of genetic algorithms for the 2-page crossing number problem
Genetic algorithms have been applied to solve the 2-page crossing number problem
successfully, but since they work with one global population, the search time and
space are limited. Parallelisation provides an attractive prospect to improve the efficiency
and solution quality of genetic algorithms. This paper investigates the complexity
of parallel genetic algorithms (PGAs) based on two evaluation measures:
Computation-time to Communication-time and Population-size to Chromosomesize.
Moreover, the paper unifies the framework of PGA models with the function
PGA (subpopulation size; cluster size, migration period; topology), and explores
the performance of PGAs for the 2-page crossing number problem
- …