15 research outputs found
Slices, slabs, and sections of the unit hypercube
Using combinatorial methods, we derive several formulas for the volume of
convex bodies obtained by intersecting a unit hypercube with a halfspace, or
with a hyperplane of codimension 1, or with a flat defined by two parallel
hyperplanes. We also describe some of the history of these problems, dating to
Polya's Ph.D. thesis, and we discuss several applications of these formulas.Comment: 11 pages; minor corrections to reference
Highly parallel, high-precision numerical integration
This paper describes a scheme for rapidly computing numerical values of definite integrals to very high accuracy, ranging from ordinary machine precision to hundreds or thousands of digits, even for functions with singularities or infinite derivatives at endpoints. Such a scheme is of interest not only in computational physics and computational chemistry, but also in experimental mathematics, where high-precision numerical values of definite integrals can be used to numerically discover new identities. This paper discusses techniques for a parallel implementation of this scheme, then presents performance results for 1-D and 2-D test suites. Results are also given for a certain problem from mathematical physics, which features a difficult singularity, confirming a conjecture to 20,000 digit accuracy. The performance rate for this latter calculation on 1024 CPUs is 690 Gflop/s. We believe that this and one other 20,000-digit integral evaluation that we report are the highest-precision non-trivial numerical integrations performed to date