|
|
|
@ -197,3 +197,233 @@ INPevaluate(char **line, int *error, int gobble) |
|
|
|
return (sign * mantis * |
|
|
|
pow(10.0, (double) (expo1 + expsgn * expo2))); |
|
|
|
} |
|
|
|
|
|
|
|
|
|
|
|
/* In addition to fcn INPevaluate() above, allow values like 4k7, |
|
|
|
similar to the RKM code (used bei inp2c and inp2r) */ |
|
|
|
double |
|
|
|
INPevaluateRKM(char** line, int* error, int gobble) |
|
|
|
/* gobble: non-zero to gobble rest of token, zero to leave it alone */ |
|
|
|
{ |
|
|
|
char* token; |
|
|
|
char* here; |
|
|
|
double mantis; |
|
|
|
double deci; |
|
|
|
int expo1; |
|
|
|
int expo2; |
|
|
|
int expo3; |
|
|
|
int sign; |
|
|
|
int expsgn; |
|
|
|
char* tmpline; |
|
|
|
bool hasmulti = FALSE; |
|
|
|
|
|
|
|
/* setup */ |
|
|
|
tmpline = *line; |
|
|
|
|
|
|
|
if (gobble) { |
|
|
|
/* MW. INPgetUTok should be called with gobble=0 or it make |
|
|
|
* errors in v(1,2) exp */ |
|
|
|
*error = INPgetUTok(line, &token, 0); |
|
|
|
if (*error) |
|
|
|
return (0.0); |
|
|
|
} |
|
|
|
else { |
|
|
|
token = *line; |
|
|
|
*error = 0; |
|
|
|
} |
|
|
|
|
|
|
|
mantis = 0; |
|
|
|
deci = 0; |
|
|
|
expo1 = 0; |
|
|
|
expo2 = 0; |
|
|
|
expo3 = 0; |
|
|
|
sign = 1; |
|
|
|
expsgn = 1; |
|
|
|
|
|
|
|
/* loop through all of the input token */ |
|
|
|
here = token; |
|
|
|
|
|
|
|
if (*here == '+') |
|
|
|
here++; /* plus, so do nothing except skip it */ |
|
|
|
else if (*here == '-') { /* minus, so skip it, and change sign */ |
|
|
|
here++; |
|
|
|
sign = -1; |
|
|
|
} |
|
|
|
|
|
|
|
if ((*here == '\0') || ((!(isdigit_c(*here))) && (*here != '.'))) { |
|
|
|
/* number looks like just a sign! */ |
|
|
|
*error = 1; |
|
|
|
if (gobble) { |
|
|
|
FREE(token); |
|
|
|
/* back out the 'gettok' operation */ |
|
|
|
*line = tmpline; |
|
|
|
} |
|
|
|
return (0); |
|
|
|
} |
|
|
|
|
|
|
|
while (isdigit_c(*here)) { |
|
|
|
/* digit, so accumulate it. */ |
|
|
|
mantis = 10 * mantis + *here - '0'; |
|
|
|
here++; |
|
|
|
} |
|
|
|
|
|
|
|
if (*here == '\0') { |
|
|
|
/* reached the end of token - done. */ |
|
|
|
if (gobble) { |
|
|
|
FREE(token); |
|
|
|
} |
|
|
|
else { |
|
|
|
*line = here; |
|
|
|
} |
|
|
|
return ((double)mantis * sign); |
|
|
|
} |
|
|
|
|
|
|
|
if (*here == ':') { |
|
|
|
/* ':' is no longer used for subcircuit node numbering |
|
|
|
but is part of ternary function a?b:c |
|
|
|
FIXME : subcircuit models still use ':' for model numbering |
|
|
|
Will this hurt somewhere? */ |
|
|
|
if (gobble) { |
|
|
|
FREE(token); |
|
|
|
} |
|
|
|
else { |
|
|
|
*line = here; |
|
|
|
} |
|
|
|
return ((double)mantis * sign); |
|
|
|
} |
|
|
|
|
|
|
|
/* after decimal point! */ |
|
|
|
if (*here == '.') { |
|
|
|
/* found a decimal point! */ |
|
|
|
here++; /* skip to next character */ |
|
|
|
|
|
|
|
if (*here == '\0') { |
|
|
|
/* number ends in the decimal point */ |
|
|
|
if (gobble) { |
|
|
|
FREE(token); |
|
|
|
} |
|
|
|
else { |
|
|
|
*line = here; |
|
|
|
} |
|
|
|
return ((double)mantis * sign); |
|
|
|
} |
|
|
|
|
|
|
|
while (isdigit_c(*here)) { |
|
|
|
/* digit, so accumulate it. */ |
|
|
|
mantis = 10 * mantis + *here - '0'; |
|
|
|
expo1 = expo1 - 1; |
|
|
|
here++; |
|
|
|
} |
|
|
|
} |
|
|
|
|
|
|
|
/* now look for "E","e",etc to indicate an exponent */ |
|
|
|
if ((*here == 'E') || (*here == 'e') || (*here == 'D') || (*here == 'd')) { |
|
|
|
|
|
|
|
/* have an exponent, so skip the e */ |
|
|
|
here++; |
|
|
|
|
|
|
|
/* now look for exponent sign */ |
|
|
|
if (*here == '+') |
|
|
|
here++; /* just skip + */ |
|
|
|
else if (*here == '-') { |
|
|
|
here++; /* skip over minus sign */ |
|
|
|
expsgn = -1; /* and make a negative exponent */ |
|
|
|
/* now look for the digits of the exponent */ |
|
|
|
} |
|
|
|
|
|
|
|
while (isdigit_c(*here)) { |
|
|
|
expo2 = 10 * expo2 + *here - '0'; |
|
|
|
here++; |
|
|
|
} |
|
|
|
} |
|
|
|
|
|
|
|
/* now we have all of the numeric part of the number, time to |
|
|
|
* look for the scale factor (alphabetic) |
|
|
|
*/ |
|
|
|
switch (*here) { |
|
|
|
case 't': |
|
|
|
case 'T': |
|
|
|
expo1 = expo1 + 12; |
|
|
|
hasmulti = TRUE; |
|
|
|
break; |
|
|
|
case 'g': |
|
|
|
case 'G': |
|
|
|
expo1 = expo1 + 9; |
|
|
|
hasmulti = TRUE; |
|
|
|
break; |
|
|
|
case 'k': |
|
|
|
case 'K': |
|
|
|
expo1 = expo1 + 3; |
|
|
|
hasmulti = TRUE; |
|
|
|
break; |
|
|
|
case 'u': |
|
|
|
case 'U': |
|
|
|
expo1 = expo1 - 6; |
|
|
|
hasmulti = TRUE; |
|
|
|
break; |
|
|
|
case 'r': |
|
|
|
case 'R': |
|
|
|
expo1 = expo1; |
|
|
|
hasmulti = TRUE; |
|
|
|
break; |
|
|
|
case 'n': |
|
|
|
case 'N': |
|
|
|
expo1 = expo1 - 9; |
|
|
|
hasmulti = TRUE; |
|
|
|
break; |
|
|
|
case 'p': |
|
|
|
case 'P': |
|
|
|
expo1 = expo1 - 12; |
|
|
|
hasmulti = TRUE; |
|
|
|
break; |
|
|
|
case 'f': |
|
|
|
case 'F': |
|
|
|
expo1 = expo1 - 15; |
|
|
|
hasmulti = TRUE; |
|
|
|
break; |
|
|
|
case 'm': |
|
|
|
case 'M': |
|
|
|
if (((here[1] == 'E') || (here[1] == 'e')) && |
|
|
|
((here[2] == 'G') || (here[2] == 'g'))) |
|
|
|
{ |
|
|
|
expo1 = expo1 + 6; /* Meg */ |
|
|
|
here += 2; |
|
|
|
hasmulti = TRUE; |
|
|
|
} |
|
|
|
else if (((here[1] == 'I') || (here[1] == 'i')) && |
|
|
|
((here[2] == 'L') || (here[2] == 'l'))) |
|
|
|
{ |
|
|
|
expo1 = expo1 - 6; |
|
|
|
mantis *= 25.4; /* Mil */ |
|
|
|
} |
|
|
|
else { |
|
|
|
expo1 = expo1 - 3; /* m, milli */ |
|
|
|
hasmulti = TRUE; |
|
|
|
} |
|
|
|
break; |
|
|
|
default: |
|
|
|
break; |
|
|
|
} |
|
|
|
|
|
|
|
/* read a digit after multiplier */ |
|
|
|
if (hasmulti) { |
|
|
|
here++; |
|
|
|
while (isdigit_c(*here)) { |
|
|
|
deci = 10 * deci + *here - '0'; |
|
|
|
expo3 = expo3 - 1; |
|
|
|
here++; |
|
|
|
} |
|
|
|
mantis = mantis + deci * pow(10.0, (double)expo3); |
|
|
|
} |
|
|
|
|
|
|
|
if (gobble) { |
|
|
|
FREE(token); |
|
|
|
} |
|
|
|
else { |
|
|
|
*line = here; |
|
|
|
} |
|
|
|
|
|
|
|
return (sign * mantis * |
|
|
|
pow(10.0, (double)(expo1 + expsgn * expo2))); |
|
|
|
} |
|
|
|
|