Searched refs:wk1i (Results 1 – 1 of 1) sorted by relevance
591 double ew, wn4r, wk1r, wk1i, wk2r, wk2i, wk3r, wk3i; in M_cft1st() local636 wk1i = sin(ew * kr); in M_cft1st()637 wk2r = 1 - 2 * wk1i * wk1i; in M_cft1st()638 wk2i = 2 * wk1i * wk1r; in M_cft1st()639 wk3r = wk1r - 2 * wk2i * wk1i; in M_cft1st()640 wk3i = 2 * wk2i * wk1r - wk1i; in M_cft1st()657 a[j + 2] = wk1r * x0r - wk1i * x0i; in M_cft1st()658 a[j + 3] = wk1r * x0i + wk1i * x0r; in M_cft1st()663 x0r = wn4r * (wk1r - wk1i); in M_cft1st()664 wk1i = wn4r * (wk1r + wk1i); in M_cft1st()[all …]