home *** CD-ROM | disk | FTP | other *** search
/ Power Programmierung 2 / Power-Programmierung CD 2 (Tewi)(1994).iso / c / library / dos / math / cephes / remes / remesg.c < prev    next >
Encoding:
C/C++ Source or Header  |  1992-11-17  |  1.7 KB  |  38 lines
  1. /* remesg.c */
  2. /* Global symbol declarations for Remes program */
  3.  
  4. #include "remes.h"
  5.  
  6. int IPS[P] = {0,};         /* simq() work vector */
  7. double AA[P*P] = {0.0,};   /* coefficient matrix */
  8. double BB[P] = {0.0,};     /* right hand side vector */
  9. double param[P] = {0.0,};  /* solution vector */
  10. double xx[P] = {0.0,};     /* points in approximation interval */
  11. double mm[P] = {0.0};      /* points of maximum deviation */
  12. double yy[P] = {0.0};      /* the error maxima */
  13. double step[P] = {0.0};    /* step sizes for searching */
  14. double qx = 0.0;           /* function argument */
  15. double qy = 0.0;           /* accurate function value */
  16. double qyaprx = 0.0;       /* approximate function value */
  17. double dev = 0.0;          /* deviation at solution points */
  18. double apstrt = 0.0;       /* start of approximation interval */
  19. double apwidt = 0.0;       /* width of interval */
  20. double apend = 0.0;        /* end of interval */
  21. double xm = 0.0;           /* variables for search of maximum */
  22. double xn = 0.0;
  23. double ym = 0.0;
  24. double yn = 0.0;
  25. double delta = 0.0;
  26. double eclose = 0.0;
  27. double farther = 0.0;
  28. double spread = 0.0;      /* error spread from last iteration */
  29. int esign = 0;            /* sign of error returned by geterr */
  30. int n = 0;                /* degree of numerator polynomial */
  31. int d = 0;                /* degree of denominator polynomial */
  32. int nd1 = 0;              /* n + d + 1 */
  33. int neq = 0;              /* number of equations to solve */
  34. int relerr = 0;           /* relative error criterion flag */
  35. int search = 0;           /* flag for automatic search */
  36. int iter = 0;             /* search-solve iteration number */
  37. int askitr = 0;           /* iteration on which to stop */
  38.