| 1 | #include "slalib.h"
|
|---|
| 2 | #include "slamac.h"
|
|---|
| 3 | #include <string.h>
|
|---|
| 4 | #include <limits.h>
|
|---|
| 5 |
|
|---|
| 6 | static int idchi ( int, char*, int*, double* );
|
|---|
| 7 |
|
|---|
| 8 | void slaIntin ( char *string, int *nstrt, long *ireslt, int *jflag )
|
|---|
| 9 | /*
|
|---|
| 10 | ** - - - - - - - - -
|
|---|
| 11 | ** s l a I n t i n
|
|---|
| 12 | ** - - - - - - - - -
|
|---|
| 13 | **
|
|---|
| 14 | ** Convert free-format input into a long integer.
|
|---|
| 15 | **
|
|---|
| 16 | ** Given:
|
|---|
| 17 | ** string char* string containing number to be decoded
|
|---|
| 18 | ** nstrt int* where to start decode (1st = 1)
|
|---|
| 19 | ** ireslt long* current value of result
|
|---|
| 20 | **
|
|---|
| 21 | ** Returned:
|
|---|
| 22 | ** nstrt int* advanced to next number
|
|---|
| 23 | ** ireslt long* result
|
|---|
| 24 | ** jflag int* status: -1 = -OK, 0 = +OK, 1 = null, 2 = error
|
|---|
| 25 | **
|
|---|
| 26 | ** Called: idchi
|
|---|
| 27 | **
|
|---|
| 28 | ** Notes:
|
|---|
| 29 | **
|
|---|
| 30 | ** 1 The reason slaIntin has separate OK status values for +
|
|---|
| 31 | ** and - is to enable minus zero to be detected. This is
|
|---|
| 32 | ** of crucial importance when decoding mixed-radix numbers.
|
|---|
| 33 | ** For example, an angle expressed as deg, arcmin, arcsec
|
|---|
| 34 | ** may have a leading minus sign but a zero degrees field.
|
|---|
| 35 | **
|
|---|
| 36 | ** 2 A TAB is interpreted as a space.
|
|---|
| 37 | **
|
|---|
| 38 | ** 3 The basic format is the sequence of fields #^, where
|
|---|
| 39 | ** # is a sign character + or -, and ^ means a string of
|
|---|
| 40 | ** decimal digits.
|
|---|
| 41 | **
|
|---|
| 42 | ** 4 Spaces:
|
|---|
| 43 | **
|
|---|
| 44 | ** . Leading spaces are ignored.
|
|---|
| 45 | **
|
|---|
| 46 | ** . Spaces between the sign and the number are allowed.
|
|---|
| 47 | **
|
|---|
| 48 | ** . Trailing spaces are ignored; the first signifies
|
|---|
| 49 | ** end of decoding and subsequent ones are skipped.
|
|---|
| 50 | **
|
|---|
| 51 | ** 5 Delimiters:
|
|---|
| 52 | **
|
|---|
| 53 | ** . Any character other than +,-,0-9 or space may be
|
|---|
| 54 | ** used to signal the end of the number and terminate
|
|---|
| 55 | ** decoding.
|
|---|
| 56 | **
|
|---|
| 57 | ** . Comma is recognized by slaIntin as a special case; it
|
|---|
| 58 | ** is skipped, leaving the pointer on the next character.
|
|---|
| 59 | ** See 9, below.
|
|---|
| 60 | **
|
|---|
| 61 | ** 6 The sign is optional. The default is +.
|
|---|
| 62 | **
|
|---|
| 63 | ** 7 A "null result" occurs when the string of characters being
|
|---|
| 64 | ** decoded does not begin with +,- or 0-9, or consists
|
|---|
| 65 | ** entirely of spaces. When this condition is detected, jflag
|
|---|
| 66 | ** is set to 1 and ireslt is left untouched.
|
|---|
| 67 | **
|
|---|
| 68 | ** 8 nstrt = 1 for the first character in the string.
|
|---|
| 69 | **
|
|---|
| 70 | ** 9 On return from slaIntin, nstrt is set ready for the next
|
|---|
| 71 | ** decode - following trailing blanks and any comma. If a
|
|---|
| 72 | ** delimiter other than comma is being used, nstrt must be
|
|---|
| 73 | ** incremented before the next call to slaIntin, otherwise
|
|---|
| 74 | ** all subsequent calls will return a null result.
|
|---|
| 75 | **
|
|---|
| 76 | ** 10 Errors (jflag=2) occur when:
|
|---|
| 77 | **
|
|---|
| 78 | ** . there is a + or - but no number; or
|
|---|
| 79 | **
|
|---|
| 80 | ** . the number is larger than LONG_MAX.
|
|---|
| 81 | **
|
|---|
| 82 | ** 11 When an error has been detected, nstrt is left
|
|---|
| 83 | ** pointing to the character following the last
|
|---|
| 84 | ** one used before the error came to light.
|
|---|
| 85 | **
|
|---|
| 86 | ** 12 See also slaFlotin and slaDfltin.
|
|---|
| 87 | **
|
|---|
| 88 | ** Last revision: 6 November 1999
|
|---|
| 89 | **
|
|---|
| 90 | ** Copyright P.T.Wallace. All rights reserved.
|
|---|
| 91 | */
|
|---|
| 92 |
|
|---|
| 93 | /* Definitions shared between slaIntin and idchi */
|
|---|
| 94 | #define NUMBER 0
|
|---|
| 95 | #define SPACE 1
|
|---|
| 96 | #define PLUS 2
|
|---|
| 97 | #define MINUS 3
|
|---|
| 98 | #define COMMA 4
|
|---|
| 99 | #define OTHER 5
|
|---|
| 100 | #define END 6
|
|---|
| 101 |
|
|---|
| 102 | {
|
|---|
| 103 | int l_string, nptr;
|
|---|
| 104 | double digit;
|
|---|
| 105 |
|
|---|
| 106 | /* Current state of the decode and the values it can take */
|
|---|
| 107 |
|
|---|
| 108 | int state;
|
|---|
| 109 |
|
|---|
| 110 | #define seek_sign 100
|
|---|
| 111 | #define neg 200
|
|---|
| 112 | #define seek_1st_digit 300
|
|---|
| 113 | #define accept_digit 400
|
|---|
| 114 | #define seek_digit 410
|
|---|
| 115 | #define end_of_field 1600
|
|---|
| 116 | #define build_result 1610
|
|---|
| 117 | #define seeking_end_of_field 1630
|
|---|
| 118 | #define next_field_OK 1720
|
|---|
| 119 | #define next_field_default 9100
|
|---|
| 120 | #define null_field 9110
|
|---|
| 121 | #define next_field_error 9200
|
|---|
| 122 | #define error 9210
|
|---|
| 123 | #define done 9900
|
|---|
| 124 |
|
|---|
| 125 |
|
|---|
| 126 | int j;
|
|---|
| 127 | double dres;
|
|---|
| 128 |
|
|---|
| 129 |
|
|---|
| 130 | /* Find string length */
|
|---|
| 131 | l_string = strlen ( string );
|
|---|
| 132 |
|
|---|
| 133 | /* Current character index (1st = 0) */
|
|---|
| 134 | nptr = *nstrt - 1;
|
|---|
| 135 |
|
|---|
| 136 | /* Set defaults: result & sign */
|
|---|
| 137 | dres = 0.0;
|
|---|
| 138 | j = 0;
|
|---|
| 139 |
|
|---|
| 140 | /* Initialize state to "looking for sign" */
|
|---|
| 141 | state = seek_sign;
|
|---|
| 142 |
|
|---|
| 143 | /* Loop until decode is complete */
|
|---|
| 144 | while ( state != done ) {
|
|---|
| 145 | switch ( state ) {
|
|---|
| 146 |
|
|---|
| 147 | case seek_sign :
|
|---|
| 148 |
|
|---|
| 149 | /* Look for sign */
|
|---|
| 150 | switch ( idchi ( l_string, string, &nptr, &digit ) ) {
|
|---|
| 151 | case NUMBER :
|
|---|
| 152 | state = accept_digit;
|
|---|
| 153 | break;
|
|---|
| 154 | case SPACE :
|
|---|
| 155 | state = seek_sign;
|
|---|
| 156 | break;
|
|---|
| 157 | case PLUS :
|
|---|
| 158 | state = seek_1st_digit;
|
|---|
| 159 | break;
|
|---|
| 160 | case MINUS :
|
|---|
| 161 | state = neg;
|
|---|
| 162 | break;
|
|---|
| 163 | case OTHER :
|
|---|
| 164 | state = next_field_default;
|
|---|
| 165 | break;
|
|---|
| 166 | case COMMA :
|
|---|
| 167 | case END :
|
|---|
| 168 | state = null_field;
|
|---|
| 169 | break;
|
|---|
| 170 | default :
|
|---|
| 171 | state = error;
|
|---|
| 172 | }
|
|---|
| 173 | break;
|
|---|
| 174 |
|
|---|
| 175 | case neg :
|
|---|
| 176 |
|
|---|
| 177 | /* Negative result */
|
|---|
| 178 | j = -1;
|
|---|
| 179 |
|
|---|
| 180 | case seek_1st_digit :
|
|---|
| 181 |
|
|---|
| 182 | /* Look for first leading decimal */
|
|---|
| 183 | switch ( idchi ( l_string, string, &nptr, &digit ) ) {
|
|---|
| 184 | case NUMBER :
|
|---|
| 185 | state = accept_digit;
|
|---|
| 186 | break;
|
|---|
| 187 | case SPACE :
|
|---|
| 188 | state = seek_1st_digit;
|
|---|
| 189 | break;
|
|---|
| 190 | case PLUS :
|
|---|
| 191 | case MINUS :
|
|---|
| 192 | case COMMA :
|
|---|
| 193 | case OTHER :
|
|---|
| 194 | state = next_field_error;
|
|---|
| 195 | break;
|
|---|
| 196 | case END :
|
|---|
| 197 | default :
|
|---|
| 198 | state = error;
|
|---|
| 199 | }
|
|---|
| 200 | break;
|
|---|
| 201 |
|
|---|
| 202 | case accept_digit :
|
|---|
| 203 |
|
|---|
| 204 | /* Accept decimals */
|
|---|
| 205 | dres = dres * 1e1 + digit;
|
|---|
| 206 | state = ( fabs ( dres ) <= LONG_MAX) ?
|
|---|
| 207 | seek_digit : next_field_error;
|
|---|
| 208 | break;
|
|---|
| 209 |
|
|---|
| 210 | case seek_digit :
|
|---|
| 211 |
|
|---|
| 212 | /* Look for next decimal */
|
|---|
| 213 | switch ( idchi ( l_string, string, &nptr, &digit ) ) {
|
|---|
| 214 | case NUMBER :
|
|---|
| 215 | state = accept_digit;
|
|---|
| 216 | break;
|
|---|
| 217 | case SPACE :
|
|---|
| 218 | state = build_result;
|
|---|
| 219 | break;
|
|---|
| 220 | case PLUS :
|
|---|
| 221 | case MINUS :
|
|---|
| 222 | case COMMA :
|
|---|
| 223 | case OTHER :
|
|---|
| 224 | state = end_of_field;
|
|---|
| 225 | break;
|
|---|
| 226 | case END :
|
|---|
| 227 | state = build_result;
|
|---|
| 228 | break;
|
|---|
| 229 | default :
|
|---|
| 230 | state = error;
|
|---|
| 231 | }
|
|---|
| 232 | break;
|
|---|
| 233 |
|
|---|
| 234 | case end_of_field :
|
|---|
| 235 |
|
|---|
| 236 | /* Off the end of the field: move pointer back */
|
|---|
| 237 | nptr--;
|
|---|
| 238 |
|
|---|
| 239 | case build_result :
|
|---|
| 240 |
|
|---|
| 241 | /* Make the result */
|
|---|
| 242 | if ( j ) dres = - dres;
|
|---|
| 243 | *ireslt = (long) ( dnint ( dres ) );
|
|---|
| 244 |
|
|---|
| 245 | case seeking_end_of_field :
|
|---|
| 246 |
|
|---|
| 247 | /* Skip to end of field */
|
|---|
| 248 | switch ( idchi ( l_string, string, &nptr, &digit ) ) {
|
|---|
| 249 | case SPACE :
|
|---|
| 250 | state = seeking_end_of_field;
|
|---|
| 251 | break;
|
|---|
| 252 | case NUMBER :
|
|---|
| 253 | case PLUS :
|
|---|
| 254 | case MINUS :
|
|---|
| 255 | case OTHER :
|
|---|
| 256 | state = next_field_OK;
|
|---|
| 257 | break;
|
|---|
| 258 | case COMMA :
|
|---|
| 259 | case END :
|
|---|
| 260 | state = done;
|
|---|
| 261 | break;
|
|---|
| 262 | default :
|
|---|
| 263 | state = error;
|
|---|
| 264 | }
|
|---|
| 265 | break;
|
|---|
| 266 |
|
|---|
| 267 | case next_field_OK :
|
|---|
| 268 |
|
|---|
| 269 | /* Next field terminates successful decode */
|
|---|
| 270 | nptr--;
|
|---|
| 271 | state = done;
|
|---|
| 272 | break;
|
|---|
| 273 |
|
|---|
| 274 | case next_field_default :
|
|---|
| 275 |
|
|---|
| 276 | /* Next field terminates null decode */
|
|---|
| 277 | nptr--;
|
|---|
| 278 |
|
|---|
| 279 | case null_field :
|
|---|
| 280 |
|
|---|
| 281 | /* Null decode */
|
|---|
| 282 | j = 1;
|
|---|
| 283 | state = done;
|
|---|
| 284 | break;
|
|---|
| 285 |
|
|---|
| 286 | case next_field_error :
|
|---|
| 287 |
|
|---|
| 288 | /* Next field detected prematurely */
|
|---|
| 289 | nptr--;
|
|---|
| 290 |
|
|---|
| 291 | case error :
|
|---|
| 292 |
|
|---|
| 293 | /* Decode has failed: set bad status */
|
|---|
| 294 | j = 2;
|
|---|
| 295 | state = done;
|
|---|
| 296 | break;
|
|---|
| 297 |
|
|---|
| 298 | default :
|
|---|
| 299 | state = error;
|
|---|
| 300 | }
|
|---|
| 301 | }
|
|---|
| 302 |
|
|---|
| 303 | /* Finished: return updated pointer and the status */
|
|---|
| 304 | *nstrt = nptr + 1;
|
|---|
| 305 | *jflag = j;
|
|---|
| 306 | }
|
|---|
| 307 |
|
|---|
| 308 | static int idchi ( int l_string, char *string, int *nptr, double *digit )
|
|---|
| 309 | /*
|
|---|
| 310 | ** - - - - -
|
|---|
| 311 | ** i d c h i
|
|---|
| 312 | ** - - - - -
|
|---|
| 313 | **
|
|---|
| 314 | ** Internal routine used by slaIntin:
|
|---|
| 315 | **
|
|---|
| 316 | ** identify next character in string.
|
|---|
| 317 | **
|
|---|
| 318 | ** Given:
|
|---|
| 319 | ** l_string int length of string
|
|---|
| 320 | ** string char* string
|
|---|
| 321 | ** nptr int* character to be identified (1st = 0)
|
|---|
| 322 | **
|
|---|
| 323 | ** Returned:
|
|---|
| 324 | ** nptr int* incremented unless end of field
|
|---|
| 325 | ** digit double* 0.0 - 9.0 if character was a numeral
|
|---|
| 326 | **
|
|---|
| 327 | ** Returned (function value):
|
|---|
| 328 | ** int vector for identified character:
|
|---|
| 329 | **
|
|---|
| 330 | ** value meaning
|
|---|
| 331 | **
|
|---|
| 332 | ** NUMBER 0-9
|
|---|
| 333 | ** SPACE space or tab
|
|---|
| 334 | ** PLUS +
|
|---|
| 335 | ** MINUS -
|
|---|
| 336 | ** COMMA ,
|
|---|
| 337 | ** OTHER else
|
|---|
| 338 | ** END outside field
|
|---|
| 339 | **
|
|---|
| 340 | ** Last revision: 24 June 1996
|
|---|
| 341 | **
|
|---|
| 342 | ** Copyright P.T.Wallace. All rights reserved.
|
|---|
| 343 | */
|
|---|
| 344 | {
|
|---|
| 345 | int ivec, ictab;
|
|---|
| 346 | char c;
|
|---|
| 347 |
|
|---|
| 348 | /* Character/vector tables */
|
|---|
| 349 |
|
|---|
| 350 | #define NCREC (15)
|
|---|
| 351 | static char kctab[NCREC] = { '0','1','2','3','4','5',
|
|---|
| 352 | '6','7','8','9',
|
|---|
| 353 | ' ','\t',
|
|---|
| 354 | '+',
|
|---|
| 355 | '-',
|
|---|
| 356 | ',' };
|
|---|
| 357 |
|
|---|
| 358 | static int kvtab[NCREC] = { NUMBER, NUMBER, NUMBER, NUMBER, NUMBER,
|
|---|
| 359 | NUMBER, NUMBER, NUMBER, NUMBER, NUMBER,
|
|---|
| 360 | SPACE, SPACE,
|
|---|
| 361 | PLUS,
|
|---|
| 362 | MINUS,
|
|---|
| 363 | COMMA };
|
|---|
| 364 |
|
|---|
| 365 |
|
|---|
| 366 | /* Initialize returned value */
|
|---|
| 367 | ivec = OTHER;
|
|---|
| 368 |
|
|---|
| 369 | /* Pointer outside field? */
|
|---|
| 370 | if ( *nptr < 0 || *nptr >= l_string ) {
|
|---|
| 371 |
|
|---|
| 372 | /* Yes: prepare returned value */
|
|---|
| 373 | ivec = END;
|
|---|
| 374 |
|
|---|
| 375 | } else {
|
|---|
| 376 |
|
|---|
| 377 | /* Not end of field: identify character */
|
|---|
| 378 | c = string [ *nptr ];
|
|---|
| 379 | for ( ictab = 0; ictab < NCREC; ictab++ ) {
|
|---|
| 380 | if ( kctab [ ictab ] == c ) {
|
|---|
| 381 |
|
|---|
| 382 | /* Recognized */
|
|---|
| 383 | ivec = kvtab [ ictab ];
|
|---|
| 384 |
|
|---|
| 385 | /* Allow for numerals */
|
|---|
| 386 | *digit = (double) ictab;
|
|---|
| 387 |
|
|---|
| 388 | /* Quit the loop */
|
|---|
| 389 | break;
|
|---|
| 390 | }
|
|---|
| 391 | }
|
|---|
| 392 |
|
|---|
| 393 | /* Increment pointer */
|
|---|
| 394 | ( *nptr )++;
|
|---|
| 395 | }
|
|---|
| 396 |
|
|---|
| 397 | /* Return the value identifying the character */
|
|---|
| 398 | return ivec;
|
|---|
| 399 | }
|
|---|