Integrals along lines passing through arrays of integer data are used in many applications. Existing algorithms employ parametric methods and floating point calculations to determine the data values that contribute to the line integrals, along with some form of interpolation to weight these data values. A new, nonparametric integer‐arithmetic noninterpolating algorithm (NI0) and an extension using first‐order interpolation (NI1) are presented here. These algorithms are compared for accuracy and speed with both Siddon’s [Med. Phys. 12, 252 (1985)] parametric floating point algorithm using no interpolation (PF0) and our extension using first‐order interpolation (PF1). NI1 gives line integral values significantly closer to those of PF1 than does PF0 and runs nine times faster on a VAXstation 2000. NI0 gives less accurate line integral values than does PF0 but runs 19 times faster.
Skip Nav Destination
Article navigation
Research Article|
March 01 1990
A comparison of line integral algorithms
J. Ross Mitchell;
J. Ross Mitchell
Allan Blair Memorial Clinic, Regina, Saskatchewan S4T 7T1, Canada
Department of Computer Science, University of Regina, Regina, Saskatchewan S4S 0A2, Canada
Search for other works by this author on:
Peter Dickof;
Peter Dickof
Allan Blair Memorial Clinic, Regina, Saskatchewan S4T 7T1, Canada
Department of Computer Science, University of Regina, Regina, Saskatchewan S4S 0A2, Canada
Search for other works by this author on:
Alan G. Law
Alan G. Law
Allan Blair Memorial Clinic, Regina, Saskatchewan S4T 7T1, Canada
Department of Computer Science, University of Regina, Regina, Saskatchewan S4S 0A2, Canada
Search for other works by this author on:
Comput. Phys. 4, 166–172 (1990)
Article history
Received:
March 29 1989
Accepted:
August 08 1989
Citation
J. Ross Mitchell, Peter Dickof, Alan G. Law; A comparison of line integral algorithms. Comput. Phys. 1 March 1990; 4 (2): 166–172. https://doi.org/10.1063/1.168381
Download citation file:
Citing articles via
Related Content
Global sensitivity analysis of a nonlinear system using animated basins of attraction
Comput. Phys. (November 1992)
Using multiple RISC CPUs in parallel to study charm quarks
Comput. Phys. (July 1992)