Lines Matching refs:k2
275 int j, jc, jf, jj, k, k1, k2, k3, k4, kk, kt, nn, ns, nt; in FFTRADIX() local
439 k2 = kk + kspan; in FFTRADIX()
440 ak = Re [k2]; in FFTRADIX()
441 bk = Im [k2]; in FFTRADIX()
442 Re [k2] = Re [kk] - ak; in FFTRADIX()
443 Im [k2] = Im [kk] - bk; in FFTRADIX()
446 kk = k2 + kspan; in FFTRADIX()
458 k2 = kk + kspan; in FFTRADIX()
459 ak = Re [kk] - Re [k2]; in FFTRADIX()
460 bk = Im [kk] - Im [k2]; in FFTRADIX()
461 Re [kk] += Re [k2]; in FFTRADIX()
462 Im [kk] += Im [k2]; in FFTRADIX()
463 Re [k2] = c1 * ak - s1 * bk; in FFTRADIX()
464 Im [k2] = s1 * ak + c1 * bk; in FFTRADIX()
465 kk = k2 + kspan; in FFTRADIX()
467 k2 = kk - nt; in FFTRADIX()
469 kk = k1 - k2; in FFTRADIX()
470 } while (kk > k2); in FFTRADIX()
477 } while (kk < k2); in FFTRADIX()
493 k2 = k1 + kspan; in FFTRADIX()
494 k3 = k2 + kspan; in FFTRADIX()
495 akp = Re [kk] + Re [k2]; in FFTRADIX()
496 akm = Re [kk] - Re [k2]; in FFTRADIX()
499 bkp = Im [kk] + Im [k2]; in FFTRADIX()
500 bkm = Im [kk] - Im [k2]; in FFTRADIX()
521 Re [k2] = ajp; in FFTRADIX()
524 Im [k2] = bjp; in FFTRADIX()
528 Re [k2] = ajp * c2 - bjp * s2; in FFTRADIX()
531 Im [k2] = ajp * s2 + bjp * c2; in FFTRADIX()
570 k2 = k1 + kspan; in FFTRADIX()
573 aj = Re [k1] + Re [k2]; in FFTRADIX()
574 bj = Im [k1] + Im [k2]; in FFTRADIX()
579 aj = (Re [k1] - Re [k2]) * s60; in FFTRADIX()
580 bj = (Im [k1] - Im [k2]) * s60; in FFTRADIX()
582 Re [k2] = ak + bj; in FFTRADIX()
584 Im [k2] = bk - aj; in FFTRADIX()
585 kk = k2 + kspan; in FFTRADIX()
597 k2 = k1 + kspan; in FFTRADIX()
598 k3 = k2 + kspan; in FFTRADIX()
604 ajp = Re [k2] + Re [k3]; in FFTRADIX()
605 ajm = Re [k2] - Re [k3]; in FFTRADIX()
606 bjp = Im [k2] + Im [k3]; in FFTRADIX()
607 bjm = Im [k2] - Im [k3]; in FFTRADIX()
624 Re [k2] = ak - bj; in FFTRADIX()
626 Im [k2] = bk + aj; in FFTRADIX()
658 k2 = kk + ispan; in FFTRADIX()
664 k2 -= kspan; in FFTRADIX()
666 Rtmp [j - 1] = Re [k1] + Re [k2]; in FFTRADIX()
668 Itmp [j - 1] = Im [k1] + Im [k2]; in FFTRADIX()
671 Rtmp [j - 1] = Re [k1] - Re [k2]; in FFTRADIX()
672 Itmp [j - 1] = Im [k1] - Im [k2]; in FFTRADIX()
674 } while (k1 < k2); in FFTRADIX()
678 k2 = kk + ispan; in FFTRADIX()
682 k2 -= kspan; in FFTRADIX()
704 Re [k2] = ak + bj; in FFTRADIX()
705 Im [k2] = bk - aj; in FFTRADIX()
771 k2 = kspan; in FFTRADIX()
781 ak = Re [kk]; Re [kk] = Re [k2]; Re [k2] = ak; in FFTRADIX()
782 bk = Im [kk]; Im [kk] = Im [k2]; Im [k2] = bk; in FFTRADIX()
784 k2 += inc; in FFTRADIX()
787 k2 += ns - jc; in FFTRADIX()
789 k2 = k2 - nt + kspan; in FFTRADIX()
791 } while (k2 < (ns-1)); in FFTRADIX()
794 k2 -= fftstate->Perm [j - 1]; in FFTRADIX()
796 k2 = fftstate->Perm [j] + k2; in FFTRADIX()
797 } while (k2 > fftstate->Perm [j - 1]); in FFTRADIX()
800 if (kk < (k2-1)) in FFTRADIX()
803 k2 += kspan; in FFTRADIX()
804 } while (k2 < (ns-1)); in FFTRADIX()
811 ak = Re [kk]; Re [kk] = Re [k2]; Re [k2] = ak; in FFTRADIX()
812 bk = Im [kk]; Im [kk] = Im [k2]; Im [k2] = bk; in FFTRADIX()
814 k2 += kspan; in FFTRADIX()
815 } while (k2 < (ns-1)); in FFTRADIX()
818 k2 -= fftstate->Perm [j - 1]; in FFTRADIX()
820 k2 = fftstate->Perm [j] + k2; in FFTRADIX()
821 } while (k2 >= fftstate->Perm [j - 1]); in FFTRADIX()
824 if (kk < k2) in FFTRADIX()
827 k2 += kspan; in FFTRADIX()
828 } while (k2 < (ns-1)); in FFTRADIX()
852 k2 = fftstate->factor [kt - 1]; in FFTRADIX()
858 while (jj >= k2) { in FFTRADIX()
859 jj -= k2; in FFTRADIX()
860 k2 = kk; in FFTRADIX()
909 k2 = 0; in FFTRADIX()
911 k2++; in FFTRADIX()
912 Rtmp [k2 - 1] = Re [k1]; in FFTRADIX()
913 Itmp [k2 - 1] = Im [k1]; in FFTRADIX()
918 k2 = k1 - jc * (k + fftstate->Perm [k - 1]); in FFTRADIX()
921 Re [k1] = Re [k2]; in FFTRADIX()
922 Im [k1] = Im [k2]; in FFTRADIX()
924 k2 -= inc; in FFTRADIX()
926 kk = k2 + 1; in FFTRADIX()
929 k2 = 0; in FFTRADIX()
931 k2++; in FFTRADIX()
932 Re [k1] = Rtmp [k2 - 1]; in FFTRADIX()
933 Im [k1] = Itmp [k2 - 1]; in FFTRADIX()