![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to math.c CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM_contrib2 / asir2000 / builtin|
Return to math.c CVS log | Up to [local] / OpenXM_contrib2 / asir2000 / builtin |
1.1 ! noro 1: /* $OpenXM: OpenXM/src/asir99/builtin/math.c,v 1.1.1.1 1999/11/10 08:12:25 noro Exp $ */
! 2: #include "ca.h"
! 3: #include <math.h>
! 4: #include "parse.h"
! 5:
! 6: void Pdsqrt(),Pdsin(),Pdcos(),Pdtan(),Pdasin(),Pdacos(),Pdatan(),Pdlog(),Pdexp();
! 7: void Pabs(),Pdfloor(),Pdceil(),Pdrint();
! 8:
! 9: struct ftab math_tab[] = {
! 10: {"dsqrt",Pdsqrt,1},
! 11: {"dabs",Pabs,1},
! 12: {"dsin",Pdsin,1},
! 13: {"dcos",Pdcos,1},
! 14: {"dtan",Pdtan,1},
! 15: {"dlog",Pdlog,1},
! 16: {"dexp",Pdexp,1},
! 17: {"dasin",Pdasin,1},
! 18: {"dacos",Pdacos,1},
! 19: {"datan",Pdatan,1},
! 20: {"dfloor",Pdfloor,1},
! 21: {"dceil",Pdceil,1},
! 22: {"drint",Pdrint,1},
! 23: {0,0,0},
! 24: };
! 25:
! 26: void Pabs(arg,rp)
! 27: NODE arg;
! 28: Real *rp;
! 29: {
! 30: double s;
! 31:
! 32: s = fabs(ToReal(ARG0(arg)));
! 33: MKReal(s,*rp);
! 34: }
! 35:
! 36: void Pdsqrt(arg,rp)
! 37: NODE arg;
! 38: Real *rp;
! 39: {
! 40: double s;
! 41:
! 42: s = sqrt(ToReal(ARG0(arg)));
! 43: MKReal(s,*rp);
! 44: }
! 45:
! 46: void Pdsin(arg,rp)
! 47: NODE arg;
! 48: Real *rp;
! 49: {
! 50: double s;
! 51:
! 52: s = sin(ToReal(ARG0(arg)));
! 53: MKReal(s,*rp);
! 54: }
! 55:
! 56: void Pdcos(arg,rp)
! 57: NODE arg;
! 58: Real *rp;
! 59: {
! 60: double s;
! 61:
! 62: s = cos(ToReal(ARG0(arg)));
! 63: MKReal(s,*rp);
! 64: }
! 65:
! 66: void Pdtan(arg,rp)
! 67: NODE arg;
! 68: Real *rp;
! 69: {
! 70: double s;
! 71:
! 72: s = tan(ToReal(ARG0(arg)));
! 73: MKReal(s,*rp);
! 74: }
! 75:
! 76: void Pdasin(arg,rp)
! 77: NODE arg;
! 78: Real *rp;
! 79: {
! 80: double s;
! 81:
! 82: s = asin(ToReal(ARG0(arg)));
! 83: MKReal(s,*rp);
! 84: }
! 85:
! 86: void Pdacos(arg,rp)
! 87: NODE arg;
! 88: Real *rp;
! 89: {
! 90: double s;
! 91:
! 92: s = acos(ToReal(ARG0(arg)));
! 93: MKReal(s,*rp);
! 94: }
! 95:
! 96: void Pdatan(arg,rp)
! 97: NODE arg;
! 98: Real *rp;
! 99: {
! 100: double s;
! 101:
! 102: s = atan(ToReal(ARG0(arg)));
! 103: MKReal(s,*rp);
! 104: }
! 105:
! 106: void Pdlog(arg,rp)
! 107: NODE arg;
! 108: Real *rp;
! 109: {
! 110: double s;
! 111:
! 112: s = log(ToReal(ARG0(arg)));
! 113: MKReal(s,*rp);
! 114: }
! 115:
! 116: void Pdexp(arg,rp)
! 117: NODE arg;
! 118: Real *rp;
! 119: {
! 120: double s;
! 121:
! 122: s = exp(ToReal(ARG0(arg)));
! 123: MKReal(s,*rp);
! 124: }
! 125:
! 126: void Pdfloor(arg,rp)
! 127: NODE arg;
! 128: Q *rp;
! 129: {
! 130: L a;
! 131: unsigned int au,al;
! 132: int sgn;
! 133: Q q;
! 134: double d;
! 135:
! 136: if ( !ARG0(arg) ) {
! 137: *rp = 0;
! 138: return;
! 139: }
! 140: d = floor(ToReal(ARG0(arg)));
! 141: if ( d < -9.223372036854775808e18 || d >= 9.223372036854775808e18 )
! 142: error("dfloor : OverFlow");
! 143: a = (L)d;
! 144: if ( a < 0 ) {
! 145: sgn = -1;
! 146: a = -a;
! 147: } else
! 148: sgn = 1;
! 149: #if defined(i386) || defined(__alpha) || defined(VISUAL)
! 150: au = ((unsigned int *)&a)[1];
! 151: al = ((unsigned int *)&a)[0];
! 152: #else
! 153: al = ((unsigned int *)&a)[1];
! 154: au = ((unsigned int *)&a)[0];
! 155: #endif
! 156: if ( au ) {
! 157: NEWQ(q); SGN(q) = sgn; NM(q)=NALLOC(2); DN(q)=0;
! 158: PL(NM(q))=2; BD(NM(q))[0]=al; BD(NM(q))[1] = au;
! 159: } else {
! 160: UTOQ(al,q); SGN(q) = sgn;
! 161: }
! 162: *rp = q;
! 163: }
! 164:
! 165: void Pdceil(arg,rp)
! 166: NODE arg;
! 167: Q *rp;
! 168: {
! 169: L a;
! 170: unsigned int au,al;
! 171: int sgn;
! 172: Q q;
! 173: double d;
! 174:
! 175: if ( !ARG0(arg) ) {
! 176: *rp = 0;
! 177: return;
! 178: }
! 179: d = ceil(ToReal(ARG0(arg)));
! 180: if ( d < -9.223372036854775808e18 || d >= 9.223372036854775808e18 )
! 181: error("dceil : OverFlow");
! 182: a = (L)d;
! 183: if ( a < 0 ) {
! 184: sgn = -1;
! 185: a = -a;
! 186: } else
! 187: sgn = 1;
! 188: #if defined(i386) || defined(__alpha) || defined(VISUAL)
! 189: au = ((unsigned int *)&a)[1];
! 190: al = ((unsigned int *)&a)[0];
! 191: #else
! 192: al = ((unsigned int *)&a)[1];
! 193: au = ((unsigned int *)&a)[0];
! 194: #endif
! 195: if ( au ) {
! 196: NEWQ(q); SGN(q) = sgn; NM(q)=NALLOC(2); DN(q)=0;
! 197: PL(NM(q))=2; BD(NM(q))[0]=al; BD(NM(q))[1] = au;
! 198: } else {
! 199: UTOQ(al,q); SGN(q) = sgn;
! 200: }
! 201: *rp = q;
! 202: }
! 203:
! 204: void Pdrint(arg,rp)
! 205: NODE arg;
! 206: Q *rp;
! 207: {
! 208: L a;
! 209: unsigned int au,al;
! 210: int sgn;
! 211: Q q;
! 212: double d;
! 213:
! 214: if ( !ARG0(arg) ) {
! 215: *rp = 0;
! 216: return;
! 217: }
! 218: #if defined(VISUAL)
! 219: d = ToReal(ARG0(arg));
! 220: if ( d > 0 )
! 221: d = floor(d+0.5);
! 222: else
! 223: d = ceil(d-0.5);
! 224: #else
! 225: d = rint(ToReal(ARG0(arg)));
! 226: #endif
! 227: if ( d < -9.223372036854775808e18 || d >= 9.223372036854775808e18 )
! 228: error("drint : OverFlow");
! 229: a = (L)d;
! 230: if ( a < 0 ) {
! 231: sgn = -1;
! 232: a = -a;
! 233: } else
! 234: sgn = 1;
! 235: #if defined(i386) || defined(__alpha) || defined(VISUAL)
! 236: au = ((unsigned int *)&a)[1];
! 237: al = ((unsigned int *)&a)[0];
! 238: #else
! 239: al = ((unsigned int *)&a)[1];
! 240: au = ((unsigned int *)&a)[0];
! 241: #endif
! 242: if ( au ) {
! 243: NEWQ(q); SGN(q) = sgn; NM(q)=NALLOC(2); DN(q)=0;
! 244: PL(NM(q))=2; BD(NM(q))[0]=al; BD(NM(q))[1] = au;
! 245: } else {
! 246: UTOQ(al,q); SGN(q) = sgn;
! 247: }
! 248: *rp = q;
! 249: }