Documentation |
On this page… |
---|
Prepare Your Algorithm for Code Acceleration or Code Generation Check for Fixed-Point Support for Functions Used in Your Algorithm Manage Data Types and Control Bit Growth |
A best practice for structuring your code is to separate your core algorithm from other code that you use to test and verify the results. Create a test file to call your original MATLAB^{®} algorithm and fixed-point versions of the algorithm. For example, as shown in the following table, you might set up some input data to feed into your algorithm, and then, after you process that data, create some plots to verify the results. Since you need to convert only the algorithmic portion to fixed-point, it is more efficient to structure your code so that you have a test file, in which you create your inputs, call your algorithm, and plot the results, and one (or more) algorithmic files, in which you do the core processing.
Original code | Best Practice | Modified code |
---|---|---|
% TEST INPUT x = randn(100,1); % ALGORITHM y = zeros(size(x)); y(1) = x(1); for n=2:length(x) y(n)=y(n-1) + x(n); end % VERIFY RESULTS yExpected=cumsum(x); plot(y-yExpected) title('Error') | Issue Generation of test input and verification of results are intermingled with the algorithm code. Fix Create a test file that is separate from your algorithm. Put the algorithm in its own function. | Test file % TEST INPUT x = randn(100,1); % ALGORITHM y = cumulative_sum(x); % VERIFY RESULTS yExpected = cumsum(x); plot(y-yExpected) title('Error') Algorithm in its own function function y = cumulative_sum(x) y = zeros(size(x)); y(1) = x(1); for n=2:length(x) y(n) = y(n-1) + x(n); end end |
You can use the test file to:
Verify that your floating-point algorithm behaves as you expect before you convert it to fixed point. The floating-point algorithm behavior is the baseline against which you compare the behavior of the fixed-point versions of your algorithm.
Propose fixed-point data types.
Compare the behavior of the fixed-point versions of your algorithm to the floating-point baseline.
Help you determine initial values for static ranges.
By default, the Fixed-Point Converter app shows code coverage results. Your test files should exercise the algorithm over its full operating range so that the simulation ranges are accurate. For example, for a filter, realistic inputs are impulses, sums of sinusoids, and chirp signals. With these inputs, using linear theory, you can verify that the outputs are correct. Signals that produce maximum output are useful for verifying that your system does not overflow. The quality of the proposed fixed-point data types depends on how well the test files cover the operating range of the algorithm with the accuracy that you want. Reviewing code coverage results help you verify that your test file is exercising the algorithm adequately. Review code flagged with a red code coverage bar because this code is not executed. If the code coverage is inadequate, modify the test file or add more test files to increase coverage. For more information see Code Coverage.
The automated conversion process instruments your code and provides data type proposals to help you convert your algorithm to fixed point.
MATLAB algorithms that you want to convert to fixed point automatically must comply with code generation requirements and rules. To view the subset of the MATLAB language that is supported for code generation, see Functions and Objects Supported for C and C++ Code Generation — Alphabetical List.
To help you identify unsupported functions or constructs in your MATLAB code, add the %#codegen pragma to the top of your MATLAB file. The MATLAB code analyzer flags functions and constructs that are not available in the subset of the MATLAB language supported for code generation. This advice appears in real-time as you edit your code in the MATLAB editor. For more information, see Check Code Using the MATLAB Code Analyzer. The software provides a link to a report that identifies calls to functions and the use of data types that are not supported for code generation. For more information, see Check Code Using the Code Generation Readiness Tool.
The Fixed-Point Converter app flags unsupported function calls found in your algorithm on the Function Replacements tab. For example, if you use the fft, which is not supported for fixed point. The tool adds an entry to the table on this tab and indicates that you need to specify a replacement function to use for fixed-point operations.
You can specify additional replacement functions. For example, functions like sin, cos,and sqrt may support fixed point, but for better efficiency, you may want to consider an alternative implementation like a lookup table or CORDIC-based algorithm. The app provides an option to generate lookup table approximations for continuous and stateless single-input, single-output functions in your original MATLAB® code. See Replacing Functions Using Lookup Table Approximations.
The automated fixed-point conversion process automatically manages data types and controls bit growth. It controls bit growth by using subscripted assignments, that is, assignments that use the colon (:) operator, in the generated code. When you use subscripted assignments, MATLAB overwrites the value of the left-hand side argument but retains the existing data type and array size. In addition to preventing bit growth, subscripted assignment reduces the number of casts in the generated fixed-point code and makes the code more readable.
Before you start the conversion, consider your goals for converting to fixed point. Are you implementing your algorithm in C or HDL? What are your target constraints? The answers to these questions determine many fixed-point properties such as the available word length, fraction length, and math modes, as well as available math libraries.
To set up these properties, use the Advanced settings.
For more information, see Specify Type Proposal Options.
Create a test file to validate that the floating-point algorithm works as expected before converting it to fixed point. You can use the same test file to propose fixed-point data types, and to compare fixed-point results to the floating-point baseline after the conversion. For more information, see Running a Simulation and Log Data for Histogram .
To fine-tune fixed-point type settings, use the histogram. To log data for histograms, in the app, click the arrow to the right of Run and Detect Ranges and select Log data for histogram.
After simulation and static analysis:
To view the histogram for a variable, on the Variables tab, click the Proposed Type field for that variable.
You can view the effect of changing the proposed data types by dragging the edges of the bounding box in the histogram window to change the proposed data type and selecting or clearing the Signed option.
If the values overflow and the range cannot fit the proposed type, the table shows proposed types in red.
When the tool applies data types, it generates an html report that provides overflow information and highlights overflows in red. Review the proposed data types.
fimath properties define the rules for performing arithmetic operations on fi objects, including math, rounding, and overflow properties. You can use the fimath ProductMode and SumMode properties to retain optimal data types for C or HDL. HDL can have arbitrary word length types in the generated HDL code whereas C requires container types (uint8, uint16, uint32). Use the Advanced settings, see Specify Type Proposal Options.
C. The KeepLSB setting for ProductMode and SumMode models the behavior of integer operations in the C language, while KeepMSB models the behavior of many DSP devices. Different rounding methods require different amounts of overhead code. Setting the RoundingMethod property to Floor, which is equivalent to two's complement truncation, provides the most efficient rounding implementation. Similarly, the standard method for handling overflows is to wrap using modulo arithmetic. Other overflow handling methods create costly logic. Whenever possible, set OverflowAction to Wrap.
MATLAB Code | Best Practice | Generated C Code | |
---|---|---|---|
Code being compiled function y = adder(a,b) y = a + b; end | Issue With the default word length set to 16 and the default fimath settings, additional code is generated to implement saturation overflow, nearest rounding, and full-precision arithmetic. | int adder(short a, short b) { int y; int i0; int i1; int i2; int i3; i0 = a; i1 = b; if ((i0 & 65536) != 0) { i2 = i0 | -65536; } else { i2 = i0 & 65535; } if ((i1 & 65536) != 0) { i3 = i1 | -65536; } else { i3 = i1 & 65535; } i0 = i2 + i3; if ((i0 & 65536) != 0) { y = i0 | -65536; } else { y = i0 & 65535; } return y; } | |
Fix To make the generated C code more efficient, choose fixed-point math settings that match your processor types. To customize fixed-point type proposals, use the Fixed-Point Converter app Advanced settings. Select fimath and then set: | int adder(short a, short b) { return a + b; } | ||
Rounding method | Floor | ||
Overflow action | Wrap | ||
Product mode | KeepLSB | ||
Sum mode | KeepLSB | ||
Product word length | 32 | ||
Sum word length | 32 |
HDL. For HDL code generation, set:
ProductMode and SumMode to FullPrecision
Overflow action to Wrap
Rounding method to Floor
Some MATLAB built-in functions can be made more efficient for fixed-point implementation. For example, you can replace a built-in function with a Lookup table implementation, or a CORDIC implementation, which requires only iterative shift-add operations. For more information, see Function Replacements.
Often, division is not fully supported by hardware and can result in slow processing. When your algorithm requires a division, consider replacing it with one of the following options:
Use bit shifting when the denominator is a power of two. For example, bitsra(x,3) instead of x/8.
Multiply by the inverse when the denominator is constant. For example, x*0.2 instead of x/5.
If the divisor is not constant, use a temporary variable for the division. Doing so results in a more efficient data type proposal and, if overflows occur, makes it easier to see which expression is overflowing.
For more efficient code, the automated fixed-point conversion process eliminates floating-point variables. The one exception to this is loop indices because they usually become integer types. It is good practice to inspect the fixed-point code after conversion to verify that there are no floating-point variables in the generated fixed-point code.
For the automated workflow, do not use explicit double or single casts in your MATLAB algorithm to insulate functions that do not support fixed-point data types. The automated conversion tool does not support these casts.
Instead of using casts, supply a replacement function. For more information, see Function Replacements.