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

Create a vector thdnm from THD calculation, e.g. thd11 

New vector fournosave: if set, suppress creation of
vectors thd11 and fourier, usefull when writing a raw file.
This commit is contained in:
Holger Vogt 2025-03-10 12:53:40 +01:00
parent 83d3890490
commit 16fbe0fb8b
1 changed files with 41 additions and 24 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

65
src/frontend/fourier.c
View File

@ -50,6 +50,7 @@ fourier(wordlist *wl, struct plot *current_plot)
char xbuf[20]; char xbuf[20];
int shift; int shift;
int rv = 1; int rv = 1;
bool foursave = TRUE;
struct dvec *n; struct dvec *n;
int newveccount = 1; int newveccount = 1;
@ -71,6 +72,8 @@ fourier(wordlist *wl, struct plot *current_plot)
polydegree = 1; polydegree = 1;
if (!cp_getvar("fourgridsize", CP_NUM, &fourgridsize, 0) || fourgridsize < 1) if (!cp_getvar("fourgridsize", CP_NUM, &fourgridsize, 0) || fourgridsize < 1)
fourgridsize = DEF_FOURGRIDSIZE; fourgridsize = DEF_FOURGRIDSIZE;
if (cp_getvar("fournosave", CP_BOOL, NULL, 0))
foursave = FALSE;
time = current_plot->pl_scale; time = current_plot->pl_scale;
if (!isreal(time)) { if (!isreal(time)) {
@ -188,34 +191,48 @@ fourier(wordlist *wl, struct plot *current_plot)
} }
fputs("\n", cp_out); fputs("\n", cp_out);
/* create and assign a new vector n */ if (foursave) {
/* with size 3 * nfreqs in current plot */ /* create a vector for THD */
/* generate name for new vector, using vec->name */ n = dvec_alloc(tprintf("thd%d%d", callstof, newveccount),
n = dvec_alloc(tprintf("fourier%d%d", callstof, newveccount), SV_NOTYPE,
SV_NOTYPE, VF_REAL | VF_PERMANENT,
VF_REAL | VF_PERMANENT, 1, NULL);
3 * nfreqs, NULL);
n->v_numdims = 2; n->v_numdims = 1;
n->v_dims[0] = 3;
n->v_dims[1] = nfreqs;
vec_new(n); vec_new(n);
/* store data in vector: freq, mag, phase */ n->v_realdata[0] = thd;
for (i = 0; i < nfreqs; i++) {
n->v_realdata[i] = freq[i]; /* create and assign a new vector n */
n->v_realdata[i + nfreqs] = mag[i]; /* with size 3 * nfreqs in current plot */
n->v_realdata[i + 2 * nfreqs] = phase[i]; /* generate name for new vector, using vec->name */
n = dvec_alloc(tprintf("fourier%d%d", callstof, newveccount),
SV_NOTYPE,
VF_REAL | VF_PERMANENT,
3 * nfreqs, NULL);
n->v_numdims = 2;
n->v_dims[0] = 3;
n->v_dims[1] = nfreqs;
vec_new(n);
/* store data in vector: freq, mag, phase */
for (i = 0; i < nfreqs; i++) {
n->v_realdata[i] = freq[i];
n->v_realdata[i + nfreqs] = mag[i];
n->v_realdata[i + 2 * nfreqs] = phase[i];
}
newveccount++;
if (polydegree) {
tfree(timescale);
tfree(data);
}
timescale = NULL;
data = NULL;
} }
newveccount++;
if (polydegree) {
tfree(timescale);
tfree(data);
}
timescale = NULL;
data = NULL;
} }
} }