Search · four archives
Search · four archives
14 papers · ranked by Valyu relevance
Daichi Mukunoki, Katsuhisa Ozaki
—To obtain accurate results in numerical computation, high-precision arithmetic is a straightforward approach. However, most processors lack hardware support for floatingpoint formats beyond double precision (FP64). Double-word arithmetic (Dekker 1971) extends precision by using standard floating-point operations to…
Stephen M. Watt
The usual formulation of efficient division uses Newton iteration to compute an inverse in a related domain where multiplicative inverses exist. On one hand, Newton iteration allows quotients to be calculated using an efficient multiplication method. On the other hand, working in another domain is not always desirable…
Jeffrey Hurchalla
This paper presents an algorithm for the integer multiplicative inverse (mod 2 w ) which completes in the fewest cycles known for modern microprocessors, when using the native bit width w for the modulus 2 w . The algorithm is a modification of a method by Dumas, and for computers it slightly increases generality and…
Jérémy Berthomieu, Stef Graillat, Dimitri Lesnoff, Théo Mary
This article is concerned with the efficient computation of modular matrix multiplication C = AB mod p, a key kernel in computer algebra. We focus on floating-point arithmetic, which allows for using efficient matrix multiplication libraries. However, the existing approach is limited to primes p with bitsize at most…
Eric B. Olsen
Residue Number Systems (RNS) offer efficient modular arithmetic and natural parallelism, but direct integer division in RNS remains a difficult and comparatively underdeveloped operation. This paper builds on the type-II division algorithm of Szabo and Tanaka and reformulates it for more efficient hardware…
Fábio Lourenço Romano
numbers, using floating-point arithmetic Authors: ['Fábio Lourenço Romano'] In this paper, an optimized version of classical Bombelli's algorithm for computing integer square roots is presented. In particular, floating-point arithmetic is used to compute the initial guess of each digit of the root, following similar…
Carlos E. Valencia, Ralihe R. Villagrán
In [5] was given an algorithm that computes arithmetical structures on matrices. We use some of the ideas contained there to get an algorithm that computes arithmetical structures over dominated polynomials. A dominated polynomial is an integer multivariate polynomial such that contains a monomial which is divided by…
Piotr Miska, Nadir Murru, Giuliano Romeo
The problem of developing an arithmetic for continued fractions (in order to perform, e.g., sums and products) does not have a straightforward solution and has been addressed by several authors. In 1972, Gosper provided an algorithm to solve this problem. In this paper, we extend this approach in order to develop an…
Mayer Goldberg
This work presents and extends a known spigot-algorithm for computing square-roots, digit-by-digit, that is suitable for calculation by hand or an abacus, using only addition and subtraction. We offer an elementary proof of correctness for the original algorithm, then present a corresponding spigot-algorithm for…
Ishan Banerjee, Amites Sarkar
On Christmas Day 1640, Pierre de Fermat stated his famous "two-square" theorem in a letter to Marin Mersenne. This theorem, which was first stated by Girard in 1625, and first proved by Euler in 1749, states that an odd prime p is the sum of two squares if (and only if) p ≡ 1 (mod 4). Among the many admirers of the…
José I. Liberati
In the Chakravala method, discovered by Bhaskara II in the 12th century, the core idea is that given a triple (a, b, k) (which satisfies a 2 − db,2 = k), we can compose it with the trivial triple (m, 1, m2 − d) (by setting l = 1 in (1.2)) to obtain a new triple (am + db, a + b, m, k(m2 − d)) which can be scaled down by…
Dominic van der Zypen
Generalizing this to infinite bit-strings we get a binary operation on P(N), the power-set of N (which we identify with the collection of infinite bit-strings). We show that this operation is "group-like" in that it has a neutral element, inverses, but it is not associative. There are a lot of questions left, which the…
Mayank Deora, Pinakpani Pal
diophantine equations Authors: ['Mayank Deora' 'Pinakpani Pal'] Solving two variable linear diophantine equations has applications in many cryptographic protocols such as RSA and Elliptic curve cryptography. Extended euclid's algorithm is the most widely used algorithm to solve these equations. We revisit two…
M. Janos Uray
In this paper we analyze the computational cost of various operations performed symbolically in real algebraic number fields where the elements are represented as polynomials of a primitive element of the field. We give bounds on the costs in terms of several parameters, including the degree of the field and the…