Eigen  3.4.90 (git rev 5a9f66fb35d03a4da9ef8976e67a61b30aa16dcf)
 
Loading...
Searching...
No Matches
ZVector/Complex.h
1// This file is part of Eigen, a lightweight C++ template library
2// for linear algebra.
3//
4// Copyright (C) 2010 Gael Guennebaud <[email protected]>
5// Copyright (C) 2016 Konstantinos Margaritis <[email protected]>
6//
7// This Source Code Form is subject to the terms of the Mozilla
8// Public License v. 2.0. If a copy of the MPL was not distributed
9// with this file, You can obtain one at http://mozilla.org/MPL/2.0/.
10
11#ifndef EIGEN_COMPLEX32_ZVECTOR_H
12#define EIGEN_COMPLEX32_ZVECTOR_H
13
14// IWYU pragma: private
15#include "../../InternalHeaderCheck.h"
16
17namespace Eigen {
18
19namespace internal {
20
21#if !defined(__ARCH__) || (defined(__ARCH__) && __ARCH__ >= 12)
22static Packet4ui p4ui_CONJ_XOR = {0x00000000, 0x80000000, 0x00000000,
23 0x80000000}; // vec_mergeh((Packet4ui)p4i_ZERO, (Packet4ui)p4f_MZERO);
24#endif
25
26static Packet2ul p2ul_CONJ_XOR1 =
27 (Packet2ul)vec_sld((Packet4ui)p2d_ZERO_, (Packet4ui)p2l_ZERO, 8); //{ 0x8000000000000000, 0x0000000000000000 };
28static Packet2ul p2ul_CONJ_XOR2 =
29 (Packet2ul)vec_sld((Packet4ui)p2l_ZERO, (Packet4ui)p2d_ZERO_, 8); //{ 0x8000000000000000, 0x0000000000000000 };
30
31struct Packet1cd {
32 EIGEN_STRONG_INLINE Packet1cd() {}
33 EIGEN_STRONG_INLINE explicit Packet1cd(const Packet2d& a) : v(a) {}
34 Packet2d v;
35};
36
37struct Packet2cf {
38 EIGEN_STRONG_INLINE Packet2cf() {}
39 EIGEN_STRONG_INLINE explicit Packet2cf(const Packet4f& a) : v(a) {}
40#if !defined(__ARCH__) || (defined(__ARCH__) && __ARCH__ < 12)
41 union {
42 Packet4f v;
43 Packet1cd cd[2];
44 };
45#else
46 Packet4f v;
47#endif
48};
49
50template <>
51struct packet_traits<std::complex<float> > : default_packet_traits {
52 typedef Packet2cf type;
53 typedef Packet2cf half;
54 enum {
55 Vectorizable = 1,
56 AlignedOnScalar = 1,
57 size = 2,
58
59 HasAdd = 1,
60 HasSub = 1,
61 HasMul = 1,
62 HasDiv = 1,
63 HasLog = 1,
64 HasExp = 1,
65 HasNegate = 1,
66 HasAbs = 0,
67 HasAbs2 = 0,
68 HasMin = 0,
69 HasMax = 0,
70 HasBlend = 1,
71 HasSetLinear = 0
72 };
73};
74
75template <>
76struct packet_traits<std::complex<double> > : default_packet_traits {
77 typedef Packet1cd type;
78 typedef Packet1cd half;
79 enum {
80 Vectorizable = 1,
81 AlignedOnScalar = 1,
82 size = 1,
83
84 HasAdd = 1,
85 HasSub = 1,
86 HasMul = 1,
87 HasDiv = 1,
88 HasLog = 1,
89 HasNegate = 1,
90 HasAbs = 0,
91 HasAbs2 = 0,
92 HasMin = 0,
93 HasMax = 0,
94 HasSetLinear = 0
95 };
96};
97
98template <>
99struct unpacket_traits<Packet2cf> {
100 typedef std::complex<float> type;
101 enum {
102 size = 2,
103 alignment = Aligned16,
104 vectorizable = true,
105 masked_load_available = false,
106 masked_store_available = false
107 };
108 typedef Packet2cf half;
109 typedef Packet4f as_real;
110};
111template <>
112struct unpacket_traits<Packet1cd> {
113 typedef std::complex<double> type;
114 enum {
115 size = 1,
116 alignment = Aligned16,
117 vectorizable = true,
118 masked_load_available = false,
119 masked_store_available = false
120 };
121 typedef Packet1cd half;
122 typedef Packet2d as_real;
123};
124
125/* Forward declaration */
126EIGEN_STRONG_INLINE void ptranspose(PacketBlock<Packet2cf, 2>& kernel);
127
128/* complex<double> first */
129template <>
130EIGEN_STRONG_INLINE Packet1cd pload<Packet1cd>(const std::complex<double>* from) {
131 EIGEN_DEBUG_ALIGNED_LOAD return Packet1cd(pload<Packet2d>((const double*)from));
132}
133template <>
134EIGEN_STRONG_INLINE Packet1cd ploadu<Packet1cd>(const std::complex<double>* from) {
135 EIGEN_DEBUG_UNALIGNED_LOAD return Packet1cd(ploadu<Packet2d>((const double*)from));
136}
137template <>
138EIGEN_STRONG_INLINE void pstore<std::complex<double> >(std::complex<double>* to, const Packet1cd& from) {
139 EIGEN_DEBUG_ALIGNED_STORE pstore((double*)to, from.v);
140}
141template <>
142EIGEN_STRONG_INLINE void pstoreu<std::complex<double> >(std::complex<double>* to, const Packet1cd& from) {
143 EIGEN_DEBUG_UNALIGNED_STORE pstoreu((double*)to, from.v);
144}
145
146template <>
147EIGEN_STRONG_INLINE Packet1cd
148pset1<Packet1cd>(const std::complex<double>& from) { /* here we really have to use unaligned loads :( */
149 return ploadu<Packet1cd>(&from);
150}
151
152template <>
153EIGEN_DEVICE_FUNC inline Packet1cd pgather<std::complex<double>, Packet1cd>(const std::complex<double>* from,
154 Index stride EIGEN_UNUSED) {
155 return pload<Packet1cd>(from);
156}
157template <>
158EIGEN_DEVICE_FUNC inline void pscatter<std::complex<double>, Packet1cd>(std::complex<double>* to, const Packet1cd& from,
159 Index stride EIGEN_UNUSED) {
160 pstore<std::complex<double> >(to, from);
161}
162template <>
163EIGEN_STRONG_INLINE Packet1cd padd<Packet1cd>(const Packet1cd& a, const Packet1cd& b) {
164 return Packet1cd(a.v + b.v);
165}
166template <>
167EIGEN_STRONG_INLINE Packet1cd psub<Packet1cd>(const Packet1cd& a, const Packet1cd& b) {
168 return Packet1cd(a.v - b.v);
169}
170template <>
171EIGEN_STRONG_INLINE Packet1cd pnegate(const Packet1cd& a) {
172 return Packet1cd(pnegate(Packet2d(a.v)));
173}
174template <>
175EIGEN_STRONG_INLINE Packet1cd pconj(const Packet1cd& a) {
176 return Packet1cd((Packet2d)vec_xor((Packet2d)a.v, (Packet2d)p2ul_CONJ_XOR2));
177}
178template <>
179EIGEN_STRONG_INLINE Packet1cd pmul<Packet1cd>(const Packet1cd& a, const Packet1cd& b) {
180 Packet2d a_re, a_im, v1, v2;
181
182 // Permute and multiply the real parts of a and b
183 a_re = vec_perm(a.v, a.v, p16uc_PSET64_HI);
184 // Get the imaginary parts of a
185 a_im = vec_perm(a.v, a.v, p16uc_PSET64_LO);
186 // multiply a_re * b
187 v1 = vec_madd(a_re, b.v, p2d_ZERO);
188 // multiply a_im * b and get the conjugate result
189 v2 = vec_madd(a_im, b.v, p2d_ZERO);
190 v2 = (Packet2d)vec_sld((Packet4ui)v2, (Packet4ui)v2, 8);
191 v2 = (Packet2d)vec_xor((Packet2d)v2, (Packet2d)p2ul_CONJ_XOR1);
192
193 return Packet1cd(v1 + v2);
194}
195template <>
196EIGEN_STRONG_INLINE Packet1cd pand<Packet1cd>(const Packet1cd& a, const Packet1cd& b) {
197 return Packet1cd(vec_and(a.v, b.v));
198}
199template <>
200EIGEN_STRONG_INLINE Packet1cd por<Packet1cd>(const Packet1cd& a, const Packet1cd& b) {
201 return Packet1cd(vec_or(a.v, b.v));
202}
203template <>
204EIGEN_STRONG_INLINE Packet1cd pxor<Packet1cd>(const Packet1cd& a, const Packet1cd& b) {
205 return Packet1cd(vec_xor(a.v, b.v));
206}
207template <>
208EIGEN_STRONG_INLINE Packet1cd pandnot<Packet1cd>(const Packet1cd& a, const Packet1cd& b) {
209 return Packet1cd(vec_and(a.v, vec_nor(b.v, b.v)));
210}
211template <>
212EIGEN_STRONG_INLINE Packet1cd ploaddup<Packet1cd>(const std::complex<double>* from) {
213 return pset1<Packet1cd>(*from);
214}
215template <>
216EIGEN_STRONG_INLINE Packet1cd pcmp_eq(const Packet1cd& a, const Packet1cd& b) {
217 Packet2d eq = vec_cmpeq(a.v, b.v);
218 Packet2d tmp = {eq[1], eq[0]};
219 return (Packet1cd)pand<Packet2d>(eq, tmp);
220}
221
222template <>
223EIGEN_STRONG_INLINE void prefetch<std::complex<double> >(const std::complex<double>* addr) {
224 EIGEN_ZVECTOR_PREFETCH(addr);
225}
226
227template <>
228EIGEN_STRONG_INLINE std::complex<double> pfirst<Packet1cd>(const Packet1cd& a) {
229 EIGEN_ALIGN16 std::complex<double> res;
230 pstore<std::complex<double> >(&res, a);
231
232 return res;
233}
234
235template <>
236EIGEN_STRONG_INLINE Packet1cd preverse(const Packet1cd& a) {
237 return a;
238}
239template <>
240EIGEN_STRONG_INLINE std::complex<double> predux<Packet1cd>(const Packet1cd& a) {
241 return pfirst(a);
242}
243template <>
244EIGEN_STRONG_INLINE std::complex<double> predux_mul<Packet1cd>(const Packet1cd& a) {
245 return pfirst(a);
246}
247EIGEN_MAKE_CONJ_HELPER_CPLX_REAL(Packet1cd, Packet2d)
248
249template <>
250EIGEN_STRONG_INLINE Packet1cd pdiv<Packet1cd>(const Packet1cd& a, const Packet1cd& b) {
251 return pdiv_complex(a, b);
252}
253
254template <>
255EIGEN_STRONG_INLINE Packet1cd plog<Packet1cd>(const Packet1cd& a, const Packet1cd& b) {
256 return plog_complex(a, b);
257}
258
259EIGEN_STRONG_INLINE Packet1cd pcplxflip /*<Packet1cd>*/ (const Packet1cd& x) {
260 return Packet1cd(preverse(Packet2d(x.v)));
261}
262
263EIGEN_STRONG_INLINE void ptranspose(PacketBlock<Packet1cd, 2>& kernel) {
264 Packet2d tmp = vec_perm(kernel.packet[0].v, kernel.packet[1].v, p16uc_TRANSPOSE64_HI);
265 kernel.packet[1].v = vec_perm(kernel.packet[0].v, kernel.packet[1].v, p16uc_TRANSPOSE64_LO);
266 kernel.packet[0].v = tmp;
267}
268
269/* complex<float> follows */
270template <>
271EIGEN_STRONG_INLINE Packet2cf pload<Packet2cf>(const std::complex<float>* from) {
272 EIGEN_DEBUG_ALIGNED_LOAD return Packet2cf(pload<Packet4f>((const float*)from));
273}
274template <>
275EIGEN_STRONG_INLINE Packet2cf ploadu<Packet2cf>(const std::complex<float>* from) {
276 EIGEN_DEBUG_UNALIGNED_LOAD return Packet2cf(ploadu<Packet4f>((const float*)from));
277}
278template <>
279EIGEN_STRONG_INLINE void pstore<std::complex<float> >(std::complex<float>* to, const Packet2cf& from) {
280 EIGEN_DEBUG_ALIGNED_STORE pstore((float*)to, from.v);
281}
282template <>
283EIGEN_STRONG_INLINE void pstoreu<std::complex<float> >(std::complex<float>* to, const Packet2cf& from) {
284 EIGEN_DEBUG_UNALIGNED_STORE pstoreu((float*)to, from.v);
285}
286
287template <>
288EIGEN_STRONG_INLINE std::complex<float> pfirst<Packet2cf>(const Packet2cf& a) {
289 EIGEN_ALIGN16 std::complex<float> res[2];
290 pstore<std::complex<float> >(res, a);
291
292 return res[0];
293}
294
295#if !defined(__ARCH__) || (defined(__ARCH__) && __ARCH__ < 12)
296template <>
297EIGEN_STRONG_INLINE Packet2cf pset1<Packet2cf>(const std::complex<float>& from) {
298 Packet2cf res;
299 res.cd[0] = Packet1cd(vec_ld2f((const float*)&from));
300 res.cd[1] = res.cd[0];
301 return res;
302}
303#else
304template <>
305EIGEN_STRONG_INLINE Packet2cf pset1<Packet2cf>(const std::complex<float>& from) {
306 Packet2cf res;
307 if ((std::ptrdiff_t(&from) % 16) == 0)
308 res.v = pload<Packet4f>((const float*)&from);
309 else
310 res.v = ploadu<Packet4f>((const float*)&from);
311 res.v = vec_perm(res.v, res.v, p16uc_PSET64_HI);
312 return res;
313}
314#endif
315
316template <>
317EIGEN_DEVICE_FUNC inline Packet2cf pgather<std::complex<float>, Packet2cf>(const std::complex<float>* from,
318 Index stride) {
319 EIGEN_ALIGN16 std::complex<float> af[2];
320 af[0] = from[0 * stride];
321 af[1] = from[1 * stride];
322 return pload<Packet2cf>(af);
323}
324template <>
325EIGEN_DEVICE_FUNC inline void pscatter<std::complex<float>, Packet2cf>(std::complex<float>* to, const Packet2cf& from,
326 Index stride) {
327 EIGEN_ALIGN16 std::complex<float> af[2];
328 pstore<std::complex<float> >((std::complex<float>*)af, from);
329 to[0 * stride] = af[0];
330 to[1 * stride] = af[1];
331}
332
333template <>
334EIGEN_STRONG_INLINE Packet2cf padd<Packet2cf>(const Packet2cf& a, const Packet2cf& b) {
335 return Packet2cf(padd<Packet4f>(a.v, b.v));
336}
337template <>
338EIGEN_STRONG_INLINE Packet2cf psub<Packet2cf>(const Packet2cf& a, const Packet2cf& b) {
339 return Packet2cf(psub<Packet4f>(a.v, b.v));
340}
341template <>
342EIGEN_STRONG_INLINE Packet2cf pnegate(const Packet2cf& a) {
343 return Packet2cf(pnegate(Packet4f(a.v)));
344}
345
346template <>
347EIGEN_STRONG_INLINE Packet2cf pand<Packet2cf>(const Packet2cf& a, const Packet2cf& b) {
348 return Packet2cf(pand<Packet4f>(a.v, b.v));
349}
350template <>
351EIGEN_STRONG_INLINE Packet2cf por<Packet2cf>(const Packet2cf& a, const Packet2cf& b) {
352 return Packet2cf(por<Packet4f>(a.v, b.v));
353}
354template <>
355EIGEN_STRONG_INLINE Packet2cf pxor<Packet2cf>(const Packet2cf& a, const Packet2cf& b) {
356 return Packet2cf(pxor<Packet4f>(a.v, b.v));
357}
358template <>
359EIGEN_STRONG_INLINE Packet2cf pandnot<Packet2cf>(const Packet2cf& a, const Packet2cf& b) {
360 return Packet2cf(pandnot<Packet4f>(a.v, b.v));
361}
362
363template <>
364EIGEN_STRONG_INLINE Packet2cf ploaddup<Packet2cf>(const std::complex<float>* from) {
365 return pset1<Packet2cf>(*from);
366}
367
368template <>
369EIGEN_STRONG_INLINE void prefetch<std::complex<float> >(const std::complex<float>* addr) {
370 EIGEN_ZVECTOR_PREFETCH(addr);
371}
372
373#if !defined(__ARCH__) || (defined(__ARCH__) && __ARCH__ < 12)
374
375template <>
376EIGEN_STRONG_INLINE Packet2cf pcmp_eq(const Packet2cf& a, const Packet2cf& b) {
377 Packet4f eq = pcmp_eq<Packet4f>(a.v, b.v);
378 Packet2cf res;
379 Packet2d tmp1 = {eq.v4f[0][1], eq.v4f[0][0]};
380 Packet2d tmp2 = {eq.v4f[1][1], eq.v4f[1][0]};
381 res.v.v4f[0] = pand<Packet2d>(eq.v4f[0], tmp1);
382 res.v.v4f[1] = pand<Packet2d>(eq.v4f[1], tmp2);
383 return res;
384}
385
386template <>
387EIGEN_STRONG_INLINE Packet2cf pconj(const Packet2cf& a) {
388 Packet2cf res;
389 res.v.v4f[0] = pconj(Packet1cd(reinterpret_cast<Packet2d>(a.v.v4f[0]))).v;
390 res.v.v4f[1] = pconj(Packet1cd(reinterpret_cast<Packet2d>(a.v.v4f[1]))).v;
391 return res;
392}
393
394template <>
395EIGEN_STRONG_INLINE Packet2cf pmul<Packet2cf>(const Packet2cf& a, const Packet2cf& b) {
396 Packet2cf res;
397 res.v.v4f[0] =
398 pmul(Packet1cd(reinterpret_cast<Packet2d>(a.v.v4f[0])), Packet1cd(reinterpret_cast<Packet2d>(b.v.v4f[0]))).v;
399 res.v.v4f[1] =
400 pmul(Packet1cd(reinterpret_cast<Packet2d>(a.v.v4f[1])), Packet1cd(reinterpret_cast<Packet2d>(b.v.v4f[1]))).v;
401 return res;
402}
403
404template <>
405EIGEN_STRONG_INLINE Packet2cf preverse(const Packet2cf& a) {
406 Packet2cf res;
407 res.cd[0] = a.cd[1];
408 res.cd[1] = a.cd[0];
409 return res;
410}
411
412template <>
413EIGEN_STRONG_INLINE std::complex<float> predux<Packet2cf>(const Packet2cf& a) {
414 std::complex<float> res;
415 Packet1cd b = padd<Packet1cd>(a.cd[0], a.cd[1]);
416 vec_st2f(b.v, (float*)&res);
417 return res;
418}
419
420template <>
421EIGEN_STRONG_INLINE std::complex<float> predux_mul<Packet2cf>(const Packet2cf& a) {
422 std::complex<float> res;
423 Packet1cd b = pmul<Packet1cd>(a.cd[0], a.cd[1]);
424 vec_st2f(b.v, (float*)&res);
425 return res;
426}
427
428EIGEN_MAKE_CONJ_HELPER_CPLX_REAL(Packet2cf, Packet4f)
429
430template <>
431EIGEN_STRONG_INLINE Packet2cf pdiv<Packet2cf>(const Packet2cf& a, const Packet2cf& b) {
432 return pdiv_complex(a, b);
433}
434
435template <>
436EIGEN_STRONG_INLINE Packet2cf plog<Packet2cf>(const Packet2cf& a, const Packet2cf& b) {
437 return plog_complex(a, b);
438}
439
440template <>
441EIGEN_STRONG_INLINE Packet2cf pexp<Packet2cf>(const Packet2cf& a, const Packet2cf& b) {
442 return pexp_complex(a, b);
443}
444
445EIGEN_STRONG_INLINE Packet2cf pcplxflip /*<Packet2cf>*/ (const Packet2cf& x) {
446 Packet2cf res;
447 res.cd[0] = pcplxflip(x.cd[0]);
448 res.cd[1] = pcplxflip(x.cd[1]);
449 return res;
450}
451
452EIGEN_STRONG_INLINE void ptranspose(PacketBlock<Packet2cf, 2>& kernel) {
453 Packet1cd tmp = kernel.packet[0].cd[1];
454 kernel.packet[0].cd[1] = kernel.packet[1].cd[0];
455 kernel.packet[1].cd[0] = tmp;
456}
457
458template <>
459EIGEN_STRONG_INLINE Packet2cf pblend(const Selector<2>& ifPacket, const Packet2cf& thenPacket,
460 const Packet2cf& elsePacket) {
461 Packet2cf result;
462 const Selector<4> ifPacket4 = {ifPacket.select[0], ifPacket.select[0], ifPacket.select[1], ifPacket.select[1]};
463 result.v = pblend<Packet4f>(ifPacket4, thenPacket.v, elsePacket.v);
464 return result;
465}
466#else
467template <>
468EIGEN_STRONG_INLINE Packet2cf pcmp_eq(const Packet2cf& a, const Packet2cf& b) {
469 Packet4f eq = vec_cmpeq(a.v, b.v);
470 Packet4f tmp = {eq[1], eq[0], eq[3], eq[2]};
471 return (Packet2cf)pand<Packet4f>(eq, tmp);
472}
473template <>
474EIGEN_STRONG_INLINE Packet2cf pconj(const Packet2cf& a) {
475 return Packet2cf(pxor<Packet4f>(a.v, reinterpret_cast<Packet4f>(p4ui_CONJ_XOR)));
476}
477template <>
478EIGEN_STRONG_INLINE Packet2cf pmul<Packet2cf>(const Packet2cf& a, const Packet2cf& b) {
479 Packet4f a_re, a_im, prod, prod_im;
480
481 // Permute and multiply the real parts of a and b
482 a_re = vec_perm(a.v, a.v, p16uc_PSET32_WODD);
483
484 // Get the imaginary parts of a
485 a_im = vec_perm(a.v, a.v, p16uc_PSET32_WEVEN);
486
487 // multiply a_im * b and get the conjugate result
488 prod_im = a_im * b.v;
489 prod_im = pxor<Packet4f>(prod_im, reinterpret_cast<Packet4f>(p4ui_CONJ_XOR));
490 // permute back to a proper order
491 prod_im = vec_perm(prod_im, prod_im, p16uc_COMPLEX32_REV);
492
493 // multiply a_re * b, add prod_im
494 prod = pmadd<Packet4f>(a_re, b.v, prod_im);
495
496 return Packet2cf(prod);
497}
498
499template <>
500EIGEN_STRONG_INLINE Packet2cf preverse(const Packet2cf& a) {
501 Packet4f rev_a;
502 rev_a = vec_perm(a.v, a.v, p16uc_COMPLEX32_REV2);
503 return Packet2cf(rev_a);
504}
505
506template <>
507EIGEN_STRONG_INLINE std::complex<float> predux<Packet2cf>(const Packet2cf& a) {
508 Packet4f b;
509 b = vec_sld(a.v, a.v, 8);
510 b = padd<Packet4f>(a.v, b);
511 return pfirst<Packet2cf>(Packet2cf(b));
512}
513
514template <>
515EIGEN_STRONG_INLINE std::complex<float> predux_mul<Packet2cf>(const Packet2cf& a) {
516 Packet4f b;
517 Packet2cf prod;
518 b = vec_sld(a.v, a.v, 8);
519 prod = pmul<Packet2cf>(a, Packet2cf(b));
520
521 return pfirst<Packet2cf>(prod);
522}
523
524EIGEN_MAKE_CONJ_HELPER_CPLX_REAL(Packet2cf, Packet4f)
525
526template <>
527EIGEN_STRONG_INLINE Packet2cf pdiv<Packet2cf>(const Packet2cf& a, const Packet2cf& b) {
528 return pdiv_complex(a, b);
529}
530
531template <>
532EIGEN_STRONG_INLINE Packet2cf pcplxflip<Packet2cf>(const Packet2cf& x) {
533 return Packet2cf(vec_perm(x.v, x.v, p16uc_COMPLEX32_REV));
534}
535
536EIGEN_STRONG_INLINE void ptranspose(PacketBlock<Packet2cf, 2>& kernel) {
537 Packet4f tmp = vec_perm(kernel.packet[0].v, kernel.packet[1].v, p16uc_TRANSPOSE64_HI);
538 kernel.packet[1].v = vec_perm(kernel.packet[0].v, kernel.packet[1].v, p16uc_TRANSPOSE64_LO);
539 kernel.packet[0].v = tmp;
540}
541
542template <>
543EIGEN_STRONG_INLINE Packet2cf pblend(const Selector<2>& ifPacket, const Packet2cf& thenPacket,
544 const Packet2cf& elsePacket) {
545 Packet2cf result;
546 result.v = reinterpret_cast<Packet4f>(
547 pblend<Packet2d>(ifPacket, reinterpret_cast<Packet2d>(thenPacket.v), reinterpret_cast<Packet2d>(elsePacket.v)));
548 return result;
549}
550#endif
551
552} // end namespace internal
553
554} // end namespace Eigen
555
556#endif // EIGEN_COMPLEX32_ZVECTOR_H
Eigen::Aligned16
@ Aligned16
Definition Constants.h:237
Eigen
Namespace containing all symbols from the Eigen library.
Definition Core:137