|
72
|
1 module dmd.ComplexExp;
|
|
|
2
|
|
114
|
3 import dmd.common;
|
|
72
|
4 import dmd.Expression;
|
|
|
5 import dmd.InterState;
|
|
|
6 import dmd.Type;
|
|
|
7 import dmd.OutBuffer;
|
|
|
8 import dmd.Loc;
|
|
|
9 import dmd.Scope;
|
|
|
10 import dmd.IRState;
|
|
|
11 import dmd.HdrGenState;
|
|
0
|
12 import dmd.Type;
|
|
|
13 import dmd.TOK;
|
|
|
14 import dmd.TY;
|
|
|
15 import dmd.Port;
|
|
|
16 import dmd.Complex;
|
|
129
|
17 import dmd.expression.Util;
|
|
72
|
18
|
|
|
19 import dmd.backend.dt_t;
|
|
0
|
20 import dmd.backend.elem;
|
|
|
21 import dmd.backend.Util;
|
|
|
22 import dmd.backend.TYM;
|
|
|
23 import dmd.backend.mTY;
|
|
72
|
24
|
|
0
|
25 class ComplexExp : Expression
|
|
|
26 {
|
|
|
27 Complex!(real) value;
|
|
|
28
|
|
|
29 this(Loc loc, Complex!(real) value, Type type)
|
|
72
|
30 {
|
|
0
|
31 super(loc, TOK.TOKcomplex80, ComplexExp.sizeof);
|
|
|
32 this.value = value;
|
|
|
33 this.type = type;
|
|
|
34 //printf("ComplexExp.ComplexExp(%s)\n", toChars());
|
|
|
35 }
|
|
|
36
|
|
72
|
37 override bool equals(Object o)
|
|
0
|
38 {
|
|
|
39 assert(false);
|
|
|
40 }
|
|
|
41
|
|
72
|
42 override Expression semantic(Scope sc)
|
|
0
|
43 {
|
|
|
44 if (!type)
|
|
|
45 type = Type.tcomplex80;
|
|
|
46 else
|
|
|
47 type = type.semantic(loc, sc);
|
|
|
48 return this;
|
|
|
49 }
|
|
|
50
|
|
72
|
51 override Expression interpret(InterState istate)
|
|
0
|
52 {
|
|
|
53 assert(false);
|
|
|
54 }
|
|
|
55
|
|
72
|
56 override string toChars()
|
|
0
|
57 {
|
|
|
58 assert(false);
|
|
|
59 }
|
|
|
60
|
|
72
|
61 override ulong toInteger()
|
|
0
|
62 {
|
|
|
63 return cast(ulong) toReal();
|
|
|
64 }
|
|
|
65
|
|
72
|
66 override ulong toUInteger()
|
|
0
|
67 {
|
|
|
68 return cast(long) toReal();
|
|
|
69 }
|
|
|
70
|
|
72
|
71 override real toReal()
|
|
0
|
72 {
|
|
|
73 return value.re;
|
|
|
74 }
|
|
|
75
|
|
72
|
76 override real toImaginary()
|
|
0
|
77 {
|
|
|
78 return value.im;
|
|
|
79 }
|
|
|
80
|
|
72
|
81 override Complex!(real) toComplex()
|
|
0
|
82 {
|
|
|
83 return value;
|
|
|
84 }
|
|
|
85
|
|
72
|
86 override Expression castTo(Scope sc, Type t)
|
|
0
|
87 {
|
|
|
88 Expression e = this;
|
|
|
89 if (type != t)
|
|
|
90 {
|
|
|
91 if (type.iscomplex() && t.iscomplex())
|
|
|
92 {
|
|
|
93 e = copy();
|
|
|
94 e.type = t;
|
|
|
95 }
|
|
|
96 else
|
|
|
97 e = Expression.castTo(sc, t);
|
|
|
98 }
|
|
|
99 return e;
|
|
|
100 }
|
|
|
101
|
|
72
|
102 override int isConst()
|
|
0
|
103 {
|
|
53
|
104 return 1;
|
|
0
|
105 }
|
|
|
106
|
|
72
|
107 override bool isBool(bool result)
|
|
0
|
108 {
|
|
123
|
109 if (result)
|
|
|
110 return value != Complex!(real).zero;
|
|
|
111 else
|
|
|
112 return value == Complex!(real).zero;
|
|
0
|
113 }
|
|
|
114
|
|
72
|
115 override void toCBuffer(OutBuffer buf, HdrGenState* hgs)
|
|
0
|
116 {
|
|
|
117 assert(false);
|
|
|
118 }
|
|
|
119
|
|
72
|
120 override void toMangleBuffer(OutBuffer buf)
|
|
0
|
121 {
|
|
129
|
122 buf.writeByte('c');
|
|
|
123 real r = toReal();
|
|
|
124 realToMangleBuffer(buf, r);
|
|
|
125 buf.writeByte('c'); // separate the two
|
|
|
126 r = toImaginary();
|
|
|
127 realToMangleBuffer(buf, r);
|
|
0
|
128 }
|
|
129
|
129
|
|
0
|
130 version (_DH) {
|
|
|
131 OutBuffer hexp;
|
|
|
132 }
|
|
|
133
|
|
72
|
134 override elem* toElem(IRState* irs)
|
|
0
|
135 {
|
|
|
136 eve c;
|
|
|
137 tym_t ty;
|
|
|
138
|
|
|
139 //printf("ComplexExp.toElem(%p) %s\n", this, toChars());
|
|
|
140
|
|
|
141 ///memset(&c, 0, c.sizeof);
|
|
|
142
|
|
|
143 ty = type.totym();
|
|
|
144 switch (tybasic(ty))
|
|
|
145 {
|
|
|
146 case TYcfloat:
|
|
|
147 {
|
|
|
148 c.Vcfloat.re = cast(float) value.re;
|
|
|
149 if (Port.isSignallingNan(value.re)) {
|
|
|
150 (cast(uint*)&c.Vcfloat.re)[0] &= 0xFFBFFFFFL;
|
|
|
151 std.stdio.writeln("float.re is snan");
|
|
|
152 }
|
|
|
153 c.Vcfloat.im = cast(float) value.im;
|
|
|
154 if (Port.isSignallingNan(value.im)) {
|
|
|
155 (cast(uint*)&c.Vcfloat.im)[0] &= 0xFFBFFFFFL;
|
|
|
156 std.stdio.writeln("float.im is snan");
|
|
|
157 }
|
|
|
158 break;
|
|
|
159 }
|
|
|
160
|
|
|
161 case TYcdouble:
|
|
|
162 {
|
|
|
163 c.Vcdouble.re = cast(double) value.re;
|
|
|
164 if (Port.isSignallingNan(value.re)) {
|
|
|
165 std.stdio.writeln("double.re is snan");
|
|
|
166 (cast(uint*)&c.Vcdouble.re)[1] &= 0xFFF7FFFFL;
|
|
|
167 }
|
|
|
168 c.Vcdouble.im = cast(double) value.im;
|
|
|
169 if (Port.isSignallingNan(value.im)) {
|
|
|
170 (cast(uint*)&c.Vcdouble.im)[1] &= 0xFFF7FFFFL;
|
|
|
171 std.stdio.writeln("double.im is snan");
|
|
|
172 }
|
|
|
173 break;
|
|
|
174 }
|
|
|
175
|
|
|
176 case TYcldouble:
|
|
|
177 {
|
|
|
178 static if (true) {
|
|
|
179 c.Vcldouble = value;
|
|
|
180 } else {
|
|
|
181 {
|
|
|
182 ushort* p = cast(ushort*)&c.Vcldouble;
|
|
|
183 for (int i = 0; i < (LNGDBLSIZE*2)/2; i++) printf("%04x ", p[i]);
|
|
|
184 printf("\n");
|
|
|
185 }
|
|
|
186 c.Vcldouble.im = im;
|
|
|
187 {
|
|
|
188 ushort* p = cast(ushort*)&c.Vcldouble;
|
|
|
189 for (int i = 0; i < (LNGDBLSIZE*2)/2; i++) printf("%04x ", p[i]);
|
|
|
190 printf("\n");
|
|
|
191 }
|
|
|
192 c.Vcldouble.re = re;
|
|
|
193 {
|
|
|
194 ushort* p = cast(ushort*)&c.Vcldouble;
|
|
|
195 for (int i = 0; i < (LNGDBLSIZE*2)/2; i++) printf("%04x ", p[i]);
|
|
|
196 printf("\n");
|
|
|
197 }
|
|
|
198 }
|
|
|
199 break;
|
|
|
200 }
|
|
|
201
|
|
|
202 default:
|
|
|
203 assert(0);
|
|
|
204 }
|
|
|
205 return el_const(ty, &c);
|
|
|
206 }
|
|
|
207
|
|
|
208 static private char[6] zeropad;
|
|
|
209
|
|
72
|
210 override dt_t** toDt(dt_t** pdt)
|
|
0
|
211 {
|
|
|
212 //printf("ComplexExp.toDt() '%s'\n", toChars());
|
|
|
213 float fvalue;
|
|
|
214 double dvalue;
|
|
|
215 real evalue;
|
|
|
216
|
|
|
217 switch (type.toBasetype().ty)
|
|
|
218 {
|
|
|
219 case Tcomplex32:
|
|
|
220 fvalue = value.re;
|
|
|
221 pdt = dtnbytes(pdt,4,cast(char*)&fvalue);
|
|
|
222 fvalue = value.im;
|
|
|
223 pdt = dtnbytes(pdt,4,cast(char*)&fvalue);
|
|
|
224 break;
|
|
|
225
|
|
|
226 case Tcomplex64:
|
|
|
227 dvalue = value.re;
|
|
|
228 pdt = dtnbytes(pdt,8,cast(char*)&dvalue);
|
|
|
229 dvalue = value.im;
|
|
|
230 pdt = dtnbytes(pdt,8,cast(char*)&dvalue);
|
|
|
231 break;
|
|
|
232
|
|
|
233 case Tcomplex80:
|
|
|
234 evalue = value.re;
|
|
|
235 pdt = dtnbytes(pdt,REALSIZE - REALPAD,cast(char*)&evalue);
|
|
|
236 pdt = dtnbytes(pdt,REALPAD,zeropad.ptr);
|
|
|
237 evalue = value.im;
|
|
|
238 pdt = dtnbytes(pdt,REALSIZE - REALPAD, cast(char*)&evalue);
|
|
|
239 pdt = dtnbytes(pdt,REALPAD,zeropad.ptr);
|
|
|
240 break;
|
|
|
241
|
|
|
242 default:
|
|
|
243 assert(0);
|
|
|
244 break;
|
|
|
245 }
|
|
|
246 return pdt;
|
|
|
247 }
|
|
|
248 }
|
|
|
249
|
