COMPUTATIONALLY EFFICIENT MIXED PRECISION FLOATING POINT WAVEFORM GENERATION

Abstract

Computationally efficient mixed precision floating point waveform generation takes advantage of the high-speed generation of waveforms with single-precision floating point numbers while reducing the generally unacceptable loss of precision of pure single-precision floating point to generate any waveform that repeats in 2. This approaches computes a reference phase in double precision as the modulus of the phase with 2 and then computes offsets to that value in single precision. The double precision reference phase is recomputed as needed depending on how quickly the phase grows and how large a machine epsilon is desired.

Claims

1. A method of generating a waveform, comprising: for each sample number n and a reference sample number r, using a single-precision floating point unit (FPU) to compute a single-precision phase as an offset from a single-precision reference phase by a phase kernel indexed by the reference sample number r; using the single-precision FPU to compute a single-precision value for a waveform that repeats in 2 as a function of the single-precision phase for the current sample number n; incrementing the sample number n and the reference sample number r; if the reference sample number r does not meet a threshold, returning control to the single-precision FPU to compute the single-precision phase for the next sample number n; and if the threshold is met, using a double-precision FPU to compute a double-precision reference phase as a modulus of 2 of the phase kernel indexed by the sample number n, converting the double-precision reference phase to single-precision to update the single-precision reference phase and resetting the reference sample number r, and returning control to the single-precision FPU to compute the single-precision phase for the next sample number n.

2. The method of claim 1, wherein the phase kernel indexed by the reference sample number r is $=2**f_c\left(\frac{r}{f.Math.s}\right),$ the phase kernel indexed by the sample number n is $=2**f_c\left(\frac{n}{f.Math.s}\right),$ and the waveform is a sinusoid.

3. The method of claim 1, wherein the threshold is set to ensure that a machine epsilon of the single-precision phase value does not exceed a specified maximum.

4. The method of claim 3, wherein the threshold $r_{m.Math..Math.\mathrm{ax}}=\frac{2^{.Math.l.Math..Math.{\mathrm{og}}_2\left(2^{p-1}.Math.{\mathrm{.Math.}}_{m.Math..Math.\mathrm{ax}}\right).Math.-2.Math..Math.}}{}$ where .sub.max is the specified machine epsilon maximum for the single-precision phase value and p is the number of bits in the significand of the single precision phase value.

5. The method of claim 1, wherein the threshold is at least 1,000.

6. The method of claim 1, wherein all subsequent processing of the waveform is performed with single-precision floating point.

7. An apparatus for generating a waveform, comprising: a single-precision floating point unit (FPU) to compute a single-precision phase based as an offset from a single-precision reference phase by a phase kernel indexed by a reference sample number r and compute single-precision values for a waveform that repeats in 2 as a function of the single-precision phase for a sample number n; an arithmetic logic unit (ALU) configured to increment the sample number n and the reference sample number r and determine whether the reference sample number r meets a threshold, if not the ALU returns control to the single-precision FPU to compute the next phase and waveform value; and a double-precision FPU to compute a double-precision reference phase as a modulus of 2 of the phase kernel indexed by the sample number n and convert the double-precision reference phase to single-precision to update the single-precision reference phase, said ALU configured to reset the reference sample number r, and return control to the single-precision FPU to compute the next phase and waveform value.

8. The apparatus of claim 7, wherein the phase kernel indexed by the reference sample number r is $=2**f_c\left(\frac{r}{f.Math.s}\right),$ the phase kernel indexed by the sample number n is $=2**f_c\left(\frac{n}{f.Math.s}\right),$ and the waveform is a sinusoid.

9. The apparatus of claim 7, wherein the threshold is set to ensure that a machine epsilon of the single-precision phase value does not exceed a specified maximum.

10. The apparatus of claim 9, wherein the threshold $r_{m.Math..Math.\mathrm{ax}}=\frac{2^{.Math.l.Math..Math.{\mathrm{og}}_2\left(2^{p-1}.Math.{\mathrm{.Math.}}_{m.Math..Math.\mathrm{ax}}\right).Math.-2.Math..Math.}}{}$ where .sub.max is the specified machine epsilon maximum for the single-precision phase value and p is the number of bits in the significand of the single precision phase value.

11. The apparatus of claim 7, wherein the threshold is at least 1,000.

12. The apparatus of claim 7, further comprising additional single-precision FPU that perform all subsequent processing of the waveform in single-precision floating point.

13. A method of generating a waveform, comprising: for each sample number n and a reference sample number r, using an L-bit floating point unit (FPU) to compute an L-bit phase as an offset from an L-bit reference phase by a phase kernel indexed by the reference sample number r; using the an L-bit FPU to compute an L-bit value for a waveform that repeats in 2 as a function of the L-bit phase for the current sample number n; incrementing the sample number n and the reference sample number r; if the reference sample number r does not meet a threshold, returning control to the L-bit FPU to compute the L-bit phase for the next sample number n; and if the threshold is met, using an M-bit FPU, where M>L, to compute an M-bit reference phase as a modulus of 2 of the phase kernel indexed by the sample number n, converting the M-bit reference phase to L-bits to update the L-bit reference phase and resetting the reference sample number r, and returning control to the L-bit FPU to compute the L-bit phase for the next sample number n.

14. The method of claim 13, wherein the threshold is set to ensure that a machine epsilon of the L-bit phase value does not exceed a specified maximum.

15. The method of claim 13, wherein M is at least twice L.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0015] FIGS. 1a and 1b, as described above, illustrate the use of N-bit floating point units to generate a sinusoidal waveform;

[0016] FIGS. 2a and 2b illustrate the use of a mix of single-precision and double-precision floating point units to generate single-precision values of a waveform that repeats in 2;

[0017] FIG. 3 is a plot of a linear phase for a waveform and the linear phase modulus 2;

[0018] FIG. 4 is a diagram plots illustrating the computation of waveform values computing the reference phase modulus 2 plus the offset;

[0019] FIG. 5 is an exemplary plot and a highly zoomed portion of that plot of machine epsilon for conventional single-precision and double precision floating point and mixed-precision floating point in accordance with the present invention;

[0020] FIG. 6 is a plot of threshold values and computational efficiency versus maximum machine epsilon; and

[0021] FIG. 7 is a block diagram for mixed-precision waveform generation.

DETAILED DESCRIPTION OF THE INVENTION

[0022] The problem with known approaches for waveform generation is that the linear phase grows unbounded for waveforms that repeat in 2. Every time the phase value doubles, the machine epsilon doubles. In many applications, with single-precision floating point the machine epsilon is unacceptably large. Although the machine epsilon continues to grow with phase for double-precision, the number of bits allocated to the significand, 53 as compared to 24 for single-precision, keeps the machine epsilon acceptably small for almost all applications. However, the system must bear the burden of generating the waveform with double-precision and performing all subsequent processing with double-precision.

[0023] Computationally efficient mixed precision floating point waveform generation takes advantage of the high-speed generation of waveforms with single-precision floating point numbers while reducing the generally unacceptable loss of precision of pure single-precision floating point to generate any waveform that repeats in 2. This approaches computes a reference phase in double precision as the modulus of the phase with 2 and then computes offsets to that value in single precision. The double precision reference phase is recomputed as needed depending on how quickly the phase grows and how large a machine epsilon is desired. The waveform value is represented with single-precision. Thus, all downstream process can be performed with single-precision resulting in considerable computational savings. This approach can be generalized to use an L-bit FPU to compute the phase and waveform values and an M-bit FPU, where M>L, to compute the reference phase.

[0024] Without loss of generality the mixed-precision approach for waveform generation will be described using the IEEE standard definitions of single and double-precision. This mixed-precision technique is valid for any waveform that repeats in 2, which includes any sinusoidal signal. Therefore, any waveform that repeats in 2 may be represented as one or more sinusoidal signals to represent, for example, continuous wave (CW), Amplitude Modulation (AM), Single-Sideband Modulation (SSB), Frequency Modulation (FM), Pulse Modulation (PM), and Quadrature Amplitude Modulation (QAM). Any signal of interest may be represented using the mixed-precision technique. For simplicity, the exemplary embodiment generates a single sinusoid.

[0025] As shown in FIGS. 2a-2b, a mixed-precision waveform generator 200 uses a combination of single-precision and double-precision floating-point units 202 and 204, respectively, to generate the samples of the waveform 205 with a single-precision floating-point representation. In this example, the double-precision floating point unit 202 and an ALU 206 are part of a processor 208 and the single-precision floating point units 204 and an ALU 210 are part of another processor 212. Note, a single processor may support both the single and double-precision FPUs or a single FPU can function as either a 64 or 32 bit FPU. Both processors are driven by a clocking signal 214. All subsequent calculations are then computed in single-precision floating point. The use of double-precision floating point to calculate a reference phase for generating the samples of the transmitter signals establishes an initial higher level of accuracy. The resulting waveform values are more accurate than they would otherwise be using single-precision throughout. Computations are saved both in the initial generation of the waveform and in the subsequent processing of single-precision floating point numbers.

[0026] For each sample number n and a reference sample number r and initial ref=0, the single-precision floating point unit 204 computes a single-precision phase

[00003]$=2**f_c\left(\frac{r}{f.Math.s}\right)+.Math..Math.\mathrm{ref}$

(step 250) as an offset of a phase kernel

[00004]$2**f_c\left(\frac{r}{f.Math.s}\right)$

indexed by the reference sample number r from a single-precision reference phase ref where f.sub.c is the signal frequency and f.sub.s is the sampling frequency and then computes the waveform value 205 S=A*cos(B) or S=A*e.sup.i (step 252) based on that phase. An arithmetic logic unit (ALU) increments the sample number n and the reference sample number r (step 254) and if the reference sample number does not meet a threshold (step 256) returns control to the single-precision floating point unit 204 to compute a next phase (step 250). If the threshold is met, the double-precision floating point unit 202 updates the reference phase

[00005]$.Math.\mathrm{ref}=\mathrm{mod}\left(2**f_c\left(\frac{n}{f.Math.s}\right),.Math.2.Math.\right)$

(step 258) by computing the phase kernel indexed by the incremented sample number n and taking the modulus with 2. The ALU resets the reference sample number r=0 and FPU 204 casts the double-precision reference phase to single-precision (step 260) and returns control to the single-precision floating point unit 204 to compute the phase for the next sample number (step 250). The threshold is determined by a maximum machine epsilon to keep the phase value bounded.

[0027] It is worth noting that if FPUs of equal precision were employed to perform all of the calculations of FIG. 2b, the flow diagram degenerates to that shown in FIG. 1b. It is the combination of using different precision FPUs with the methodology depicted in FIG. 2b; high precision FPU to compute the reference phase and low precision FPU to compute the offset that achieves the improvement in accuracy and computational efficiency. It is critical that the calculation of the reference phase is performed with greater precision than the computation of the phase and waveform values. Computational efficiency for waveform generation is achieved by the fact that the inner loop of computing the single-precision phase is performed far more often than the outer loop of computing the double-precision reference phase.

[0028] The approach of computing a reference phase in double precision as the modulus of the phase with 2 and then computing offsets to that value in single precision until the reference phase needs to be updated is illustrated in FIGS. 3 and 4. As shown in FIG. 3, a linear phase 300 of a waveform that repeats in 2 increases with time without bound. The linear phase modulus 2 provides a reference phase 302 that is bounded by 0 and 2. Therefore the range, hence machine epsilon of the phase is bounded. When this is coupled with a double-precision floating point computation of the reference phase 302 the machine epsilon of the phase remains small, approximately that of a pure double-precision floating point calculation.

[0029] In an exemplary calculation, as shown in FIG. 4, the threshold is set to 5. The single-precision FPU calculates the phase and waveform value for the first five samples based on an initial reference phase of zero as depicted by inner loop 350. The double-precision FPU than calculates the reference phase 352. The single-precision FPU calculates the phase and waveform value for the next five samples as offsets from the reference phase 352 as depicted by outer loop 354 and the process repeats.

[0030] As previously described, for any floating point representation of a number, the smallest difference between any two numbers is known as the machine epsilon. The machine epsilon is dependent on the base of the number system, the magnitude of the value, and the size of the significand. The machine epsilon of a base-2 floating point number () of value n with significand size p is

[00006]$\mathrm{.Math.}=\frac{2^{.Math.l.Math..Math.{\mathrm{og}}_2.Math.n.Math.}}{2^{p-1}}.$

Solving for a desired machine epsilon .sub.max with precision p, the upper limit for phase value is .sub.max=2.sup.log.sup.2.sup.2.sup.p-1.sup..sup.max.sup.). Assuming that the starting phase of the signal is <2, and increments by no more than at every iteration (keeping at or below Nyquist frequency), the maximum threshold for updating the single precision reference phase is:

[00007]$\begin{array}{cc}{r}_{m.Math..Math.\mathrm{ax}}=\frac{2^{.Math.l.Math..Math.{\mathrm{og}}_2\left(2^{p-1}.Math.{\mathrm{.Math.}}_{m.Math..Math.\mathrm{ax}}\right).Math.-2.Math..Math.}}{}& \left(2\right)\end{array}$

[0031] The enhanced accuracy of mixed precision floating point waveform generation is illustrated in comparison to straight single-precision and double precision in FIG. 5. In this example, the input value (or range) is set at 6,000,000 and the machine epsilon is calculated for that value. A double-precision floating point representation 400 closely approximates the real number for the phase value and the mixed-precision floating point representation 402 is virtually indistinguishable. The machine epsilon for both is very small. As shown in the inset to FIG. 5b, the error 404 of the mixed-precision value from the double-precision value is extremely small. By contrast, the single-precision floating point representation 406 is a coarse stair-case approximate to the real phase value because the machine epsilon is relatively large. The near double-precision accuracy of the mixed-precision calculation of the phase carries over to an improved accuracy of the waveform values.

[0032] As shown in FIG. 6, the threshold value r.sub.max 450 follows an increasing stair-case function with increasing maximum machine epsilon .sub.max. The larger the allowable maximum machine epsilon the greater the number of single-precision calculations of the offset phase that can be made before the double-precision reference phase must be updated. The computational efficiency 452 as measured by the percent of iterations requiring double-precision is inversely proportional to the maximum machine epsilon and the threshold. As the threshold increases, computational efficiency improves. For a given application, a system designer must trade the desire for mixed-precision representation that approximates double-precision performance against a computational efficiency that approximates single-precision performance. The benefit of the mixed-precision approach is that it can achieve most of the performance benefits of double-precision and most of the computational efficiency benefits of single-precision in waveform generation. For example, in some cases the threshold value may be at least 1,000. The percentage of double-precisions calculations is very small yet the accuracy closely approximates double-precision.

[0033] As shown in FIG. 7, most systems and applications include both the generation of a waveform in a floating point representation (step 500) and subsequent downstream processing of that waveform in its native floating point representation (step 502). In addition to the initial benefits in improved accuracy and computational efficiency to generate the waveform, a major benefit of the mixed-precision floating point computation are the improved accuracy that propagates down stream and the reduction in the computational burden. Accuracy is improved because the waveform itself approximates the double-precision floating point representation. Computational efficiency is improved because all subsequent processing is performing using single precision floating point rather than double precision.

[0034] While several illustrative embodiments of the invention have been shown and described, numerous variations and alternate embodiments will occur to those skilled in the art. Such variations and alternate embodiments are contemplated, and can be made without departing from the spirit and scope of the invention as defined in the appended claims.