luke
/
ngspice
mirror of https://git.code.sf.net/p/ngspice/ngspice
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

Merge branch 'pre-master' into bt_dev

This commit is contained in:
Brian Taylor 2023-07-03 13:31:00 -07:00
parent fe26418374 4c4a9c78ec
commit a2ca37401b
1 changed files with 5 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
src/maths/cmaths/cmath4.c
View File

@ -933,13 +933,17 @@ cx_ifft(void *data, short int type, int length, int *newlength, short int *newty
printf("IFFT: Frequency span: %g Hz, input length: %d\n", 1/span*length, length); printf("IFFT: Frequency span: %g Hz, input length: %d\n", 1/span*length, length);
printf("IFFT: Time resolution: %g s, output length: %d\n", span/(tpts-1), tpts); printf("IFFT: Time resolution: %g s, output length: %d\n", span/(tpts-1), tpts);
in = fftw_malloc(sizeof(fftw_complex) * (unsigned int) length); in = fftw_malloc(sizeof(fftw_complex) * (unsigned int) tpts);
out = fftw_malloc(sizeof(fftw_complex) * (unsigned int) tpts); out = fftw_malloc(sizeof(fftw_complex) * (unsigned int) tpts);
for (i = 0; i < length; i++) { for (i = 0; i < length; i++) {
in[i][0] = indata[i].cx_real; in[i][0] = indata[i].cx_real;
in[i][1] = indata[i].cx_imag; in[i][1] = indata[i].cx_imag;
} }
for (i = length; i < tpts; i++) {
in[i][0] = 0.0;
in[i][1] = 0.0;
}
plan_backward = fftw_plan_dft_1d(tpts, in, out, FFTW_BACKWARD, FFTW_ESTIMATE); plan_backward = fftw_plan_dft_1d(tpts, in, out, FFTW_BACKWARD, FFTW_ESTIMATE);