Annotation of OpenXM_contrib2/asir2000/builtin/math.c, Revision 1.7
1.2 noro 1: /*
2: * Copyright (c) 1994-2000 FUJITSU LABORATORIES LIMITED
3: * All rights reserved.
4: *
5: * FUJITSU LABORATORIES LIMITED ("FLL") hereby grants you a limited,
6: * non-exclusive and royalty-free license to use, copy, modify and
7: * redistribute, solely for non-commercial and non-profit purposes, the
8: * computer program, "Risa/Asir" ("SOFTWARE"), subject to the terms and
9: * conditions of this Agreement. For the avoidance of doubt, you acquire
10: * only a limited right to use the SOFTWARE hereunder, and FLL or any
11: * third party developer retains all rights, including but not limited to
12: * copyrights, in and to the SOFTWARE.
13: *
14: * (1) FLL does not grant you a license in any way for commercial
15: * purposes. You may use the SOFTWARE only for non-commercial and
16: * non-profit purposes only, such as academic, research and internal
17: * business use.
18: * (2) The SOFTWARE is protected by the Copyright Law of Japan and
19: * international copyright treaties. If you make copies of the SOFTWARE,
20: * with or without modification, as permitted hereunder, you shall affix
21: * to all such copies of the SOFTWARE the above copyright notice.
22: * (3) An explicit reference to this SOFTWARE and its copyright owner
23: * shall be made on your publication or presentation in any form of the
24: * results obtained by use of the SOFTWARE.
25: * (4) In the event that you modify the SOFTWARE, you shall notify FLL by
1.3 noro 26: * e-mail at risa-admin@sec.flab.fujitsu.co.jp of the detailed specification
1.2 noro 27: * for such modification or the source code of the modified part of the
28: * SOFTWARE.
29: *
30: * THE SOFTWARE IS PROVIDED AS IS WITHOUT ANY WARRANTY OF ANY KIND. FLL
31: * MAKES ABSOLUTELY NO WARRANTIES, EXPRESSED, IMPLIED OR STATUTORY, AND
32: * EXPRESSLY DISCLAIMS ANY IMPLIED WARRANTY OF MERCHANTABILITY, FITNESS
33: * FOR A PARTICULAR PURPOSE OR NONINFRINGEMENT OF THIRD PARTIES'
34: * RIGHTS. NO FLL DEALER, AGENT, EMPLOYEES IS AUTHORIZED TO MAKE ANY
35: * MODIFICATIONS, EXTENSIONS, OR ADDITIONS TO THIS WARRANTY.
36: * UNDER NO CIRCUMSTANCES AND UNDER NO LEGAL THEORY, TORT, CONTRACT,
37: * OR OTHERWISE, SHALL FLL BE LIABLE TO YOU OR ANY OTHER PERSON FOR ANY
38: * DIRECT, INDIRECT, SPECIAL, INCIDENTAL, PUNITIVE OR CONSEQUENTIAL
39: * DAMAGES OF ANY CHARACTER, INCLUDING, WITHOUT LIMITATION, DAMAGES
40: * ARISING OUT OF OR RELATING TO THE SOFTWARE OR THIS AGREEMENT, DAMAGES
41: * FOR LOSS OF GOODWILL, WORK STOPPAGE, OR LOSS OF DATA, OR FOR ANY
42: * DAMAGES, EVEN IF FLL SHALL HAVE BEEN INFORMED OF THE POSSIBILITY OF
43: * SUCH DAMAGES, OR FOR ANY CLAIM BY ANY OTHER PARTY. EVEN IF A PART
44: * OF THE SOFTWARE HAS BEEN DEVELOPED BY A THIRD PARTY, THE THIRD PARTY
45: * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,
46: * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.
47: *
1.7 ! noro 48: * $OpenXM: OpenXM_contrib2/asir2000/builtin/math.c,v 1.6 2003/12/11 00:19:24 noro Exp $
1.2 noro 49: */
1.1 noro 50: #include "ca.h"
51: #include <math.h>
52: #include "parse.h"
53:
54: void Pdsqrt(),Pdsin(),Pdcos(),Pdtan(),Pdasin(),Pdacos(),Pdatan(),Pdlog(),Pdexp();
55: void Pabs(),Pdfloor(),Pdceil(),Pdrint();
56:
57: struct ftab math_tab[] = {
58: {"dsqrt",Pdsqrt,1},
59: {"dabs",Pabs,1},
60: {"dsin",Pdsin,1},
61: {"dcos",Pdcos,1},
62: {"dtan",Pdtan,1},
63: {"dlog",Pdlog,1},
64: {"dexp",Pdexp,1},
65: {"dasin",Pdasin,1},
66: {"dacos",Pdacos,1},
67: {"datan",Pdatan,1},
1.5 noro 68: {"floor",Pdfloor,1},
1.1 noro 69: {"dfloor",Pdfloor,1},
1.5 noro 70: {"ceil",Pdceil,1},
1.1 noro 71: {"dceil",Pdceil,1},
1.5 noro 72: {"rint",Pdrint,1},
1.1 noro 73: {"drint",Pdrint,1},
74: {0,0,0},
75: };
76:
77: void Pabs(arg,rp)
78: NODE arg;
79: Real *rp;
80: {
81: double s;
82:
83: s = fabs(ToReal(ARG0(arg)));
84: MKReal(s,*rp);
85: }
86:
87: void Pdsqrt(arg,rp)
88: NODE arg;
89: Real *rp;
90: {
91: double s;
92:
93: s = sqrt(ToReal(ARG0(arg)));
94: MKReal(s,*rp);
95: }
96:
97: void Pdsin(arg,rp)
98: NODE arg;
99: Real *rp;
100: {
101: double s;
102:
103: s = sin(ToReal(ARG0(arg)));
104: MKReal(s,*rp);
105: }
106:
107: void Pdcos(arg,rp)
108: NODE arg;
109: Real *rp;
110: {
111: double s;
112:
113: s = cos(ToReal(ARG0(arg)));
114: MKReal(s,*rp);
115: }
116:
117: void Pdtan(arg,rp)
118: NODE arg;
119: Real *rp;
120: {
121: double s;
122:
123: s = tan(ToReal(ARG0(arg)));
124: MKReal(s,*rp);
125: }
126:
127: void Pdasin(arg,rp)
128: NODE arg;
129: Real *rp;
130: {
131: double s;
132:
133: s = asin(ToReal(ARG0(arg)));
134: MKReal(s,*rp);
135: }
136:
137: void Pdacos(arg,rp)
138: NODE arg;
139: Real *rp;
140: {
141: double s;
142:
143: s = acos(ToReal(ARG0(arg)));
144: MKReal(s,*rp);
145: }
146:
147: void Pdatan(arg,rp)
148: NODE arg;
149: Real *rp;
150: {
151: double s;
152:
153: s = atan(ToReal(ARG0(arg)));
154: MKReal(s,*rp);
155: }
156:
157: void Pdlog(arg,rp)
158: NODE arg;
159: Real *rp;
160: {
161: double s;
162:
163: s = log(ToReal(ARG0(arg)));
164: MKReal(s,*rp);
165: }
166:
167: void Pdexp(arg,rp)
168: NODE arg;
169: Real *rp;
170: {
171: double s;
172:
173: s = exp(ToReal(ARG0(arg)));
174: MKReal(s,*rp);
175: }
176:
177: void Pdfloor(arg,rp)
178: NODE arg;
179: Q *rp;
180: {
181: L a;
182: unsigned int au,al;
183: int sgn;
184: Q q;
185: double d;
186:
187: if ( !ARG0(arg) ) {
188: *rp = 0;
189: return;
190: }
191: d = floor(ToReal(ARG0(arg)));
192: if ( d < -9.223372036854775808e18 || d >= 9.223372036854775808e18 )
193: error("dfloor : OverFlow");
1.4 noro 194: if ( !d ) {
195: *rp = 0;
196: return;
197: }
1.1 noro 198: a = (L)d;
199: if ( a < 0 ) {
200: sgn = -1;
201: a = -a;
202: } else
203: sgn = 1;
1.7 ! noro 204: #if defined(i386) || defined(__alpha) || defined(VISUAL) || defined(__x86_64)
1.1 noro 205: au = ((unsigned int *)&a)[1];
206: al = ((unsigned int *)&a)[0];
207: #else
208: al = ((unsigned int *)&a)[1];
209: au = ((unsigned int *)&a)[0];
210: #endif
211: if ( au ) {
212: NEWQ(q); SGN(q) = sgn; NM(q)=NALLOC(2); DN(q)=0;
213: PL(NM(q))=2; BD(NM(q))[0]=al; BD(NM(q))[1] = au;
214: } else {
215: UTOQ(al,q); SGN(q) = sgn;
216: }
217: *rp = q;
218: }
219:
220: void Pdceil(arg,rp)
221: NODE arg;
222: Q *rp;
223: {
224: L a;
225: unsigned int au,al;
226: int sgn;
227: Q q;
228: double d;
229:
230: if ( !ARG0(arg) ) {
231: *rp = 0;
232: return;
233: }
234: d = ceil(ToReal(ARG0(arg)));
235: if ( d < -9.223372036854775808e18 || d >= 9.223372036854775808e18 )
236: error("dceil : OverFlow");
1.4 noro 237: if ( !d ) {
238: *rp = 0;
239: return;
240: }
1.1 noro 241: a = (L)d;
242: if ( a < 0 ) {
243: sgn = -1;
244: a = -a;
245: } else
246: sgn = 1;
1.7 ! noro 247: #if defined(i386) || defined(__alpha) || defined(VISUAL) || defined(__x86_64)
1.1 noro 248: au = ((unsigned int *)&a)[1];
249: al = ((unsigned int *)&a)[0];
250: #else
251: al = ((unsigned int *)&a)[1];
252: au = ((unsigned int *)&a)[0];
253: #endif
254: if ( au ) {
255: NEWQ(q); SGN(q) = sgn; NM(q)=NALLOC(2); DN(q)=0;
256: PL(NM(q))=2; BD(NM(q))[0]=al; BD(NM(q))[1] = au;
257: } else {
258: UTOQ(al,q); SGN(q) = sgn;
259: }
260: *rp = q;
261: }
262:
263: void Pdrint(arg,rp)
264: NODE arg;
265: Q *rp;
266: {
267: L a;
268: unsigned int au,al;
269: int sgn;
270: Q q;
271: double d;
272:
273: if ( !ARG0(arg) ) {
274: *rp = 0;
275: return;
276: }
277: #if defined(VISUAL)
278: d = ToReal(ARG0(arg));
279: if ( d > 0 )
280: d = floor(d+0.5);
281: else
282: d = ceil(d-0.5);
283: #else
284: d = rint(ToReal(ARG0(arg)));
285: #endif
286: if ( d < -9.223372036854775808e18 || d >= 9.223372036854775808e18 )
287: error("drint : OverFlow");
288: a = (L)d;
289: if ( a < 0 ) {
290: sgn = -1;
291: a = -a;
292: } else
293: sgn = 1;
1.7 ! noro 294: #if defined(i386) || defined(__alpha) || defined(VISUAL) || defined(__x86_64)
1.1 noro 295: au = ((unsigned int *)&a)[1];
296: al = ((unsigned int *)&a)[0];
297: #else
298: al = ((unsigned int *)&a)[1];
299: au = ((unsigned int *)&a)[0];
300: #endif
301: if ( au ) {
302: NEWQ(q); SGN(q) = sgn; NM(q)=NALLOC(2); DN(q)=0;
303: PL(NM(q))=2; BD(NM(q))[0]=al; BD(NM(q))[1] = au;
1.6 noro 304: } else if ( al ) {
1.1 noro 305: UTOQ(al,q); SGN(q) = sgn;
1.6 noro 306: } else
307: q = 0;
1.1 noro 308: *rp = q;
309: }
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>
![[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