vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 1 | /*- |
| 2 | * Copyright (c) 2004 Lev Walkin <vlm@lionet.info>. All rights reserved. |
| 3 | * Redistribution and modifications are permitted subject to BSD license. |
| 4 | */ |
vlm | 6c59384 | 2004-10-26 09:03:31 +0000 | [diff] [blame] | 5 | #if defined(__alpha) |
vlm | b8ce6f0 | 2004-10-26 09:35:25 +0000 | [diff] [blame] | 6 | #define _ISOC99_SOURCE /* For quiet NAN, through bits/nan.h */ |
| 7 | #define _BSD_SOURCE /* To reintroduce finite(3) */ |
| 8 | #include <sys/resource.h> /* For INFINITY */ |
vlm | 6c59384 | 2004-10-26 09:03:31 +0000 | [diff] [blame] | 9 | #endif |
vlm | 39ba4c4 | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 10 | #include <asn_internal.h> |
vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 11 | #include <stdlib.h> /* for strtod(3) */ |
| 12 | #include <math.h> |
| 13 | #include <errno.h> |
vlm | 6c59384 | 2004-10-26 09:03:31 +0000 | [diff] [blame] | 14 | #include <REAL.h> |
vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 15 | |
| 16 | #undef INT_MAX |
| 17 | #define INT_MAX ((int)(((unsigned int)-1) >> 1)) |
| 18 | |
vlm | b8ce6f0 | 2004-10-26 09:35:25 +0000 | [diff] [blame] | 19 | #if !(defined(NAN) || defined(INFINITY)) |
| 20 | static volatile double real_zero __attribute__ ((unused)) = 0.0; |
| 21 | #endif |
vlm | 4005630 | 2004-09-27 22:18:34 +0000 | [diff] [blame] | 22 | #ifndef NAN |
vlm | b1f8d61 | 2004-09-27 22:16:18 +0000 | [diff] [blame] | 23 | #define NAN (real_zero/real_zero) |
vlm | 1d3ed28 | 2004-09-14 13:40:42 +0000 | [diff] [blame] | 24 | #endif |
vlm | b8ce6f0 | 2004-10-26 09:35:25 +0000 | [diff] [blame] | 25 | #ifndef INFINITY |
| 26 | #define INFINITY (1.0/real_zero) |
| 27 | #endif |
vlm | 1d3ed28 | 2004-09-14 13:40:42 +0000 | [diff] [blame] | 28 | |
vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 29 | /* |
| 30 | * REAL basic type description. |
| 31 | */ |
vlm | ef6355b | 2004-09-29 13:26:15 +0000 | [diff] [blame] | 32 | static ber_tlv_tag_t asn_DEF_REAL_tags[] = { |
vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 33 | (ASN_TAG_CLASS_UNIVERSAL | (9 << 2)) |
| 34 | }; |
vlm | ef6355b | 2004-09-29 13:26:15 +0000 | [diff] [blame] | 35 | asn_TYPE_descriptor_t asn_DEF_REAL = { |
vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 36 | "REAL", |
vlm | 9de248e | 2004-10-20 15:50:55 +0000 | [diff] [blame] | 37 | "REAL", |
vlm | 6678cb1 | 2004-09-26 13:10:40 +0000 | [diff] [blame] | 38 | ASN__PRIMITIVE_TYPE_free, |
vlm | 39ba4c4 | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 39 | REAL_print, |
vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 40 | asn_generic_no_constraint, |
vlm | 6678cb1 | 2004-09-26 13:10:40 +0000 | [diff] [blame] | 41 | ber_decode_primitive, |
| 42 | der_encode_primitive, |
vlm | cce2ddb | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 43 | REAL_decode_xer, |
vlm | 39ba4c4 | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 44 | REAL_encode_xer, |
vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 45 | 0, /* Use generic outmost tag fetcher */ |
vlm | ef6355b | 2004-09-29 13:26:15 +0000 | [diff] [blame] | 46 | asn_DEF_REAL_tags, |
| 47 | sizeof(asn_DEF_REAL_tags) / sizeof(asn_DEF_REAL_tags[0]), |
| 48 | asn_DEF_REAL_tags, /* Same as above */ |
| 49 | sizeof(asn_DEF_REAL_tags) / sizeof(asn_DEF_REAL_tags[0]), |
vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 50 | 0, 0, /* No members */ |
| 51 | 0 /* No specifics */ |
| 52 | }; |
| 53 | |
vlm | cce2ddb | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 54 | typedef enum specialRealValue { |
| 55 | SRV__NOT_A_NUMBER, |
| 56 | SRV__MINUS_INFINITY, |
| 57 | SRV__PLUS_INFINITY |
| 58 | } specialRealValue_e; |
| 59 | static struct specialRealValue_s { |
| 60 | char *string; |
vlm | b848b2a | 2004-10-21 14:02:19 +0000 | [diff] [blame] | 61 | size_t length; |
vlm | b8ce6f0 | 2004-10-26 09:35:25 +0000 | [diff] [blame] | 62 | long dv; |
vlm | cce2ddb | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 63 | } specialRealValue[] = { |
| 64 | #define SRV_SET(foo, val) { foo, sizeof(foo) - 1, val } |
vlm | b8ce6f0 | 2004-10-26 09:35:25 +0000 | [diff] [blame] | 65 | SRV_SET("<NOT-A-NUMBER/>", 0), |
| 66 | SRV_SET("<MINUS-INFINITY/>", -1), |
| 67 | SRV_SET("<PLUS-INFINITY/>", 1), |
vlm | cce2ddb | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 68 | #undef SRV_SET |
| 69 | }; |
| 70 | |
vlm | 39ba4c4 | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 71 | ssize_t |
| 72 | REAL__dump(double d, int canonical, asn_app_consume_bytes_f *cb, void *app_key) { |
vlm | 0a1738c | 2004-09-27 20:52:18 +0000 | [diff] [blame] | 73 | char local_buf[64]; |
vlm | 39ba4c4 | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 74 | char *buf = local_buf; |
| 75 | ssize_t buflen = sizeof(local_buf); |
vlm | 0a1738c | 2004-09-27 20:52:18 +0000 | [diff] [blame] | 76 | const char *fmt = canonical?"%.15E":"%.15f"; |
vlm | 39ba4c4 | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 77 | ssize_t ret; |
| 78 | |
vlm | 6678cb1 | 2004-09-26 13:10:40 +0000 | [diff] [blame] | 79 | /* |
| 80 | * Check whether it is a special value. |
| 81 | */ |
vlm | b1f8d61 | 2004-09-27 22:16:18 +0000 | [diff] [blame] | 82 | /* fpclassify(3) is not portable yet */ |
| 83 | if(isnan(d)) { |
vlm | cce2ddb | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 84 | buf = specialRealValue[SRV__NOT_A_NUMBER].string; |
| 85 | buflen = specialRealValue[SRV__NOT_A_NUMBER].length; |
vlm | b1f8d61 | 2004-09-27 22:16:18 +0000 | [diff] [blame] | 86 | return (cb(buf, buflen, app_key) < 0) ? -1 : buflen; |
vlm | 6534a8d | 2004-10-20 15:40:04 +0000 | [diff] [blame] | 87 | } else if(!finite(d)) { |
vlm | b1f8d61 | 2004-09-27 22:16:18 +0000 | [diff] [blame] | 88 | if(copysign(1.0, d) < 0.0) { |
vlm | cce2ddb | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 89 | buf = specialRealValue[SRV__MINUS_INFINITY].string; |
| 90 | buflen = specialRealValue[SRV__MINUS_INFINITY].length; |
vlm | b1f8d61 | 2004-09-27 22:16:18 +0000 | [diff] [blame] | 91 | } else { |
vlm | cce2ddb | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 92 | buf = specialRealValue[SRV__PLUS_INFINITY].string; |
| 93 | buflen = specialRealValue[SRV__PLUS_INFINITY].length; |
vlm | b1f8d61 | 2004-09-27 22:16:18 +0000 | [diff] [blame] | 94 | } |
| 95 | return (cb(buf, buflen, app_key) < 0) ? -1 : buflen; |
| 96 | } else if(ilogb(d) <= -INT_MAX) { |
| 97 | if(copysign(1.0, d) < 0.0) { |
| 98 | buf = "-0"; |
| 99 | buflen = 2; |
| 100 | } else { |
| 101 | buf = "0"; |
| 102 | buflen = 1; |
vlm | 6678cb1 | 2004-09-26 13:10:40 +0000 | [diff] [blame] | 103 | } |
| 104 | return (cb(buf, buflen, app_key) < 0) ? -1 : buflen; |
| 105 | } |
| 106 | |
| 107 | /* |
| 108 | * Use the libc's double printing, hopefully they got it right. |
| 109 | */ |
vlm | 39ba4c4 | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 110 | do { |
| 111 | ret = snprintf(buf, buflen, fmt, d); |
| 112 | if(ret < 0) { |
| 113 | buflen <<= 1; |
| 114 | } else if(ret >= buflen) { |
| 115 | buflen = ret + 1; |
| 116 | } else { |
| 117 | buflen = ret; |
| 118 | break; |
| 119 | } |
| 120 | if(buf != local_buf) free(buf); |
vlm | 6678cb1 | 2004-09-26 13:10:40 +0000 | [diff] [blame] | 121 | buf = (char *)MALLOC(buflen); |
vlm | 39ba4c4 | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 122 | if(!buf) return -1; |
| 123 | } while(1); |
| 124 | |
vlm | 39ba4c4 | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 125 | if(canonical) { |
vlm | 0a1738c | 2004-09-27 20:52:18 +0000 | [diff] [blame] | 126 | /* |
| 127 | * Transform the "[-]d.dddE+-dd" output into "[-]d.dddE[-]d" |
vlm | cce2ddb | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 128 | * Check that snprintf() constructed the output correctly. |
vlm | 0a1738c | 2004-09-27 20:52:18 +0000 | [diff] [blame] | 129 | */ |
vlm | 39ba4c4 | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 130 | char *dot, *E; |
| 131 | char *end = buf + buflen; |
vlm | 0a1738c | 2004-09-27 20:52:18 +0000 | [diff] [blame] | 132 | char *last_zero; |
vlm | 39ba4c4 | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 133 | |
vlm | cce2ddb | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 134 | dot = (buf[0] == 0x2d /* '-' */) ? (buf + 2) : (buf + 1); |
vlm | 39ba4c4 | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 135 | if(*dot >= 0x30) { |
| 136 | errno = EINVAL; |
| 137 | return -1; /* Not a dot, really */ |
| 138 | } |
vlm | 0a1738c | 2004-09-27 20:52:18 +0000 | [diff] [blame] | 139 | *dot = 0x2e; /* Replace possible comma */ |
vlm | 39ba4c4 | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 140 | |
vlm | 0a1738c | 2004-09-27 20:52:18 +0000 | [diff] [blame] | 141 | for(last_zero = dot + 2, E = dot; dot < end; E++) { |
| 142 | if(*E == 0x45) { |
| 143 | char *expptr = ++E; |
| 144 | char *s = expptr; |
| 145 | int sign; |
vlm | cce2ddb | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 146 | if(*expptr == 0x2b /* '+' */) { |
vlm | 0a1738c | 2004-09-27 20:52:18 +0000 | [diff] [blame] | 147 | /* Skip the "+" */ |
vlm | 39ba4c4 | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 148 | buflen -= 1; |
vlm | 0a1738c | 2004-09-27 20:52:18 +0000 | [diff] [blame] | 149 | sign = 0; |
| 150 | } else { |
| 151 | sign = 1; |
vlm | 39ba4c4 | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 152 | s++; |
| 153 | } |
vlm | 0a1738c | 2004-09-27 20:52:18 +0000 | [diff] [blame] | 154 | expptr++; |
| 155 | if(expptr > end) { |
vlm | 39ba4c4 | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 156 | errno = EINVAL; |
| 157 | return -1; |
| 158 | } |
vlm | 0a1738c | 2004-09-27 20:52:18 +0000 | [diff] [blame] | 159 | if(*expptr == 0x30) { |
| 160 | buflen--; |
| 161 | expptr++; |
| 162 | } |
| 163 | if(*last_zero == 0x30) { |
| 164 | *last_zero = 0x45; /* E */ |
vlm | cce2ddb | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 165 | buflen -= s - (last_zero + 1); |
vlm | 0a1738c | 2004-09-27 20:52:18 +0000 | [diff] [blame] | 166 | s = last_zero + 1; |
vlm | cce2ddb | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 167 | if(sign) { |
| 168 | *s++ = 0x2d /* '-' */; |
| 169 | buflen++; |
| 170 | } |
vlm | 0a1738c | 2004-09-27 20:52:18 +0000 | [diff] [blame] | 171 | } |
| 172 | for(; expptr <= end; s++, expptr++) |
| 173 | *s = *expptr; |
| 174 | break; |
| 175 | } else if(*E == 0x30) { |
| 176 | if(*last_zero != 0x30) |
| 177 | last_zero = E; |
vlm | 39ba4c4 | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 178 | } |
| 179 | } |
| 180 | if(E == end) { |
| 181 | errno = EINVAL; |
| 182 | return -1; /* No promised E */ |
| 183 | } |
vlm | 0a1738c | 2004-09-27 20:52:18 +0000 | [diff] [blame] | 184 | } else { |
| 185 | /* |
| 186 | * Remove trailing zeros. |
| 187 | */ |
| 188 | char *end = buf + buflen; |
| 189 | char *last_zero = end; |
vlm | 0919fcc | 2005-04-25 21:08:25 +0000 | [diff] [blame] | 190 | int stoplooking = 0; |
vlm | 0a1738c | 2004-09-27 20:52:18 +0000 | [diff] [blame] | 191 | char *z; |
| 192 | for(z = end - 1; z > buf; z--) { |
| 193 | switch(*z) { |
vlm | 0919fcc | 2005-04-25 21:08:25 +0000 | [diff] [blame] | 194 | case 0x30: |
| 195 | if(!stoplooking) |
| 196 | last_zero = z; |
| 197 | continue; |
vlm | 0a1738c | 2004-09-27 20:52:18 +0000 | [diff] [blame] | 198 | case 0x31: case 0x32: case 0x33: case 0x34: |
| 199 | case 0x35: case 0x36: case 0x37: case 0x38: case 0x39: |
vlm | 0919fcc | 2005-04-25 21:08:25 +0000 | [diff] [blame] | 200 | stoplooking = 1; |
vlm | 0a1738c | 2004-09-27 20:52:18 +0000 | [diff] [blame] | 201 | continue; |
| 202 | default: /* Catch dot and other separators */ |
vlm | 0919fcc | 2005-04-25 21:08:25 +0000 | [diff] [blame] | 203 | /* |
| 204 | * Replace possible comma (which may even |
| 205 | * be not a comma at all: locale-defined). |
| 206 | */ |
| 207 | *z = 0x2e; |
vlm | 0a1738c | 2004-09-27 20:52:18 +0000 | [diff] [blame] | 208 | if(last_zero == z + 1) { /* leave x.0 */ |
| 209 | last_zero++; |
| 210 | } |
| 211 | buflen = last_zero - buf; |
| 212 | *last_zero = '\0'; |
| 213 | break; |
| 214 | } |
| 215 | break; |
| 216 | } |
vlm | 39ba4c4 | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 217 | } |
| 218 | |
| 219 | ret = cb(buf, buflen, app_key); |
| 220 | if(buf != local_buf) free(buf); |
| 221 | return (ret < 0) ? -1 : buflen; |
| 222 | } |
| 223 | |
vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 224 | int |
vlm | ef6355b | 2004-09-29 13:26:15 +0000 | [diff] [blame] | 225 | REAL_print(asn_TYPE_descriptor_t *td, const void *sptr, int ilevel, |
vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 226 | asn_app_consume_bytes_f *cb, void *app_key) { |
| 227 | const REAL_t *st = (const REAL_t *)sptr; |
vlm | 6678cb1 | 2004-09-26 13:10:40 +0000 | [diff] [blame] | 228 | ssize_t ret; |
vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 229 | double d; |
vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 230 | |
| 231 | (void)td; /* Unused argument */ |
| 232 | (void)ilevel; /* Unused argument */ |
| 233 | |
vlm | 39ba4c4 | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 234 | if(!st || !st->buf) |
vlm | 6678cb1 | 2004-09-26 13:10:40 +0000 | [diff] [blame] | 235 | ret = cb("<absent>", 8, app_key); |
vlm | ef6355b | 2004-09-29 13:26:15 +0000 | [diff] [blame] | 236 | else if(asn_REAL2double(st, &d)) |
vlm | 6678cb1 | 2004-09-26 13:10:40 +0000 | [diff] [blame] | 237 | ret = cb("<error>", 7, app_key); |
| 238 | else |
| 239 | ret = REAL__dump(d, 0, cb, app_key); |
vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 240 | |
vlm | 6678cb1 | 2004-09-26 13:10:40 +0000 | [diff] [blame] | 241 | return (ret < 0) ? -1 : 0; |
vlm | 39ba4c4 | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 242 | } |
vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 243 | |
vlm | 39ba4c4 | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 244 | asn_enc_rval_t |
vlm | ef6355b | 2004-09-29 13:26:15 +0000 | [diff] [blame] | 245 | REAL_encode_xer(asn_TYPE_descriptor_t *td, void *sptr, |
vlm | 39ba4c4 | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 246 | int ilevel, enum xer_encoder_flags_e flags, |
| 247 | asn_app_consume_bytes_f *cb, void *app_key) { |
| 248 | REAL_t *st = (REAL_t *)sptr; |
| 249 | asn_enc_rval_t er; |
| 250 | double d; |
| 251 | |
| 252 | (void)ilevel; |
| 253 | |
vlm | ef6355b | 2004-09-29 13:26:15 +0000 | [diff] [blame] | 254 | if(!st || !st->buf || asn_REAL2double(st, &d)) |
vlm | 39ba4c4 | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 255 | _ASN_ENCODE_FAILED; |
| 256 | |
| 257 | er.encoded = REAL__dump(d, flags & XER_F_CANONICAL, cb, app_key); |
| 258 | if(er.encoded < 0) _ASN_ENCODE_FAILED; |
| 259 | |
| 260 | return er; |
vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 261 | } |
| 262 | |
vlm | cce2ddb | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 263 | |
| 264 | /* |
| 265 | * Decode the chunk of XML text encoding REAL. |
| 266 | */ |
vlm | 4df9cc1 | 2005-03-09 22:19:25 +0000 | [diff] [blame] | 267 | static enum xer_pbd_rval |
| 268 | REAL__xer_body_decode(asn_TYPE_descriptor_t *td, void *sptr, const void *chunk_buf, size_t chunk_size) { |
vlm | b848b2a | 2004-10-21 14:02:19 +0000 | [diff] [blame] | 269 | REAL_t *st = (REAL_t *)sptr; |
vlm | cce2ddb | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 270 | double value; |
vlm | 4df9cc1 | 2005-03-09 22:19:25 +0000 | [diff] [blame] | 271 | const char *xerdata = (const char *)chunk_buf; |
vlm | cce2ddb | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 272 | char *endptr = 0; |
| 273 | char *b; |
| 274 | |
vlm | 80a4859 | 2005-02-25 12:10:27 +0000 | [diff] [blame] | 275 | (void)td; |
| 276 | |
vlm | 4df9cc1 | 2005-03-09 22:19:25 +0000 | [diff] [blame] | 277 | if(!chunk_size) return XPBD_BROKEN_ENCODING; |
vlm | cce2ddb | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 278 | |
| 279 | /* |
| 280 | * Decode an XMLSpecialRealValue: <MINUS-INFINITY>, etc. |
| 281 | */ |
| 282 | if(xerdata[0] == 0x3c /* '<' */) { |
| 283 | size_t i; |
| 284 | for(i = 0; i < sizeof(specialRealValue) |
| 285 | / sizeof(specialRealValue[0]); i++) { |
| 286 | struct specialRealValue_s *srv = &specialRealValue[i]; |
vlm | b8ce6f0 | 2004-10-26 09:35:25 +0000 | [diff] [blame] | 287 | double dv; |
| 288 | |
vlm | cce2ddb | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 289 | if(srv->length != chunk_size |
| 290 | || memcmp(srv->string, chunk_buf, chunk_size)) |
| 291 | continue; |
| 292 | |
vlm | b8ce6f0 | 2004-10-26 09:35:25 +0000 | [diff] [blame] | 293 | /* |
| 294 | * It could've been done using |
| 295 | * (double)srv->dv / real_zero, |
| 296 | * but it summons fp exception on some platforms. |
| 297 | */ |
| 298 | switch(srv->dv) { |
| 299 | case -1: dv = - INFINITY; break; |
| 300 | case 0: dv = NAN; break; |
| 301 | case 1: dv = INFINITY; break; |
vlm | 4df9cc1 | 2005-03-09 22:19:25 +0000 | [diff] [blame] | 302 | default: return XPBD_SYSTEM_FAILURE; |
vlm | b8ce6f0 | 2004-10-26 09:35:25 +0000 | [diff] [blame] | 303 | } |
| 304 | |
vlm | 4df9cc1 | 2005-03-09 22:19:25 +0000 | [diff] [blame] | 305 | if(asn_double2REAL(st, dv)) |
| 306 | return XPBD_SYSTEM_FAILURE; |
vlm | cce2ddb | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 307 | |
vlm | 4df9cc1 | 2005-03-09 22:19:25 +0000 | [diff] [blame] | 308 | return XPBD_BODY_CONSUMED; |
vlm | cce2ddb | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 309 | } |
| 310 | ASN_DEBUG("Unknown XMLSpecialRealValue"); |
vlm | 4df9cc1 | 2005-03-09 22:19:25 +0000 | [diff] [blame] | 311 | return XPBD_BROKEN_ENCODING; |
vlm | cce2ddb | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 312 | } |
| 313 | |
| 314 | /* |
| 315 | * Copy chunk into the nul-terminated string, and run strtod. |
| 316 | */ |
vlm | b848b2a | 2004-10-21 14:02:19 +0000 | [diff] [blame] | 317 | b = (char *)MALLOC(chunk_size + 1); |
vlm | 4df9cc1 | 2005-03-09 22:19:25 +0000 | [diff] [blame] | 318 | if(!b) return XPBD_SYSTEM_FAILURE; |
vlm | cce2ddb | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 319 | memcpy(b, chunk_buf, chunk_size); |
vlm | b848b2a | 2004-10-21 14:02:19 +0000 | [diff] [blame] | 320 | b[chunk_size] = 0; /* nul-terminate */ |
vlm | cce2ddb | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 321 | |
| 322 | value = strtod(b, &endptr); |
| 323 | free(b); |
vlm | 4df9cc1 | 2005-03-09 22:19:25 +0000 | [diff] [blame] | 324 | if(endptr == b) return XPBD_BROKEN_ENCODING; |
vlm | cce2ddb | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 325 | |
| 326 | if(asn_double2REAL(st, value)) |
vlm | 4df9cc1 | 2005-03-09 22:19:25 +0000 | [diff] [blame] | 327 | return XPBD_SYSTEM_FAILURE; |
vlm | cce2ddb | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 328 | |
vlm | 4df9cc1 | 2005-03-09 22:19:25 +0000 | [diff] [blame] | 329 | return XPBD_BODY_CONSUMED; |
vlm | cce2ddb | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 330 | } |
| 331 | |
| 332 | asn_dec_rval_t |
| 333 | REAL_decode_xer(asn_codec_ctx_t *opt_codec_ctx, |
| 334 | asn_TYPE_descriptor_t *td, void **sptr, const char *opt_mname, |
vlm | b02dcc6 | 2005-03-10 18:52:02 +0000 | [diff] [blame] | 335 | const void *buf_ptr, size_t size) { |
vlm | cce2ddb | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 336 | |
| 337 | return xer_decode_primitive(opt_codec_ctx, td, |
vlm | b848b2a | 2004-10-21 14:02:19 +0000 | [diff] [blame] | 338 | sptr, sizeof(REAL_t), opt_mname, |
vlm | cce2ddb | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 339 | buf_ptr, size, REAL__xer_body_decode); |
| 340 | } |
| 341 | |
| 342 | |
vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 343 | int |
vlm | ef6355b | 2004-09-29 13:26:15 +0000 | [diff] [blame] | 344 | asn_REAL2double(const REAL_t *st, double *dbl_value) { |
vlm | 39ba4c4 | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 345 | unsigned int octv; |
vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 346 | |
| 347 | if(!st || !st->buf) { |
| 348 | errno = EINVAL; |
| 349 | return -1; |
| 350 | } |
| 351 | |
| 352 | if(st->size == 0) { |
| 353 | *dbl_value = 0; |
| 354 | return 0; |
| 355 | } |
| 356 | |
| 357 | octv = st->buf[0]; /* unsigned byte */ |
| 358 | |
| 359 | switch(octv & 0xC0) { |
| 360 | case 0x40: /* X.690: 8.5.8 */ |
| 361 | /* "SpecialRealValue" */ |
| 362 | |
| 363 | /* Be liberal in what you accept... |
| 364 | if(st->size != 1) ... |
| 365 | */ |
| 366 | |
| 367 | switch(st->buf[0]) { |
| 368 | case 0x40: /* 01000000: PLUS-INFINITY */ |
vlm | b8ce6f0 | 2004-10-26 09:35:25 +0000 | [diff] [blame] | 369 | *dbl_value = INFINITY; |
vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 370 | return 0; |
| 371 | case 0x41: /* 01000001: MINUS-INFINITY */ |
vlm | b8ce6f0 | 2004-10-26 09:35:25 +0000 | [diff] [blame] | 372 | *dbl_value = - INFINITY; |
vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 373 | return 0; |
| 374 | /* |
| 375 | * The following cases are defined by |
| 376 | * X.690 Amendment 1 (10/03) |
| 377 | */ |
| 378 | case 0x42: /* 01000010: NOT-A-NUMBER */ |
| 379 | *dbl_value = NAN; |
| 380 | return 0; |
| 381 | case 0x43: /* 01000011: minus zero */ |
vlm | 8f6ded0 | 2004-09-27 21:36:59 +0000 | [diff] [blame] | 382 | *dbl_value = -0.0; |
vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 383 | return 0; |
| 384 | } |
| 385 | |
| 386 | errno = EINVAL; |
| 387 | return -1; |
| 388 | case 0x00: { /* X.690: 8.5.6 */ |
| 389 | /* |
| 390 | * Decimal. NR{1,2,3} format. |
| 391 | */ |
| 392 | double d; |
| 393 | |
| 394 | assert(st->buf[st->size - 1] == 0); /* Security, vashu mat' */ |
| 395 | |
| 396 | d = strtod((char *)st->buf, 0); |
| 397 | if(finite(d)) { |
| 398 | *dbl_value = d; |
| 399 | return 0; |
| 400 | } else { |
| 401 | errno = ERANGE; |
| 402 | return 0; |
| 403 | } |
| 404 | } |
| 405 | } |
| 406 | |
| 407 | /* |
| 408 | * Binary representation. |
| 409 | */ |
| 410 | { |
| 411 | double m; |
| 412 | int expval; /* exponent value */ |
| 413 | unsigned int elen; /* exponent value length, in octets */ |
| 414 | unsigned int scaleF; |
| 415 | unsigned int baseF; |
| 416 | uint8_t *ptr; |
| 417 | uint8_t *end; |
| 418 | int sign; |
| 419 | |
| 420 | switch((octv & 0x30) >> 4) { |
| 421 | case 0x00: baseF = 1; break; /* base 2 */ |
| 422 | case 0x01: baseF = 3; break; /* base 8 */ |
| 423 | case 0x02: baseF = 4; break; /* base 16 */ |
| 424 | default: |
| 425 | /* Reserved field, can't parse now. */ |
| 426 | errno = EINVAL; |
| 427 | return -1; |
| 428 | } |
| 429 | |
| 430 | sign = (octv & 0x40); /* bit 7 */ |
| 431 | scaleF = (octv & 0x0C) >> 2; /* bits 4 to 3 */ |
| 432 | |
vlm | 39ba4c4 | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 433 | if(st->size <= (int)(1 + (octv & 0x03))) { |
vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 434 | errno = EINVAL; |
| 435 | return -1; |
| 436 | } |
| 437 | |
| 438 | if((octv & 0x03) == 0x11) { |
| 439 | /* 8.5.6.4, case d) */ |
| 440 | elen = st->buf[1]; /* unsigned binary number */ |
vlm | 39ba4c4 | 2004-09-22 16:06:28 +0000 | [diff] [blame] | 441 | if(elen == 0 || st->size <= (int)(2 + elen)) { |
vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 442 | errno = EINVAL; |
| 443 | return -1; |
| 444 | } |
| 445 | ptr = &st->buf[2]; |
| 446 | } else { |
| 447 | elen = (octv & 0x03); |
| 448 | ptr = &st->buf[1]; |
| 449 | } |
| 450 | |
| 451 | /* Fetch the multibyte exponent */ |
| 452 | expval = (int)(*(int8_t *)ptr); |
| 453 | end = ptr + elen + 1; |
| 454 | for(ptr++; ptr < end; ptr++) |
| 455 | expval = (expval * 256) + *ptr; |
| 456 | |
| 457 | m = 0.0; /* Initial mantissa value */ |
| 458 | |
| 459 | /* Okay, the exponent is here. Now, what about mantissa? */ |
| 460 | end = st->buf + st->size; |
| 461 | if(ptr < end) { |
| 462 | for(; ptr < end; ptr++) |
vlm | e4e13e5 | 2005-04-28 22:56:36 +0000 | [diff] [blame] | 463 | m = ldexp(m, 8) + *ptr; |
vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 464 | } |
| 465 | |
vlm | cce2ddb | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 466 | if(0) |
vlm | e4e13e5 | 2005-04-28 22:56:36 +0000 | [diff] [blame] | 467 | ASN_DEBUG("m=%.10f, scF=%d, bF=%d, expval=%d, ldexp()=%f, ldexp()=%f", |
vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 468 | m, scaleF, baseF, expval, |
| 469 | ldexp(m, expval * baseF + scaleF), |
vlm | e4e13e5 | 2005-04-28 22:56:36 +0000 | [diff] [blame] | 470 | ldexp(m, scaleF) * pow(pow(2, baseF), expval) |
vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 471 | ); |
| 472 | |
| 473 | /* |
| 474 | * (S * N * 2^F) * B^E |
| 475 | * Essentially: |
vlm | e4e13e5 | 2005-04-28 22:56:36 +0000 | [diff] [blame] | 476 | m = ldexp(m, scaleF) * pow(pow(2, base), expval); |
vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 477 | */ |
| 478 | m = ldexp(m, expval * baseF + scaleF); |
| 479 | if(finite(m)) { |
| 480 | *dbl_value = sign ? -m : m; |
| 481 | } else { |
| 482 | errno = ERANGE; |
| 483 | return -1; |
| 484 | } |
| 485 | |
| 486 | } /* if(binary_format) */ |
| 487 | |
| 488 | return 0; |
| 489 | } |
| 490 | |
| 491 | /* |
| 492 | * Assume IEEE 754 floating point: standard 64 bit double. |
| 493 | * [1 bit sign] [11 bits exponent] [52 bits mantissa] |
| 494 | */ |
| 495 | int |
vlm | ef6355b | 2004-09-29 13:26:15 +0000 | [diff] [blame] | 496 | asn_double2REAL(REAL_t *st, double dbl_value) { |
vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 497 | #ifdef WORDS_BIGENDIAN /* Known to be big-endian */ |
| 498 | int littleEndian = 0; |
| 499 | #else /* need to test: have no explicit information */ |
| 500 | unsigned int LE = 1; |
| 501 | int littleEndian = *(unsigned char *)&LE; |
| 502 | #endif |
| 503 | uint8_t buf[16]; /* More than enough for 8-byte dbl_value */ |
| 504 | uint8_t dscr[sizeof(dbl_value)]; /* double value scratch pad */ |
| 505 | /* Assertion guards: won't even compile, if unexpected double size */ |
| 506 | char assertion_buffer1[9 - sizeof(dbl_value)] __attribute__((unused)); |
| 507 | char assertion_buffer2[sizeof(dbl_value) - 7] __attribute__((unused)); |
| 508 | uint8_t *ptr = buf; |
| 509 | uint8_t *mstop; /* Last byte of mantissa */ |
| 510 | unsigned int mval; /* Value of the last byte of mantissa */ |
| 511 | unsigned int bmsign; /* binary mask with sign */ |
| 512 | unsigned int buflen; |
| 513 | unsigned int accum; |
| 514 | int expval; |
| 515 | |
| 516 | if(!st) { |
| 517 | errno = EINVAL; |
| 518 | return -1; |
| 519 | } |
| 520 | |
vlm | 4636198 | 2004-09-14 13:58:10 +0000 | [diff] [blame] | 521 | /* |
| 522 | * ilogb(+-0) returns -INT_MAX or INT_MIN (platform-dependent) |
vlm | cb4aaa7 | 2005-07-01 08:28:18 +0000 | [diff] [blame] | 523 | * ilogb(+-inf) returns INT_MAX, logb(+-inf) returns +inf |
vlm | 8f6ded0 | 2004-09-27 21:36:59 +0000 | [diff] [blame] | 524 | * ilogb(NaN) returns INT_MIN or INT_MAX (platform-dependent) |
vlm | 4636198 | 2004-09-14 13:58:10 +0000 | [diff] [blame] | 525 | */ |
vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 526 | expval = ilogb(dbl_value); |
vlm | 8f6ded0 | 2004-09-27 21:36:59 +0000 | [diff] [blame] | 527 | if(expval <= -INT_MAX /* Also catches +-0 and maybe isnan() */ |
| 528 | || expval == INT_MAX /* catches isfin() and maybe isnan() */ |
vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 529 | ) { |
| 530 | if(!st->buf || st->size < 2) { |
vlm | 6678cb1 | 2004-09-26 13:10:40 +0000 | [diff] [blame] | 531 | ptr = (uint8_t *)MALLOC(2); |
vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 532 | if(!ptr) return -1; |
| 533 | st->buf = ptr; |
| 534 | } |
| 535 | /* fpclassify(3) is not portable yet */ |
vlm | 8f6ded0 | 2004-09-27 21:36:59 +0000 | [diff] [blame] | 536 | if(isnan(dbl_value)) { |
| 537 | st->buf[0] = 0x42; /* NaN */ |
| 538 | st->buf[1] = 0; |
| 539 | st->size = 1; |
vlm | 6534a8d | 2004-10-20 15:40:04 +0000 | [diff] [blame] | 540 | } else if(!finite(dbl_value)) { |
vlm | b02ccb2 | 2004-09-14 13:50:21 +0000 | [diff] [blame] | 541 | if(copysign(1.0, dbl_value) < 0.0) { |
vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 542 | st->buf[0] = 0x41; /* MINUS-INFINITY */ |
| 543 | } else { |
| 544 | st->buf[0] = 0x40; /* PLUS-INFINITY */ |
| 545 | } |
| 546 | st->buf[1] = 0; |
| 547 | st->size = 1; |
vlm | b1f8d61 | 2004-09-27 22:16:18 +0000 | [diff] [blame] | 548 | } else { |
| 549 | if(copysign(1.0, dbl_value) < 0.0) { |
| 550 | st->buf[0] = 0x80 | 0x40; |
| 551 | st->buf[1] = 0; |
| 552 | st->size = 2; |
| 553 | } else { |
| 554 | /* no content octets: positive zero */ |
| 555 | st->buf[0] = 0; /* JIC */ |
| 556 | st->size = 0; |
| 557 | } |
vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 558 | } |
| 559 | return 0; |
| 560 | } |
| 561 | |
| 562 | if(littleEndian) { |
| 563 | uint8_t *s = ((uint8_t *)&dbl_value) + sizeof(dbl_value) - 2; |
vlm | 4636198 | 2004-09-14 13:58:10 +0000 | [diff] [blame] | 564 | uint8_t *start = ((uint8_t *)&dbl_value); |
vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 565 | uint8_t *d; |
| 566 | |
| 567 | bmsign = 0x80 | ((s[1] >> 1) & 0x40); /* binary mask & - */ |
vlm | 4636198 | 2004-09-14 13:58:10 +0000 | [diff] [blame] | 568 | for(mstop = d = dscr; s >= start; d++, s--) { |
vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 569 | *d = *s; |
| 570 | if(*d) mstop = d; |
| 571 | } |
| 572 | } else { |
| 573 | uint8_t *s = ((uint8_t *)&dbl_value) + 1; |
| 574 | uint8_t *end = ((uint8_t *)&dbl_value) + sizeof(double); |
| 575 | uint8_t *d; |
| 576 | |
| 577 | bmsign = 0x80 | ((s[-1] >> 1) & 0x40); /* binary mask & - */ |
| 578 | for(mstop = d = dscr; s < end; d++, s++) { |
| 579 | *d = *s; |
| 580 | if(*d) mstop = d; |
| 581 | } |
| 582 | } |
| 583 | |
| 584 | /* Remove parts of the exponent, leave mantissa and explicit 1. */ |
| 585 | dscr[0] = 0x10 | (dscr[0] & 0x0f); |
| 586 | |
| 587 | /* Adjust exponent in a very unobvious way */ |
| 588 | expval -= 8 * ((mstop - dscr) + 1) - 4; |
| 589 | |
| 590 | /* This loop ensures DER conformance by forcing mantissa odd: 11.3.1 */ |
| 591 | mval = *mstop; |
| 592 | if(mval && !(mval & 1)) { |
| 593 | unsigned int shift_count = 1; |
| 594 | unsigned int ishift; |
| 595 | uint8_t *mptr; |
| 596 | |
| 597 | /* |
| 598 | * Figure out what needs to be done to make mantissa odd. |
| 599 | */ |
| 600 | if(!(mval & 0x0f)) /* Speed-up a little */ |
| 601 | shift_count = 4; |
| 602 | while(((mval >> shift_count) & 1) == 0) |
| 603 | shift_count++; |
| 604 | |
| 605 | ishift = 8 - shift_count; |
| 606 | accum = 0; |
| 607 | |
| 608 | /* Go over the buffer, shifting it shift_count bits right. */ |
| 609 | for(mptr = dscr; mptr <= mstop; mptr++) { |
| 610 | mval = *mptr; |
| 611 | *mptr = accum | (mval >> shift_count); |
| 612 | accum = mval << ishift; |
| 613 | } |
| 614 | |
| 615 | /* Adjust mantissa appropriately. */ |
| 616 | expval += shift_count; |
| 617 | } |
| 618 | |
| 619 | if(expval < 0) { |
| 620 | if((expval >> 7) == -1) { |
| 621 | *ptr++ = bmsign | 0x00; |
| 622 | *ptr++ = expval; |
| 623 | } else if((expval >> 15) == -1) { |
| 624 | *ptr++ = bmsign | 0x01; |
| 625 | *ptr++ = expval >> 8; |
| 626 | *ptr++ = expval; |
| 627 | } else { |
vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 628 | *ptr++ = bmsign | 0x02; |
| 629 | *ptr++ = expval >> 16; |
| 630 | *ptr++ = expval >> 8; |
| 631 | *ptr++ = expval; |
| 632 | } |
| 633 | } else if(expval <= 0x7f) { |
| 634 | *ptr++ = bmsign | 0x00; |
| 635 | *ptr++ = expval; |
| 636 | } else if(expval <= 0x7fff) { |
| 637 | *ptr++ = bmsign | 0x01; |
| 638 | *ptr++ = expval >> 8; |
| 639 | *ptr++ = expval; |
| 640 | } else { |
| 641 | assert(expval <= 0x7fffff); |
| 642 | *ptr++ = bmsign | 0x02; |
| 643 | *ptr++ = expval >> 16; |
| 644 | *ptr++ = expval >> 8; |
| 645 | *ptr++ = expval; |
| 646 | } |
| 647 | |
| 648 | buflen = (mstop - dscr) + 1; |
| 649 | memcpy(ptr, dscr, buflen); |
| 650 | ptr += buflen; |
| 651 | buflen = ptr - buf; |
| 652 | |
vlm | a8f4a4c | 2005-01-17 14:32:45 +0000 | [diff] [blame] | 653 | ptr = (uint8_t *)MALLOC(buflen + 1); |
vlm | 785435b | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 654 | if(!ptr) return -1; |
| 655 | |
| 656 | memcpy(ptr, buf, buflen); |
| 657 | buf[buflen] = 0; /* JIC */ |
| 658 | |
| 659 | if(st->buf) FREEMEM(st->buf); |
| 660 | st->buf = ptr; |
| 661 | st->size = buflen; |
| 662 | |
| 663 | return 0; |
| 664 | } |