Home | History | Annotate | Download | only in source

Lines Matching refs:kk

276   int j, jc, jf, jj, k, k1, k2, k3, k4, kk, kt, nn, ns, nt;
430 kk = 0;
440 k2 = kk + kspan;
443 Re [k2] = Re [kk] - ak;
444 Im [k2] = Im [kk] - bk;
445 Re [kk] += ak;
446 Im [kk] += bk;
447 kk = k2 + kspan;
448 } while (kk < nn);
449 kk -= nn;
450 } while (kk < jc);
451 if (kk >= kspan)
459 k2 = kk + kspan;
460 ak = Re [kk] - Re [k2];
461 bk = Im [kk] - Im [k2];
462 Re [kk] += Re [k2];
463 Im [kk] += Im [k2];
466 kk = k2 + kspan;
467 } while (kk < (nt-1));
468 k2 = kk - nt;
470 kk = k1 - k2;
471 } while (kk > k2);
477 kk += jc;
478 } while (kk < k2);
480 kk = (k1 - kspan + 1) / 2 + jc - 1;
481 } while (kk < (jc + jc));
493 k1 = kk + kspan;
496 akp = Re [kk] + Re [k2];
497 akm = Re [kk] - Re [k2];
500 bkp = Im [kk] + Im [k2];
501 bkm = Im [kk] - Im [k2];
504 Re [kk] = akp + ajp;
505 Im [kk] = bkp + bjp;
535 kk = k3 + kspan;
536 } while (kk < nt);
548 kk = kk - nt + jc;
549 } while (kk < kspan);
550 kk = kk - kspan + inc;
551 } while (kk < jc);
570 k1 = kk + kspan;
572 ak = Re [kk];
573 bk = Im [kk];
576 Re [kk] = ak + aj;
577 Im [kk] = bk + bj;
586 kk = k2 + kspan;
587 } while (kk < (nn - 1));
588 kk -= nn;
589 } while (kk < kspan);
597 k1 = kk + kspan;
609 aa = Re [kk];
610 bb = Im [kk];
611 Re [kk] = aa + akp + ajp;
612 Im [kk] = bb + bkp + bjp;
629 kk = k4 + kspan;
630 } while (kk < (nn-1));
631 kk -= nn;
632 } while (kk < kspan);
658 k1 = kk;
659 k2 = kk + ispan;
660 ak = aa = Re [kk];
661 bk = bb = Im [kk];
676 Re [kk] = ak;
677 Im [kk] = bk;
678 k1 = kk;
679 k2 = kk + ispan;
709 kk += ispan;
710 } while (kk < nn);
711 kk -= nn;
712 } while (kk < kspan);
719 kk = jc;
726 kk += kspan;
729 ak = Re [kk];
730 Re [kk] = c2 * ak - s2 * Im [kk];
731 Im [kk] = s2 * ak + c2 * Im [kk];
732 kk += ispan;
733 } while (kk < nt);
737 kk = kk - nt + kspan;
738 } while (kk < ispan);
744 kk = kk - ispan + jc;
745 } while (kk < kspan);
746 kk = kk - kspan + jc + inc;
747 } while (kk < (jc + jc));
771 kk = jc;
779 k = kk + jc;
781 /* swap Re [kk] <> Re [k2], Im [kk] <> Im [k2] */
782 ak = Re [kk]; Re [kk] = Re [k2]; Re [k2] = ak;
783 bk = Im [kk]; Im [kk] = Im [k2]; Im [k2] = bk;
784 kk += inc;
786 } while (kk < (k-1));
787 kk += ns - jc;
789 } while (kk < (nt-1));
791 kk = kk - nt + jc;
801 if (kk < (k2-1))
803 kk += jc;
806 } while (kk < (ns-1));
811 /* swap Re [kk] <> Re [k2], Im [kk] <> Im [k2] */
812 ak = Re [kk]; Re [kk] = Re [k2]; Re [k2] = ak;
813 bk = Im [kk]; Im [kk] = Im [k2]; Im [k2] = bk;
814 kk += inc;
825 if (kk < k2)
827 kk += inc;
830 } while (kk < (ns-1));
854 kk = fftstate->factor [k - 1];
858 jj += kk;
861 k2 = kk;
863 kk = fftstate->factor [k - 1];
864 jj += kk;
873 kk = fftstate->Perm [j - 1];
874 } while (kk < 0);
875 if (kk != j) {
877 k = kk;
878 kk = fftstate->Perm [k - 1];
879 fftstate->Perm [k - 1] = -kk;
880 } while (kk != j);
881 k3 = kk;
908 kk = jc * k + ii + jj;
909 k1 = kk + kspan - 1;
916 } while (k1 != (kk-1));
918 k1 = kk + kspan - 1;
926 } while (k1 != (kk-1));
927 kk = k2 + 1;
929 k1 = kk + kspan - 1;
936 } while (k1 != (kk-1));