1,881 research outputs found
Introducing Molly: Distributed Memory Parallelization with LLVM
Programming for distributed memory machines has always been a tedious task,
but necessary because compilers have not been sufficiently able to optimize for
such machines themselves. Molly is an extension to the LLVM compiler toolchain
that is able to distribute and reorganize workload and data if the program is
organized in statically determined loop control-flows. These are represented as
polyhedral integer-point sets that allow program transformations applied on
them. Memory distribution and layout can be declared by the programmer as
needed and the necessary asynchronous MPI communication is generated
automatically. The primary motivation is to run Lattice QCD simulations on IBM
Blue Gene/Q supercomputers, but since the implementation is not yet completed,
this paper shows the capabilities on Conway's Game of Life
Computing the vertices of tropical polyhedra using directed hypergraphs
We establish a characterization of the vertices of a tropical polyhedron
defined as the intersection of finitely many half-spaces. We show that a point
is a vertex if, and only if, a directed hypergraph, constructed from the
subdifferentials of the active constraints at this point, admits a unique
strongly connected component that is maximal with respect to the reachability
relation (all the other strongly connected components have access to it). This
property can be checked in almost linear-time. This allows us to develop a
tropical analogue of the classical double description method, which computes a
minimal internal representation (in terms of vertices) of a polyhedron defined
externally (by half-spaces or hyperplanes). We provide theoretical worst case
complexity bounds and report extensive experimental tests performed using the
library TPLib, showing that this method outperforms the other existing
approaches.Comment: 29 pages (A4), 10 figures, 1 table; v2: Improved algorithm in section
5 (using directed hypergraphs), detailed appendix; v3: major revision of the
article (adding tropical hyperplanes, alternative method by arrangements,
etc); v4: minor revisio
Probabilistic structural analysis by extremum methods
The objective is to demonstrate discrete extremum methods of structural analysis as a tool for structural system reliability evaluation. Specifically, linear and multiobjective linear programming models for analysis of rigid plastic frames under proportional and multiparametric loadings, respectively, are considered. Kinematic and static approaches for analysis form a primal-dual pair in each of these models and have a polyhedral format. Duality relations link extreme points and hyperplanes of these polyhedra and lead naturally to dual methods for system reliability evaluation
Numerical Algebraic Geometry: A New Perspective on String and Gauge Theories
The interplay rich between algebraic geometry and string and gauge theories
has recently been immensely aided by advances in computational algebra.
However, these symbolic (Gr\"{o}bner) methods are severely limited by
algorithmic issues such as exponential space complexity and being highly
sequential. In this paper, we introduce a novel paradigm of numerical algebraic
geometry which in a plethora of situations overcomes these short-comings. Its
so-called 'embarrassing parallelizability' allows us to solve many problems and
extract physical information which elude the symbolic methods. We describe the
method and then use it to solve various problems arising from physics which
could not be otherwise solved.Comment: 36 page
Efficient solution of 3D electromagnetic eddy-current problems within the finite volume framework of OpenFOAM
Eddy-current problems occur in a wide range of industrial and metallurgical
applications where conducting material is processed inductively. Motivated by
realising coupled multi-physics simulations, we present a new method for the
solution of such problems in the finite volume framework of foam-extend, an
extended version of the very popular OpenFOAM software. The numerical procedure
involves a semi-coupled multi-mesh approach to solve Maxwell's equations for
non-magnetic materials by means of the Coulomb gauged magnetic vector potential
and the electric scalar potential. The concept is further extended on the basis
of the impressed and reduced magnetic vector potential and its usage in
accordance with Biot-Savart's law to achieve a very efficient overall modelling
even for complex three-dimensional geometries. Moreover, we present a special
discretisation scheme to account for possible discontinuities in the electrical
conductivity. To complement our numerical method, an extensive validation is
completing the paper, which provides insight into the behaviour and the
potential of our approach.Comment: 47 pages, improved figures, updated references, fixed typos, reverse
phase shift, consistent use of inner produc
- …