1/*
2 *  Licensed to the Apache Software Foundation (ASF) under one or more
3 *  contributor license agreements.  See the NOTICE file distributed with
4 *  this work for additional information regarding copyright ownership.
5 *  The ASF licenses this file to You under the Apache License, Version 2.0
6 *  (the "License"); you may not use this file except in compliance with
7 *  the License.  You may obtain a copy of the License at
8 *
9 *     http://www.apache.org/licenses/LICENSE-2.0
10 *
11 *  Unless required by applicable law or agreed to in writing, software
12 *  distributed under the License is distributed on an "AS IS" BASIS,
13 *  WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
14 *  See the License for the specific language governing permissions and
15 *  limitations under the License.
16 */
17
18#if !defined(cbigint_h)
19#define cbigint_h
20
21#include "JNIHelp.h"
22
23#include <sys/types.h>
24#include <sys/param.h>
25#include <stdint.h>
26
27/* IEEE floats consist of: sign bit, exponent field, significand field
28    single:  31 = sign bit, 30..23 = exponent (8 bits), 22..0 = significand (23 bits)
29    double:  63 = sign bit, 62..52 = exponent (11 bits), 51..0 = significand (52 bits)
30    inf                ==    (all exponent bits set) and (all mantissa bits clear)
31    nan                ==    (all exponent bits set) and (at least one mantissa bit set)
32    finite             ==    (at least one exponent bit clear)
33    zero               ==    (all exponent bits clear) and (all mantissa bits clear)
34    denormal           ==    (all exponent bits clear) and (at least one mantissa bit set)
35    positive           ==    sign bit clear
36    negative           ==    sign bit set
37*/
38#if __BYTE_ORDER == __LITTLE_ENDIAN
39#define DOUBLE_LO_OFFSET        0
40#define DOUBLE_HI_OFFSET        1
41#define LONG_LO_OFFSET          0
42#define LONG_HI_OFFSET          1
43#else
44#define DOUBLE_LO_OFFSET        1
45#define DOUBLE_HI_OFFSET        0
46#define LONG_LO_OFFSET          1
47#define LONG_HI_OFFSET          0
48#endif
49
50#define DOUBLE_EXPONENT_MASK_HI 0x7FF00000
51#define DOUBLE_MANTISSA_MASK_HI 0x000FFFFF
52
53union U64U32DBL {
54    uint64_t    u64val;
55    uint32_t    u32val[2];
56    int32_t    i32val[2];
57    double  dval;
58};
59
60#define DOUBLE_TO_LONGBITS(dbl) (*(reinterpret_cast<uint64_t*>(&dbl)))
61#define FLOAT_TO_INTBITS(flt) (*(reinterpret_cast<uint32_t*>(&flt)))
62#define INTBITS_TO_FLOAT(bits) (*(reinterpret_cast<float*>(&bits)))
63
64/* Replace P_FLOAT_HI and P_FLOAT_LOW */
65/* These macros are used to access the high and low 32-bit parts of a double (64-bit) value. */
66#define LOW_U32_FROM_DBL_PTR(dblptr) ((reinterpret_cast<U64U32DBL*>(dblptr))->u32val[DOUBLE_LO_OFFSET])
67#define HIGH_U32_FROM_DBL_PTR(dblptr) ((reinterpret_cast<U64U32DBL*>(dblptr))->u32val[DOUBLE_HI_OFFSET])
68#define LOW_I32_FROM_DBL_PTR(dblptr) ((reinterpret_cast<U64U32DBL*>(dblptr))->i32val[DOUBLE_LO_OFFSET])
69#define HIGH_I32_FROM_DBL_PTR(dblptr) ((reinterpret_cast<U64U32DBL*>(dblptr))->i32val[DOUBLE_HI_OFFSET])
70#define LOW_U32_FROM_DBL(dbl) LOW_U32_FROM_DBL_PTR(&(dbl))
71#define HIGH_U32_FROM_DBL(dbl) HIGH_U32_FROM_DBL_PTR(&(dbl))
72#define LOW_U32_FROM_LONG64_PTR(long64ptr) ((reinterpret_cast<U64U32DBL*>(long64ptr))->u32val[LONG_LO_OFFSET])
73#define HIGH_U32_FROM_LONG64_PTR(long64ptr) ((reinterpret_cast<U64U32DBL*>(long64ptr))->u32val[LONG_HI_OFFSET])
74#define LOW_I32_FROM_LONG64_PTR(long64ptr) ((reinterpret_cast<U64U32DBL*>(long64ptr))->i32val[LONG_LO_OFFSET])
75#define HIGH_I32_FROM_LONG64_PTR(long64ptr) ((reinterpret_cast<U64U32DBL*>(long64ptr))->i32val[LONG_HI_OFFSET])
76#define LOW_U32_FROM_LONG64(long64) LOW_U32_FROM_LONG64_PTR(&(long64))
77#define HIGH_U32_FROM_LONG64(long64) HIGH_U32_FROM_LONG64_PTR(&(long64))
78#define LOW_I32_FROM_LONG64(long64) LOW_I32_FROM_LONG64_PTR(&(long64))
79#define HIGH_I32_FROM_LONG64(long64) HIGH_I32_FROM_LONG64_PTR(&(long64))
80#define IS_DENORMAL_DBL_PTR(dblptr) (((HIGH_U32_FROM_DBL_PTR(dblptr) & DOUBLE_EXPONENT_MASK_HI) == 0) && ((HIGH_U32_FROM_DBL_PTR(dblptr) & DOUBLE_MANTISSA_MASK_HI) != 0 || (LOW_U32_FROM_DBL_PTR(dblptr) != 0)))
81#define IS_DENORMAL_DBL(dbl) IS_DENORMAL_DBL_PTR(&(dbl))
82
83#define LOW_U32_FROM_VAR(u64)     LOW_U32_FROM_LONG64(u64)
84#define LOW_U32_FROM_PTR(u64ptr)  LOW_U32_FROM_LONG64_PTR(u64ptr)
85#define HIGH_U32_FROM_VAR(u64)    HIGH_U32_FROM_LONG64(u64)
86#define HIGH_U32_FROM_PTR(u64ptr) HIGH_U32_FROM_LONG64_PTR(u64ptr)
87
88void multiplyHighPrecision(uint64_t* arg1, int32_t length1, uint64_t* arg2, int32_t length2,
89        uint64_t* result, int32_t length);
90uint32_t simpleAppendDecimalDigitHighPrecision(uint64_t* arg1, int32_t length, uint64_t digit);
91jdouble toDoubleHighPrecision(uint64_t* arg, int32_t length);
92uint64_t doubleMantissa(jdouble z);
93int32_t compareHighPrecision(uint64_t* arg1, int32_t length1, uint64_t* arg2, int32_t length2);
94int32_t highestSetBitHighPrecision(uint64_t* arg, int32_t length);
95void subtractHighPrecision(uint64_t* arg1, int32_t length1, uint64_t* arg2, int32_t length2);
96int32_t doubleExponent(jdouble z);
97int32_t addHighPrecision(uint64_t* arg1, int32_t length1, uint64_t* arg2, int32_t length2);
98int32_t lowestSetBit(uint64_t* y);
99int32_t timesTenToTheEHighPrecision(uint64_t* result, int32_t length, jint e);
100int32_t highestSetBit(uint64_t* y);
101int32_t lowestSetBitHighPrecision(uint64_t* arg, int32_t length);
102void simpleShiftLeftHighPrecision(uint64_t* arg1, int32_t length, int32_t arg2);
103uint32_t floatMantissa(jfloat z);
104int32_t simpleAddHighPrecision(uint64_t* arg1, int32_t length, uint64_t arg2);
105int32_t floatExponent(jfloat z);
106
107#endif                          /* cbigint_h */
108