BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Fredrik Johansson

Fredrik Johansson appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
Symbolic ComputationMathematical SoftwareNumerical Analysismath.COmath.NAmath.NTastro-phmath-phmath.CVmath.MPmath.PRNeurons and Cognition

Catalog footprint

What is connected

15works
12topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2022arXiv

Computing elementary functions using multi-prime argument reduction

We describe an algorithm for arbitrary-precision computation of the elementary functions (exp, log, sin, atan, etc.) which, after a cheap precomputation, gives roughly a factor-two speedup over previous state-of-the-art algorithms at precision from a few thousand bits up to millions of bits. Following an idea of Sch{ö}nhage, we perform argument reduction using Diophantine combinations of logarithms of primes; our contribution is to use a large set of primes instead of a single pair, aided by a fast algorithm to solve the associated integer relation problem. We also list new, optimized Machin-like formulas for the necessary logarithm and arctangent precomputations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Passage-of-time Model of the Cerebellar Purkinje Cell

The cerebellar Purkinje cell controlling eyeblinks can learn, remember and reproduce the interstimulus interval in a classical conditioning paradigm. Given temporally separated inputs, the cerebellar Purkinje cell learns to pause its tonic inhibition of a motor pathway with high temporal precision so that an overt blink occurs at the right time. Most models relegate the Purkinje cell's passage-of-time representation to afferent granule cells, a subpopulation of which is supposedly selected for synaptic depression in order to make the Purkinje cell pause. However, granule cell models have recently faced two crucial challenges: 1) bypassing granule cells and directly stimulating the Purkinje cell's pre-synaptic fibers during training still produces a well-timed pause, and 2) the Purkinje cell can reproduce the learned pause, invariant to the temporal structure of probe stimulation. Here, we present a passage-of-time model which is internal to the Purkinje cell and is invariant to probe structure. The model accurately simulates the Purkinje cell learning mechanism and makes testable electrophysiological predictions. Importantly, the model is a numerical proof-of-principle for a biological learning mechanism which does not rely on changes of synaptic strength.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Arb: Efficient Arbitrary-Precision Midpoint-Radius Interval Arithmetic

Arb is a C library for arbitrary-precision interval arithmetic using the midpoint-radius representation, also known as ball arithmetic. It supports real and complex numbers, polynomials, power series, matrices, and evaluation of many special functions. The core number types are designed for versatility and speed in a range of scenarios, allowing performance that is competitive with non-interval arbitrary-precision types such as MPFR and MPC floating-point numbers. We discuss the low-level number representation, strategies for precision and error bounds, and the implementation of efficient polynomial arithmetic with interval coefficients.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Computing hypergeometric functions rigorously

We present an efficient implementation of hypergeometric functions in arbitrary-precision interval arithmetic. The functions ${}_0F_1$, ${}_1F_1$, ${}_2F_1$ and ${}_2F_0$ (or the Kummer $U$-function) are supported for unrestricted complex parameters and argument, and by extension, we cover exponential and trigonometric integrals, error functions, Fresnel integrals, incomplete gamma and beta functions, Bessel functions, Airy functions, Legendre functions, Jacobi polynomials, complete elliptic integrals, and other special functions. The output can be used directly for interval computations or to generate provably correct floating-point approximations in any format. Performance is competitive with earlier arbitrary-precision software, and sometimes orders of magnitude faster. We also partially cover the generalized hypergeometric function ${}_pF_q$ and computation of high-order parameter derivatives.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Efficient implementation of elementary functions in the medium-precision range

We describe a new implementation of the elementary transcendental functions exp, sin, cos, log and atan for variable precision up to approximately 4096 bits. Compared to the MPFR library, we achieve a maximum speedup ranging from a factor 3 for cos to 30 for atan. Our implementation uses table-based argument reduction together with rectangular splitting to evaluate Taylor series. We collect denominators to reduce the number of divisions in the Taylor series, and avoid overhead by doing all multiprecision arithmetic using the mpn layer of the GMP library. Our implementation provides rigorous error bounds.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A bound for the error term in the Brent-McMillan algorithm

The Brent-McMillan algorithm B3 (1980), when implemented with binary splitting, is the fastest known algorithm for high-precision computation of Euler's constant. However, no rigorous error bound for the algorithm has ever been published. We provide such a bound and justify the empirical observations of Brent and McMillan. We also give bounds on the error in the asymptotic expansions of functions related to modified Bessel functions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A fast algorithm for reversion of power series

We give an algorithm for reversion of formal power series, based on an efficient way to implement the Lagrange inversion formula. Our algorithm requires $O(n^{1/2}(M(n) + MM(n^{1/2})))$ operations where $M(n)$ and $MM(n)$ are the costs of polynomial and matrix multiplication respectively. This matches the asymptotic complexity of an algorithm of Brent and Kung, but we achieve a constant factor speedup whose magnitude depends on the polynomial and matrix multiplication algorithms used. Benchmarks confirm that the algorithm performs well in practice.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Evaluating parametric holonomic sequences using rectangular splitting

We adapt the rectangular splitting technique of Paterson and Stockmeyer to the problem of evaluating terms in holonomic sequences that depend on a parameter. This approach allows computing the $n$-th term in a recurrent sequence of suitable type using $O(n^{1/2})$ "expensive" operations at the cost of an increased number of "cheap" operations. Rectangular splitting has little overhead and can perform better than either naive evaluation or asymptotically faster algorithms for ranges of $n$ encountered in applications. As an example, fast numerical evaluation of the gamma function is investigated. Our work generalizes two previous algorithms of Smith.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Finding Hyperexponential Solutions of Linear ODEs by Numerical Evaluation

We present a new algorithm for computing hyperexponential solutions of ordinary linear differential equations with polynomial coefficients. The algorithm relies on interpreting formal series solutions at the singular points as analytic functions and evaluating them numerically at some common ordinary point. The numerical data is used to determine a small number of combinations of the formal series that may give rise to hyperexponential solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Rigorous high-precision computation of the Hurwitz zeta function and its derivatives

We study the use of the Euler-Maclaurin formula to numerically evaluate the Hurwitz zeta function $ζ(s,a)$ for $s, a \in \mathbb{C}$, along with an arbitrary number of derivatives with respect to $s$, to arbitrary precision with rigorous error bounds. Techniques that lead to a fast implementation are discussed. We present new record computations of Stieltjes constants, Keiper-Li coefficients and the first nontrivial zero of the Riemann zeta function, obtained using an open source implementation of the algorithms described in this paper.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Using functional equations to enumerate 1324-avoiding permutations

We consider the problem of enumerating permutations with exactly r occurrences of the pattern 1324 and derive functional equations for this general case as well as for the pattern avoidance (r=0) case. The functional equations lead to a new algorithm for enumerating length n permutations that avoid 1324. This approach is used to enumerate the 1324-avoiders up to n=31. We also extend those functional equations to account for the number of inversions and derive analogous algorithms.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

No observational constraints from hypothetical collisions of hypothetical dark halo primordial black holes with galactic objects

It was suggested by several authors that hypothetical primordial black holes (PBHs) may contribute to the dark matter in our Galaxy. There are strong constraints based on the Hawking evaporation that practically exclude PBHs with masses m~1e15-1e16g and smaller as significant contributors to the Galactic dark matter. Similarly, PBHs with masses greater than about 1e26g are practically excluded by the gravitational lensing observation. The mass range between 10e16g<m<10e26g is unconstrained. In this paper, we examine possible observational signatures in the unexplored mass range, investigating hypothetical collisions of PBHs with main sequence stars, red giants, white dwarfs, and neutron stars in our Galaxy. This has previously been discussed as possibly leading to an observable photon eruption due to shock production during the encounter. We find that such collisions are either too rare to be observed (if the PBH masses are typically larger than about 1e20g), or produce too little power to be detected (if the masses are smaller than about 1e20g).

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Scaling limits of anisotropic Hastings-Levitov clusters

We consider a variation of the standard Hastings-Levitov model HL(0), in which growth is anisotropic. Two natural scaling limits are established and we give precise descriptions of the effects of the anisotropy. We show that the limit shapes can be realised as Loewner hulls and that the evolution of harmonic measure on the cluster boundary can be described by the solution to a deterministic ordinary differential equation related to the Loewner equation. We also characterise the stochastic fluctuations around the deterministic limit flow.

0 citations0 reviewsNo code linkedOpen access