Mohsen Samadani
/
Algorithm-testing
Code for measuring the signal with a specified length and sampling rate, and saving it on a SD card.
Envelope/ifft.cpp@1:a514e4de034d, 2017-07-14 (annotated)
- Committer:
- msamadani
- Date:
- Fri Jul 14 14:41:10 2017 +0000
- Revision:
- 1:a514e4de034d
envSignal function works.
Who changed what in which revision?
User | Revision | Line number | New contents of line |
---|---|---|---|
msamadani | 1:a514e4de034d | 1 | // |
msamadani | 1:a514e4de034d | 2 | // Academic License - for use in teaching, academic research, and meeting |
msamadani | 1:a514e4de034d | 3 | // course requirements at degree granting institutions only. Not for |
msamadani | 1:a514e4de034d | 4 | // government, commercial, or other organizational use. |
msamadani | 1:a514e4de034d | 5 | // File: ifft.cpp |
msamadani | 1:a514e4de034d | 6 | // |
msamadani | 1:a514e4de034d | 7 | // MATLAB Coder version : 3.3 |
msamadani | 1:a514e4de034d | 8 | // C/C++ source code generated on : 13-Jul-2017 15:47:47 |
msamadani | 1:a514e4de034d | 9 | // |
msamadani | 1:a514e4de034d | 10 | |
msamadani | 1:a514e4de034d | 11 | // Include Files |
msamadani | 1:a514e4de034d | 12 | #include "envSignal.h" |
msamadani | 1:a514e4de034d | 13 | #include "ifft.h" |
msamadani | 1:a514e4de034d | 14 | #include "envSignal_emxutil.h" |
msamadani | 1:a514e4de034d | 15 | |
msamadani | 1:a514e4de034d | 16 | // Function Definitions |
msamadani | 1:a514e4de034d | 17 | |
msamadani | 1:a514e4de034d | 18 | // |
msamadani | 1:a514e4de034d | 19 | // Arguments : const emxArray_creal32_T *x |
msamadani | 1:a514e4de034d | 20 | // int n1_unsigned |
msamadani | 1:a514e4de034d | 21 | // const emxArray_real32_T *costab |
msamadani | 1:a514e4de034d | 22 | // const emxArray_real32_T *sintab |
msamadani | 1:a514e4de034d | 23 | // emxArray_creal32_T *y |
msamadani | 1:a514e4de034d | 24 | // Return Type : void |
msamadani | 1:a514e4de034d | 25 | // |
msamadani | 1:a514e4de034d | 26 | void b_r2br_r2dit_trig(const emxArray_creal32_T *x, int n1_unsigned, const |
msamadani | 1:a514e4de034d | 27 | emxArray_real32_T *costab, const emxArray_real32_T *sintab, emxArray_creal32_T |
msamadani | 1:a514e4de034d | 28 | *y) |
msamadani | 1:a514e4de034d | 29 | { |
msamadani | 1:a514e4de034d | 30 | int j; |
msamadani | 1:a514e4de034d | 31 | int nRowsD2; |
msamadani | 1:a514e4de034d | 32 | int nRowsD4; |
msamadani | 1:a514e4de034d | 33 | int iDelta2; |
msamadani | 1:a514e4de034d | 34 | int iy; |
msamadani | 1:a514e4de034d | 35 | int ix; |
msamadani | 1:a514e4de034d | 36 | int ju; |
msamadani | 1:a514e4de034d | 37 | int i; |
msamadani | 1:a514e4de034d | 38 | boolean_T tst; |
msamadani | 1:a514e4de034d | 39 | float temp_re; |
msamadani | 1:a514e4de034d | 40 | float temp_im; |
msamadani | 1:a514e4de034d | 41 | float r; |
msamadani | 1:a514e4de034d | 42 | float twid_im; |
msamadani | 1:a514e4de034d | 43 | int ihi; |
msamadani | 1:a514e4de034d | 44 | j = x->size[0]; |
msamadani | 1:a514e4de034d | 45 | if (!(j < n1_unsigned)) { |
msamadani | 1:a514e4de034d | 46 | j = n1_unsigned; |
msamadani | 1:a514e4de034d | 47 | } |
msamadani | 1:a514e4de034d | 48 | |
msamadani | 1:a514e4de034d | 49 | nRowsD2 = n1_unsigned / 2; |
msamadani | 1:a514e4de034d | 50 | nRowsD4 = nRowsD2 / 2; |
msamadani | 1:a514e4de034d | 51 | iDelta2 = y->size[0]; |
msamadani | 1:a514e4de034d | 52 | y->size[0] = n1_unsigned; |
msamadani | 1:a514e4de034d | 53 | emxEnsureCapacity((emxArray__common *)y, iDelta2, sizeof(creal32_T)); |
msamadani | 1:a514e4de034d | 54 | if (n1_unsigned > x->size[0]) { |
msamadani | 1:a514e4de034d | 55 | iy = y->size[0]; |
msamadani | 1:a514e4de034d | 56 | iDelta2 = y->size[0]; |
msamadani | 1:a514e4de034d | 57 | y->size[0] = iy; |
msamadani | 1:a514e4de034d | 58 | emxEnsureCapacity((emxArray__common *)y, iDelta2, sizeof(creal32_T)); |
msamadani | 1:a514e4de034d | 59 | for (iDelta2 = 0; iDelta2 < iy; iDelta2++) { |
msamadani | 1:a514e4de034d | 60 | y->data[iDelta2].re = 0.0F; |
msamadani | 1:a514e4de034d | 61 | y->data[iDelta2].im = 0.0F; |
msamadani | 1:a514e4de034d | 62 | } |
msamadani | 1:a514e4de034d | 63 | } |
msamadani | 1:a514e4de034d | 64 | |
msamadani | 1:a514e4de034d | 65 | ix = 0; |
msamadani | 1:a514e4de034d | 66 | ju = 0; |
msamadani | 1:a514e4de034d | 67 | iy = 0; |
msamadani | 1:a514e4de034d | 68 | for (i = 1; i < j; i++) { |
msamadani | 1:a514e4de034d | 69 | y->data[iy] = x->data[ix]; |
msamadani | 1:a514e4de034d | 70 | iDelta2 = n1_unsigned; |
msamadani | 1:a514e4de034d | 71 | tst = true; |
msamadani | 1:a514e4de034d | 72 | while (tst) { |
msamadani | 1:a514e4de034d | 73 | iDelta2 >>= 1; |
msamadani | 1:a514e4de034d | 74 | ju ^= iDelta2; |
msamadani | 1:a514e4de034d | 75 | tst = ((ju & iDelta2) == 0); |
msamadani | 1:a514e4de034d | 76 | } |
msamadani | 1:a514e4de034d | 77 | |
msamadani | 1:a514e4de034d | 78 | iy = ju; |
msamadani | 1:a514e4de034d | 79 | ix++; |
msamadani | 1:a514e4de034d | 80 | } |
msamadani | 1:a514e4de034d | 81 | |
msamadani | 1:a514e4de034d | 82 | y->data[iy] = x->data[ix]; |
msamadani | 1:a514e4de034d | 83 | if (n1_unsigned > 1) { |
msamadani | 1:a514e4de034d | 84 | for (i = 0; i <= n1_unsigned - 2; i += 2) { |
msamadani | 1:a514e4de034d | 85 | temp_re = y->data[i + 1].re; |
msamadani | 1:a514e4de034d | 86 | temp_im = y->data[i + 1].im; |
msamadani | 1:a514e4de034d | 87 | y->data[i + 1].re = y->data[i].re - y->data[i + 1].re; |
msamadani | 1:a514e4de034d | 88 | y->data[i + 1].im = y->data[i].im - y->data[i + 1].im; |
msamadani | 1:a514e4de034d | 89 | y->data[i].re += temp_re; |
msamadani | 1:a514e4de034d | 90 | y->data[i].im += temp_im; |
msamadani | 1:a514e4de034d | 91 | } |
msamadani | 1:a514e4de034d | 92 | } |
msamadani | 1:a514e4de034d | 93 | |
msamadani | 1:a514e4de034d | 94 | iy = 2; |
msamadani | 1:a514e4de034d | 95 | iDelta2 = 4; |
msamadani | 1:a514e4de034d | 96 | ix = 1 + ((nRowsD4 - 1) << 2); |
msamadani | 1:a514e4de034d | 97 | while (nRowsD4 > 0) { |
msamadani | 1:a514e4de034d | 98 | for (i = 0; i < ix; i += iDelta2) { |
msamadani | 1:a514e4de034d | 99 | temp_re = y->data[i + iy].re; |
msamadani | 1:a514e4de034d | 100 | temp_im = y->data[i + iy].im; |
msamadani | 1:a514e4de034d | 101 | y->data[i + iy].re = y->data[i].re - temp_re; |
msamadani | 1:a514e4de034d | 102 | y->data[i + iy].im = y->data[i].im - temp_im; |
msamadani | 1:a514e4de034d | 103 | y->data[i].re += temp_re; |
msamadani | 1:a514e4de034d | 104 | y->data[i].im += temp_im; |
msamadani | 1:a514e4de034d | 105 | } |
msamadani | 1:a514e4de034d | 106 | |
msamadani | 1:a514e4de034d | 107 | ju = 1; |
msamadani | 1:a514e4de034d | 108 | for (j = nRowsD4; j < nRowsD2; j += nRowsD4) { |
msamadani | 1:a514e4de034d | 109 | r = costab->data[j]; |
msamadani | 1:a514e4de034d | 110 | twid_im = sintab->data[j]; |
msamadani | 1:a514e4de034d | 111 | i = ju; |
msamadani | 1:a514e4de034d | 112 | ihi = ju + ix; |
msamadani | 1:a514e4de034d | 113 | while (i < ihi) { |
msamadani | 1:a514e4de034d | 114 | temp_re = r * y->data[i + iy].re - twid_im * y->data[i + iy].im; |
msamadani | 1:a514e4de034d | 115 | temp_im = r * y->data[i + iy].im + twid_im * y->data[i + iy].re; |
msamadani | 1:a514e4de034d | 116 | y->data[i + iy].re = y->data[i].re - temp_re; |
msamadani | 1:a514e4de034d | 117 | y->data[i + iy].im = y->data[i].im - temp_im; |
msamadani | 1:a514e4de034d | 118 | y->data[i].re += temp_re; |
msamadani | 1:a514e4de034d | 119 | y->data[i].im += temp_im; |
msamadani | 1:a514e4de034d | 120 | i += iDelta2; |
msamadani | 1:a514e4de034d | 121 | } |
msamadani | 1:a514e4de034d | 122 | |
msamadani | 1:a514e4de034d | 123 | ju++; |
msamadani | 1:a514e4de034d | 124 | } |
msamadani | 1:a514e4de034d | 125 | |
msamadani | 1:a514e4de034d | 126 | nRowsD4 /= 2; |
msamadani | 1:a514e4de034d | 127 | iy = iDelta2; |
msamadani | 1:a514e4de034d | 128 | iDelta2 += iDelta2; |
msamadani | 1:a514e4de034d | 129 | ix -= iy; |
msamadani | 1:a514e4de034d | 130 | } |
msamadani | 1:a514e4de034d | 131 | |
msamadani | 1:a514e4de034d | 132 | if (y->size[0] > 1) { |
msamadani | 1:a514e4de034d | 133 | r = 1.0F / (float)y->size[0]; |
msamadani | 1:a514e4de034d | 134 | iDelta2 = y->size[0]; |
msamadani | 1:a514e4de034d | 135 | emxEnsureCapacity((emxArray__common *)y, iDelta2, sizeof(creal32_T)); |
msamadani | 1:a514e4de034d | 136 | iy = y->size[0]; |
msamadani | 1:a514e4de034d | 137 | for (iDelta2 = 0; iDelta2 < iy; iDelta2++) { |
msamadani | 1:a514e4de034d | 138 | y->data[iDelta2].re *= r; |
msamadani | 1:a514e4de034d | 139 | y->data[iDelta2].im *= r; |
msamadani | 1:a514e4de034d | 140 | } |
msamadani | 1:a514e4de034d | 141 | } |
msamadani | 1:a514e4de034d | 142 | } |
msamadani | 1:a514e4de034d | 143 | |
msamadani | 1:a514e4de034d | 144 | // |
msamadani | 1:a514e4de034d | 145 | // Arguments : int n1 |
msamadani | 1:a514e4de034d | 146 | // boolean_T useRadix2 |
msamadani | 1:a514e4de034d | 147 | // int *N2blue |
msamadani | 1:a514e4de034d | 148 | // int *nRows |
msamadani | 1:a514e4de034d | 149 | // Return Type : void |
msamadani | 1:a514e4de034d | 150 | // |
msamadani | 1:a514e4de034d | 151 | void get_algo_sizes(int n1, boolean_T useRadix2, int *N2blue, int *nRows) |
msamadani | 1:a514e4de034d | 152 | { |
msamadani | 1:a514e4de034d | 153 | int nn1m1; |
msamadani | 1:a514e4de034d | 154 | int pmax; |
msamadani | 1:a514e4de034d | 155 | int pmin; |
msamadani | 1:a514e4de034d | 156 | boolean_T exitg1; |
msamadani | 1:a514e4de034d | 157 | int p; |
msamadani | 1:a514e4de034d | 158 | int pow2p; |
msamadani | 1:a514e4de034d | 159 | *N2blue = 1; |
msamadani | 1:a514e4de034d | 160 | if (useRadix2) { |
msamadani | 1:a514e4de034d | 161 | *nRows = n1; |
msamadani | 1:a514e4de034d | 162 | } else { |
msamadani | 1:a514e4de034d | 163 | if (n1 > 0) { |
msamadani | 1:a514e4de034d | 164 | nn1m1 = (n1 + n1) - 1; |
msamadani | 1:a514e4de034d | 165 | pmax = 31; |
msamadani | 1:a514e4de034d | 166 | if (nn1m1 > MIN_int32_T) { |
msamadani | 1:a514e4de034d | 167 | if (nn1m1 < 0) { |
msamadani | 1:a514e4de034d | 168 | nn1m1 = -nn1m1; |
msamadani | 1:a514e4de034d | 169 | } |
msamadani | 1:a514e4de034d | 170 | |
msamadani | 1:a514e4de034d | 171 | if (nn1m1 <= 1) { |
msamadani | 1:a514e4de034d | 172 | pmax = 0; |
msamadani | 1:a514e4de034d | 173 | } else { |
msamadani | 1:a514e4de034d | 174 | pmin = 0; |
msamadani | 1:a514e4de034d | 175 | exitg1 = false; |
msamadani | 1:a514e4de034d | 176 | while ((!exitg1) && (pmax - pmin > 1)) { |
msamadani | 1:a514e4de034d | 177 | p = (pmin + pmax) >> 1; |
msamadani | 1:a514e4de034d | 178 | pow2p = 1 << p; |
msamadani | 1:a514e4de034d | 179 | if (pow2p == nn1m1) { |
msamadani | 1:a514e4de034d | 180 | pmax = p; |
msamadani | 1:a514e4de034d | 181 | exitg1 = true; |
msamadani | 1:a514e4de034d | 182 | } else if (pow2p > nn1m1) { |
msamadani | 1:a514e4de034d | 183 | pmax = p; |
msamadani | 1:a514e4de034d | 184 | } else { |
msamadani | 1:a514e4de034d | 185 | pmin = p; |
msamadani | 1:a514e4de034d | 186 | } |
msamadani | 1:a514e4de034d | 187 | } |
msamadani | 1:a514e4de034d | 188 | } |
msamadani | 1:a514e4de034d | 189 | } |
msamadani | 1:a514e4de034d | 190 | |
msamadani | 1:a514e4de034d | 191 | *N2blue = 1 << pmax; |
msamadani | 1:a514e4de034d | 192 | } |
msamadani | 1:a514e4de034d | 193 | |
msamadani | 1:a514e4de034d | 194 | *nRows = *N2blue; |
msamadani | 1:a514e4de034d | 195 | } |
msamadani | 1:a514e4de034d | 196 | } |
msamadani | 1:a514e4de034d | 197 | |
msamadani | 1:a514e4de034d | 198 | // |
msamadani | 1:a514e4de034d | 199 | // Arguments : const emxArray_creal32_T *x |
msamadani | 1:a514e4de034d | 200 | // emxArray_creal32_T *y |
msamadani | 1:a514e4de034d | 201 | // Return Type : void |
msamadani | 1:a514e4de034d | 202 | // |
msamadani | 1:a514e4de034d | 203 | void ifft(const emxArray_creal32_T *x, emxArray_creal32_T *y) |
msamadani | 1:a514e4de034d | 204 | { |
msamadani | 1:a514e4de034d | 205 | int n1; |
msamadani | 1:a514e4de034d | 206 | emxArray_real32_T *costab1q; |
msamadani | 1:a514e4de034d | 207 | int nInt2; |
msamadani | 1:a514e4de034d | 208 | boolean_T useRadix2; |
msamadani | 1:a514e4de034d | 209 | int N2blue; |
msamadani | 1:a514e4de034d | 210 | int nd2; |
msamadani | 1:a514e4de034d | 211 | float e; |
msamadani | 1:a514e4de034d | 212 | int nRowsD4; |
msamadani | 1:a514e4de034d | 213 | int k; |
msamadani | 1:a514e4de034d | 214 | emxArray_real32_T *costab; |
msamadani | 1:a514e4de034d | 215 | emxArray_real32_T *sintab; |
msamadani | 1:a514e4de034d | 216 | emxArray_real32_T *sintabinv; |
msamadani | 1:a514e4de034d | 217 | emxArray_creal32_T *wwc; |
msamadani | 1:a514e4de034d | 218 | int nInt2m1; |
msamadani | 1:a514e4de034d | 219 | int b_y; |
msamadani | 1:a514e4de034d | 220 | float denom_re; |
msamadani | 1:a514e4de034d | 221 | emxArray_creal32_T *fy; |
msamadani | 1:a514e4de034d | 222 | float wwc_im; |
msamadani | 1:a514e4de034d | 223 | emxArray_creal32_T *fv; |
msamadani | 1:a514e4de034d | 224 | float wwc_re; |
msamadani | 1:a514e4de034d | 225 | float fv_im; |
msamadani | 1:a514e4de034d | 226 | float b_wwc_im; |
msamadani | 1:a514e4de034d | 227 | float fv_re; |
msamadani | 1:a514e4de034d | 228 | float b_fv_im; |
msamadani | 1:a514e4de034d | 229 | float b_fv_re; |
msamadani | 1:a514e4de034d | 230 | n1 = x->size[0]; |
msamadani | 1:a514e4de034d | 231 | if (x->size[0] == 0) { |
msamadani | 1:a514e4de034d | 232 | nInt2 = y->size[0]; |
msamadani | 1:a514e4de034d | 233 | y->size[0] = 0; |
msamadani | 1:a514e4de034d | 234 | emxEnsureCapacity((emxArray__common *)y, nInt2, sizeof(creal32_T)); |
msamadani | 1:a514e4de034d | 235 | } else { |
msamadani | 1:a514e4de034d | 236 | emxInit_real32_T(&costab1q, 2); |
msamadani | 1:a514e4de034d | 237 | useRadix2 = ((x->size[0] & (x->size[0] - 1)) == 0); |
msamadani | 1:a514e4de034d | 238 | get_algo_sizes(x->size[0], useRadix2, &N2blue, &nd2); |
msamadani | 1:a514e4de034d | 239 | e = 6.28318548F / (float)nd2; |
msamadani | 1:a514e4de034d | 240 | nRowsD4 = nd2 / 2 / 2; |
msamadani | 1:a514e4de034d | 241 | nInt2 = costab1q->size[0] * costab1q->size[1]; |
msamadani | 1:a514e4de034d | 242 | costab1q->size[0] = 1; |
msamadani | 1:a514e4de034d | 243 | costab1q->size[1] = nRowsD4 + 1; |
msamadani | 1:a514e4de034d | 244 | emxEnsureCapacity((emxArray__common *)costab1q, nInt2, sizeof(float)); |
msamadani | 1:a514e4de034d | 245 | costab1q->data[0] = 1.0F; |
msamadani | 1:a514e4de034d | 246 | nd2 = nRowsD4 / 2; |
msamadani | 1:a514e4de034d | 247 | for (k = 1; k <= nd2; k++) { |
msamadani | 1:a514e4de034d | 248 | costab1q->data[k] = std::cos(e * (float)k); |
msamadani | 1:a514e4de034d | 249 | } |
msamadani | 1:a514e4de034d | 250 | |
msamadani | 1:a514e4de034d | 251 | for (k = nd2 + 1; k < nRowsD4; k++) { |
msamadani | 1:a514e4de034d | 252 | costab1q->data[k] = std::sin(e * (float)(nRowsD4 - k)); |
msamadani | 1:a514e4de034d | 253 | } |
msamadani | 1:a514e4de034d | 254 | |
msamadani | 1:a514e4de034d | 255 | costab1q->data[nRowsD4] = 0.0F; |
msamadani | 1:a514e4de034d | 256 | emxInit_real32_T(&costab, 2); |
msamadani | 1:a514e4de034d | 257 | emxInit_real32_T(&sintab, 2); |
msamadani | 1:a514e4de034d | 258 | emxInit_real32_T(&sintabinv, 2); |
msamadani | 1:a514e4de034d | 259 | if (!useRadix2) { |
msamadani | 1:a514e4de034d | 260 | nd2 = costab1q->size[1] - 1; |
msamadani | 1:a514e4de034d | 261 | nRowsD4 = (costab1q->size[1] - 1) << 1; |
msamadani | 1:a514e4de034d | 262 | nInt2 = costab->size[0] * costab->size[1]; |
msamadani | 1:a514e4de034d | 263 | costab->size[0] = 1; |
msamadani | 1:a514e4de034d | 264 | costab->size[1] = nRowsD4 + 1; |
msamadani | 1:a514e4de034d | 265 | emxEnsureCapacity((emxArray__common *)costab, nInt2, sizeof(float)); |
msamadani | 1:a514e4de034d | 266 | nInt2 = sintab->size[0] * sintab->size[1]; |
msamadani | 1:a514e4de034d | 267 | sintab->size[0] = 1; |
msamadani | 1:a514e4de034d | 268 | sintab->size[1] = nRowsD4 + 1; |
msamadani | 1:a514e4de034d | 269 | emxEnsureCapacity((emxArray__common *)sintab, nInt2, sizeof(float)); |
msamadani | 1:a514e4de034d | 270 | costab->data[0] = 1.0F; |
msamadani | 1:a514e4de034d | 271 | sintab->data[0] = 0.0F; |
msamadani | 1:a514e4de034d | 272 | nInt2 = sintabinv->size[0] * sintabinv->size[1]; |
msamadani | 1:a514e4de034d | 273 | sintabinv->size[0] = 1; |
msamadani | 1:a514e4de034d | 274 | sintabinv->size[1] = nRowsD4 + 1; |
msamadani | 1:a514e4de034d | 275 | emxEnsureCapacity((emxArray__common *)sintabinv, nInt2, sizeof(float)); |
msamadani | 1:a514e4de034d | 276 | for (k = 1; k <= nd2; k++) { |
msamadani | 1:a514e4de034d | 277 | sintabinv->data[k] = costab1q->data[nd2 - k]; |
msamadani | 1:a514e4de034d | 278 | } |
msamadani | 1:a514e4de034d | 279 | |
msamadani | 1:a514e4de034d | 280 | for (k = costab1q->size[1]; k <= nRowsD4; k++) { |
msamadani | 1:a514e4de034d | 281 | sintabinv->data[k] = costab1q->data[k - nd2]; |
msamadani | 1:a514e4de034d | 282 | } |
msamadani | 1:a514e4de034d | 283 | |
msamadani | 1:a514e4de034d | 284 | for (k = 1; k <= nd2; k++) { |
msamadani | 1:a514e4de034d | 285 | costab->data[k] = costab1q->data[k]; |
msamadani | 1:a514e4de034d | 286 | sintab->data[k] = -costab1q->data[nd2 - k]; |
msamadani | 1:a514e4de034d | 287 | } |
msamadani | 1:a514e4de034d | 288 | |
msamadani | 1:a514e4de034d | 289 | for (k = costab1q->size[1]; k <= nRowsD4; k++) { |
msamadani | 1:a514e4de034d | 290 | costab->data[k] = -costab1q->data[nRowsD4 - k]; |
msamadani | 1:a514e4de034d | 291 | sintab->data[k] = -costab1q->data[k - nd2]; |
msamadani | 1:a514e4de034d | 292 | } |
msamadani | 1:a514e4de034d | 293 | } else { |
msamadani | 1:a514e4de034d | 294 | nd2 = costab1q->size[1] - 1; |
msamadani | 1:a514e4de034d | 295 | nRowsD4 = (costab1q->size[1] - 1) << 1; |
msamadani | 1:a514e4de034d | 296 | nInt2 = costab->size[0] * costab->size[1]; |
msamadani | 1:a514e4de034d | 297 | costab->size[0] = 1; |
msamadani | 1:a514e4de034d | 298 | costab->size[1] = nRowsD4 + 1; |
msamadani | 1:a514e4de034d | 299 | emxEnsureCapacity((emxArray__common *)costab, nInt2, sizeof(float)); |
msamadani | 1:a514e4de034d | 300 | nInt2 = sintab->size[0] * sintab->size[1]; |
msamadani | 1:a514e4de034d | 301 | sintab->size[0] = 1; |
msamadani | 1:a514e4de034d | 302 | sintab->size[1] = nRowsD4 + 1; |
msamadani | 1:a514e4de034d | 303 | emxEnsureCapacity((emxArray__common *)sintab, nInt2, sizeof(float)); |
msamadani | 1:a514e4de034d | 304 | costab->data[0] = 1.0F; |
msamadani | 1:a514e4de034d | 305 | sintab->data[0] = 0.0F; |
msamadani | 1:a514e4de034d | 306 | for (k = 1; k <= nd2; k++) { |
msamadani | 1:a514e4de034d | 307 | costab->data[k] = costab1q->data[k]; |
msamadani | 1:a514e4de034d | 308 | sintab->data[k] = costab1q->data[nd2 - k]; |
msamadani | 1:a514e4de034d | 309 | } |
msamadani | 1:a514e4de034d | 310 | |
msamadani | 1:a514e4de034d | 311 | for (k = costab1q->size[1]; k <= nRowsD4; k++) { |
msamadani | 1:a514e4de034d | 312 | costab->data[k] = -costab1q->data[nRowsD4 - k]; |
msamadani | 1:a514e4de034d | 313 | sintab->data[k] = costab1q->data[k - nd2]; |
msamadani | 1:a514e4de034d | 314 | } |
msamadani | 1:a514e4de034d | 315 | |
msamadani | 1:a514e4de034d | 316 | nInt2 = sintabinv->size[0] * sintabinv->size[1]; |
msamadani | 1:a514e4de034d | 317 | sintabinv->size[0] = 1; |
msamadani | 1:a514e4de034d | 318 | sintabinv->size[1] = 0; |
msamadani | 1:a514e4de034d | 319 | emxEnsureCapacity((emxArray__common *)sintabinv, nInt2, sizeof(float)); |
msamadani | 1:a514e4de034d | 320 | } |
msamadani | 1:a514e4de034d | 321 | |
msamadani | 1:a514e4de034d | 322 | emxFree_real32_T(&costab1q); |
msamadani | 1:a514e4de034d | 323 | if (useRadix2) { |
msamadani | 1:a514e4de034d | 324 | b_r2br_r2dit_trig(x, x->size[0], costab, sintab, y); |
msamadani | 1:a514e4de034d | 325 | } else { |
msamadani | 1:a514e4de034d | 326 | emxInit_creal32_T(&wwc, 1); |
msamadani | 1:a514e4de034d | 327 | nInt2m1 = (x->size[0] + x->size[0]) - 1; |
msamadani | 1:a514e4de034d | 328 | nInt2 = wwc->size[0]; |
msamadani | 1:a514e4de034d | 329 | wwc->size[0] = nInt2m1; |
msamadani | 1:a514e4de034d | 330 | emxEnsureCapacity((emxArray__common *)wwc, nInt2, sizeof(creal32_T)); |
msamadani | 1:a514e4de034d | 331 | nd2 = x->size[0]; |
msamadani | 1:a514e4de034d | 332 | nRowsD4 = 0; |
msamadani | 1:a514e4de034d | 333 | wwc->data[x->size[0] - 1].re = 1.0F; |
msamadani | 1:a514e4de034d | 334 | wwc->data[x->size[0] - 1].im = 0.0F; |
msamadani | 1:a514e4de034d | 335 | nInt2 = x->size[0] << 1; |
msamadani | 1:a514e4de034d | 336 | for (k = 1; k < n1; k++) { |
msamadani | 1:a514e4de034d | 337 | b_y = (k << 1) - 1; |
msamadani | 1:a514e4de034d | 338 | if (nInt2 - nRowsD4 <= b_y) { |
msamadani | 1:a514e4de034d | 339 | nRowsD4 += b_y - nInt2; |
msamadani | 1:a514e4de034d | 340 | } else { |
msamadani | 1:a514e4de034d | 341 | nRowsD4 += b_y; |
msamadani | 1:a514e4de034d | 342 | } |
msamadani | 1:a514e4de034d | 343 | |
msamadani | 1:a514e4de034d | 344 | e = 3.14159274F * (float)nRowsD4 / (float)x->size[0]; |
msamadani | 1:a514e4de034d | 345 | if (e == 0.0F) { |
msamadani | 1:a514e4de034d | 346 | denom_re = 1.0F; |
msamadani | 1:a514e4de034d | 347 | e = 0.0F; |
msamadani | 1:a514e4de034d | 348 | } else { |
msamadani | 1:a514e4de034d | 349 | denom_re = std::cos(e); |
msamadani | 1:a514e4de034d | 350 | e = std::sin(e); |
msamadani | 1:a514e4de034d | 351 | } |
msamadani | 1:a514e4de034d | 352 | |
msamadani | 1:a514e4de034d | 353 | wwc->data[nd2 - 2].re = denom_re; |
msamadani | 1:a514e4de034d | 354 | wwc->data[nd2 - 2].im = -e; |
msamadani | 1:a514e4de034d | 355 | nd2--; |
msamadani | 1:a514e4de034d | 356 | } |
msamadani | 1:a514e4de034d | 357 | |
msamadani | 1:a514e4de034d | 358 | nd2 = 0; |
msamadani | 1:a514e4de034d | 359 | for (k = nInt2m1 - 1; k >= n1; k--) { |
msamadani | 1:a514e4de034d | 360 | wwc->data[k] = wwc->data[nd2]; |
msamadani | 1:a514e4de034d | 361 | nd2++; |
msamadani | 1:a514e4de034d | 362 | } |
msamadani | 1:a514e4de034d | 363 | |
msamadani | 1:a514e4de034d | 364 | nRowsD4 = x->size[0]; |
msamadani | 1:a514e4de034d | 365 | nd2 = x->size[0]; |
msamadani | 1:a514e4de034d | 366 | nInt2 = y->size[0]; |
msamadani | 1:a514e4de034d | 367 | y->size[0] = nd2; |
msamadani | 1:a514e4de034d | 368 | emxEnsureCapacity((emxArray__common *)y, nInt2, sizeof(creal32_T)); |
msamadani | 1:a514e4de034d | 369 | nd2 = 0; |
msamadani | 1:a514e4de034d | 370 | for (k = 0; k + 1 <= nRowsD4; k++) { |
msamadani | 1:a514e4de034d | 371 | denom_re = wwc->data[(n1 + k) - 1].re; |
msamadani | 1:a514e4de034d | 372 | e = wwc->data[(n1 + k) - 1].im; |
msamadani | 1:a514e4de034d | 373 | wwc_im = x->data[nd2].re; |
msamadani | 1:a514e4de034d | 374 | wwc_re = x->data[nd2].im; |
msamadani | 1:a514e4de034d | 375 | fv_im = x->data[nd2].im; |
msamadani | 1:a514e4de034d | 376 | b_wwc_im = x->data[nd2].re; |
msamadani | 1:a514e4de034d | 377 | y->data[k].re = denom_re * wwc_im + e * wwc_re; |
msamadani | 1:a514e4de034d | 378 | y->data[k].im = denom_re * fv_im - e * b_wwc_im; |
msamadani | 1:a514e4de034d | 379 | nd2++; |
msamadani | 1:a514e4de034d | 380 | } |
msamadani | 1:a514e4de034d | 381 | |
msamadani | 1:a514e4de034d | 382 | while (nRowsD4 + 1 <= n1) { |
msamadani | 1:a514e4de034d | 383 | y->data[nRowsD4].re = 0.0F; |
msamadani | 1:a514e4de034d | 384 | y->data[nRowsD4].im = 0.0F; |
msamadani | 1:a514e4de034d | 385 | nRowsD4++; |
msamadani | 1:a514e4de034d | 386 | } |
msamadani | 1:a514e4de034d | 387 | |
msamadani | 1:a514e4de034d | 388 | emxInit_creal32_T(&fy, 1); |
msamadani | 1:a514e4de034d | 389 | emxInit_creal32_T(&fv, 1); |
msamadani | 1:a514e4de034d | 390 | r2br_r2dit_trig_impl(y, 0, N2blue, costab, sintab, fy); |
msamadani | 1:a514e4de034d | 391 | r2br_r2dit_trig(wwc, N2blue, costab, sintab, fv); |
msamadani | 1:a514e4de034d | 392 | nInt2 = fy->size[0]; |
msamadani | 1:a514e4de034d | 393 | emxEnsureCapacity((emxArray__common *)fy, nInt2, sizeof(creal32_T)); |
msamadani | 1:a514e4de034d | 394 | nd2 = fy->size[0]; |
msamadani | 1:a514e4de034d | 395 | for (nInt2 = 0; nInt2 < nd2; nInt2++) { |
msamadani | 1:a514e4de034d | 396 | e = fy->data[nInt2].re; |
msamadani | 1:a514e4de034d | 397 | wwc_im = fy->data[nInt2].im; |
msamadani | 1:a514e4de034d | 398 | fv_re = fv->data[nInt2].re; |
msamadani | 1:a514e4de034d | 399 | b_fv_im = fv->data[nInt2].im; |
msamadani | 1:a514e4de034d | 400 | fy->data[nInt2].re = e * fv_re - wwc_im * b_fv_im; |
msamadani | 1:a514e4de034d | 401 | fy->data[nInt2].im = e * b_fv_im + wwc_im * fv_re; |
msamadani | 1:a514e4de034d | 402 | } |
msamadani | 1:a514e4de034d | 403 | |
msamadani | 1:a514e4de034d | 404 | b_r2br_r2dit_trig(fy, N2blue, costab, sintabinv, fv); |
msamadani | 1:a514e4de034d | 405 | nd2 = 0; |
msamadani | 1:a514e4de034d | 406 | denom_re = (float)x->size[0]; |
msamadani | 1:a514e4de034d | 407 | k = (int)(float)x->size[0] - 1; |
msamadani | 1:a514e4de034d | 408 | emxFree_creal32_T(&fy); |
msamadani | 1:a514e4de034d | 409 | while (k + 1 <= wwc->size[0]) { |
msamadani | 1:a514e4de034d | 410 | e = wwc->data[k].re; |
msamadani | 1:a514e4de034d | 411 | fv_re = fv->data[k].re; |
msamadani | 1:a514e4de034d | 412 | wwc_im = wwc->data[k].im; |
msamadani | 1:a514e4de034d | 413 | b_fv_im = fv->data[k].im; |
msamadani | 1:a514e4de034d | 414 | wwc_re = wwc->data[k].re; |
msamadani | 1:a514e4de034d | 415 | fv_im = fv->data[k].im; |
msamadani | 1:a514e4de034d | 416 | b_wwc_im = wwc->data[k].im; |
msamadani | 1:a514e4de034d | 417 | b_fv_re = fv->data[k].re; |
msamadani | 1:a514e4de034d | 418 | y->data[nd2].re = e * fv_re + wwc_im * b_fv_im; |
msamadani | 1:a514e4de034d | 419 | y->data[nd2].im = wwc_re * fv_im - b_wwc_im * b_fv_re; |
msamadani | 1:a514e4de034d | 420 | e = wwc->data[k].re; |
msamadani | 1:a514e4de034d | 421 | fv_re = fv->data[k].re; |
msamadani | 1:a514e4de034d | 422 | wwc_im = wwc->data[k].im; |
msamadani | 1:a514e4de034d | 423 | b_fv_im = fv->data[k].im; |
msamadani | 1:a514e4de034d | 424 | wwc_re = wwc->data[k].re; |
msamadani | 1:a514e4de034d | 425 | fv_im = fv->data[k].im; |
msamadani | 1:a514e4de034d | 426 | b_wwc_im = wwc->data[k].im; |
msamadani | 1:a514e4de034d | 427 | b_fv_re = fv->data[k].re; |
msamadani | 1:a514e4de034d | 428 | y->data[nd2].re = e * fv_re + wwc_im * b_fv_im; |
msamadani | 1:a514e4de034d | 429 | y->data[nd2].im = wwc_re * fv_im - b_wwc_im * b_fv_re; |
msamadani | 1:a514e4de034d | 430 | e = y->data[nd2].re; |
msamadani | 1:a514e4de034d | 431 | wwc_im = y->data[nd2].im; |
msamadani | 1:a514e4de034d | 432 | if (wwc_im == 0.0F) { |
msamadani | 1:a514e4de034d | 433 | y->data[nd2].re = e / denom_re; |
msamadani | 1:a514e4de034d | 434 | y->data[nd2].im = 0.0F; |
msamadani | 1:a514e4de034d | 435 | } else if (e == 0.0F) { |
msamadani | 1:a514e4de034d | 436 | y->data[nd2].re = 0.0F; |
msamadani | 1:a514e4de034d | 437 | y->data[nd2].im = wwc_im / denom_re; |
msamadani | 1:a514e4de034d | 438 | } else { |
msamadani | 1:a514e4de034d | 439 | y->data[nd2].re = e / denom_re; |
msamadani | 1:a514e4de034d | 440 | y->data[nd2].im = wwc_im / denom_re; |
msamadani | 1:a514e4de034d | 441 | } |
msamadani | 1:a514e4de034d | 442 | |
msamadani | 1:a514e4de034d | 443 | nd2++; |
msamadani | 1:a514e4de034d | 444 | k++; |
msamadani | 1:a514e4de034d | 445 | } |
msamadani | 1:a514e4de034d | 446 | |
msamadani | 1:a514e4de034d | 447 | emxFree_creal32_T(&fv); |
msamadani | 1:a514e4de034d | 448 | emxFree_creal32_T(&wwc); |
msamadani | 1:a514e4de034d | 449 | } |
msamadani | 1:a514e4de034d | 450 | |
msamadani | 1:a514e4de034d | 451 | emxFree_real32_T(&sintabinv); |
msamadani | 1:a514e4de034d | 452 | emxFree_real32_T(&sintab); |
msamadani | 1:a514e4de034d | 453 | emxFree_real32_T(&costab); |
msamadani | 1:a514e4de034d | 454 | } |
msamadani | 1:a514e4de034d | 455 | } |
msamadani | 1:a514e4de034d | 456 | |
msamadani | 1:a514e4de034d | 457 | // |
msamadani | 1:a514e4de034d | 458 | // Arguments : const emxArray_creal32_T *x |
msamadani | 1:a514e4de034d | 459 | // int n1_unsigned |
msamadani | 1:a514e4de034d | 460 | // const emxArray_real32_T *costab |
msamadani | 1:a514e4de034d | 461 | // const emxArray_real32_T *sintab |
msamadani | 1:a514e4de034d | 462 | // emxArray_creal32_T *y |
msamadani | 1:a514e4de034d | 463 | // Return Type : void |
msamadani | 1:a514e4de034d | 464 | // |
msamadani | 1:a514e4de034d | 465 | void r2br_r2dit_trig(const emxArray_creal32_T *x, int n1_unsigned, const |
msamadani | 1:a514e4de034d | 466 | emxArray_real32_T *costab, const emxArray_real32_T *sintab, |
msamadani | 1:a514e4de034d | 467 | emxArray_creal32_T *y) |
msamadani | 1:a514e4de034d | 468 | { |
msamadani | 1:a514e4de034d | 469 | int j; |
msamadani | 1:a514e4de034d | 470 | int nRowsD2; |
msamadani | 1:a514e4de034d | 471 | int nRowsD4; |
msamadani | 1:a514e4de034d | 472 | int iy; |
msamadani | 1:a514e4de034d | 473 | int iDelta; |
msamadani | 1:a514e4de034d | 474 | int ix; |
msamadani | 1:a514e4de034d | 475 | int ju; |
msamadani | 1:a514e4de034d | 476 | int i; |
msamadani | 1:a514e4de034d | 477 | boolean_T tst; |
msamadani | 1:a514e4de034d | 478 | float temp_re; |
msamadani | 1:a514e4de034d | 479 | float temp_im; |
msamadani | 1:a514e4de034d | 480 | float twid_re; |
msamadani | 1:a514e4de034d | 481 | float twid_im; |
msamadani | 1:a514e4de034d | 482 | int ihi; |
msamadani | 1:a514e4de034d | 483 | j = x->size[0]; |
msamadani | 1:a514e4de034d | 484 | if (!(j < n1_unsigned)) { |
msamadani | 1:a514e4de034d | 485 | j = n1_unsigned; |
msamadani | 1:a514e4de034d | 486 | } |
msamadani | 1:a514e4de034d | 487 | |
msamadani | 1:a514e4de034d | 488 | nRowsD2 = n1_unsigned / 2; |
msamadani | 1:a514e4de034d | 489 | nRowsD4 = nRowsD2 / 2; |
msamadani | 1:a514e4de034d | 490 | iy = y->size[0]; |
msamadani | 1:a514e4de034d | 491 | y->size[0] = n1_unsigned; |
msamadani | 1:a514e4de034d | 492 | emxEnsureCapacity((emxArray__common *)y, iy, sizeof(creal32_T)); |
msamadani | 1:a514e4de034d | 493 | if (n1_unsigned > x->size[0]) { |
msamadani | 1:a514e4de034d | 494 | iDelta = y->size[0]; |
msamadani | 1:a514e4de034d | 495 | iy = y->size[0]; |
msamadani | 1:a514e4de034d | 496 | y->size[0] = iDelta; |
msamadani | 1:a514e4de034d | 497 | emxEnsureCapacity((emxArray__common *)y, iy, sizeof(creal32_T)); |
msamadani | 1:a514e4de034d | 498 | for (iy = 0; iy < iDelta; iy++) { |
msamadani | 1:a514e4de034d | 499 | y->data[iy].re = 0.0F; |
msamadani | 1:a514e4de034d | 500 | y->data[iy].im = 0.0F; |
msamadani | 1:a514e4de034d | 501 | } |
msamadani | 1:a514e4de034d | 502 | } |
msamadani | 1:a514e4de034d | 503 | |
msamadani | 1:a514e4de034d | 504 | ix = 0; |
msamadani | 1:a514e4de034d | 505 | ju = 0; |
msamadani | 1:a514e4de034d | 506 | iy = 0; |
msamadani | 1:a514e4de034d | 507 | for (i = 1; i < j; i++) { |
msamadani | 1:a514e4de034d | 508 | y->data[iy] = x->data[ix]; |
msamadani | 1:a514e4de034d | 509 | iDelta = n1_unsigned; |
msamadani | 1:a514e4de034d | 510 | tst = true; |
msamadani | 1:a514e4de034d | 511 | while (tst) { |
msamadani | 1:a514e4de034d | 512 | iDelta >>= 1; |
msamadani | 1:a514e4de034d | 513 | ju ^= iDelta; |
msamadani | 1:a514e4de034d | 514 | tst = ((ju & iDelta) == 0); |
msamadani | 1:a514e4de034d | 515 | } |
msamadani | 1:a514e4de034d | 516 | |
msamadani | 1:a514e4de034d | 517 | iy = ju; |
msamadani | 1:a514e4de034d | 518 | ix++; |
msamadani | 1:a514e4de034d | 519 | } |
msamadani | 1:a514e4de034d | 520 | |
msamadani | 1:a514e4de034d | 521 | y->data[iy] = x->data[ix]; |
msamadani | 1:a514e4de034d | 522 | if (n1_unsigned > 1) { |
msamadani | 1:a514e4de034d | 523 | for (i = 0; i <= n1_unsigned - 2; i += 2) { |
msamadani | 1:a514e4de034d | 524 | temp_re = y->data[i + 1].re; |
msamadani | 1:a514e4de034d | 525 | temp_im = y->data[i + 1].im; |
msamadani | 1:a514e4de034d | 526 | y->data[i + 1].re = y->data[i].re - y->data[i + 1].re; |
msamadani | 1:a514e4de034d | 527 | y->data[i + 1].im = y->data[i].im - y->data[i + 1].im; |
msamadani | 1:a514e4de034d | 528 | y->data[i].re += temp_re; |
msamadani | 1:a514e4de034d | 529 | y->data[i].im += temp_im; |
msamadani | 1:a514e4de034d | 530 | } |
msamadani | 1:a514e4de034d | 531 | } |
msamadani | 1:a514e4de034d | 532 | |
msamadani | 1:a514e4de034d | 533 | iDelta = 2; |
msamadani | 1:a514e4de034d | 534 | iy = 4; |
msamadani | 1:a514e4de034d | 535 | ix = 1 + ((nRowsD4 - 1) << 2); |
msamadani | 1:a514e4de034d | 536 | while (nRowsD4 > 0) { |
msamadani | 1:a514e4de034d | 537 | for (i = 0; i < ix; i += iy) { |
msamadani | 1:a514e4de034d | 538 | temp_re = y->data[i + iDelta].re; |
msamadani | 1:a514e4de034d | 539 | temp_im = y->data[i + iDelta].im; |
msamadani | 1:a514e4de034d | 540 | y->data[i + iDelta].re = y->data[i].re - temp_re; |
msamadani | 1:a514e4de034d | 541 | y->data[i + iDelta].im = y->data[i].im - temp_im; |
msamadani | 1:a514e4de034d | 542 | y->data[i].re += temp_re; |
msamadani | 1:a514e4de034d | 543 | y->data[i].im += temp_im; |
msamadani | 1:a514e4de034d | 544 | } |
msamadani | 1:a514e4de034d | 545 | |
msamadani | 1:a514e4de034d | 546 | ju = 1; |
msamadani | 1:a514e4de034d | 547 | for (j = nRowsD4; j < nRowsD2; j += nRowsD4) { |
msamadani | 1:a514e4de034d | 548 | twid_re = costab->data[j]; |
msamadani | 1:a514e4de034d | 549 | twid_im = sintab->data[j]; |
msamadani | 1:a514e4de034d | 550 | i = ju; |
msamadani | 1:a514e4de034d | 551 | ihi = ju + ix; |
msamadani | 1:a514e4de034d | 552 | while (i < ihi) { |
msamadani | 1:a514e4de034d | 553 | temp_re = twid_re * y->data[i + iDelta].re - twid_im * y->data[i + |
msamadani | 1:a514e4de034d | 554 | iDelta].im; |
msamadani | 1:a514e4de034d | 555 | temp_im = twid_re * y->data[i + iDelta].im + twid_im * y->data[i + |
msamadani | 1:a514e4de034d | 556 | iDelta].re; |
msamadani | 1:a514e4de034d | 557 | y->data[i + iDelta].re = y->data[i].re - temp_re; |
msamadani | 1:a514e4de034d | 558 | y->data[i + iDelta].im = y->data[i].im - temp_im; |
msamadani | 1:a514e4de034d | 559 | y->data[i].re += temp_re; |
msamadani | 1:a514e4de034d | 560 | y->data[i].im += temp_im; |
msamadani | 1:a514e4de034d | 561 | i += iy; |
msamadani | 1:a514e4de034d | 562 | } |
msamadani | 1:a514e4de034d | 563 | |
msamadani | 1:a514e4de034d | 564 | ju++; |
msamadani | 1:a514e4de034d | 565 | } |
msamadani | 1:a514e4de034d | 566 | |
msamadani | 1:a514e4de034d | 567 | nRowsD4 /= 2; |
msamadani | 1:a514e4de034d | 568 | iDelta = iy; |
msamadani | 1:a514e4de034d | 569 | iy += iy; |
msamadani | 1:a514e4de034d | 570 | ix -= iDelta; |
msamadani | 1:a514e4de034d | 571 | } |
msamadani | 1:a514e4de034d | 572 | } |
msamadani | 1:a514e4de034d | 573 | |
msamadani | 1:a514e4de034d | 574 | // |
msamadani | 1:a514e4de034d | 575 | // Arguments : const emxArray_creal32_T *x |
msamadani | 1:a514e4de034d | 576 | // int xoffInit |
msamadani | 1:a514e4de034d | 577 | // int unsigned_nRows |
msamadani | 1:a514e4de034d | 578 | // const emxArray_real32_T *costab |
msamadani | 1:a514e4de034d | 579 | // const emxArray_real32_T *sintab |
msamadani | 1:a514e4de034d | 580 | // emxArray_creal32_T *y |
msamadani | 1:a514e4de034d | 581 | // Return Type : void |
msamadani | 1:a514e4de034d | 582 | // |
msamadani | 1:a514e4de034d | 583 | void r2br_r2dit_trig_impl(const emxArray_creal32_T *x, int xoffInit, int |
msamadani | 1:a514e4de034d | 584 | unsigned_nRows, const emxArray_real32_T *costab, const emxArray_real32_T |
msamadani | 1:a514e4de034d | 585 | *sintab, emxArray_creal32_T *y) |
msamadani | 1:a514e4de034d | 586 | { |
msamadani | 1:a514e4de034d | 587 | int j; |
msamadani | 1:a514e4de034d | 588 | int nRowsD2; |
msamadani | 1:a514e4de034d | 589 | int nRowsD4; |
msamadani | 1:a514e4de034d | 590 | int iy; |
msamadani | 1:a514e4de034d | 591 | int iDelta; |
msamadani | 1:a514e4de034d | 592 | int ix; |
msamadani | 1:a514e4de034d | 593 | int ju; |
msamadani | 1:a514e4de034d | 594 | int i; |
msamadani | 1:a514e4de034d | 595 | boolean_T tst; |
msamadani | 1:a514e4de034d | 596 | float temp_re; |
msamadani | 1:a514e4de034d | 597 | float temp_im; |
msamadani | 1:a514e4de034d | 598 | float twid_re; |
msamadani | 1:a514e4de034d | 599 | float twid_im; |
msamadani | 1:a514e4de034d | 600 | int ihi; |
msamadani | 1:a514e4de034d | 601 | j = x->size[0]; |
msamadani | 1:a514e4de034d | 602 | if (!(j < unsigned_nRows)) { |
msamadani | 1:a514e4de034d | 603 | j = unsigned_nRows; |
msamadani | 1:a514e4de034d | 604 | } |
msamadani | 1:a514e4de034d | 605 | |
msamadani | 1:a514e4de034d | 606 | nRowsD2 = unsigned_nRows / 2; |
msamadani | 1:a514e4de034d | 607 | nRowsD4 = nRowsD2 / 2; |
msamadani | 1:a514e4de034d | 608 | iy = y->size[0]; |
msamadani | 1:a514e4de034d | 609 | y->size[0] = unsigned_nRows; |
msamadani | 1:a514e4de034d | 610 | emxEnsureCapacity((emxArray__common *)y, iy, sizeof(creal32_T)); |
msamadani | 1:a514e4de034d | 611 | if (unsigned_nRows > x->size[0]) { |
msamadani | 1:a514e4de034d | 612 | iDelta = y->size[0]; |
msamadani | 1:a514e4de034d | 613 | iy = y->size[0]; |
msamadani | 1:a514e4de034d | 614 | y->size[0] = iDelta; |
msamadani | 1:a514e4de034d | 615 | emxEnsureCapacity((emxArray__common *)y, iy, sizeof(creal32_T)); |
msamadani | 1:a514e4de034d | 616 | for (iy = 0; iy < iDelta; iy++) { |
msamadani | 1:a514e4de034d | 617 | y->data[iy].re = 0.0F; |
msamadani | 1:a514e4de034d | 618 | y->data[iy].im = 0.0F; |
msamadani | 1:a514e4de034d | 619 | } |
msamadani | 1:a514e4de034d | 620 | } |
msamadani | 1:a514e4de034d | 621 | |
msamadani | 1:a514e4de034d | 622 | ix = xoffInit; |
msamadani | 1:a514e4de034d | 623 | ju = 0; |
msamadani | 1:a514e4de034d | 624 | iy = 0; |
msamadani | 1:a514e4de034d | 625 | for (i = 1; i < j; i++) { |
msamadani | 1:a514e4de034d | 626 | y->data[iy] = x->data[ix]; |
msamadani | 1:a514e4de034d | 627 | iDelta = unsigned_nRows; |
msamadani | 1:a514e4de034d | 628 | tst = true; |
msamadani | 1:a514e4de034d | 629 | while (tst) { |
msamadani | 1:a514e4de034d | 630 | iDelta >>= 1; |
msamadani | 1:a514e4de034d | 631 | ju ^= iDelta; |
msamadani | 1:a514e4de034d | 632 | tst = ((ju & iDelta) == 0); |
msamadani | 1:a514e4de034d | 633 | } |
msamadani | 1:a514e4de034d | 634 | |
msamadani | 1:a514e4de034d | 635 | iy = ju; |
msamadani | 1:a514e4de034d | 636 | ix++; |
msamadani | 1:a514e4de034d | 637 | } |
msamadani | 1:a514e4de034d | 638 | |
msamadani | 1:a514e4de034d | 639 | y->data[iy] = x->data[ix]; |
msamadani | 1:a514e4de034d | 640 | if (unsigned_nRows > 1) { |
msamadani | 1:a514e4de034d | 641 | for (i = 0; i <= unsigned_nRows - 2; i += 2) { |
msamadani | 1:a514e4de034d | 642 | temp_re = y->data[i + 1].re; |
msamadani | 1:a514e4de034d | 643 | temp_im = y->data[i + 1].im; |
msamadani | 1:a514e4de034d | 644 | y->data[i + 1].re = y->data[i].re - y->data[i + 1].re; |
msamadani | 1:a514e4de034d | 645 | y->data[i + 1].im = y->data[i].im - y->data[i + 1].im; |
msamadani | 1:a514e4de034d | 646 | y->data[i].re += temp_re; |
msamadani | 1:a514e4de034d | 647 | y->data[i].im += temp_im; |
msamadani | 1:a514e4de034d | 648 | } |
msamadani | 1:a514e4de034d | 649 | } |
msamadani | 1:a514e4de034d | 650 | |
msamadani | 1:a514e4de034d | 651 | iDelta = 2; |
msamadani | 1:a514e4de034d | 652 | iy = 4; |
msamadani | 1:a514e4de034d | 653 | ix = 1 + ((nRowsD4 - 1) << 2); |
msamadani | 1:a514e4de034d | 654 | while (nRowsD4 > 0) { |
msamadani | 1:a514e4de034d | 655 | for (i = 0; i < ix; i += iy) { |
msamadani | 1:a514e4de034d | 656 | temp_re = y->data[i + iDelta].re; |
msamadani | 1:a514e4de034d | 657 | temp_im = y->data[i + iDelta].im; |
msamadani | 1:a514e4de034d | 658 | y->data[i + iDelta].re = y->data[i].re - temp_re; |
msamadani | 1:a514e4de034d | 659 | y->data[i + iDelta].im = y->data[i].im - temp_im; |
msamadani | 1:a514e4de034d | 660 | y->data[i].re += temp_re; |
msamadani | 1:a514e4de034d | 661 | y->data[i].im += temp_im; |
msamadani | 1:a514e4de034d | 662 | } |
msamadani | 1:a514e4de034d | 663 | |
msamadani | 1:a514e4de034d | 664 | ju = 1; |
msamadani | 1:a514e4de034d | 665 | for (j = nRowsD4; j < nRowsD2; j += nRowsD4) { |
msamadani | 1:a514e4de034d | 666 | twid_re = costab->data[j]; |
msamadani | 1:a514e4de034d | 667 | twid_im = sintab->data[j]; |
msamadani | 1:a514e4de034d | 668 | i = ju; |
msamadani | 1:a514e4de034d | 669 | ihi = ju + ix; |
msamadani | 1:a514e4de034d | 670 | while (i < ihi) { |
msamadani | 1:a514e4de034d | 671 | temp_re = twid_re * y->data[i + iDelta].re - twid_im * y->data[i + |
msamadani | 1:a514e4de034d | 672 | iDelta].im; |
msamadani | 1:a514e4de034d | 673 | temp_im = twid_re * y->data[i + iDelta].im + twid_im * y->data[i + |
msamadani | 1:a514e4de034d | 674 | iDelta].re; |
msamadani | 1:a514e4de034d | 675 | y->data[i + iDelta].re = y->data[i].re - temp_re; |
msamadani | 1:a514e4de034d | 676 | y->data[i + iDelta].im = y->data[i].im - temp_im; |
msamadani | 1:a514e4de034d | 677 | y->data[i].re += temp_re; |
msamadani | 1:a514e4de034d | 678 | y->data[i].im += temp_im; |
msamadani | 1:a514e4de034d | 679 | i += iy; |
msamadani | 1:a514e4de034d | 680 | } |
msamadani | 1:a514e4de034d | 681 | |
msamadani | 1:a514e4de034d | 682 | ju++; |
msamadani | 1:a514e4de034d | 683 | } |
msamadani | 1:a514e4de034d | 684 | |
msamadani | 1:a514e4de034d | 685 | nRowsD4 /= 2; |
msamadani | 1:a514e4de034d | 686 | iDelta = iy; |
msamadani | 1:a514e4de034d | 687 | iy += iy; |
msamadani | 1:a514e4de034d | 688 | ix -= iDelta; |
msamadani | 1:a514e4de034d | 689 | } |
msamadani | 1:a514e4de034d | 690 | } |
msamadani | 1:a514e4de034d | 691 | |
msamadani | 1:a514e4de034d | 692 | // |
msamadani | 1:a514e4de034d | 693 | // File trailer for ifft.cpp |
msamadani | 1:a514e4de034d | 694 | // |
msamadani | 1:a514e4de034d | 695 | // [EOF] |
msamadani | 1:a514e4de034d | 696 | // |
msamadani | 1:a514e4de034d | 697 |