BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Fredrik Johansson

Fredrik Johansson contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
Symbolic Computationmath.COmath.NTMathematical SoftwareNumerical Analysisastro-phmath-phmath.CVmath.MPmath.NAmath.PR

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
11topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

11 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
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 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