| 598 | A_EQ_CB_MUL_C (CH(i,k,u2),WA(u2-1,i),db) \ |
| 599 | } |
| 600 | NOINLINE static void pass5f (size_t ido, size_t l1, const cmplx * restrict cc, |
| 601 | cmplx * restrict ch, const cmplx * restrict wa) |
| 602 | { |
| 603 | const size_t cdim=5; |
| 604 | const double tw1r= 0.3090169943749474241, |
| 605 | tw1i= -0.95105651629515357212, |
| 606 | tw2r= -0.8090169943749474241, |
| 607 | tw2i= -0.58778525229247312917; |
| 608 | |
| 609 | if (ido==1) |
| 610 | for (size_t k=0; k<l1; ++k) |
| 611 | { |
| 612 | PREP5(0) |
| 613 | PARTSTEP5a(1,4,tw1r,tw2r,+tw1i,+tw2i) |
| 614 | PARTSTEP5a(2,3,tw2r,tw1r,+tw2i,-tw1i) |
| 615 | } |
| 616 | else |
| 617 | for (size_t k=0; k<l1; ++k) |
| 618 | { |
| 619 | { |
| 620 | PREP5(0) |
| 621 | PARTSTEP5a(1,4,tw1r,tw2r,+tw1i,+tw2i) |
| 622 | PARTSTEP5a(2,3,tw2r,tw1r,+tw2i,-tw1i) |
| 623 | } |
| 624 | for (size_t i=1; i<ido; ++i) |
| 625 | { |
| 626 | PREP5(i) |
| 627 | PARTSTEP5f(1,4,tw1r,tw2r,+tw1i,+tw2i) |
| 628 | PARTSTEP5f(2,3,tw2r,tw1r,+tw2i,-tw1i) |
| 629 | } |
| 630 | } |
| 631 | } |
| 632 | |
| 633 | #define PREP7(idx) \ |
| 634 | cmplx t1 = CC(idx,0,k), t2, t3, t4, t5, t6, t7; \ |