Lines Matching refs:k2
276 int j, jc, jf, jj, k, k1, k2, k3, k4, kk, kt, nn, ns, nt;
444 k2 = kk + kspan;
445 ak = Re [k2];
446 bk = Im [k2];
447 Re [k2] = Re [kk] - ak;
448 Im [k2] = Im [kk] - bk;
451 kk = k2 + kspan;
463 k2 = kk + kspan;
464 ak = Re [kk] - Re [k2];
465 bk = Im [kk] - Im [k2];
466 Re [kk] += Re [k2];
467 Im [kk] += Im [k2];
468 Re [k2] = c1 * ak - s1 * bk;
469 Im [k2] = s1 * ak + c1 * bk;
470 kk = k2 + kspan;
472 k2 = kk - nt;
474 kk = k1 - k2;
475 } while (kk > k2);
482 } while (kk < k2);
498 k2 = k1 + kspan;
499 k3 = k2 + kspan;
500 akp = Re [kk] + Re [k2];
501 akm = Re [kk] - Re [k2];
504 bkp = Im [kk] + Im [k2];
505 bkm = Im [kk] - Im [k2];
526 Re [k2] = ajp;
529 Im [k2] = bjp;
533 Re [k2] = ajp * c2 - bjp * s2;
536 Im [k2] = ajp * s2 + bjp * c2;
575 k2 = k1 + kspan;
578 aj = Re [k1] + Re [k2];
579 bj = Im [k1] + Im [k2];
584 aj = (Re [k1] - Re [k2]) * s60;
585 bj = (Im [k1] - Im [k2]) * s60;
587 Re [k2] = ak + bj;
589 Im [k2] = bk - aj;
590 kk = k2 + kspan;
602 k2 = k1 + kspan;
603 k3 = k2 + kspan;
609 ajp = Re [k2] + Re [k3];
610 ajm = Re [k2] - Re [k3];
611 bjp = Im [k2] + Im [k3];
612 bjm = Im [k2] - Im [k3];
629 Re [k2] = ak - bj;
631 Im [k2] = bk + aj;
663 k2 = kk + ispan;
669 k2 -= kspan;
671 Rtmp [j - 1] = Re [k1] + Re [k2];
673 Itmp [j - 1] = Im [k1] + Im [k2];
676 Rtmp [j - 1] = Re [k1] - Re [k2];
677 Itmp [j - 1] = Im [k1] - Im [k2];
679 } while (k1 < k2);
683 k2 = kk + ispan;
687 k2 -= kspan;
709 Re [k2] = ak + bj;
710 Im [k2] = bk - aj;
776 k2 = kspan;
785 /* swap Re [kk] <> Re [k2], Im [kk] <> Im [k2] */
786 ak = Re [kk]; Re [kk] = Re [k2]; Re [k2] = ak;
787 bk = Im [kk]; Im [kk] = Im [k2]; Im [k2] = bk;
789 k2 += inc;
792 k2 += ns - jc;
794 k2 = k2 - nt + kspan;
796 } while (k2 < (ns-1));
799 k2 -= fftstate->Perm [j - 1];
801 k2 = fftstate->Perm [j] + k2;
802 } while (k2 > fftstate->Perm [j - 1]);
805 if (kk < (k2-1))
808 k2 += kspan;
809 } while (k2 < (ns-1));
815 /* swap Re [kk] <> Re [k2], Im [kk] <> Im [k2] */
816 ak = Re [kk]; Re [kk] = Re [k2]; Re [k2] = ak;
817 bk = Im [kk]; Im [kk] = Im [k2]; Im [k2] = bk;
819 k2 += kspan;
820 } while (k2 < (ns-1));
823 k2 -= fftstate->Perm [j - 1];
825 k2 = fftstate->Perm [j] + k2;
826 } while (k2 >= fftstate->Perm [j - 1]);
829 if (kk < k2)
832 k2 += kspan;
833 } while (k2 < (ns-1));
857 k2 = fftstate->factor [kt - 1];
863 while (jj >= k2) {
864 jj -= k2;
865 k2 = kk;
914 k2 = 0;
916 k2++;
917 Rtmp [k2 - 1] = Re [k1];
918 Itmp [k2 - 1] = Im [k1];
923 k2 = k1 - jc * (k + fftstate->Perm [k - 1]);
926 Re [k1] = Re [k2];
927 Im [k1] = Im [k2];
929 k2 -= inc;
931 kk = k2 + 1;
934 k2 = 0;
936 k2++;
937 Re [k1] = Rtmp [k2 - 1];
938 Im [k1] = Itmp [k2 - 1];