Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 1 | /*- |
Lev Walkin | a5972be | 2017-09-29 23:15:58 -0700 | [diff] [blame] | 2 | * Copyright (c) 2004-2017 Lev Walkin <vlm@lionet.info>. All rights reserved. |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 3 | * Redistribution and modifications are permitted subject to BSD license. |
| 4 | */ |
Lev Walkin | e9360e6 | 2013-03-28 00:21:57 -0700 | [diff] [blame] | 5 | #define _ISOC99_SOURCE /* For ilogb() and quiet NAN */ |
Vasil Velichkov | adf6594 | 2016-08-01 16:52:02 +0300 | [diff] [blame] | 6 | #ifndef _BSD_SOURCE |
Lev Walkin | 1aea698 | 2004-10-26 09:35:25 +0000 | [diff] [blame] | 7 | #define _BSD_SOURCE /* To reintroduce finite(3) */ |
Vasil Velichkov | adf6594 | 2016-08-01 16:52:02 +0300 | [diff] [blame] | 8 | #endif |
Lev Walkin | 73823e1 | 2013-03-28 00:24:48 -0700 | [diff] [blame] | 9 | #if defined(__alpha) |
Lev Walkin | 1aea698 | 2004-10-26 09:35:25 +0000 | [diff] [blame] | 10 | #include <sys/resource.h> /* For INFINITY */ |
Lev Walkin | 3370016 | 2004-10-26 09:03:31 +0000 | [diff] [blame] | 11 | #endif |
Lev Walkin | a9cc46e | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 12 | #include <asn_internal.h> |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 13 | #include <stdlib.h> /* for strtod(3) */ |
| 14 | #include <math.h> |
| 15 | #include <errno.h> |
Lev Walkin | 3370016 | 2004-10-26 09:03:31 +0000 | [diff] [blame] | 16 | #include <REAL.h> |
Lev Walkin | 725883b | 2006-10-09 12:07:58 +0000 | [diff] [blame] | 17 | #include <OCTET_STRING.h> |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 18 | |
| 19 | #undef INT_MAX |
| 20 | #define INT_MAX ((int)(((unsigned int)-1) >> 1)) |
| 21 | |
Lev Walkin | 1aea698 | 2004-10-26 09:35:25 +0000 | [diff] [blame] | 22 | #if !(defined(NAN) || defined(INFINITY)) |
Lev Walkin | 3f99563 | 2017-09-26 18:27:32 -0700 | [diff] [blame] | 23 | static volatile double real_zero CC_NOTUSED = 0.0; |
Lev Walkin | 1aea698 | 2004-10-26 09:35:25 +0000 | [diff] [blame] | 24 | #endif |
Lev Walkin | ae1bce9 | 2004-09-27 22:18:34 +0000 | [diff] [blame] | 25 | #ifndef NAN |
Lev Walkin | c51e7d6 | 2004-09-27 22:16:18 +0000 | [diff] [blame] | 26 | #define NAN (real_zero/real_zero) |
Lev Walkin | 40319a1 | 2004-09-14 13:40:42 +0000 | [diff] [blame] | 27 | #endif |
Lev Walkin | 1aea698 | 2004-10-26 09:35:25 +0000 | [diff] [blame] | 28 | #ifndef INFINITY |
| 29 | #define INFINITY (1.0/real_zero) |
| 30 | #endif |
Lev Walkin | 40319a1 | 2004-09-14 13:40:42 +0000 | [diff] [blame] | 31 | |
Lev Walkin | eff98a5 | 2017-09-13 22:24:35 +0000 | [diff] [blame] | 32 | #if defined(__clang__) |
| 33 | /* |
| 34 | * isnan() is defined using generic selections and won't compile in |
| 35 | * strict C89 mode because of too fancy system's standard library. |
| 36 | * However, prior to C11 the math had a perfectly working isnan() |
| 37 | * in the math library. |
| 38 | * Disable generic selection warning so we can test C89 mode with newer libc. |
| 39 | */ |
| 40 | #pragma clang diagnostic push |
| 41 | #pragma clang diagnostic ignored "-Wc11-extensions" |
| 42 | static int asn_isnan(double d) { |
| 43 | return isnan(d); |
| 44 | } |
| 45 | static int asn_isfinite(double d) { |
Lev Walkin | daeb216 | 2014-01-13 23:08:35 -0800 | [diff] [blame] | 46 | #ifdef isfinite |
Lev Walkin | eff98a5 | 2017-09-13 22:24:35 +0000 | [diff] [blame] | 47 | return isfinite(d); /* ISO C99 */ |
Lev Walkin | daeb216 | 2014-01-13 23:08:35 -0800 | [diff] [blame] | 48 | #else |
Lev Walkin | eff98a5 | 2017-09-13 22:24:35 +0000 | [diff] [blame] | 49 | return finite(d); /* Deprecated on Mac OS X 10.9 */ |
Lev Walkin | daeb216 | 2014-01-13 23:08:35 -0800 | [diff] [blame] | 50 | #endif |
Lev Walkin | eff98a5 | 2017-09-13 22:24:35 +0000 | [diff] [blame] | 51 | } |
| 52 | #pragma clang diagnostic pop |
| 53 | #else /* !clang */ |
| 54 | #define asn_isnan(v) isnan(v) |
| 55 | #ifdef isfinite |
| 56 | #define asn_isfinite(d) isfinite(d) /* ISO C99 */ |
| 57 | #else |
| 58 | #define asn_isfinite(d) finite(d) /* Deprecated on Mac OS X 10.9 */ |
| 59 | #endif |
| 60 | #endif /* clang */ |
| 61 | |
Lev Walkin | daeb216 | 2014-01-13 23:08:35 -0800 | [diff] [blame] | 62 | |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 63 | /* |
| 64 | * REAL basic type description. |
| 65 | */ |
Wim Lewis | 18c2ec9 | 2014-07-29 11:30:10 -0700 | [diff] [blame] | 66 | static const ber_tlv_tag_t asn_DEF_REAL_tags[] = { |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 67 | (ASN_TAG_CLASS_UNIVERSAL | (9 << 2)) |
| 68 | }; |
Bi-Ruei, Chiu | 1f87ac0 | 2017-08-20 01:25:45 +0800 | [diff] [blame] | 69 | asn_TYPE_operation_t asn_OP_REAL = { |
| 70 | ASN__PRIMITIVE_TYPE_free, |
| 71 | REAL_print, |
| 72 | REAL_compare, |
Bi-Ruei, Chiu | 1f87ac0 | 2017-08-20 01:25:45 +0800 | [diff] [blame] | 73 | ber_decode_primitive, |
| 74 | der_encode_primitive, |
| 75 | REAL_decode_xer, |
| 76 | REAL_encode_xer, |
Lev Walkin | cc15947 | 2017-07-06 08:26:36 -0700 | [diff] [blame] | 77 | #ifdef ASN_DISABLE_OER_SUPPORT |
Bi-Ruei, Chiu | 1f87ac0 | 2017-08-20 01:25:45 +0800 | [diff] [blame] | 78 | 0, |
| 79 | 0, |
Lev Walkin | cc15947 | 2017-07-06 08:26:36 -0700 | [diff] [blame] | 80 | #else |
Bi-Ruei, Chiu | 1f87ac0 | 2017-08-20 01:25:45 +0800 | [diff] [blame] | 81 | 0, |
| 82 | 0, |
Lev Walkin | cc15947 | 2017-07-06 08:26:36 -0700 | [diff] [blame] | 83 | #endif /* ASN_DISABLE_OER_SUPPORT */ |
Lev Walkin | b33425f | 2017-07-14 14:59:52 +0400 | [diff] [blame] | 84 | #ifdef ASN_DISABLE_PER_SUPPORT |
Bi-Ruei, Chiu | 1f87ac0 | 2017-08-20 01:25:45 +0800 | [diff] [blame] | 85 | 0, |
| 86 | 0, |
Lev Walkin | b33425f | 2017-07-14 14:59:52 +0400 | [diff] [blame] | 87 | #else |
Bi-Ruei, Chiu | 1f87ac0 | 2017-08-20 01:25:45 +0800 | [diff] [blame] | 88 | REAL_decode_uper, |
| 89 | REAL_encode_uper, |
Lev Walkin | b33425f | 2017-07-14 14:59:52 +0400 | [diff] [blame] | 90 | #endif /* ASN_DISABLE_PER_SUPPORT */ |
Lev Walkin | a5972be | 2017-09-29 23:15:58 -0700 | [diff] [blame] | 91 | REAL_random_fill, |
Bi-Ruei, Chiu | 1f87ac0 | 2017-08-20 01:25:45 +0800 | [diff] [blame] | 92 | 0 /* Use generic outmost tag fetcher */ |
| 93 | }; |
| 94 | asn_TYPE_descriptor_t asn_DEF_REAL = { |
| 95 | "REAL", |
| 96 | "REAL", |
| 97 | &asn_OP_REAL, |
Bi-Ruei, Chiu | 1f87ac0 | 2017-08-20 01:25:45 +0800 | [diff] [blame] | 98 | asn_DEF_REAL_tags, |
| 99 | sizeof(asn_DEF_REAL_tags) / sizeof(asn_DEF_REAL_tags[0]), |
| 100 | asn_DEF_REAL_tags, /* Same as above */ |
| 101 | sizeof(asn_DEF_REAL_tags) / sizeof(asn_DEF_REAL_tags[0]), |
Lev Walkin | a5972be | 2017-09-29 23:15:58 -0700 | [diff] [blame] | 102 | { 0, 0, asn_generic_no_constraint }, |
Bi-Ruei, Chiu | 1f87ac0 | 2017-08-20 01:25:45 +0800 | [diff] [blame] | 103 | 0, |
| 104 | 0, /* No members */ |
| 105 | 0 /* No specifics */ |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 106 | }; |
| 107 | |
Lev Walkin | 5f56091 | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 108 | typedef enum specialRealValue { |
| 109 | SRV__NOT_A_NUMBER, |
| 110 | SRV__MINUS_INFINITY, |
| 111 | SRV__PLUS_INFINITY |
| 112 | } specialRealValue_e; |
| 113 | static struct specialRealValue_s { |
| 114 | char *string; |
Lev Walkin | 8471cec | 2004-10-21 14:02:19 +0000 | [diff] [blame] | 115 | size_t length; |
Lev Walkin | 1aea698 | 2004-10-26 09:35:25 +0000 | [diff] [blame] | 116 | long dv; |
Lev Walkin | 5f56091 | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 117 | } specialRealValue[] = { |
| 118 | #define SRV_SET(foo, val) { foo, sizeof(foo) - 1, val } |
Lev Walkin | 1aea698 | 2004-10-26 09:35:25 +0000 | [diff] [blame] | 119 | SRV_SET("<NOT-A-NUMBER/>", 0), |
| 120 | SRV_SET("<MINUS-INFINITY/>", -1), |
| 121 | SRV_SET("<PLUS-INFINITY/>", 1), |
Lev Walkin | 5f56091 | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 122 | #undef SRV_SET |
| 123 | }; |
| 124 | |
Lev Walkin | a9cc46e | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 125 | ssize_t |
| 126 | REAL__dump(double d, int canonical, asn_app_consume_bytes_f *cb, void *app_key) { |
Lev Walkin | 92b35d2 | 2004-09-27 20:52:18 +0000 | [diff] [blame] | 127 | char local_buf[64]; |
Lev Walkin | a9cc46e | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 128 | char *buf = local_buf; |
| 129 | ssize_t buflen = sizeof(local_buf); |
Lev Walkin | 92b35d2 | 2004-09-27 20:52:18 +0000 | [diff] [blame] | 130 | const char *fmt = canonical?"%.15E":"%.15f"; |
Lev Walkin | a9cc46e | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 131 | ssize_t ret; |
| 132 | |
Lev Walkin | 8e8078a | 2004-09-26 13:10:40 +0000 | [diff] [blame] | 133 | /* |
| 134 | * Check whether it is a special value. |
| 135 | */ |
Lev Walkin | c51e7d6 | 2004-09-27 22:16:18 +0000 | [diff] [blame] | 136 | /* fpclassify(3) is not portable yet */ |
Lev Walkin | eff98a5 | 2017-09-13 22:24:35 +0000 | [diff] [blame] | 137 | if(asn_isnan(d)) { |
Lev Walkin | 5f56091 | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 138 | buf = specialRealValue[SRV__NOT_A_NUMBER].string; |
| 139 | buflen = specialRealValue[SRV__NOT_A_NUMBER].length; |
Lev Walkin | c51e7d6 | 2004-09-27 22:16:18 +0000 | [diff] [blame] | 140 | return (cb(buf, buflen, app_key) < 0) ? -1 : buflen; |
Lev Walkin | eff98a5 | 2017-09-13 22:24:35 +0000 | [diff] [blame] | 141 | } else if(!asn_isfinite(d)) { |
Lev Walkin | c51e7d6 | 2004-09-27 22:16:18 +0000 | [diff] [blame] | 142 | if(copysign(1.0, d) < 0.0) { |
Lev Walkin | 5f56091 | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 143 | buf = specialRealValue[SRV__MINUS_INFINITY].string; |
| 144 | buflen = specialRealValue[SRV__MINUS_INFINITY].length; |
Lev Walkin | c51e7d6 | 2004-09-27 22:16:18 +0000 | [diff] [blame] | 145 | } else { |
Lev Walkin | 5f56091 | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 146 | buf = specialRealValue[SRV__PLUS_INFINITY].string; |
| 147 | buflen = specialRealValue[SRV__PLUS_INFINITY].length; |
Lev Walkin | c51e7d6 | 2004-09-27 22:16:18 +0000 | [diff] [blame] | 148 | } |
| 149 | return (cb(buf, buflen, app_key) < 0) ? -1 : buflen; |
| 150 | } else if(ilogb(d) <= -INT_MAX) { |
| 151 | if(copysign(1.0, d) < 0.0) { |
| 152 | buf = "-0"; |
| 153 | buflen = 2; |
| 154 | } else { |
| 155 | buf = "0"; |
| 156 | buflen = 1; |
Lev Walkin | 8e8078a | 2004-09-26 13:10:40 +0000 | [diff] [blame] | 157 | } |
| 158 | return (cb(buf, buflen, app_key) < 0) ? -1 : buflen; |
| 159 | } |
| 160 | |
| 161 | /* |
| 162 | * Use the libc's double printing, hopefully they got it right. |
| 163 | */ |
Lev Walkin | a9cc46e | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 164 | do { |
| 165 | ret = snprintf(buf, buflen, fmt, d); |
| 166 | if(ret < 0) { |
| 167 | buflen <<= 1; |
| 168 | } else if(ret >= buflen) { |
| 169 | buflen = ret + 1; |
| 170 | } else { |
| 171 | buflen = ret; |
| 172 | break; |
| 173 | } |
Lev Walkin | 419f675 | 2006-09-13 04:02:00 +0000 | [diff] [blame] | 174 | if(buf != local_buf) FREEMEM(buf); |
Lev Walkin | 8e8078a | 2004-09-26 13:10:40 +0000 | [diff] [blame] | 175 | buf = (char *)MALLOC(buflen); |
Lev Walkin | a9cc46e | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 176 | if(!buf) return -1; |
| 177 | } while(1); |
| 178 | |
Lev Walkin | a9cc46e | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 179 | if(canonical) { |
Lev Walkin | 92b35d2 | 2004-09-27 20:52:18 +0000 | [diff] [blame] | 180 | /* |
| 181 | * Transform the "[-]d.dddE+-dd" output into "[-]d.dddE[-]d" |
Lev Walkin | 5f56091 | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 182 | * Check that snprintf() constructed the output correctly. |
Lev Walkin | 92b35d2 | 2004-09-27 20:52:18 +0000 | [diff] [blame] | 183 | */ |
Lev Walkin | 8ca13c8 | 2016-01-24 22:40:00 -0800 | [diff] [blame] | 184 | char *dot; |
Lev Walkin | a9cc46e | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 185 | char *end = buf + buflen; |
Lev Walkin | 92b35d2 | 2004-09-27 20:52:18 +0000 | [diff] [blame] | 186 | char *last_zero; |
Lev Walkin | 8ca13c8 | 2016-01-24 22:40:00 -0800 | [diff] [blame] | 187 | char *prev_zero; |
| 188 | char *s; |
| 189 | |
| 190 | enum { |
| 191 | LZSTATE_NOTHING, |
| 192 | LZSTATE_SEEN_ZERO |
| 193 | } lz_state = LZSTATE_NOTHING; |
Lev Walkin | a9cc46e | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 194 | |
Lev Walkin | 5f56091 | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 195 | dot = (buf[0] == 0x2d /* '-' */) ? (buf + 2) : (buf + 1); |
Lev Walkin | a9cc46e | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 196 | if(*dot >= 0x30) { |
Lev Walkin | 78d917d | 2012-01-07 15:29:08 -0800 | [diff] [blame] | 197 | if(buf != local_buf) FREEMEM(buf); |
Lev Walkin | a9cc46e | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 198 | errno = EINVAL; |
| 199 | return -1; /* Not a dot, really */ |
| 200 | } |
Lev Walkin | 92b35d2 | 2004-09-27 20:52:18 +0000 | [diff] [blame] | 201 | *dot = 0x2e; /* Replace possible comma */ |
Lev Walkin | a9cc46e | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 202 | |
Lev Walkin | 8ca13c8 | 2016-01-24 22:40:00 -0800 | [diff] [blame] | 203 | for(prev_zero = last_zero = s = dot + 2; s < end; s++) { |
| 204 | switch(*s) { |
| 205 | case 0x45: /* 'E' */ |
| 206 | if(lz_state == LZSTATE_SEEN_ZERO) |
| 207 | last_zero = prev_zero; |
| 208 | break; |
| 209 | case 0x30: /* '0' */ |
| 210 | if(lz_state == LZSTATE_NOTHING) |
| 211 | prev_zero = s; |
| 212 | lz_state = LZSTATE_SEEN_ZERO; |
| 213 | continue; |
| 214 | default: |
| 215 | lz_state = LZSTATE_NOTHING; |
| 216 | continue; |
| 217 | } |
| 218 | break; |
| 219 | } |
| 220 | |
| 221 | if(s == end) { |
Lev Walkin | 78d917d | 2012-01-07 15:29:08 -0800 | [diff] [blame] | 222 | if(buf != local_buf) FREEMEM(buf); |
Lev Walkin | a9cc46e | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 223 | errno = EINVAL; |
| 224 | return -1; /* No promised E */ |
| 225 | } |
Lev Walkin | 8ca13c8 | 2016-01-24 22:40:00 -0800 | [diff] [blame] | 226 | |
| 227 | assert(*s == 0x45); |
| 228 | { |
| 229 | char *E = s; |
| 230 | char *expptr = ++E; |
| 231 | char *s = expptr; |
| 232 | int sign; |
| 233 | |
| 234 | if(*expptr == 0x2b /* '+' */) { |
| 235 | /* Skip the "+" */ |
| 236 | buflen -= 1; |
| 237 | sign = 0; |
| 238 | } else { |
| 239 | sign = 1; |
| 240 | s++; |
| 241 | } |
| 242 | expptr++; |
| 243 | if(expptr > end) { |
| 244 | if(buf != local_buf) FREEMEM(buf); |
| 245 | errno = EINVAL; |
| 246 | return -1; |
| 247 | } |
| 248 | if(*expptr == 0x30) { |
| 249 | buflen--; |
| 250 | expptr++; |
| 251 | } |
| 252 | if(*last_zero == 0x30) { |
| 253 | *last_zero = 0x45; /* E */ |
| 254 | buflen -= s - (last_zero + 1); |
| 255 | s = last_zero + 1; |
| 256 | if(sign) { |
| 257 | *s++ = 0x2d /* '-' */; |
| 258 | buflen++; |
| 259 | } |
| 260 | } |
| 261 | for(; expptr <= end; s++, expptr++) |
| 262 | *s = *expptr; |
| 263 | } |
Lev Walkin | 92b35d2 | 2004-09-27 20:52:18 +0000 | [diff] [blame] | 264 | } else { |
| 265 | /* |
| 266 | * Remove trailing zeros. |
| 267 | */ |
| 268 | char *end = buf + buflen; |
| 269 | char *last_zero = end; |
Lev Walkin | f0b808d | 2005-04-25 21:08:25 +0000 | [diff] [blame] | 270 | int stoplooking = 0; |
Lev Walkin | 92b35d2 | 2004-09-27 20:52:18 +0000 | [diff] [blame] | 271 | char *z; |
| 272 | for(z = end - 1; z > buf; z--) { |
| 273 | switch(*z) { |
Lev Walkin | f0b808d | 2005-04-25 21:08:25 +0000 | [diff] [blame] | 274 | case 0x30: |
| 275 | if(!stoplooking) |
| 276 | last_zero = z; |
| 277 | continue; |
Lev Walkin | 92b35d2 | 2004-09-27 20:52:18 +0000 | [diff] [blame] | 278 | case 0x31: case 0x32: case 0x33: case 0x34: |
| 279 | case 0x35: case 0x36: case 0x37: case 0x38: case 0x39: |
Lev Walkin | f0b808d | 2005-04-25 21:08:25 +0000 | [diff] [blame] | 280 | stoplooking = 1; |
Lev Walkin | 92b35d2 | 2004-09-27 20:52:18 +0000 | [diff] [blame] | 281 | continue; |
| 282 | default: /* Catch dot and other separators */ |
Lev Walkin | f0b808d | 2005-04-25 21:08:25 +0000 | [diff] [blame] | 283 | /* |
| 284 | * Replace possible comma (which may even |
| 285 | * be not a comma at all: locale-defined). |
| 286 | */ |
| 287 | *z = 0x2e; |
Lev Walkin | 92b35d2 | 2004-09-27 20:52:18 +0000 | [diff] [blame] | 288 | if(last_zero == z + 1) { /* leave x.0 */ |
| 289 | last_zero++; |
| 290 | } |
| 291 | buflen = last_zero - buf; |
| 292 | *last_zero = '\0'; |
| 293 | break; |
| 294 | } |
| 295 | break; |
| 296 | } |
Lev Walkin | a9cc46e | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 297 | } |
| 298 | |
| 299 | ret = cb(buf, buflen, app_key); |
Lev Walkin | 419f675 | 2006-09-13 04:02:00 +0000 | [diff] [blame] | 300 | if(buf != local_buf) FREEMEM(buf); |
Lev Walkin | a9cc46e | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 301 | return (ret < 0) ? -1 : buflen; |
| 302 | } |
| 303 | |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 304 | int |
Lev Walkin | 5e03376 | 2004-09-29 13:26:15 +0000 | [diff] [blame] | 305 | REAL_print(asn_TYPE_descriptor_t *td, const void *sptr, int ilevel, |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 306 | asn_app_consume_bytes_f *cb, void *app_key) { |
| 307 | const REAL_t *st = (const REAL_t *)sptr; |
Lev Walkin | 8e8078a | 2004-09-26 13:10:40 +0000 | [diff] [blame] | 308 | ssize_t ret; |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 309 | double d; |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 310 | |
| 311 | (void)td; /* Unused argument */ |
| 312 | (void)ilevel; /* Unused argument */ |
| 313 | |
Lev Walkin | a9cc46e | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 314 | if(!st || !st->buf) |
Lev Walkin | 8e8078a | 2004-09-26 13:10:40 +0000 | [diff] [blame] | 315 | ret = cb("<absent>", 8, app_key); |
Lev Walkin | 5e03376 | 2004-09-29 13:26:15 +0000 | [diff] [blame] | 316 | else if(asn_REAL2double(st, &d)) |
Lev Walkin | 8e8078a | 2004-09-26 13:10:40 +0000 | [diff] [blame] | 317 | ret = cb("<error>", 7, app_key); |
| 318 | else |
| 319 | ret = REAL__dump(d, 0, cb, app_key); |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 320 | |
Lev Walkin | 8e8078a | 2004-09-26 13:10:40 +0000 | [diff] [blame] | 321 | return (ret < 0) ? -1 : 0; |
Lev Walkin | a9cc46e | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 322 | } |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 323 | |
Lev Walkin | cd2f48e | 2017-08-10 02:14:59 -0700 | [diff] [blame] | 324 | int |
| 325 | REAL_compare(const asn_TYPE_descriptor_t *td, const void *aptr, |
| 326 | const void *bptr) { |
| 327 | const REAL_t *a = aptr; |
| 328 | const REAL_t *b = bptr; |
| 329 | |
| 330 | (void)td; |
| 331 | |
| 332 | if(a && b) { |
| 333 | double adbl, bdbl; |
| 334 | int ra, rb; |
| 335 | ra = asn_REAL2double(a, &adbl); |
| 336 | rb = asn_REAL2double(b, &bdbl); |
| 337 | if(ra == 0 && rb == 0) { |
Lev Walkin | eff98a5 | 2017-09-13 22:24:35 +0000 | [diff] [blame] | 338 | if(asn_isnan(adbl)) { |
| 339 | if(asn_isnan(bdbl)) { |
Lev Walkin | cd2f48e | 2017-08-10 02:14:59 -0700 | [diff] [blame] | 340 | return 0; |
| 341 | } else { |
| 342 | return -1; |
| 343 | } |
Lev Walkin | eff98a5 | 2017-09-13 22:24:35 +0000 | [diff] [blame] | 344 | } else if(asn_isnan(bdbl)) { |
Lev Walkin | cd2f48e | 2017-08-10 02:14:59 -0700 | [diff] [blame] | 345 | return 1; |
| 346 | } |
| 347 | /* Value comparison. */ |
| 348 | if(adbl < bdbl) { |
| 349 | return -1; |
| 350 | } else if(adbl > bdbl) { |
| 351 | return 1; |
| 352 | } else { |
| 353 | return 0; |
| 354 | } |
| 355 | } else if(ra) { |
| 356 | return -1; |
| 357 | } else { |
| 358 | return 1; |
| 359 | } |
| 360 | } else if(!a) { |
| 361 | return -1; |
| 362 | } else { |
| 363 | return 1; |
| 364 | } |
| 365 | } |
| 366 | |
Lev Walkin | a9cc46e | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 367 | asn_enc_rval_t |
Lev Walkin | 5e03376 | 2004-09-29 13:26:15 +0000 | [diff] [blame] | 368 | REAL_encode_xer(asn_TYPE_descriptor_t *td, void *sptr, |
Lev Walkin | a9cc46e | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 369 | int ilevel, enum xer_encoder_flags_e flags, |
| 370 | asn_app_consume_bytes_f *cb, void *app_key) { |
| 371 | REAL_t *st = (REAL_t *)sptr; |
| 372 | asn_enc_rval_t er; |
| 373 | double d; |
| 374 | |
| 375 | (void)ilevel; |
| 376 | |
Lev Walkin | 5e03376 | 2004-09-29 13:26:15 +0000 | [diff] [blame] | 377 | if(!st || !st->buf || asn_REAL2double(st, &d)) |
Lev Walkin | 7c1dc05 | 2016-03-14 03:08:15 -0700 | [diff] [blame] | 378 | ASN__ENCODE_FAILED; |
Lev Walkin | a9cc46e | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 379 | |
| 380 | er.encoded = REAL__dump(d, flags & XER_F_CANONICAL, cb, app_key); |
Lev Walkin | 7c1dc05 | 2016-03-14 03:08:15 -0700 | [diff] [blame] | 381 | if(er.encoded < 0) ASN__ENCODE_FAILED; |
Lev Walkin | a9cc46e | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 382 | |
Lev Walkin | 7c1dc05 | 2016-03-14 03:08:15 -0700 | [diff] [blame] | 383 | ASN__ENCODED_OK(er); |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 384 | } |
| 385 | |
Lev Walkin | 5f56091 | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 386 | |
| 387 | /* |
| 388 | * Decode the chunk of XML text encoding REAL. |
| 389 | */ |
Lev Walkin | 0fab1a6 | 2005-03-09 22:19:25 +0000 | [diff] [blame] | 390 | static enum xer_pbd_rval |
| 391 | REAL__xer_body_decode(asn_TYPE_descriptor_t *td, void *sptr, const void *chunk_buf, size_t chunk_size) { |
Lev Walkin | 8471cec | 2004-10-21 14:02:19 +0000 | [diff] [blame] | 392 | REAL_t *st = (REAL_t *)sptr; |
Lev Walkin | 5f56091 | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 393 | double value; |
Lev Walkin | 0fab1a6 | 2005-03-09 22:19:25 +0000 | [diff] [blame] | 394 | const char *xerdata = (const char *)chunk_buf; |
Lev Walkin | 5f56091 | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 395 | char *endptr = 0; |
| 396 | char *b; |
| 397 | |
Lev Walkin | e0b56e0 | 2005-02-25 12:10:27 +0000 | [diff] [blame] | 398 | (void)td; |
| 399 | |
Lev Walkin | 0fab1a6 | 2005-03-09 22:19:25 +0000 | [diff] [blame] | 400 | if(!chunk_size) return XPBD_BROKEN_ENCODING; |
Lev Walkin | 5f56091 | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 401 | |
| 402 | /* |
| 403 | * Decode an XMLSpecialRealValue: <MINUS-INFINITY>, etc. |
| 404 | */ |
| 405 | if(xerdata[0] == 0x3c /* '<' */) { |
| 406 | size_t i; |
| 407 | for(i = 0; i < sizeof(specialRealValue) |
| 408 | / sizeof(specialRealValue[0]); i++) { |
| 409 | struct specialRealValue_s *srv = &specialRealValue[i]; |
Lev Walkin | 1aea698 | 2004-10-26 09:35:25 +0000 | [diff] [blame] | 410 | double dv; |
| 411 | |
Lev Walkin | 5f56091 | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 412 | if(srv->length != chunk_size |
| 413 | || memcmp(srv->string, chunk_buf, chunk_size)) |
| 414 | continue; |
| 415 | |
Lev Walkin | 1aea698 | 2004-10-26 09:35:25 +0000 | [diff] [blame] | 416 | /* |
| 417 | * It could've been done using |
| 418 | * (double)srv->dv / real_zero, |
| 419 | * but it summons fp exception on some platforms. |
| 420 | */ |
| 421 | switch(srv->dv) { |
| 422 | case -1: dv = - INFINITY; break; |
| 423 | case 0: dv = NAN; break; |
| 424 | case 1: dv = INFINITY; break; |
Lev Walkin | 0fab1a6 | 2005-03-09 22:19:25 +0000 | [diff] [blame] | 425 | default: return XPBD_SYSTEM_FAILURE; |
Lev Walkin | 1aea698 | 2004-10-26 09:35:25 +0000 | [diff] [blame] | 426 | } |
| 427 | |
Lev Walkin | 0fab1a6 | 2005-03-09 22:19:25 +0000 | [diff] [blame] | 428 | if(asn_double2REAL(st, dv)) |
| 429 | return XPBD_SYSTEM_FAILURE; |
Lev Walkin | 5f56091 | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 430 | |
Lev Walkin | 0fab1a6 | 2005-03-09 22:19:25 +0000 | [diff] [blame] | 431 | return XPBD_BODY_CONSUMED; |
Lev Walkin | 5f56091 | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 432 | } |
| 433 | ASN_DEBUG("Unknown XMLSpecialRealValue"); |
Lev Walkin | 0fab1a6 | 2005-03-09 22:19:25 +0000 | [diff] [blame] | 434 | return XPBD_BROKEN_ENCODING; |
Lev Walkin | 5f56091 | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 435 | } |
| 436 | |
| 437 | /* |
| 438 | * Copy chunk into the nul-terminated string, and run strtod. |
| 439 | */ |
Lev Walkin | 8471cec | 2004-10-21 14:02:19 +0000 | [diff] [blame] | 440 | b = (char *)MALLOC(chunk_size + 1); |
Lev Walkin | 0fab1a6 | 2005-03-09 22:19:25 +0000 | [diff] [blame] | 441 | if(!b) return XPBD_SYSTEM_FAILURE; |
Lev Walkin | 5f56091 | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 442 | memcpy(b, chunk_buf, chunk_size); |
Lev Walkin | 8471cec | 2004-10-21 14:02:19 +0000 | [diff] [blame] | 443 | b[chunk_size] = 0; /* nul-terminate */ |
Lev Walkin | 5f56091 | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 444 | |
| 445 | value = strtod(b, &endptr); |
Lev Walkin | 419f675 | 2006-09-13 04:02:00 +0000 | [diff] [blame] | 446 | FREEMEM(b); |
Lev Walkin | 0fab1a6 | 2005-03-09 22:19:25 +0000 | [diff] [blame] | 447 | if(endptr == b) return XPBD_BROKEN_ENCODING; |
Lev Walkin | 5f56091 | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 448 | |
| 449 | if(asn_double2REAL(st, value)) |
Lev Walkin | 0fab1a6 | 2005-03-09 22:19:25 +0000 | [diff] [blame] | 450 | return XPBD_SYSTEM_FAILURE; |
Lev Walkin | 5f56091 | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 451 | |
Lev Walkin | 0fab1a6 | 2005-03-09 22:19:25 +0000 | [diff] [blame] | 452 | return XPBD_BODY_CONSUMED; |
Lev Walkin | 5f56091 | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 453 | } |
| 454 | |
| 455 | asn_dec_rval_t |
Lev Walkin | afbf2a9 | 2017-09-12 23:30:27 -0700 | [diff] [blame] | 456 | REAL_decode_xer(const asn_codec_ctx_t *opt_codec_ctx, |
Lev Walkin | 5f56091 | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 457 | asn_TYPE_descriptor_t *td, void **sptr, const char *opt_mname, |
Lev Walkin | 8c3b854 | 2005-03-10 18:52:02 +0000 | [diff] [blame] | 458 | const void *buf_ptr, size_t size) { |
Lev Walkin | 5f56091 | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 459 | |
| 460 | return xer_decode_primitive(opt_codec_ctx, td, |
Lev Walkin | 8471cec | 2004-10-21 14:02:19 +0000 | [diff] [blame] | 461 | sptr, sizeof(REAL_t), opt_mname, |
Lev Walkin | 5f56091 | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 462 | buf_ptr, size, REAL__xer_body_decode); |
| 463 | } |
| 464 | |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 465 | int |
Lev Walkin | 5e03376 | 2004-09-29 13:26:15 +0000 | [diff] [blame] | 466 | asn_REAL2double(const REAL_t *st, double *dbl_value) { |
Lev Walkin | a9cc46e | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 467 | unsigned int octv; |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 468 | |
| 469 | if(!st || !st->buf) { |
| 470 | errno = EINVAL; |
| 471 | return -1; |
| 472 | } |
| 473 | |
| 474 | if(st->size == 0) { |
| 475 | *dbl_value = 0; |
| 476 | return 0; |
| 477 | } |
| 478 | |
| 479 | octv = st->buf[0]; /* unsigned byte */ |
| 480 | |
| 481 | switch(octv & 0xC0) { |
Lev Walkin | 5fda7d5 | 2012-01-09 03:20:19 -0800 | [diff] [blame] | 482 | case 0x40: /* X.690: 8.5.6 a) => 8.5.9 */ |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 483 | /* "SpecialRealValue" */ |
| 484 | |
| 485 | /* Be liberal in what you accept... |
Lev Walkin | 749916f | 2012-01-07 17:03:39 -0800 | [diff] [blame] | 486 | * http://en.wikipedia.org/wiki/Robustness_principle |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 487 | if(st->size != 1) ... |
| 488 | */ |
| 489 | |
| 490 | switch(st->buf[0]) { |
| 491 | case 0x40: /* 01000000: PLUS-INFINITY */ |
Lev Walkin | 1aea698 | 2004-10-26 09:35:25 +0000 | [diff] [blame] | 492 | *dbl_value = INFINITY; |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 493 | return 0; |
| 494 | case 0x41: /* 01000001: MINUS-INFINITY */ |
Lev Walkin | 1aea698 | 2004-10-26 09:35:25 +0000 | [diff] [blame] | 495 | *dbl_value = - INFINITY; |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 496 | return 0; |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 497 | case 0x42: /* 01000010: NOT-A-NUMBER */ |
| 498 | *dbl_value = NAN; |
| 499 | return 0; |
| 500 | case 0x43: /* 01000011: minus zero */ |
Lev Walkin | 2a789d9 | 2004-09-27 21:36:59 +0000 | [diff] [blame] | 501 | *dbl_value = -0.0; |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 502 | return 0; |
| 503 | } |
| 504 | |
| 505 | errno = EINVAL; |
| 506 | return -1; |
Lev Walkin | 0959ffb | 2011-06-25 15:48:52 -0700 | [diff] [blame] | 507 | case 0x00: { /* X.690: 8.5.7 */ |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 508 | /* |
Lev Walkin | 5fda7d5 | 2012-01-09 03:20:19 -0800 | [diff] [blame] | 509 | * Decimal. NR{1,2,3} format from ISO 6093. |
| 510 | * NR1: [ ]*[+-]?[0-9]+ |
| 511 | * NR2: [ ]*[+-]?([0-9]+\.[0-9]*|[0-9]*\.[0-9]+) |
| 512 | * NR3: [ ]*[+-]?([0-9]+\.[0-9]*|[0-9]*\.[0-9]+)[Ee][+-]?[0-9]+ |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 513 | */ |
| 514 | double d; |
Lev Walkin | 5fda7d5 | 2012-01-09 03:20:19 -0800 | [diff] [blame] | 515 | char *buf; |
| 516 | char *endptr; |
| 517 | int used_malloc = 0; |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 518 | |
Lev Walkin | 5fda7d5 | 2012-01-09 03:20:19 -0800 | [diff] [blame] | 519 | if(octv == 0 || (octv & 0x3C)) { |
Lev Walkin | 0959ffb | 2011-06-25 15:48:52 -0700 | [diff] [blame] | 520 | /* Remaining values of bits 6 to 1 are Reserved. */ |
| 521 | errno = EINVAL; |
| 522 | return -1; |
| 523 | } |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 524 | |
Lev Walkin | 5fda7d5 | 2012-01-09 03:20:19 -0800 | [diff] [blame] | 525 | |
Lev Walkin | 53e5ae6 | 2017-08-23 10:56:19 -0700 | [diff] [blame] | 526 | /* 1. By contract, an input buffer should be '\0'-terminated. |
Lev Walkin | 5fda7d5 | 2012-01-09 03:20:19 -0800 | [diff] [blame] | 527 | * OCTET STRING decoder ensures that, as is asn_double2REAL(). |
| 528 | * 2. ISO 6093 specifies COMMA as a possible decimal separator. |
| 529 | * However, strtod() can't always deal with COMMA. |
| 530 | * So her we fix both by reallocating, copying and fixing. |
| 531 | */ |
Lev Walkin | 53e5ae6 | 2017-08-23 10:56:19 -0700 | [diff] [blame] | 532 | if(st->buf[st->size] != '\0' || memchr(st->buf, ',', st->size)) { |
Lev Walkin | 5fda7d5 | 2012-01-09 03:20:19 -0800 | [diff] [blame] | 533 | uint8_t *p, *end; |
| 534 | char *b; |
| 535 | if(st->size > 100) { |
| 536 | /* Avoid malicious stack overflow in alloca() */ |
| 537 | buf = (char *)MALLOC(st->size); |
| 538 | if(!buf) return -1; |
| 539 | used_malloc = 1; |
| 540 | } else { |
| 541 | buf = alloca(st->size); |
| 542 | } |
| 543 | b = buf; |
| 544 | /* Copy without the first byte and with 0-termination */ |
| 545 | for(p = st->buf + 1, end = st->buf + st->size; |
| 546 | p < end; b++, p++) |
| 547 | *b = (*p == ',') ? '.' : *p; |
| 548 | *b = '\0'; |
| 549 | } else { |
| 550 | buf = (char *)&st->buf[1]; |
| 551 | } |
| 552 | |
| 553 | endptr = buf; |
| 554 | d = strtod(buf, &endptr); |
| 555 | if(*endptr != '\0') { |
| 556 | /* Format is not consistent with ISO 6093 */ |
| 557 | if(used_malloc) FREEMEM(buf); |
Lev Walkin | 0959ffb | 2011-06-25 15:48:52 -0700 | [diff] [blame] | 558 | errno = EINVAL; |
| 559 | return -1; |
| 560 | } |
Lev Walkin | 5fda7d5 | 2012-01-09 03:20:19 -0800 | [diff] [blame] | 561 | if(used_malloc) FREEMEM(buf); |
Lev Walkin | eff98a5 | 2017-09-13 22:24:35 +0000 | [diff] [blame] | 562 | if(asn_isfinite(d)) { |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 563 | *dbl_value = d; |
| 564 | return 0; |
| 565 | } else { |
| 566 | errno = ERANGE; |
Lev Walkin | 0959ffb | 2011-06-25 15:48:52 -0700 | [diff] [blame] | 567 | return -1; |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 568 | } |
| 569 | } |
| 570 | } |
| 571 | |
| 572 | /* |
| 573 | * Binary representation. |
| 574 | */ |
| 575 | { |
| 576 | double m; |
| 577 | int expval; /* exponent value */ |
| 578 | unsigned int elen; /* exponent value length, in octets */ |
Lev Walkin | 9c57f33 | 2017-09-17 22:30:49 -0700 | [diff] [blame] | 579 | int scaleF; |
| 580 | int baseF; |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 581 | uint8_t *ptr; |
| 582 | uint8_t *end; |
| 583 | int sign; |
| 584 | |
| 585 | switch((octv & 0x30) >> 4) { |
| 586 | case 0x00: baseF = 1; break; /* base 2 */ |
| 587 | case 0x01: baseF = 3; break; /* base 8 */ |
| 588 | case 0x02: baseF = 4; break; /* base 16 */ |
| 589 | default: |
| 590 | /* Reserved field, can't parse now. */ |
| 591 | errno = EINVAL; |
| 592 | return -1; |
| 593 | } |
| 594 | |
| 595 | sign = (octv & 0x40); /* bit 7 */ |
| 596 | scaleF = (octv & 0x0C) >> 2; /* bits 4 to 3 */ |
| 597 | |
Lev Walkin | 494fb70 | 2017-08-07 20:07:00 -0700 | [diff] [blame] | 598 | if(st->size <= 1 + (octv & 0x03)) { |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 599 | errno = EINVAL; |
| 600 | return -1; |
| 601 | } |
| 602 | |
Lev Walkin | e78753d | 2007-06-23 17:02:00 +0000 | [diff] [blame] | 603 | elen = (octv & 0x03); /* bits 2 to 1; 8.5.6.4 */ |
| 604 | if(elen == 0x03) { /* bits 2 to 1 = 11; 8.5.6.4, case d) */ |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 605 | elen = st->buf[1]; /* unsigned binary number */ |
Lev Walkin | 494fb70 | 2017-08-07 20:07:00 -0700 | [diff] [blame] | 606 | if(elen == 0 || st->size <= (2 + elen)) { |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 607 | errno = EINVAL; |
| 608 | return -1; |
| 609 | } |
Lev Walkin | e78753d | 2007-06-23 17:02:00 +0000 | [diff] [blame] | 610 | /* FIXME: verify constraints of case d) */ |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 611 | ptr = &st->buf[2]; |
| 612 | } else { |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 613 | ptr = &st->buf[1]; |
| 614 | } |
| 615 | |
| 616 | /* Fetch the multibyte exponent */ |
| 617 | expval = (int)(*(int8_t *)ptr); |
| 618 | end = ptr + elen + 1; |
| 619 | for(ptr++; ptr < end; ptr++) |
| 620 | expval = (expval * 256) + *ptr; |
| 621 | |
| 622 | m = 0.0; /* Initial mantissa value */ |
| 623 | |
| 624 | /* Okay, the exponent is here. Now, what about mantissa? */ |
| 625 | end = st->buf + st->size; |
Lev Walkin | 5fda7d5 | 2012-01-09 03:20:19 -0800 | [diff] [blame] | 626 | for(; ptr < end; ptr++) |
| 627 | m = ldexp(m, 8) + *ptr; |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 628 | |
Lev Walkin | 5f56091 | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 629 | if(0) |
Lev Walkin | 5fda7d5 | 2012-01-09 03:20:19 -0800 | [diff] [blame] | 630 | ASN_DEBUG("m=%.10f, scF=%d, bF=%d, expval=%d, ldexp()=%f, ldexp()=%f\n", |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 631 | m, scaleF, baseF, expval, |
| 632 | ldexp(m, expval * baseF + scaleF), |
Lev Walkin | 77aa90f | 2005-04-28 22:56:36 +0000 | [diff] [blame] | 633 | ldexp(m, scaleF) * pow(pow(2, baseF), expval) |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 634 | ); |
| 635 | |
| 636 | /* |
| 637 | * (S * N * 2^F) * B^E |
| 638 | * Essentially: |
Lev Walkin | 77aa90f | 2005-04-28 22:56:36 +0000 | [diff] [blame] | 639 | m = ldexp(m, scaleF) * pow(pow(2, base), expval); |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 640 | */ |
| 641 | m = ldexp(m, expval * baseF + scaleF); |
Lev Walkin | eff98a5 | 2017-09-13 22:24:35 +0000 | [diff] [blame] | 642 | if(asn_isfinite(m)) { |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 643 | *dbl_value = sign ? -m : m; |
| 644 | } else { |
| 645 | errno = ERANGE; |
| 646 | return -1; |
| 647 | } |
| 648 | |
| 649 | } /* if(binary_format) */ |
| 650 | |
| 651 | return 0; |
| 652 | } |
| 653 | |
| 654 | /* |
| 655 | * Assume IEEE 754 floating point: standard 64 bit double. |
| 656 | * [1 bit sign] [11 bits exponent] [52 bits mantissa] |
| 657 | */ |
| 658 | int |
Lev Walkin | 5e03376 | 2004-09-29 13:26:15 +0000 | [diff] [blame] | 659 | asn_double2REAL(REAL_t *st, double dbl_value) { |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 660 | #ifdef WORDS_BIGENDIAN /* Known to be big-endian */ |
| 661 | int littleEndian = 0; |
| 662 | #else /* need to test: have no explicit information */ |
| 663 | unsigned int LE = 1; |
| 664 | int littleEndian = *(unsigned char *)&LE; |
| 665 | #endif |
| 666 | uint8_t buf[16]; /* More than enough for 8-byte dbl_value */ |
| 667 | uint8_t dscr[sizeof(dbl_value)]; /* double value scratch pad */ |
| 668 | /* Assertion guards: won't even compile, if unexpected double size */ |
Lev Walkin | 3f99563 | 2017-09-26 18:27:32 -0700 | [diff] [blame] | 669 | char assertion_buffer1[9 - sizeof(dbl_value)] CC_NOTUSED; |
| 670 | char assertion_buffer2[sizeof(dbl_value) - 7] CC_NOTUSED; |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 671 | uint8_t *ptr = buf; |
| 672 | uint8_t *mstop; /* Last byte of mantissa */ |
| 673 | unsigned int mval; /* Value of the last byte of mantissa */ |
| 674 | unsigned int bmsign; /* binary mask with sign */ |
| 675 | unsigned int buflen; |
| 676 | unsigned int accum; |
| 677 | int expval; |
| 678 | |
| 679 | if(!st) { |
| 680 | errno = EINVAL; |
| 681 | return -1; |
| 682 | } |
| 683 | |
Lev Walkin | 057fb73 | 2004-09-14 13:58:10 +0000 | [diff] [blame] | 684 | /* |
| 685 | * ilogb(+-0) returns -INT_MAX or INT_MIN (platform-dependent) |
Lev Walkin | 2f50502 | 2005-07-01 08:28:18 +0000 | [diff] [blame] | 686 | * ilogb(+-inf) returns INT_MAX, logb(+-inf) returns +inf |
Lev Walkin | 2a789d9 | 2004-09-27 21:36:59 +0000 | [diff] [blame] | 687 | * ilogb(NaN) returns INT_MIN or INT_MAX (platform-dependent) |
Lev Walkin | 057fb73 | 2004-09-14 13:58:10 +0000 | [diff] [blame] | 688 | */ |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 689 | expval = ilogb(dbl_value); |
Lev Walkin | 2a789d9 | 2004-09-27 21:36:59 +0000 | [diff] [blame] | 690 | if(expval <= -INT_MAX /* Also catches +-0 and maybe isnan() */ |
| 691 | || expval == INT_MAX /* catches isfin() and maybe isnan() */ |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 692 | ) { |
| 693 | if(!st->buf || st->size < 2) { |
Lev Walkin | 8e8078a | 2004-09-26 13:10:40 +0000 | [diff] [blame] | 694 | ptr = (uint8_t *)MALLOC(2); |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 695 | if(!ptr) return -1; |
Lev Walkin | dcae9ce | 2017-09-18 20:13:36 -0700 | [diff] [blame] | 696 | if(st->buf) FREEMEM(st->buf); |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 697 | st->buf = ptr; |
| 698 | } |
| 699 | /* fpclassify(3) is not portable yet */ |
Lev Walkin | eff98a5 | 2017-09-13 22:24:35 +0000 | [diff] [blame] | 700 | if(asn_isnan(dbl_value)) { |
Lev Walkin | 2a789d9 | 2004-09-27 21:36:59 +0000 | [diff] [blame] | 701 | st->buf[0] = 0x42; /* NaN */ |
| 702 | st->buf[1] = 0; |
| 703 | st->size = 1; |
Lev Walkin | eff98a5 | 2017-09-13 22:24:35 +0000 | [diff] [blame] | 704 | } else if(!asn_isfinite(dbl_value)) { |
Lev Walkin | 7e03db9 | 2004-09-14 13:50:21 +0000 | [diff] [blame] | 705 | if(copysign(1.0, dbl_value) < 0.0) { |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 706 | st->buf[0] = 0x41; /* MINUS-INFINITY */ |
| 707 | } else { |
| 708 | st->buf[0] = 0x40; /* PLUS-INFINITY */ |
| 709 | } |
| 710 | st->buf[1] = 0; |
| 711 | st->size = 1; |
Lev Walkin | c51e7d6 | 2004-09-27 22:16:18 +0000 | [diff] [blame] | 712 | } else { |
Lev Walkin | 5fe72e9 | 2012-01-07 17:00:29 -0800 | [diff] [blame] | 713 | if(copysign(1.0, dbl_value) >= 0.0) { |
Lev Walkin | c51e7d6 | 2004-09-27 22:16:18 +0000 | [diff] [blame] | 714 | /* no content octets: positive zero */ |
| 715 | st->buf[0] = 0; /* JIC */ |
| 716 | st->size = 0; |
Lev Walkin | 5fe72e9 | 2012-01-07 17:00:29 -0800 | [diff] [blame] | 717 | } else { |
| 718 | /* Negative zero. #8.5.3, 8.5.9 */ |
| 719 | st->buf[0] = 0x43; |
Lev Walkin | dcae9ce | 2017-09-18 20:13:36 -0700 | [diff] [blame] | 720 | st->buf[1] = 0; |
Lev Walkin | 5fe72e9 | 2012-01-07 17:00:29 -0800 | [diff] [blame] | 721 | st->size = 1; |
Lev Walkin | c51e7d6 | 2004-09-27 22:16:18 +0000 | [diff] [blame] | 722 | } |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 723 | } |
| 724 | return 0; |
| 725 | } |
| 726 | |
| 727 | if(littleEndian) { |
| 728 | uint8_t *s = ((uint8_t *)&dbl_value) + sizeof(dbl_value) - 2; |
Lev Walkin | 057fb73 | 2004-09-14 13:58:10 +0000 | [diff] [blame] | 729 | uint8_t *start = ((uint8_t *)&dbl_value); |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 730 | uint8_t *d; |
| 731 | |
| 732 | bmsign = 0x80 | ((s[1] >> 1) & 0x40); /* binary mask & - */ |
Lev Walkin | 057fb73 | 2004-09-14 13:58:10 +0000 | [diff] [blame] | 733 | for(mstop = d = dscr; s >= start; d++, s--) { |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 734 | *d = *s; |
| 735 | if(*d) mstop = d; |
| 736 | } |
| 737 | } else { |
| 738 | uint8_t *s = ((uint8_t *)&dbl_value) + 1; |
| 739 | uint8_t *end = ((uint8_t *)&dbl_value) + sizeof(double); |
| 740 | uint8_t *d; |
| 741 | |
| 742 | bmsign = 0x80 | ((s[-1] >> 1) & 0x40); /* binary mask & - */ |
| 743 | for(mstop = d = dscr; s < end; d++, s++) { |
| 744 | *d = *s; |
| 745 | if(*d) mstop = d; |
| 746 | } |
| 747 | } |
| 748 | |
| 749 | /* Remove parts of the exponent, leave mantissa and explicit 1. */ |
| 750 | dscr[0] = 0x10 | (dscr[0] & 0x0f); |
| 751 | |
| 752 | /* Adjust exponent in a very unobvious way */ |
| 753 | expval -= 8 * ((mstop - dscr) + 1) - 4; |
| 754 | |
| 755 | /* This loop ensures DER conformance by forcing mantissa odd: 11.3.1 */ |
| 756 | mval = *mstop; |
| 757 | if(mval && !(mval & 1)) { |
Lev Walkin | 9c57f33 | 2017-09-17 22:30:49 -0700 | [diff] [blame] | 758 | int shift_count = 1; |
| 759 | int ishift; |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 760 | uint8_t *mptr; |
| 761 | |
| 762 | /* |
| 763 | * Figure out what needs to be done to make mantissa odd. |
| 764 | */ |
| 765 | if(!(mval & 0x0f)) /* Speed-up a little */ |
| 766 | shift_count = 4; |
| 767 | while(((mval >> shift_count) & 1) == 0) |
| 768 | shift_count++; |
| 769 | |
| 770 | ishift = 8 - shift_count; |
| 771 | accum = 0; |
| 772 | |
| 773 | /* Go over the buffer, shifting it shift_count bits right. */ |
| 774 | for(mptr = dscr; mptr <= mstop; mptr++) { |
| 775 | mval = *mptr; |
| 776 | *mptr = accum | (mval >> shift_count); |
| 777 | accum = mval << ishift; |
| 778 | } |
| 779 | |
Lev Walkin | b89b340 | 2012-01-09 18:26:30 -0800 | [diff] [blame] | 780 | /* Adjust exponent appropriately. */ |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 781 | expval += shift_count; |
| 782 | } |
| 783 | |
| 784 | if(expval < 0) { |
| 785 | if((expval >> 7) == -1) { |
| 786 | *ptr++ = bmsign | 0x00; |
| 787 | *ptr++ = expval; |
| 788 | } else if((expval >> 15) == -1) { |
| 789 | *ptr++ = bmsign | 0x01; |
| 790 | *ptr++ = expval >> 8; |
| 791 | *ptr++ = expval; |
| 792 | } else { |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 793 | *ptr++ = bmsign | 0x02; |
| 794 | *ptr++ = expval >> 16; |
| 795 | *ptr++ = expval >> 8; |
| 796 | *ptr++ = expval; |
| 797 | } |
| 798 | } else if(expval <= 0x7f) { |
| 799 | *ptr++ = bmsign | 0x00; |
| 800 | *ptr++ = expval; |
| 801 | } else if(expval <= 0x7fff) { |
| 802 | *ptr++ = bmsign | 0x01; |
| 803 | *ptr++ = expval >> 8; |
| 804 | *ptr++ = expval; |
| 805 | } else { |
| 806 | assert(expval <= 0x7fffff); |
| 807 | *ptr++ = bmsign | 0x02; |
| 808 | *ptr++ = expval >> 16; |
| 809 | *ptr++ = expval >> 8; |
| 810 | *ptr++ = expval; |
| 811 | } |
| 812 | |
| 813 | buflen = (mstop - dscr) + 1; |
| 814 | memcpy(ptr, dscr, buflen); |
| 815 | ptr += buflen; |
| 816 | buflen = ptr - buf; |
| 817 | |
Lev Walkin | c17d90f | 2005-01-17 14:32:45 +0000 | [diff] [blame] | 818 | ptr = (uint8_t *)MALLOC(buflen + 1); |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 819 | if(!ptr) return -1; |
| 820 | |
| 821 | memcpy(ptr, buf, buflen); |
| 822 | buf[buflen] = 0; /* JIC */ |
| 823 | |
| 824 | if(st->buf) FREEMEM(st->buf); |
| 825 | st->buf = ptr; |
| 826 | st->size = buflen; |
| 827 | |
| 828 | return 0; |
| 829 | } |
Lev Walkin | 41d972b | 2017-08-23 23:30:59 -0700 | [diff] [blame] | 830 | |
| 831 | #ifndef ASN_DISABLE_PER_SUPPORT |
| 832 | |
| 833 | asn_dec_rval_t |
Lev Walkin | afbf2a9 | 2017-09-12 23:30:27 -0700 | [diff] [blame] | 834 | REAL_decode_uper(const asn_codec_ctx_t *opt_codec_ctx, asn_TYPE_descriptor_t *td, |
Lev Walkin | 41d972b | 2017-08-23 23:30:59 -0700 | [diff] [blame] | 835 | const asn_per_constraints_t *constraints, void **sptr, |
| 836 | asn_per_data_t *pd) { |
| 837 | (void)constraints; /* No PER visible constraints */ |
| 838 | return OCTET_STRING_decode_uper(opt_codec_ctx, td, 0, sptr, pd); |
| 839 | } |
| 840 | |
| 841 | asn_enc_rval_t |
| 842 | REAL_encode_uper(asn_TYPE_descriptor_t *td, |
| 843 | const asn_per_constraints_t *constraints, void *sptr, |
| 844 | asn_per_outp_t *po) { |
| 845 | (void)constraints; /* No PER visible constraints */ |
| 846 | return OCTET_STRING_encode_uper(td, 0, sptr, po); |
| 847 | } |
| 848 | |
| 849 | #endif /* ASN_DISABLE_PER_SUPPORT */ |
Lev Walkin | a5972be | 2017-09-29 23:15:58 -0700 | [diff] [blame] | 850 | |
| 851 | |
| 852 | asn_random_fill_result_t |
| 853 | REAL_random_fill(const asn_TYPE_descriptor_t *td, void **sptr, |
| 854 | const asn_encoding_constraints_t *constraints, |
| 855 | size_t max_length) { |
| 856 | asn_random_fill_result_t result_ok = {ARFILL_OK, 1}; |
| 857 | asn_random_fill_result_t result_failed = {ARFILL_FAILED, 0}; |
| 858 | asn_random_fill_result_t result_skipped = {ARFILL_SKIPPED, 0}; |
| 859 | static const double values[] = {0, -0.0, -1, 1, INFINITY, -INFINITY, NAN}; |
| 860 | REAL_t *st; |
| 861 | double d; |
| 862 | |
| 863 | (void)constraints; |
| 864 | |
| 865 | if(max_length == 0) return result_skipped; |
| 866 | |
| 867 | d = values[asn_random_between(0, sizeof(values) / sizeof(values[0]) - 1)]; |
| 868 | |
| 869 | if(*sptr) { |
| 870 | st = *sptr; |
| 871 | } else { |
| 872 | st = (REAL_t*)(*sptr = CALLOC(1, sizeof(REAL_t))); |
| 873 | if(!st) { |
| 874 | return result_failed; |
| 875 | } |
| 876 | } |
| 877 | |
| 878 | if(asn_double2REAL(st, d)) { |
| 879 | if(st == *sptr) { |
| 880 | ASN_STRUCT_RESET(*td, st); |
| 881 | } else { |
| 882 | ASN_STRUCT_FREE(*td, st); |
| 883 | } |
| 884 | return result_failed; |
| 885 | } |
| 886 | |
| 887 | result_ok.length = st->size; |
| 888 | return result_ok; |
| 889 | } |