1 // SPDX-License-Identifier: MIT
2 /*
3 * Copyright 2023 Advanced Micro Devices, Inc.
4 *
5 * Permission is hereby granted, free of charge, to any person obtaining a
6 * copy of this software and associated documentation files (the "Software"),
7 * to deal in the Software without restriction, including without limitation
8 * the rights to use, copy, modify, merge, publish, distribute, sublicense,
9 * and/or sell copies of the Software, and to permit persons to whom the
10 * Software is furnished to do so, subject to the following conditions:
11 *
12 * The above copyright notice and this permission notice shall be included in
13 * all copies or substantial portions of the Software.
14 *
15 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
16 * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
17 * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
18 * THE COPYRIGHT HOLDER(S) OR AUTHOR(S) BE LIABLE FOR ANY CLAIM, DAMAGES OR
19 * OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE,
20 * ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR
21 * OTHER DEALINGS IN THE SOFTWARE.
22 *
23 * Authors: AMD
24 *
25 */
26 #include "dm_services.h"
27 #include "custom_float.h"
28
build_custom_float(struct fixed31_32 value,const struct custom_float_format * format,bool * negative,uint32_t * mantissa,uint32_t * exponenta)29 static bool build_custom_float(struct fixed31_32 value,
30 const struct custom_float_format *format,
31 bool *negative,
32 uint32_t *mantissa,
33 uint32_t *exponenta)
34 {
35 uint32_t exp_offset = (1 << (format->exponenta_bits - 1)) - 1;
36
37 const struct fixed31_32 mantissa_constant_plus_max_fraction =
38 dc_fixpt_from_fraction((1LL << (format->mantissa_bits + 1)) - 1,
39 1LL << format->mantissa_bits);
40
41 struct fixed31_32 mantiss;
42
43 if (dc_fixpt_eq(value, dc_fixpt_zero)) {
44 *negative = false;
45 *mantissa = 0;
46 *exponenta = 0;
47 return true;
48 }
49
50 if (dc_fixpt_lt(value, dc_fixpt_zero)) {
51 *negative = format->sign;
52 value = dc_fixpt_neg(value);
53 } else {
54 *negative = false;
55 }
56
57 if (dc_fixpt_lt(value, dc_fixpt_one)) {
58 uint32_t i = 1;
59
60 do {
61 value = dc_fixpt_shl(value, 1);
62 ++i;
63 } while (dc_fixpt_lt(value, dc_fixpt_one));
64
65 --i;
66
67 if (exp_offset <= i) {
68 *mantissa = 0;
69 *exponenta = 0;
70 return true;
71 }
72
73 *exponenta = exp_offset - i;
74 } else if (dc_fixpt_le(mantissa_constant_plus_max_fraction, value)) {
75 uint32_t i = 1;
76
77 do {
78 value = dc_fixpt_shr(value, 1);
79 ++i;
80 } while (dc_fixpt_lt(mantissa_constant_plus_max_fraction, value));
81
82 *exponenta = exp_offset + i - 1;
83 } else {
84 *exponenta = exp_offset;
85 }
86
87 mantiss = dc_fixpt_sub(value, dc_fixpt_one);
88
89 if (dc_fixpt_lt(mantiss, dc_fixpt_zero) ||
90 dc_fixpt_lt(dc_fixpt_one, mantiss))
91 mantiss = dc_fixpt_zero;
92 else
93 mantiss = dc_fixpt_shl(mantiss, format->mantissa_bits);
94
95 *mantissa = dc_fixpt_floor(mantiss);
96
97 return true;
98 }
99
setup_custom_float(const struct custom_float_format * format,bool negative,uint32_t mantissa,uint32_t exponenta,uint32_t * result)100 static bool setup_custom_float(const struct custom_float_format *format,
101 bool negative,
102 uint32_t mantissa,
103 uint32_t exponenta,
104 uint32_t *result)
105 {
106 uint32_t i = 0;
107 uint32_t j = 0;
108 uint32_t value = 0;
109
110 /* verification code:
111 * once calculation is ok we can remove it
112 */
113
114 const uint32_t mantissa_mask =
115 (1 << (format->mantissa_bits + 1)) - 1;
116
117 const uint32_t exponenta_mask =
118 (1 << (format->exponenta_bits + 1)) - 1;
119
120 if (mantissa & ~mantissa_mask) {
121 BREAK_TO_DEBUGGER();
122 mantissa = mantissa_mask;
123 }
124
125 if (exponenta & ~exponenta_mask) {
126 BREAK_TO_DEBUGGER();
127 exponenta = exponenta_mask;
128 }
129
130 /* end of verification code */
131
132 while (i < format->mantissa_bits) {
133 uint32_t mask = 1 << i;
134
135 if (mantissa & mask)
136 value |= mask;
137
138 ++i;
139 }
140
141 while (j < format->exponenta_bits) {
142 uint32_t mask = 1 << j;
143
144 if (exponenta & mask)
145 value |= mask << i;
146
147 ++j;
148 }
149
150 if (negative && format->sign)
151 value |= 1 << (i + j);
152
153 *result = value;
154
155 return true;
156 }
157
convert_to_custom_float_format(struct fixed31_32 value,const struct custom_float_format * format,uint32_t * result)158 bool convert_to_custom_float_format(struct fixed31_32 value,
159 const struct custom_float_format *format,
160 uint32_t *result)
161 {
162 uint32_t mantissa;
163 uint32_t exponenta;
164 bool negative;
165
166 return build_custom_float(value, format, &negative, &mantissa, &exponenta) &&
167 setup_custom_float(format,
168 negative,
169 mantissa,
170 exponenta,
171 result);
172 }
173
174