[BACK]Return to math.c CVS log [TXT][DIR] Up to [local] / OpenXM_contrib2 / asir2000 / builtin

Diff for /OpenXM_contrib2/asir2000/builtin/math.c between version 1.9 and 1.11

version 1.9, 2011/08/24 07:20:09 version 1.11, 2015/08/06 10:01:52
Line 45 
Line 45 
  * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,   * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,
  * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.   * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.
  *   *
  * $OpenXM: OpenXM_contrib2/asir2000/builtin/math.c,v 1.8 2003/12/26 05:47:37 noro Exp $   * $OpenXM: OpenXM_contrib2/asir2000/builtin/math.c,v 1.10 2011/12/15 16:53:26 ohara Exp $
 */  */
 #include "ca.h"  #include "ca.h"
 #include <math.h>  #include <math.h>
 #include "parse.h"  #include "parse.h"
   #if defined(VISUAL) || defined(__MINGW32__) || defined(__MINGW64__)
   #include <float.h>
   #endif
   
 void Pdsqrt(),Pdsin(),Pdcos(),Pdtan(),Pdasin(),Pdacos(),Pdatan(),Pdlog(),Pdexp();  void Pdsqrt(),Pdsin(),Pdcos(),Pdtan(),Pdasin(),Pdacos(),Pdatan(),Pdlog(),Pdexp();
 void Pabs(),Pdfloor(),Pdceil(),Pdrint(),Pdisnan();  void Pabs(),Pdfloor(),Pdceil(),Pdrint(),Pdisnan();
Line 255  Q *rp;
Line 258  Q *rp;
                 a = -a;                  a = -a;
         } else          } else
                 sgn = 1;                  sgn = 1;
 #if defined(i386) || defined(__alpha) || defined(VISUAL) || defined(__x86_64)  #if defined(i386) || defined(__alpha) || defined(VISUAL) || defined(__MINGW32__) || defined(__MINGW64__) || defined(__x86_64)
         au = ((unsigned int *)&a)[1];          au = ((unsigned int *)&a)[1];
         al = ((unsigned int *)&a)[0];          al = ((unsigned int *)&a)[0];
 #else  #else
Line 298  Q *rp;
Line 301  Q *rp;
                 a = -a;                  a = -a;
         } else          } else
                 sgn = 1;                  sgn = 1;
 #if defined(i386) || defined(__alpha) || defined(VISUAL) || defined(__x86_64)  #if defined(i386) || defined(__alpha) || defined(VISUAL) || defined(__MINGW32__) || defined(__MINGW64__) || defined(__x86_64)
         au = ((unsigned int *)&a)[1];          au = ((unsigned int *)&a)[1];
         al = ((unsigned int *)&a)[0];          al = ((unsigned int *)&a)[0];
 #else  #else
Line 328  Q *rp;
Line 331  Q *rp;
                 *rp = 0;                  *rp = 0;
                 return;                  return;
         }          }
 #if defined(VISUAL)  #if defined(VISUAL) || defined(__MINGW32__) || defined(__MINGW64__)
         d = ToReal(ARG0(arg));          d = ToReal(ARG0(arg));
         if ( d > 0 )          if ( d > 0 )
                 d = floor(d+0.5);                  d = floor(d+0.5);
Line 345  Q *rp;
Line 348  Q *rp;
                 a = -a;                  a = -a;
         } else          } else
                 sgn = 1;                  sgn = 1;
 #if defined(i386) || defined(__alpha) || defined(VISUAL) || defined(__x86_64)  #if defined(i386) || defined(__alpha) || defined(VISUAL) || defined(__MINGW32__) || defined(__MINGW64__) || defined(__x86_64)
         au = ((unsigned int *)&a)[1];          au = ((unsigned int *)&a)[1];
         al = ((unsigned int *)&a)[0];          al = ((unsigned int *)&a)[0];
 #else  #else
Line 366  void Pdisnan(NODE arg,Q *rp)
Line 369  void Pdisnan(NODE arg,Q *rp)
 {  {
         Real r;          Real r;
         double d;          double d;
   #if defined(VISUAL) || defined(__MINGW32__) || defined(__MINGW64__)
       int c;
   #endif
   
         r = (Real)ARG0(arg);          r = (Real)ARG0(arg);
         if ( !r || !NUM(r) || !REAL(r) ) {          if ( !r || !NUM(r) || !REAL(r) ) {
Line 373  void Pdisnan(NODE arg,Q *rp)
Line 379  void Pdisnan(NODE arg,Q *rp)
                 return;                  return;
         }          }
         d = ToReal(r);          d = ToReal(r);
   #if defined(VISUAL) || defined(__MINGW32__) || defined(__MINGW64__)
           c = _fpclass(d);
           if ( c == _FPCLASS_SNAN || c == _FPCLASS_QNAN ) *rp = ONE;
           else if ( c == _FPCLASS_PINF || c == _FPCLASS_NINF ) STOQ(2,*rp);
   #else
         if ( isnan(d) ) *rp = ONE;          if ( isnan(d) ) *rp = ONE;
         else if ( isinf(d) ) STOQ(2,*rp);          else if ( isinf(d) ) STOQ(2,*rp);
   #endif
         else *rp = 0;          else *rp = 0;
 }  }
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.11
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>