![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to cplxnum.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 cplxnum.c CVS log | Up to [local] / OpenXM_contrib2 / asir2000 / builtin |
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.3 ! noro 48: * $OpenXM: OpenXM_contrib2/asir2000/builtin/cplxnum.c,v 1.2 2000/08/21 08:31:19 noro Exp $
1.2 noro 49: */
1.1 noro 50: #include "ca.h"
51: #include "parse.h"
52:
53: void Pconj(),Preal(),Pimag();
54:
55: #if defined(THINK_C)
56: void cplx_conj(Obj,Obj *);
57: void cplx_real(Obj,Obj *);
58: void cplx_imag(Obj,Obj *);
59: #else
60: void cplx_conj();
61: void cplx_real();
62: void cplx_imag();
63: #endif
64:
65: struct ftab cplx_tab[] = {
66: {"conj",Pconj,1},
67: {"real",Preal,1},
68: {"imag",Pimag,1},
69: {0,0,0},
70: };
71:
72: void Pconj(arg,rp)
73: NODE arg;
74: Obj *rp;
75: {
76: cplx_conj((Obj)ARG0(arg),rp);
77: }
78:
79: void Preal(arg,rp)
80: NODE arg;
81: Obj *rp;
82: {
83: cplx_real((Obj)ARG0(arg),rp);
84: }
85:
86: void Pimag(arg,rp)
87: NODE arg;
88: Obj *rp;
89: {
90: cplx_imag((Obj)ARG0(arg),rp);
91: }
92:
93: void cplx_conj(p,r)
94: Obj p;
95: Obj *r;
96: {
97: C c;
98: DCP dc,dcr,dcr0;
99: P t;
100:
101: if ( !p )
102: *r = 0;
103: else
104: switch ( OID(p) ) {
105: case O_N:
106: if ( NID((Num)p) <= N_B )
107: *r = p;
108: else {
109: NEWC(c); c->r = ((C)p)->r; chsgnnum(((C)p)->i,&c->i);
110: *r = (Obj)c;
111: }
112: break;
113: case O_P:
114: for ( dcr0 = 0, dc = DC((P)p); dc; dc = NEXT(dc) ) {
115: NEXTDC(dcr0,dcr); cplx_conj((Obj)COEF(dc),(Obj *)&COEF(dcr));
116: DEG(dcr) = DEG(dc);
117: }
118: NEXT(dcr) = 0; MKP(VR((P)p),dcr0,t); *r = (Obj)t;
119: break;
120: default:
121: error("cplx_conj : not implemented"); break;
122: }
123: }
124:
125: void cplx_real(p,r)
126: Obj p;
127: Obj *r;
128: {
129: DCP dc,dcr,dcr0;
130: P t;
131:
132: if ( !p )
133: *r = 0;
134: else
135: switch ( OID(p) ) {
136: case O_N:
137: if ( NID((Num)p) <= N_B )
138: *r = p;
139: else
140: *r = (Obj)((C)p)->r;
141: break;
142: case O_P:
143: for ( dcr0 = 0, dc = DC((P)p); dc; dc = NEXT(dc) ) {
144: cplx_real((Obj)COEF(dc),(Obj *)&t);
145: if ( t ) {
146: NEXTDC(dcr0,dcr); DEG(dcr) = DEG(dc); COEF(dcr) = t;
147: }
148: }
149: if ( !dcr0 )
150: *r = 0;
151: else {
152: NEXT(dcr) = 0; MKP(VR((P)p),dcr0,t); *r = (Obj)t;
153: }
154: break;
155: default:
156: error("cplx_real : not implemented"); break;
157: }
158: }
159:
160: void cplx_imag(p,r)
161: Obj p;
162: Obj *r;
163: {
164: DCP dc,dcr,dcr0;
165: P t;
166:
167: if ( !p )
168: *r = 0;
169: else
170: switch ( OID(p) ) {
171: case O_N:
172: if ( NID((Num)p) <= N_B )
173: *r = 0;
174: else
175: *r = (Obj)((C)p)->i;
176: break;
177: case O_P:
178: for ( dcr0 = 0, dc = DC((P)p); dc; dc = NEXT(dc) ) {
179: cplx_imag((Obj)COEF(dc),(Obj *)&t);
180: if ( t ) {
181: NEXTDC(dcr0,dcr); DEG(dcr) = DEG(dc); COEF(dcr) = t;
182: }
183: }
184: if ( !dcr0 )
185: *r = 0;
186: else {
187: NEXT(dcr) = 0; MKP(VR((P)p),dcr0,t); *r = (Obj)t;
188: }
189: break;
190: default:
191: error("cplx_imag : not implemented"); break;
192: }
193: }