home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1992 #30 / NN_1992_30.iso / spool / comp / os / msdos / programm / 11484 < prev    next >
Encoding:
Text File  |  1992-12-17  |  2.4 KB  |  73 lines
  1. Path: sparky!uunet!spool.mu.edu!umn.edu!csus.edu!netcom.com!netcomsv!wet!mcs
  2. From: mcs@wet.UUCP (Gil Nardo)
  3. Newsgroups: comp.os.msdos.programmer
  4. Subject: Re: Math error in BC++?  Help!
  5. Message-ID: <4874@wet.UUCP>
  6. Date: 17 Dec 92 19:08:19 GMT
  7. References: <1992Dec15.211516.1664@schunix.uucp>
  8. Distribution: usa
  9. Organization: Wetware Diversions, San Francisco
  10. Lines: 61
  11.  
  12. sonix@schunix.uucp (Duane Morin) writes:
  13. > This is driving us NUTS!! Will someone please take a look at the following 
  14. > code and see if they can figure out why it does what it does?  Thank you.
  15. > 
  16. > The following code segment shows a snippet of a program that needs
  17. > to calculate
  18. > the log base 2 of 128 (which happens to be 7).  The problem is that the 
  19. > code works fine on a 386DX/40, with 387 installed, and does NOT work on 
  20. > a 486DX/33!  Help??
  21. > 
  22.  
  23. I included a shortened version of your program.  I did not completely
  24. solved it but I do have some clues as to what's causing the problem.
  25. After compiling with the assembly output option, I found that at
  26. least two methods can be used to make it return the answer of 7.
  27.  
  28. One method is to pop the answer off of the co-processor stack and
  29. pop it back on using the fstp and fld instructions prior to the
  30. call to N_FTOL (float to long, near call version). This method
  31. is equivalent to the two step processes of as shown in the code.
  32.  
  33. The other method is to substitute N_FTOL with the use of fistp
  34. (store and pop stack as integer) and followed by the use of the
  35. lower part of the of the quad word as the result.
  36.  
  37. I would guess at this point that the function N_FTOL inadvertantly
  38. modifies either the rounding or the precision bits of the
  39. co-processor control word. I will leave it to you and others to
  40. figure out why this happens on the 486 chip and not
  41. on the 386/387 chips, mainly because I need to get back to work.
  42.  
  43. ----------
  44. #include <stdio.h>
  45. #include <stdlib.h>
  46. #include <math.h>
  47. #include <float.h>
  48.  
  49. // The following constants are normally defined within the program:
  50.  
  51. #define LOG2 (log(2))
  52. #define LOGN(y)        (log ((double)y) / LOG2 )
  53. main()
  54. {
  55.      double l;
  56.     int i;
  57.  
  58.     l = LOGN(128);
  59.     i = (int)l;
  60.     printf("two step = %d\n", i); // this comes up with 7
  61.  
  62.     i = (int)(LOGN(128));
  63.     printf("one step = %d\n", i); // this comes up with 6?
  64.  
  65.     return(0);
  66. }
  67.  
  68. -- 
  69. Gil Nardo                  | mcs@wet.com
  70. Migrant Computing Services | (415)664-1032
  71. 1032 Irving Street, #435   |-----------------
  72. San Francisco, 94122       | Supernova = *!
  73.