TL;DR: In this paper, on-line algorithms for division and multiplication are developed and it is assumed that the operands as well as the result flow through the arithmetic unit in a digit-by-digit, most significant digit first fashion.
Abstract: In this paper, on-line algorithms for division and multiplication are developed. It is assumed that the operands as well as the result flow through the arithmetic unit in a digit-by-digit, most significant digit first fashion. The use of a redundant digit set, at least for the digits of the result, is required.
TL;DR: The four basic arithmetic algorithms for Hensel's code are described and their application to rational matrix computations is demonstrated by solving a system of linear equations exactly, using the Gaussian elimination procedure.
Abstract: A unique code (called Hensel's code) is derived for a rational number by truncating its infinite p-adic expansion. The four basic arithmetic algorithms for these codes are described and their application to rational matrix computations is demonstrated by solving a system of linear equations exactly, using the Gaussian elimination procedure.
TL;DR: In this paper, the authors consider the problem of computing a sequence of disjoint intervals (rectangles) of length (width) ∈ or less, each containing exactly one real (complex) zero of a polynomial A(x)∈D[x] and any positive rational number ∈, and together containing all real (composite) zeros.
Abstract: Let D be a Euclidean domain which is a subring of the field of complex numbers and in which the arithmetic operations can be algorithmically performed. Examples of D are Z , the integers, Q , the rational numbers, G , the Gaussian integers, and P , the real algebraic numbers. We are concerned with algorithms which, given any polynomial A(x)∈D[x] and any positive rational number ∈, computes a sequence of disjoint intervals (rectangles) of length (width) ∈ or less, each containing exactly one real (complex) zero of A, and together containing all real (complex) zeros. The algorithms must also compute the multiplicity of each zero. Any algorithm strictly fulfilling all of these specification may truly be described as infallible. While retaining these lofty goals we are nevertheless concerned with algorithms which are as efficient as possible, in practice as well as theory. Such infallible algorithms can be derived quite readily from several diverse mathematical theorems using “exact” arithmetic in the domain D, but these algorithms may differ significantly in their time complexities. Also, one may further improve practical efficiency by appropriate substitution of “approximate” arithmetic (e.g., interval arithmetic) for exact arithmetic in certain contexts without sacrificing infallibility (using exact arithmetic as backup). This paper surveys recent progress on this problem, starting from Heindel's 1970 implementation and analysis of Sturm's theorem, including Pinkert's 1973 application of Sturm sequences to complex zeros, the 1975 Collins-Loos algorithm based on Rolle's theorem, the 1976 Collins-Akritas modification of Uspensky's method based on Descartes' theorem, and concluding with a brief report on current research by Collins and Chou relating to use of approximate arithmetic and the “principle of argument”. Both theoretical time bounds and empirical times are presented.
TL;DR: Techniques of interval arithmetic can be used to give bounds on the global minima of unconstrained optimization problems using the design of a hypothetical chemical plant.
Abstract: We show how techniques of interval arithmetic can be used to give bounds on the global minima of unconstrained optimization problems. We illustrate the techniques using the design of a hypothetical chemical plant.
TL;DR: In this article, a real-time arithmetic such as vector conversion and display its results in a display system, by performing various arithmetic processes in the process of reading from a screen memory and displaying.
Abstract: PURPOSE: To execute real-time arithmetic such as vector conversion and display its results in a display system, by performing various arithmetic processes in the process of reading from a screen memory and displaying. COPYRIGHT: (C)1978,JPO&Japio
TL;DR: The operation and implementation of the arithmetic proof rule for the quantifier free integer arithmetic used in the PL/CV 2 program verification system is described and the general arithmetic satisfiability problem underlying the rule is shown to be NP complete.
Abstract: This paper describes the operation and implementation of the arithmetic proof rule for the quantifier free integer arithmetic used in the PL/CV 2 program verification system. The general arithmetic satisfiability problem underlying the rule is shown to be NP complete.
TL;DR: In this paper, the congruence (f ∘ g) = 0 (mod n) and the functional equation f ∘ f ∞ ∘ ∞ g ∞ = g, where ∘ is an exponential regular convolution, are studied.
Abstract: In this paper the congruence (f ∘ g)(n) = 0 (mod n) and the functional equation f ∘ f ∘ … ∘ f = g, are studied, where ∘ is an exponential regular convolution.
TL;DR: The structure for hardware realisation of Fast Walsh-Fourier Transform (FWFT) is presented and two alternative structures are proposed one based on PCM encoding and serial two's complement arithmetic, whilst the other is based on delta-sigma encoding with appropriate arithmetic operations.
Abstract: The structure for hardware realisation of Fast Walsh-Fourier Transform (FWFT) is presented in this paper. The development of the approach is based on Shanks' algorithm for FWFT and use is made in the realisation of serial storage and simple arithmetic. Two alternative structures are proposed one based on PCM encoding and serial two's complement arithmetic, whilst the other is based on delta-sigma encoding with appropriate arithmetic operations. Comparison of these alternative solutions is given in terms of hardware requirements and the mean square error produced.
TL;DR: The papers consider the spaces which occur in numerical computations on computers depending on a properly defined computer arithmetic, and give a summary of the main ideas of these papers.
Abstract: During recent years a number of papers concerning a mathematical foundation of computer arithmetic have been written. Some of these papers are still unpublished. The papers consider the spaces which occur in numerical computations on computers depending on a properly defined computer arithmetic. The following treatment gives a summary of the main ideas of these papers. Many of the proofs had to be sketched or completely omitted. In such cases the full information can be found in the references.
TL;DR: An algorithm is presented to compute inverses modulo a power of a polynomial in terms of inverse modulo that polynometric, used to make an improvement in the Kung-Tong partial fraction decomposition algorithm.
Abstract: Algorithms for general partial fraction decomposition are obtained by using modular polynomial arithmetic. An algorithm is presented to compute inverses modulo a power of a polynomial in terms of inverses modulo that polynomial. This algorithm is used to make an improvement in the Kung-Tong partial fraction decomposition algorithm.
TL;DR: A combined arithmetic unit and language support system which allows user specifications of the arithmetic and allows one set of polymorphic arithmetic operators to be defined across all combinations of user-defined operand types.
Abstract: This paper describes a combined arithmetic unit and language support system which allows user specifications of the arithmetic. Limited extensions to a high-level language, in connection with a generalized underlying arithmetic unit, allow a single skeletal unified numeric operand type to be refined into a variety of data types. The interpretation of operands by the operators is based on type descriptors, allowing one set of polymorphic arithmetic operators to be defined across all combinations of user-defined operand types. The arithmetic unit is realized in microcode to achieve efficiency.
TL;DR: It is shown how careful use of single precision, floating point arithmetic within the context of a slightly modified algorithm can make the calculation considerably faster and no less exact.
Abstract: THe Collins-Loos algorithm for computing isolating intervals for the zeros of an integer polynomial requires the evaluation of polynomials at rational points. This implies the use of arbitrary precision integer arithmetic. It is shown how careful use of single precision, floating point arithmetic within the context of a slightly modified algorithm can make the calculation considerably faster and no less exact. Typically, 95% or more of the evaluations can be done without exact arithmetic. The precise speedup depends on the relative costs of the arithmetic in a given implementation. The implementation on DEC KL-10 computer is some 5 to 10 times faster than the original Univac 1110 implementation in SAC-I.
TL;DR: Embedding, a topic discussed sporadically in the literature [l,2] most recently by APL 8nthusiasts, is the extension of the functional domain of a programming language, to includ e d a t a structures and operations not previously supported, by writin ~! subroutines only in the language being extended.
Abstract: Embedding, a topic discussed sporadically in the literature [l,2] most recently by APL 8nthusiasts [3,4], is the extension o f the functional domain of a programming language, to includ e d a t a structures and operations not previously supported, by writin ~! subroutines only in the language being extended. This latter s e r ves tO distinguish embedding from the prOCeSS Of exten d in g a lang ua CJC' by re-writing parts of the processor for that lan gu a ge .
TL;DR: It is illustrated that conflicts do arise and a way of trying to overcome them is presented to highlight the importance of integrating the new orderings to the old ones as conflicting requirements can occur.
Abstract: Sussman (1973) has shown that strategies could be reordered according to experience. In his examples no conflicts were experienced after reordering. We shall illustrate that conflicts do arise and present a way of trying to overcome them. The program written (ELS) accepts a set of partially ordered Horn clauses as input and a sequence of problems to solve. The search tree obtained is examined in order to generate a priority ordering for the clauses used. If conflicts are detected, ELS tries to determine the conditions under which the clauses involved should be applied. For example , suppose the problem 3+2=X is to be solved and we are given the clauses: assoc: subs: subz: sue : pred: Generally, the problems are solved in a breadth first manner, but priority orderings between individual clauses are respected. First all clauses that match the given subgoal are chosen. In our case both the clauses 'subs' and 'subz' match the initial goal 3+2-X, as this goal matches X1+X2=X3. Both clauses are actually tried in the search, since initially there is no priority ordering specified between them. If it had been specified that 'subz
TL;DR: This paper describes an implementation of Triplex arithmetic in single and multiple precision based on the Augment precompiler to obtain an easily used package.
Abstract: : Triplex arithmetic is a variation of interval arithmetic in which a main or rounded value is computed in addition to the end-points of the containing interval. The main value may be considered, depending on the application, to be the most probable result of the computation, with the interval bounds indicating the possible error. This paper describes an implementation of Triplex arithmetic in single and multiple precision. The implementation is based on the Augment precompiler to obtain an easily used package. (Author)
TL;DR: A fast-division algorithm, effectively implementable with available IC technology, that provides simple control provisions for the division of numbers in any format and to any fractional precision.
Abstract: A fast-division algorithm, effectively implementable with available IC technology, has been developed and implemented. Especially suited for asynchronous division units, the skip-and-set algorithm provides simple control provisions for the division of numbers in any format and to any fractional precision.
TL;DR: The specification of a truly portable implementation independent programming language requires that digital arithmetic be specified in an implementation independent manner.
Abstract: The specification of a truly portable implementation independent programming language requires that digital arithmetic be specified in an implementation independent manner.
Variable length storage is proposed for digital numbers. Linked storage and linked overflow storage are proposed so as to enable the use of variable length numbers. It is suggested that significant digit arithmetic be used for single numbers and interval arithmetic be used for interval pairs. The question is put to numerical mathematicians to specify the decision rules to be used in deciding on the number of digits to be retained when computing in either the significant digit mode or in the interval mode. The question is also raised as to what other arithmetic modes should be provided for.
The speculation is offered not for the purposes of suggesting useful hardware design techniques but rather to indicate feasibility of some design. If the implementation independent arithmetic is considered feasible, then language designers can design a language freed from present hardware restrictions.
TL;DR: A parallel computational method is described that provides a simple and fast algorithm for the evaluation of polynomials, certain rational functions and arithmetic expressions, solving a class of systems of linear equations, or performing the basic arithmetic operations in a fixed-point number representation system.
Abstract: A parallel computational method, amenable for efficient hardware-level implementation, is described. It provides a simple and fast algorithm for the evaluation of polynomials, certain rational functions and arithmetic expressions, solving a class of systems of linear equations, or performing the basic arithmetic operations in a fixed-point number representation system. The time required to perform the computation is of the order of m carry-free addition operations, m being the number of digits in the solution. In particular, the method is suitable for fast evaluation of mathematical functions in hardware.
TL;DR: This paper discusses the application of the residue number system to realizing digital signal processing elements using such arrays and advantages and disadvantages over conventional realizations are discussed.
Abstract: In the past, hardware realization of digital signal processing elements have been based upon binary arithmetic concepts. Because of the dependence between digits in binary arithmetic operations, the hardware required to construct arithmetic elements is cumbersome. In the residue number system, arithmetic operations can be performed with complete independence between digits and a corresponding reduction in hardware complexity. In fact, using current technology, arithmetic operations can be carried out using arrays of look-up tables placed in high density ROMs. This paper discusses the application of the residue number system to realizing digital signal processing elements using such arrays and advantages and disadvantages over conventional realizations are discussed. Examples are given of recursive filter and FFT butterfly element realization.