7#include "xfa/fgas/crt/cfgas_decimal.h"
15#include "core/fxcrt/fx_extension.h"
16#include "third_party/base/check.h"
18#define FXMATH_DECIMAL_SCALELIMIT 0x1c
19#define FXMATH_DECIMAL_RSHIFT32BIT(x) ((x) >> 0x10
>> 0x10
)
20#define FXMATH_DECIMAL_LSHIFT32BIT(x) ((x) << 0x10
<< 0x10
)
24inline uint8_t decimal_helper_div10(uint64_t& phi,
37inline uint8_t decimal_helper_div10_any(uint64_t nums[], uint8_t numcount) {
39 for (
int i = numcount - 1; i > 0; i--) {
44 retVal = nums[0] % 0xA;
50inline void decimal_helper_mul10(uint64_t& phi, uint64_t& pmid, uint64_t& plo) {
55 pmid = (uint32_t)pmid;
58inline void decimal_helper_mul10_any(uint64_t nums[], uint8_t numcount) {
60 for (
int i = 1; i < numcount; i++) {
62 nums[i - 1] = (uint32_t)nums[i - 1];
66inline void decimal_helper_normalize(uint64_t& phi,
70 pmid = (uint32_t)pmid;
74 pmid = (uint32_t)pmid;
77inline void decimal_helper_normalize_any(uint64_t nums[], uint8_t len) {
78 for (
int i = len - 2; i > 0; i--) {
80 nums[i] = (uint32_t)nums[i];
82 for (
int i = 0; i < len - 1; i++) {
84 nums[i] = (uint32_t)nums[i];
88inline int8_t decimal_helper_raw_compare_any(uint64_t a[],
93 for (
int i =
std::max(al - 1, bl - 1); i >= 0; i--) {
94 uint64_t l = (i >= al ? 0 : a[i]), r = (i >= bl ? 0 : b[i]);
95 retVal += (l > r ? 1 : (l < r ? -1 : 0));
102inline void decimal_helper_dec_any(uint64_t a[], uint8_t al) {
103 for (
int i = 0; i < al; i++) {
109inline void decimal_helper_inc_any(uint64_t a[], uint8_t al) {
110 for (
int i = 0; i < al; i++) {
112 if ((uint32_t)a[i] == a[i])
118inline void decimal_helper_raw_mul(uint64_t a[],
124 DCHECK(al + bl <= cl);
125 for (
int i = 0; i < cl; i++)
128 for (
int i = 0; i < al; i++) {
129 for (
int j = 0; j < bl; j++) {
130 uint64_t m = (uint64_t)a[i] * b[j];
131 c[i + j] += (uint32_t)m;
135 for (
int i = 0; i < cl - 1; i++) {
137 c[i] = (uint32_t)c[i];
139 for (
int i = 0; i < cl; i++)
140 c[i] = (uint32_t)c[i];
143inline void decimal_helper_raw_div(uint64_t a[],
149 for (
int i = 0; i < cl; i++)
152 uint64_t left[16] = {0};
153 uint64_t right[16] = {0};
155 for (
int i = 0; i < al; i++)
159 while (decimal_helper_raw_compare_any(left, al, right, al) <= 0) {
161 for (
int i = 0; i < al; i++)
162 cur[i] = left[i] + right[i];
164 for (
int i = al - 1; i >= 0; i--) {
170 decimal_helper_raw_mul(cur, al, b, bl, tmp, 16);
171 switch (decimal_helper_raw_compare_any(tmp, 16, a, al)) {
173 for (
int i = 0; i < 16; i++)
177 decimal_helper_normalize_any(left, al);
180 for (
int i = 0; i < 16; i++)
182 decimal_helper_dec_any(right, al);
185 for (
int i = 0; i <
std::min(al, cl); i++)
190 for (
int i = 0; i <
std::min(al, cl); i++)
194inline bool decimal_helper_outofrange(uint64_t a[], uint8_t al, uint8_t goal) {
195 for (
int i = goal; i < al; i++) {
202inline void decimal_helper_shrinkintorange(uint64_t a[],
206 bool bRoundUp =
false;
208 decimal_helper_outofrange(a, al, goal))) {
209 bRoundUp = decimal_helper_div10_any(a, al) >= 5;
213 decimal_helper_normalize_any(a, goal);
214 decimal_helper_inc_any(a, goal);
218inline void decimal_helper_truncate(uint64_t& phi,
222 uint8_t minscale = 0) {
223 while (scale > minscale) {
224 uint64_t thi = phi, tmid = pmid, tlo = plo;
225 if (decimal_helper_div10(thi, tmid, tlo) != 0)
241 m_uLo(
static_cast<uint32_t>(val)) {}
244 : m_uLo(
static_cast<uint32_t>(val)) {}
260 }
else if (val ==
std::numeric_limits<int32_t>::min()) {
270 float newval = fabs(val);
271 float divisor = powf(2.0, 64.0f);
272 uint64_t bottom64 =
static_cast<uint64_t>(fmodf(newval, divisor));
273 uint64_t top64 =
static_cast<uint64_t>(newval / divisor);
274 uint64_t plo = bottom64 & 0xFFFFFFFF;
275 uint64_t pmid = bottom64 >> 32;
276 uint64_t phi = top64 & 0xFFFFFFFF;
278 newval = fmodf(newval, 1.0f);
279 for (uint8_t iter = 0; iter < scale; iter++) {
280 decimal_helper_mul10(phi, pmid, plo);
282 plo +=
static_cast<uint64_t>(newval);
283 newval = fmodf(newval, 1.0f);
287 decimal_helper_normalize(phi, pmid, plo);
288 m_uHi =
static_cast<uint32_t>(phi);
289 m_uMid =
static_cast<uint32_t>(pmid);
290 m_uLo =
static_cast<uint32_t>(plo);
296 const wchar_t* str = strObj.unterminated_c_str();
297 const wchar_t* strBound = str + strObj.GetLength();
298 bool pointmet =
false;
301 while (str != strBound && *str ==
' ')
303 if (str != strBound && *str ==
'-') {
306 }
else if (str != strBound && *str ==
'+') {
318 m_uLo = m_uLo * 0xA + (*str -
'0');
329 WideString retString;
331 uint64_t phi = m_uHi;
332 uint64_t pmid = m_uMid;
333 uint64_t plo = m_uLo;
334 while (phi || pmid || plo)
335 tmpbuf
+= decimal_helper_div10(phi, pmid, plo) +
'0';
337 uint8_t outputlen = (uint8_t)tmpbuf.GetLength();
338 uint8_t scale = m_uScale;
339 while (scale >= outputlen) {
346 for (uint8_t idx = 0; idx < outputlen; idx++) {
347 if (idx == (outputlen - scale) && scale != 0)
349 retString
+= tmpbuf
[outputlen - 1 - idx
];
359 double pow = (
double)(1 << 16) * (1 << 16);
360 double base =
static_cast<
double>(m_uHi) * pow * pow +
361 static_cast<
double>(m_uMid) * pow +
static_cast<
double>(m_uLo);
362 return (m_bNeg ? -1 : 1) * base * powf(10.0f, -m_uScale);
366 uint8_t oldscale = m_uScale;
367 if (oldscale == newscale)
370 uint64_t phi = m_uHi;
371 uint64_t pmid = m_uMid;
372 uint64_t plo = m_uLo;
373 if (newscale > oldscale) {
374 for (uint8_t iter = 0; iter < newscale - oldscale; iter++)
375 decimal_helper_mul10(phi, pmid, plo);
377 m_uHi =
static_cast<uint32_t>(phi);
378 m_uMid =
static_cast<uint32_t>(pmid);
379 m_uLo =
static_cast<uint32_t>(plo);
383 uint64_t point5_hi = 0;
384 uint64_t point5_mid = 0;
385 uint64_t point5_lo = 5;
386 for (uint8_t iter = 0; iter < oldscale - newscale - 1; iter++)
387 decimal_helper_mul10(point5_hi, point5_mid, point5_lo);
392 decimal_helper_normalize(phi, pmid, plo);
393 for (uint8_t iter = 0; iter < oldscale - newscale; iter++)
394 decimal_helper_div10(phi, pmid, plo);
396 m_uHi =
static_cast<uint32_t>(phi);
397 m_uMid =
static_cast<uint32_t>(pmid);
398 m_uLo =
static_cast<uint32_t>(plo);
409 uint64_t a[3] = {m_uLo, m_uMid, m_uHi},
410 b[3] = {val.m_uLo, val.m_uMid, val.m_uHi};
412 decimal_helper_raw_mul(a, 3, b, 3, c, 6);
413 bool neg = m_bNeg ^ val.m_bNeg;
414 uint8_t scale = m_uScale + val.m_uScale;
415 decimal_helper_shrinkintorange(c, 6, 3, scale);
416 return CFGAS_Decimal(
static_cast<uint32_t>(c[0]),
static_cast<uint32_t>(c[1]),
417 static_cast<uint32_t>(c[2]), neg, scale);
424 bool neg = m_bNeg ^ val.m_bNeg;
425 uint64_t a[7] = {m_uLo, m_uMid, m_uHi},
426 b[3] = {val.m_uLo, val.m_uMid, val.m_uHi}, c[7] = {0};
428 if (m_uScale < val.m_uScale) {
429 for (
int i = val.m_uScale - m_uScale; i > 0; i--)
430 decimal_helper_mul10_any(a, 7);
432 scale = m_uScale - val.m_uScale;
435 uint8_t minscale = scale;
440 decimal_helper_mul10_any(a, 7);
444 decimal_helper_div10_any(a, 7);
446 decimal_helper_raw_div(a, 6, b, 3, c, 7);
447 decimal_helper_shrinkintorange(c, 6, 3, scale);
448 decimal_helper_truncate(c[2], c[1], c[0], scale, minscale);
449 return CFGAS_Decimal(
static_cast<uint32_t>(c[0]),
static_cast<uint32_t>(c[1]),
450 static_cast<uint32_t>(c[2]), neg, scale);
#define FXMATH_DECIMAL_SCALELIMIT
#define FXMATH_DECIMAL_RSHIFT32BIT(x)
#define FXMATH_DECIMAL_LSHIFT32BIT(x)
void SetScale(uint8_t newScale)
CFGAS_Decimal(uint64_t val)
CFGAS_Decimal(float val, uint8_t scale)
CFGAS_Decimal(WideStringView str)
CFGAS_Decimal operator*(const CFGAS_Decimal &val) const
CFGAS_Decimal(uint32_t val)
WideString ToWideString() const
CFGAS_Decimal operator/(const CFGAS_Decimal &val) const
CFGAS_Decimal(int32_t val)
WideString & operator+=(wchar_t ch)
CharType operator[](const size_t index) const
bool FXSYS_IsDecimalDigit(wchar_t c)
int FXSYS_roundf(float f)