[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.8 and 1.10

version 1.8, 2003/12/26 05:47:37 version 1.10, 2011/12/15 16:53:26
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.7 2003/12/24 08:00:38 noro Exp $   * $OpenXM: OpenXM_contrib2/asir2000/builtin/math.c,v 1.9 2011/08/24 07:20:09 noro Exp $
 */  */
 #include "ca.h"  #include "ca.h"
 #include <math.h>  #include <math.h>
 #include "parse.h"  #include "parse.h"
   #if defined(VISUAL)
   #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();  void Pabs(),Pdfloor(),Pdceil(),Pdrint(),Pdisnan();
   
 struct ftab math_tab[] = {  struct ftab math_tab[] = {
         {"dsqrt",Pdsqrt,1},          {"dsqrt",Pdsqrt,1},
Line 71  struct ftab math_tab[] = {
Line 74  struct ftab math_tab[] = {
         {"dceil",Pdceil,1},          {"dceil",Pdceil,1},
         {"rint",Pdrint,1},          {"rint",Pdrint,1},
         {"drint",Pdrint,1},          {"drint",Pdrint,1},
           {"disnan",Pdisnan,1},
         {0,0,0},          {0,0,0},
 };  };
   
Line 359  Q *rp;
Line 363  Q *rp;
         } else          } else
                 q = 0;                  q = 0;
         *rp = q;          *rp = q;
   }
   
   void Pdisnan(NODE arg,Q *rp)
   {
           Real r;
           double d;
   #if defined(VISUAL)
       int c;
   #endif
   
           r = (Real)ARG0(arg);
           if ( !r || !NUM(r) || !REAL(r) ) {
                   *rp = 0;
                   return;
           }
           d = ToReal(r);
   #if defined(VISUAL)
           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;
           else if ( isinf(d) ) STOQ(2,*rp);
   #endif
           else *rp = 0;
 }  }
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.10
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>