1,248 research outputs found
Modula-2*: An extension of Modula-2 for highly parallel programs
Parallel programs should be machine-independent, i.e., independent of properties that are likely to differ from one parallel computer to the next. Extensions are described of Modula-2 for writing highly parallel, portable programs meeting these requirements. The extensions are: synchronous and asynchronous forms of forall statement; and control of the allocation of data to processors. Sample programs written with the extensions demonstrate the clarity of parallel programs when machine-dependent details are omitted. The principles of efficiently implementing the extensions on SIMD, MIMD, and MSIMD machines are discussed. The extensions are small enough to be integrated easily into other imperative languages
Alternatives to standard puncture initial data for binary black hole evolution
Standard puncture initial data have been widely used for numerical binary
black hole evolutions despite their shortcomings, most notably the inherent
lack of gravitational radiation at the initial time that is later followed by a
burst of spurious radiation. We study the evolution of three alternative
initial data schemes. Two of the three alternatives are based on post-Newtonian
expansions that contain realistic gravitational waves. The first scheme is
based on a second-order post-Newtonian expansion in Arnowitt, Deser, and Misner
transverse-traceless (ADMTT) gauge that has been resummed to approach standard
puncture data at the black holes. The second scheme is based on asymptotic
matching of the 4-metrics of two tidally perturbed Schwarzschild solutions to a
first-order post-Newtonian expansion in ADMTT gauge away from the black holes.
The final alternative is obtained through asymptotic matching of the 4-metrics
of two tidally perturbed Schwarzschild solutions to a second-order
post-Newtonian expansion in harmonic gauge away from the black holes. When
evolved, the second scheme fails to produce quasicircular orbits (and instead
leads to a nearly head-on collision). This failure can be traced back to
inaccuracies in the extrinsic curvature due to low order matching. More
encouraging is that the latter two alternatives lead to quasicircular orbits
and show gravitational radiation from the onset of the evolution, as well as a
reduction of spurious radiation. Current deficiencies compared to standard
punctures data include more eccentric trajectories during the inspiral and
larger constraint violations, since the alternative data sets are only
approximate solutions of Einstein's equations. The eccentricity problem can be
ameliorated by adjusting the initial momentum parameters.Comment: 11 pages, 11 figures, 1 appendix, typos corrected, removed duplicate
reference, matches published versio
Quasi-equilibrium binary black hole sequences for puncture data derived from helical Killing vector conditions
We construct a sequence of binary black hole puncture data derived under the
assumptions (i) that the ADM mass of each puncture as measured in the
asymptotically flat space at the puncture stays constant along the sequence,
and (ii) that the orbits along the sequence are quasi-circular in the sense
that several necessary conditions for the existence of a helical Killing vector
are satisfied. These conditions are equality of ADM and Komar mass at infinity
and equality of the ADM and a rescaled Komar mass at each puncture. In this
paper we explicitly give results for the case of an equal mass black hole
binary without spin, but our approach can also be applied in the general case.
We find that up to numerical accuracy the apparent horizon mass also remains
constant along the sequence and that the prediction for the innermost stable
circular orbit is similar to what has been found with the effective potential
method.Comment: 6 pages, 3 figures, 1 tabl
Bosonic behavior of entangled fermions
Two bound, entangled fermions form a composite boson, which can be treated as
an elementary boson as long as the Pauli principle does not affect the behavior
of many such composite bosons. The departure of ideal bosonic behavior is
quantified by the normalization ratio of multi-composite-boson states. We
derive the two-fermion-states that extremize the normalization ratio for a
fixed single-fermion purity P, and establish general tight bounds for this
indicator. For very small purities, P<1/N^2, the upper and lower bounds
converge, which allows to quantify accurately the departure from perfectly
bosonic behavior, for any state of many composite bosons.Comment: 9 pages, 5 figures, accepted by PR
A single-domain spectral method for black hole puncture data
We calculate puncture initial data corresponding to both single and binary
black hole solutions of the constraint equations by means of a pseudo-spectral
method applied in a single spatial domain. Introducing appropriate coordinates,
these methods exhibit rapid convergence of the conformal factor and lead to
highly accurate solutions. As an application we investigate small mass ratios
of binary black holes and compare these with the corresponding test mass limit
that we obtain through a semi-analytical limiting procedure. In particular, we
compare the binding energy of puncture data in this limit with that of a test
particle in the Schwarzschild spacetime and find that it deviates by 50% from
the Schwarzschild result at the innermost stable circular orbit of
Schwarzschild, if the ADM mass at each puncture is used to define the local
black hole masses.Comment: 13 pages, 6 figures; published version with one important change, see
Fig. 4 and the corresponding changes to the tex
Improved initial data for black hole binaries by asymptotic matching of post-Newtonian and perturbed black hole solutions
We construct approximate initial data for non-spinning black hole binary
systems by asymptotically matching the 4-metrics of two tidally perturbed
Schwarzschild solutions in isotropic coordinates to a resummed post-Newtonian
4-metric in ADMTT coordinates. The specific matching procedure used here
closely follows the calculation in gr-qc/0503011, and is performed in the so
called buffer zone where both the post-Newtonian and the perturbed
Schwarzschild approximations hold. The result is that both metrics agree in the
buffer zone, up to the errors in the approximations. However, since isotropic
coordinates are very similar to ADMTT coordinates, matching yields better
results than in the previous calculation, where harmonic coordinates were used
for the post-Newtonian 4-metric. In particular, not only does matching improve
in the buffer zone, but due to the similarity between ADMTT and isotropic
coordinates the two metrics are also close to each other near the black hole
horizons. With the help of a transition function we also obtain a global smooth
4-metric which has errors on the order of the error introduced by the more
accurate of the two approximations we match. This global smoothed out 4-metric
is obtained in ADMTT coordinates which are not horizon penetrating. In
addition, we construct a further coordinate transformation that takes the
4-metric from global ADMTT coordinates to new coordinates which are similar to
Kerr-Schild coordinates near each black hole, but which remain ADMTT further
away from the black holes. These new coordinates are horizon penetrating and
lead, for example, to a lapse which is everywhere positive on the t=0 slice.
Such coordinates may be more useful in numerical simulations.Comment: 25 pages, 21 figures. Replaced with accepted versio
Efficient parallel computation on workstation clusters
We present novel hard- and software that efficiently implements
communication primitives for parallel execution on Workstation
clusters. We provide low communication latencies, minimal protocol,
zero operating system overhead, and high throughput. With this
technology, it is possible to build effective parallel systems
using off-the-shelf workstations. Our goal is to develop a standard
interfaceboard and the necessary software for interfacing any
number of computers, from a workstation to a cabinet full of
workstation-boards
Normality of numbers generated by the values of entire functions
AbstractWe show that the number generated by the q-ary integer part of an entire function of logarithmic order, where the function is evaluated over the natural numbers and the primes, respectively, is normal in base q. This is an extension of related results for polynomials over the real numbers established by Nakai and Shiokawa
Latency hiding in parallel systems: a quantitative approach
In many parallel applications, network latency causes a dramatic
loss in processor utilization. This paper examines software
pipelining as a technique for network latency hiding. It
quantifies the potential improvements with
detailed,instruction-level simulations.
The benchmarks used are the Livermore Loop kernels and BLAS Level
1.
These were parallelized and run on the instruction-level RISC
simulator DLX, extended with both a blocking and a pipelined
network. Our results show that prefetch in a pipelined network
improves performance by a factor of 2 to 9, provided the network
has sufficient bandwidth to accept at least 10 requests per
processor
- …