Lev Walkin | 4eceeba | 2007-07-23 06:48:26 +0000 | [diff] [blame] | 1 | #include <stdio.h> |
| 2 | #include <assert.h> |
| 3 | #include <math.h> |
| 4 | |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 5 | #define EMIT_ASN_DEBUG 1 |
Lev Walkin | 4eceeba | 2007-07-23 06:48:26 +0000 | [diff] [blame] | 6 | #include <REAL.h> |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 7 | |
Lev Walkin | ed46543 | 2004-09-27 20:52:36 +0000 | [diff] [blame] | 8 | static char reconstructed[2][512]; |
| 9 | static int reconstr_lens[2]; |
| 10 | |
| 11 | static int |
| 12 | callback(const void *buffer, size_t size, void *app_key) { |
| 13 | char *buf = reconstructed[app_key ? 1 : 0]; |
| 14 | int *len = &reconstr_lens[app_key ? 1 : 0]; |
| 15 | |
| 16 | if(*len + size >= sizeof(reconstructed[0])) |
| 17 | return -1; |
| 18 | |
| 19 | memcpy(buf + *len, buffer, size); |
| 20 | *len += size; |
| 21 | |
| 22 | return 0; |
| 23 | } |
| 24 | |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 25 | static void |
Lev Walkin | ed46543 | 2004-09-27 20:52:36 +0000 | [diff] [blame] | 26 | check_str_repr(double d, const char *sample, const char *canonical_sample) { |
| 27 | ssize_t s1, s2; |
| 28 | |
| 29 | reconstr_lens[1] = reconstr_lens[0] = 0; |
| 30 | |
| 31 | s1 = REAL__dump(d, 0, callback, 0); |
| 32 | assert(s1 < sizeof(reconstructed[0])); |
| 33 | assert(s1 == reconstr_lens[0]); |
| 34 | reconstructed[0][s1] = '\0'; |
| 35 | |
| 36 | s2 = REAL__dump(d, 1, callback, (void *)1); |
| 37 | assert(s2 < sizeof(reconstructed[1])); |
| 38 | assert(s2 == reconstr_lens[1]); |
| 39 | reconstructed[1][s2] = '\0'; |
| 40 | |
| 41 | if(sample) { |
Lev Walkin | f0b808d | 2005-04-25 21:08:25 +0000 | [diff] [blame] | 42 | printf("Checking %f->[%s] against [%s]%s\n", |
| 43 | d, reconstructed[0], sample, |
| 44 | canonical_sample ? " (canonical follows...)" : "" |
| 45 | ); |
Lev Walkin | ed46543 | 2004-09-27 20:52:36 +0000 | [diff] [blame] | 46 | assert(!strcmp(reconstructed[0], sample)); |
| 47 | } |
| 48 | if(canonical_sample) { |
Lev Walkin | f0b808d | 2005-04-25 21:08:25 +0000 | [diff] [blame] | 49 | printf("Checking %f->[%s] against [%s] (canonical)\n", |
| 50 | d, reconstructed[1], canonical_sample); |
Lev Walkin | ed46543 | 2004-09-27 20:52:36 +0000 | [diff] [blame] | 51 | assert(!strcmp(reconstructed[1], canonical_sample)); |
| 52 | } |
| 53 | } |
| 54 | |
| 55 | static void |
Lev Walkin | f0b808d | 2005-04-25 21:08:25 +0000 | [diff] [blame] | 56 | check_impl(REAL_t *rn, double orig_dbl, const char *sample, const char *canonical_sample, int line) { |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 57 | double val; |
| 58 | uint8_t *p, *end; |
| 59 | int ret; |
| 60 | |
Lev Walkin | f0b808d | 2005-04-25 21:08:25 +0000 | [diff] [blame] | 61 | printf("Line %d: double value %.12f [", line, orig_dbl); |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 62 | for(p = (uint8_t *)&orig_dbl, end = p + sizeof(double); p < end ; p++) |
| 63 | printf("%02x", *p); |
| 64 | printf("] (ilogb %d)\n", ilogb(orig_dbl)); |
| 65 | |
| 66 | val = frexp(orig_dbl, &ret); |
| 67 | printf("frexp(%f, %d): [", val, ret); |
| 68 | for(p = (uint8_t *)&val, end = p + sizeof(double); p < end ; p++) |
| 69 | printf("%02x", *p); |
| 70 | printf("]\n"); |
| 71 | |
Lev Walkin | 5e03376 | 2004-09-29 13:26:15 +0000 | [diff] [blame] | 72 | ret = asn_double2REAL(rn, orig_dbl); |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 73 | assert(ret == 0); |
| 74 | |
| 75 | printf("converted into ["); |
| 76 | for(p = rn->buf, end = p + rn->size; p < end; p++) |
| 77 | printf("%02x", *p); |
Lev Walkin | 2a789d9 | 2004-09-27 21:36:59 +0000 | [diff] [blame] | 78 | printf("]: %d\n", rn->size); |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 79 | |
Lev Walkin | 5e03376 | 2004-09-29 13:26:15 +0000 | [diff] [blame] | 80 | ret = asn_REAL2double(rn, &val); |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 81 | assert(ret == 0); |
| 82 | |
| 83 | printf("and back to double: ["); |
| 84 | for(p = (uint8_t *)&val, end = p + sizeof(double); p < end ; p++) |
| 85 | printf("%02x", *p); |
| 86 | printf("] (ilogb %d)\n", ilogb(val)); |
| 87 | |
Lev Walkin | 1aea698 | 2004-10-26 09:35:25 +0000 | [diff] [blame] | 88 | printf("%.12f vs %.12f\n", val, orig_dbl); |
| 89 | assert((isnan(orig_dbl) && isnan(val)) || val == orig_dbl); |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 90 | printf("OK\n"); |
Lev Walkin | ed46543 | 2004-09-27 20:52:36 +0000 | [diff] [blame] | 91 | |
| 92 | check_str_repr(val, sample, canonical_sample); |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 93 | } |
| 94 | |
| 95 | uint8_t buf_1_0[] = { 0x80, 0xcc, 0x10, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00 }; |
| 96 | uint8_t buf_1_1[] = { 0x80, 0xcc, 0x11, 0x99, 0x99, 0x99, 0x99, 0x99, 0x9a }; |
| 97 | uint8_t buf_3_14[] = { 0x80, 0xcd, 0x19, 0x1e, 0xb8, 0x51, 0xeb, 0x85, 0x1f }; |
Lev Walkin | 1aea698 | 2004-10-26 09:35:25 +0000 | [diff] [blame] | 98 | /* These ones are very interesting! They check mantissa overflow! */ |
Lev Walkin | ed46543 | 2004-09-27 20:52:36 +0000 | [diff] [blame] | 99 | uint8_t buf_mo1[] = { 0xC0, 0xc5, 0x19, 0x1e, 0xb8, 0x51, 0xeb, 0x85, 0x1f,3}; |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 100 | uint8_t buf_mo2[] = { 0x80, 0xbd, 0x19, 0x1e, 0xb8, 0x51, 0xeb, 0x85, 0x1f,3,2}; |
| 101 | |
| 102 | static void |
Lev Walkin | ed46543 | 2004-09-27 20:52:36 +0000 | [diff] [blame] | 103 | check_buf(uint8_t *buf, size_t bufsize, double verify, const char *sample, const char *canonical_sample) { |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 104 | REAL_t rn; |
| 105 | double val; |
| 106 | uint8_t *p, *end; |
| 107 | int ret; |
| 108 | |
| 109 | printf("verify double value %.12f [", verify); |
| 110 | for(p = (uint8_t *)&verify, end = p + sizeof(double); p < end ; p++) |
| 111 | printf("%02x", *p); |
| 112 | printf("] (ilogb %d)\n", ilogb(verify)); |
| 113 | |
| 114 | rn.buf = 0; |
| 115 | rn.size = 0; |
| 116 | |
Lev Walkin | 5e03376 | 2004-09-29 13:26:15 +0000 | [diff] [blame] | 117 | ret = asn_double2REAL(&rn, verify); |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 118 | assert(ret == 0); |
| 119 | |
| 120 | printf("canonical DER: ["); |
| 121 | for(p = rn.buf, end = p + rn.size; p < end; p++) |
| 122 | printf("%02x", *p); |
| 123 | printf("]\n"); |
| 124 | |
| 125 | rn.buf = buf; |
| 126 | rn.size = bufsize; |
| 127 | |
| 128 | printf("received as: ["); |
| 129 | for(p = rn.buf, end = p + rn.size; p < end; p++) |
| 130 | printf("%02x", *p); |
| 131 | printf("]\n"); |
| 132 | |
Lev Walkin | 5e03376 | 2004-09-29 13:26:15 +0000 | [diff] [blame] | 133 | ret = asn_REAL2double(&rn, &val); |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 134 | assert(ret == 0); |
| 135 | |
| 136 | printf("%.12f vs %.12f\n", verify, val); |
| 137 | |
| 138 | assert(val == verify); |
Lev Walkin | ed46543 | 2004-09-27 20:52:36 +0000 | [diff] [blame] | 139 | |
| 140 | check_str_repr(val, sample, canonical_sample); |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 141 | } |
| 142 | |
Lev Walkin | 5f56091 | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 143 | static void |
| 144 | check_xer(int fuzzy, double orig_value) { |
| 145 | asn_enc_rval_t er; |
| 146 | asn_dec_rval_t rc; |
| 147 | REAL_t st; |
| 148 | REAL_t *newst0 = 0; |
| 149 | REAL_t *newst1 = 0; |
Lev Walkin | b191938 | 2006-07-27 11:46:25 +0000 | [diff] [blame] | 150 | REAL_t **newst0p = &newst0; |
| 151 | REAL_t **newst1p = &newst1; |
Lev Walkin | 5f56091 | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 152 | double value0, value1; |
| 153 | int ret; |
| 154 | |
| 155 | memset(&st, 0, sizeof(st)); |
| 156 | ret = asn_double2REAL(&st, orig_value); |
| 157 | assert(ret == 0); |
| 158 | |
| 159 | reconstr_lens[0] = 0; |
| 160 | reconstr_lens[1] = 0; |
| 161 | er = xer_encode(&asn_DEF_REAL, &st, |
| 162 | XER_F_BASIC, callback, 0); |
| 163 | assert(er.encoded == reconstr_lens[0]); |
| 164 | er = xer_encode(&asn_DEF_REAL, &st, |
| 165 | XER_F_CANONICAL, callback, (void *)1); |
| 166 | assert(er.encoded == reconstr_lens[1]); |
| 167 | reconstructed[0][reconstr_lens[0]] = 0; |
| 168 | reconstructed[1][reconstr_lens[1]] = 0; |
| 169 | |
| 170 | printf("%f vs (%d)[%s] & (%d)%s", |
| 171 | orig_value, |
| 172 | reconstr_lens[1], reconstructed[1], |
| 173 | reconstr_lens[0], reconstructed[0] |
| 174 | ); |
| 175 | |
Lev Walkin | b191938 | 2006-07-27 11:46:25 +0000 | [diff] [blame] | 176 | rc = xer_decode(0, &asn_DEF_REAL, (void **)newst0p, |
Lev Walkin | 5f56091 | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 177 | reconstructed[0], reconstr_lens[0]); |
| 178 | assert(rc.code == RC_OK); |
| 179 | assert(rc.consumed < reconstr_lens[0]); |
| 180 | |
Lev Walkin | b191938 | 2006-07-27 11:46:25 +0000 | [diff] [blame] | 181 | rc = xer_decode(0, &asn_DEF_REAL, (void **)newst1p, |
Lev Walkin | 5f56091 | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 182 | reconstructed[1], reconstr_lens[1]); |
| 183 | assert(rc.code == RC_OK); |
| 184 | assert(rc.consumed == reconstr_lens[1]); |
| 185 | |
| 186 | ret = asn_REAL2double(newst0, &value0); |
| 187 | assert(ret == 0); |
| 188 | ret = asn_REAL2double(newst1, &value1); |
| 189 | assert(ret == 0); |
| 190 | |
Lev Walkin | 1aea698 | 2004-10-26 09:35:25 +0000 | [diff] [blame] | 191 | assert((isnan(value0) && isnan(orig_value)) |
| 192 | || value0 == orig_value |
Lev Walkin | 5f56091 | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 193 | || fuzzy); |
Lev Walkin | 1aea698 | 2004-10-26 09:35:25 +0000 | [diff] [blame] | 194 | assert((isnan(value1) && isnan(orig_value)) |
| 195 | || value1 == orig_value); |
Lev Walkin | 5f56091 | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 196 | |
| 197 | assert(newst0->size == st.size || fuzzy); |
| 198 | assert(newst1->size == st.size); |
| 199 | assert(fuzzy || memcmp(newst0->buf, st.buf, st.size) == 0); |
| 200 | assert(memcmp(newst1->buf, st.buf, st.size) == 0); |
| 201 | } |
| 202 | |
Lev Walkin | f0b808d | 2005-04-25 21:08:25 +0000 | [diff] [blame] | 203 | #define check(rn, d, str1, str2) \ |
| 204 | check_impl(rn, d, str1, str2, __LINE__) |
| 205 | |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 206 | int |
| 207 | main() { |
| 208 | REAL_t rn; |
Lev Walkin | c51e7d6 | 2004-09-27 22:16:18 +0000 | [diff] [blame] | 209 | static const double zero = 0.0; |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 210 | |
| 211 | memset(&rn, 0, sizeof(rn)); |
| 212 | |
Lev Walkin | ed46543 | 2004-09-27 20:52:36 +0000 | [diff] [blame] | 213 | check(&rn, 0.0, "0", "0"); |
| 214 | check(&rn, -0.0, "-0", "-0"); /* minus-zero */ |
Lev Walkin | c51e7d6 | 2004-09-27 22:16:18 +0000 | [diff] [blame] | 215 | check(&rn, zero/zero, "<NOT-A-NUMBER/>", "<NOT-A-NUMBER/>"); |
| 216 | check(&rn, 1.0/zero, "<PLUS-INFINITY/>", "<PLUS-INFINITY/>"); |
| 217 | check(&rn, -1.0/zero, "<MINUS-INFINITY/>", "<MINUS-INFINITY/>"); |
Lev Walkin | ed46543 | 2004-09-27 20:52:36 +0000 | [diff] [blame] | 218 | check(&rn, 1.0, "1.0", "1.0E0"); |
| 219 | check(&rn, -1.0, "-1.0", "-1.0E0"); |
Lev Walkin | ed46543 | 2004-09-27 20:52:36 +0000 | [diff] [blame] | 220 | check(&rn, 0.1, "0.1", "1.0E-1"); |
Lev Walkin | f0b808d | 2005-04-25 21:08:25 +0000 | [diff] [blame] | 221 | check(&rn, 0.01, "0.01", "1.0E-2"); |
| 222 | check(&rn, 0.02, "0.02", "2.0E-2"); |
| 223 | check(&rn, 0.09, "0.09", "9.0E-2"); |
| 224 | check(&rn, 1.5, "1.5", "1.5E0"); |
Lev Walkin | ed46543 | 2004-09-27 20:52:36 +0000 | [diff] [blame] | 225 | check(&rn, 0.33333, "0.33333", "3.3333E-1"); |
| 226 | check(&rn, 2, "2.0", "2.0E0"); |
| 227 | check(&rn, 2.1, "2.1", "2.1E0"); |
| 228 | check(&rn, 3, "3.0", "3.0E0"); |
| 229 | check(&rn, 3.1, "3.1", "3.1E0"); |
| 230 | check(&rn, 3.14, "3.14", "3.14E0"); |
| 231 | check(&rn, 3.1415, "3.1415", "3.1415E0"); |
| 232 | check(&rn, 3.141592, "3.141592", "3.141592E0"); |
| 233 | check(&rn, 3.14159265, "3.14159265", "3.14159265E0"); |
| 234 | check(&rn, -3.14159265, "-3.14159265", "-3.14159265E0"); |
| 235 | check(&rn, 14159265.0, "14159265.0", "1.4159265E7"); |
| 236 | check(&rn, -123456789123456789.0, "-123456789123456784.0", "-1.234567891234568E17"); |
Lev Walkin | f0b808d | 2005-04-25 21:08:25 +0000 | [diff] [blame] | 237 | check(&rn, 0.00000000001, "0.00000000001", "9.999999999999999E-12"); |
| 238 | check(&rn, 0.00000000002, "0.00000000002", "2.0E-11"); |
| 239 | check(&rn, 0.00000000009, "0.00000000009", "9.0E-11"); |
| 240 | check(&rn, 0.000000000002, "0.000000000002", "2.0E-12"); |
| 241 | check(&rn, 0.0000000000002, "0.0000000000002", "2.0E-13"); |
| 242 | check(&rn, 0.00000000000002, "0.00000000000002", "2.0E-14"); |
| 243 | check(&rn, 0.000000000000002, "0.000000000000002", "2.0E-15"); |
| 244 | check(&rn, 0.0000000000000002, "0.0", "2.0E-16"); |
Lev Walkin | ed46543 | 2004-09-27 20:52:36 +0000 | [diff] [blame] | 245 | check(&rn, 0.0000000000000000000001, "0.0", "1.0E-22"); |
| 246 | check(&rn, 0.000000000000000000000000000001, "0.0", "1.0E-30"); /* proved 2B a problem */ |
| 247 | check(&rn,-0.000000000000000000000000000001, "-0.0", "-1.0E-30"); /* proved 2B a problem */ |
| 248 | check(&rn, 0.0000000000010000000001000000000001, 0, 0); |
| 249 | check(&rn, 0.00000000000000000000000000000000000000000000000000000000000000000000000000000000000000000001, 0, 0); |
| 250 | check(&rn, 0.000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000001, 0, 0); |
| 251 | check(&rn,-0.000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000001, 0, 0); |
| 252 | check(&rn,-3.33333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333, 0, 0); |
| 253 | check(&rn, 0.0000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000033333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333, 0, 0); |
| 254 | check(&rn, -0.00000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000001, 0, 0); |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 255 | |
Lev Walkin | ed46543 | 2004-09-27 20:52:36 +0000 | [diff] [blame] | 256 | check_buf(buf_1_0, sizeof(buf_1_0), 1.0, "1.0", "1.0E0"); |
| 257 | check_buf(buf_1_1, sizeof(buf_1_1), 1.1, "1.1", "1.1E0"); |
| 258 | check_buf(buf_3_14, sizeof(buf_3_14), 3.14, "3.14", "3.14E0"); |
| 259 | check_buf(buf_mo1, sizeof(buf_mo1), -3.14, "-3.14", "-3.14E0"); |
| 260 | check_buf(buf_mo2, sizeof(buf_mo2), 3.14, "3.14", "3.14E0"); |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 261 | |
Lev Walkin | 5f56091 | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 262 | |
Lev Walkin | 1aea698 | 2004-10-26 09:35:25 +0000 | [diff] [blame] | 263 | #ifdef NAN |
| 264 | check_xer(0, NAN); /* "<NOT-A-NUMBER/>" */ |
| 265 | #else |
Lev Walkin | 5f56091 | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 266 | check_xer(0, zero/zero); /* "<NOT-A-NUMBER/>" */ |
Lev Walkin | 1aea698 | 2004-10-26 09:35:25 +0000 | [diff] [blame] | 267 | #endif |
| 268 | #ifdef INFINITY |
| 269 | check_xer(0, INFINITY); /* "<PLUS-INFINITY/>" */ |
| 270 | check_xer(0, -INFINITY); /* "<MINUS-INFINITY/>" */ |
| 271 | #else |
Lev Walkin | 5f56091 | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 272 | check_xer(0, 1.0/zero); /* "<PLUS-INFINITY/>" */ |
| 273 | check_xer(0, -1.0/zero); /* "<MINUS-INFINITY/>" */ |
Lev Walkin | 1aea698 | 2004-10-26 09:35:25 +0000 | [diff] [blame] | 274 | #endif |
Lev Walkin | 5f56091 | 2004-10-21 13:37:57 +0000 | [diff] [blame] | 275 | check_xer(0, 1.0); |
| 276 | check_xer(0, -1.0); |
| 277 | check_xer(0, 1.5); |
| 278 | check_xer(0, 123); |
| 279 | check_xer(1, 0.0000000000000000000001); |
| 280 | check_xer(1, -0.0000000000000000000001); |
| 281 | |
Lev Walkin | 41ba1f2 | 2004-09-14 12:46:35 +0000 | [diff] [blame] | 282 | return 0; |
| 283 | } |