home *** CD-ROM | disk | FTP | other *** search
/ Power-Programmierung / CD2.mdf / c / library / dos / math / pac / convbase.c < prev    next >
Encoding:
C/C++ Source or Header  |  1994-06-05  |  8.0 KB  |  244 lines
  1. /* convbase.c */
  2. /**********************************************************************
  3. *    File Name     : convbase.c
  4. *    Function      : use second bc pipe for conversion between radices
  5. *    Author        : Istvan Mohos, 1987
  6. ***********************************************************************/
  7.  
  8. #include "defs.h"
  9. #include "toktab.h"
  10. #define INTERMAP
  11. #include "maps.h"
  12. #undef INTERMAP
  13.  
  14. char *
  15. substivar(found, token, ib)
  16. char *token;
  17. int found, ib;
  18. { 
  19.     static char res[PIPEMAX];
  20.     static struct stk_cell *sr = &Stk[0];
  21.     static char *fid = "substivar";
  22.     char *ret = res;
  23.     int tokval;
  24.  
  25.     _TR
  26.     if (found == -1) {
  27.         if (token == ZERO) {
  28.             TR_
  29.             return(ZERO);
  30.         }
  31.         tokval = lookup(token);
  32.     }
  33.     else
  34.         tokval = found;
  35.  
  36.     switch(tokval) {
  37.         case H_:
  38.         case I_:
  39.         case J_:
  40.         case K_:
  41.         case L_:
  42.         case M_:
  43.         case N_:
  44.         case O_:
  45.         case P_:
  46.         case Q_:
  47.         case R_:
  48.         case S_:
  49.         case T_:
  50.         case U_:
  51.         case V_:
  52.         case W_:
  53.             onereg(*token - 'g');
  54.             conv_bc(Onebuf, res, 1, ib);
  55.             break;
  56.  
  57. case PI: conv_bc("3.1415926535897932384626433832795028841971693993751",
  58.                     res, 10, ib); break;
  59. case ASTRO:     conv_bc("149504200", res, 10, ib); break;
  60.                     /* km */
  61. case AMASS:     conv_bc(".00000000000000000000000165975",
  62.                     res, 10, ib); break; /* grams */
  63. case AVOGADRO:  conv_bc("602502000000000000000000",
  64.                     res, 10, ib); break; /* per g mole */
  65. case BOLTZMANN: conv_bc(".000000000000000138040", res, 10, ib); break;
  66.                     /* erg/Kelvin */
  67. case ECHARGE:   conv_bc(".000000000160202", res, 10, ib); break;
  68.                     /* esu */
  69. case CHROMA:    conv_bc("1.0594631", res, 10, ib); break;
  70. case EMASS:     conv_bc(".000000000000000000000000000910820",
  71.                     res, 10, ib); break; /* grams */
  72. case EULER:     conv_bc(".577216", res, 10, ib); break;
  73. case FARADAY:   conv_bc("96521900", res, 10, ib); break;
  74.                     /* C/kmole */
  75. case G_:        conv_bc("9.80665", res, 10, ib); break;
  76.                     /* m/s2 */
  77. case GAS:       conv_bc("83169600", res, 10, ib); break;
  78.                     /* erg/g mole Kelvin */
  79. case GRAVITY:   conv_bc(".00000000006673", res, 10, ib); break;
  80.                     /* N m2/kg2 */
  81. case HEAT:      conv_bc("4.1855", res, 10, ib); break;
  82.                     /* J/cal */
  83. case LIGHT:     conv_bc("299792.50", res, 10, ib); break;
  84.                     /* km/sec */
  85. case LIGHTYEAR: conv_bc("9460530000000", res, 10, ib); break;
  86.                     /* km */
  87. case MOONMASS:  conv_bc("73430000000000000000000", res, 10, ib); break;
  88.                     /* kg */
  89. case SUNMASS:   conv_bc("1987000000000000000000000000000",
  90.                     res, 10, ib); break; /* kg */
  91. case EARTHMASS: conv_bc("5976500000000000000000000",
  92.                     res, 10, ib); break; /* kg */
  93. case NATURAL:   conv_bc("2.7182818284590452353602874713526",
  94.                     res, 10, ib); break;
  95. case NMASS:     conv_bc(".00000000000000000000000167465",
  96.                     res, 10, ib); break; /* grams */
  97. case PARSEC:    conv_bc("30837450000000", res, 10, ib); break;
  98.                     /* km */
  99. case PARALLAX:  conv_bc("8.794", res, 10, ib); break;
  100.                     /* " */
  101. case PLANCK:    conv_bc(".00000000000000000000000000662491",
  102.                     res, 10, ib); break; /* erg sec */
  103. case PMASS:     conv_bc(".00000000000000000000000167235",
  104.                     res, 10, ib); break; /* grams */
  105. case MOONRAD:   conv_bc("1738000", res, 10, ib); break;
  106.                     /* meters */
  107. case SUNRAD:    conv_bc("696500000", res, 10, ib); break;
  108.                     /* meters */
  109. case EARTHRAD:  conv_bc("6378388", res, 10, ib); break;
  110.                     /* meters */
  111. case RYDBERG:   conv_bc("10973732.8", res, 10, ib); break;
  112.                     /* per meter */
  113. case SOUND:     conv_bc("340.505", res, 10, ib); break;
  114.                     /* meters/sec: 331.4 + (.607 * Celsius) */
  115. case STEFAN:    conv_bc(".000000056693", res, 10, ib); break;
  116.                     /* J/m2 Kelvin4 sec */
  117. case TOMOON:    conv_bc("384400", res, 10, ib); break;
  118.                     /* km */
  119. case TOSUN:     conv_bc("149500000", res, 10, ib); break;
  120.                     /* km */
  121. case WIEN:      conv_bc(".289778", res, 10, ib); break;
  122.                     /* cm Kelvin */
  123.  
  124. case MILLI:     conv_bc(".001", res, 10, ib); break;
  125. case MICRO:     conv_bc(".000001", res, 10, ib); break;
  126. case NANO:      conv_bc(".000000001", res, 10, ib); break;
  127. case PICO:      conv_bc(".000000000001", res, 10, ib); break;
  128. case FEMTO:     conv_bc(".000000000000001", res, 10, ib); break;
  129. case ATTO:      conv_bc(".000000000000000001", res, 10, ib); break;
  130. case KILO:      conv_bc("1000", res, 10, ib); break;
  131. case MEGA:      conv_bc("1000000", res, 10, ib); break;
  132. case GIGA:      conv_bc("1000000000", res, 10, ib); break;
  133. case TERA:      conv_bc("1000000000000", res, 10, ib); break;
  134. case PETA:      conv_bc("1000000000000000", res, 10, ib); break;
  135. case EXA:       conv_bc("1000000000000000000", res, 10, ib); break;
  136.  
  137.         case X_LOWER:
  138.         case X_UPPER:
  139.             sprintf(res, "%s", sixteen[Ibase]);
  140.             break;
  141.  
  142.         case BACKSLASH:
  143.             conv_bc(sr->cell, res, 1, ib);
  144.             break;
  145.  
  146.         case NOTINLIST:    /* presumably a digit string */
  147.             upcase(token);
  148.             ret = token;
  149.             break;
  150.  
  151.         default: /* token in list, but is not variable or constant */
  152.             ret = Convbuf;   /* just to warn the caller */
  153.             break;
  154.     }
  155.     TR_
  156.     return(ret);
  157. }
  158.  
  159. /* translate frombuf number string in frombase radix,
  160.    to string in tobuf in tobase radix.
  161.  
  162.    If frombase and tobase are identical no conversion is done;
  163.    buffer is simply copied.
  164.    Otherwise, conv_pipe's ibase is set to the frombase radix,
  165.    conv_pipe's obase is set to tobase, and the frombuf
  166.    string is executed through conv_pipe.
  167.  
  168.    If frombase is 0, Lastob is taken as frombase.
  169.    If frombase is -1, Ibase is taken as frombase.
  170.    If frombase is 1, frombase is deduced from the first character
  171.        of frombuf in stack string format.
  172.  
  173.    If tobuf is ZERO, result goes to Convbuf.  Else tobuf must
  174.    be at least PIPEMAX long.
  175. */
  176.  
  177. conv_bc(frombuf, tobuf, frombase, tobase)
  178. char *frombuf, *tobuf;
  179. int frombase, tobase;
  180. {
  181.     int value, maxchars;
  182.     int charcount = 0;
  183.     register char *bptr, *toasc_p;
  184.     static char *fid = "conv_bc";
  185.  
  186.     _TR
  187.     bptr = frombuf;
  188.     if (tobase == 0)
  189.         tobase = Ibase;
  190.     maxchars = strlen(frombuf);
  191.     if (tobuf == ZERO)
  192.         tobuf = Convbuf;
  193.  
  194.     switch (frombase) {
  195.         case 1:
  196.             charcount = 1;
  197.             value = *bptr & 127; /* ascii value of passed radix label */
  198.             for (toasc_p = Base_str + 16;; toasc_p--) {
  199.                 if (*toasc_p == value) {
  200.                     frombase = toasc_p - Base_str;
  201.                     break;
  202.                 }
  203.             }
  204.             while ((*++bptr & 127) == ' ') /* get to digit string */
  205.                 ++charcount;
  206.             break;
  207.  
  208.         case -1:
  209.             frombase = Ibase;
  210.             break;
  211.  
  212.         case 0:
  213.             frombase = Lastob;
  214.             break;
  215.  
  216.         default:
  217.             break;
  218.  
  219.     }
  220.  
  221.     if (frombase == tobase) {
  222.         toasc_p = tobuf;
  223.         while(charcount++ < maxchars)
  224.             if ((*bptr & 127) != ' ')
  225.                 *toasc_p++ = *bptr++ & 127;
  226.             else
  227.                 ++bptr;
  228.         *toasc_p = '\0';
  229.     }
  230.     else {
  231.         sprintf(Mop, "ibase=A;obase=%d;ibase=%d\n", tobase, frombase);
  232.         toasc_p = &Mop[strlen(Mop)];
  233.         while(charcount++ < maxchars)
  234.             *toasc_p++ = *bptr++ & 127;
  235.         *toasc_p++ = '\n';
  236.         *toasc_p = '\0';
  237.         if (write (B_write[1], Mop, strlen(Mop)) == -1)
  238.             fatal("wait_conv_pipe write");
  239.         wait_conv_pipe(tobuf);
  240.     }
  241.     TR_
  242. }
  243.  
  244.