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

Add breakpoints to achieve quadratic smoothing. 

Shift 50%-detection for breakpoints into the INIT section,
to do it only once.
This commit is contained in:
Holger Vogt 2025-10-10 20:44:21 +02:00
parent 172fb6d668
commit a1de300a82
1 changed files with 41 additions and 16 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

57
src/xspice/icm/analog/pwlts/cfunc.mod
View File

@ -13,6 +13,7 @@ AUTHORS
Pwl with time input and smoothing: pwlts
9 Sep 2022 Holger Vogt
10 Oct 2025 Holger Vogt
SUMMARY
@ -211,6 +212,46 @@ void cm_pwlts(ARGS) /* structure holding parms,
y[0] = 2. * y[1] - y[2];
y[size - 1] = 2. * y[size - 2] - y[size - 3];
}
/* See if input_domain is absolute...if so, test against */
/* breakpoint segments for violation of 50% rule... */
if (PARAM(fraction) == MIF_FALSE) {
if ( 3 < size ) {
for (i=1; i<(size-2); i++) {
/* Test for overlap...0.999999999 factor is to */
/* prevent floating point problems with comparison. */
if ( (test1 = x[i+1] - x[i]) <
(test2 = 0.999999999 * (2.0 * input_domain)) ) {
cm_message_send(limit_error);
}
}
}
}
/* add permanent breakpoints */
if (PARAM(fraction) == MIF_FALSE) {
for (i=1; i<size-1; i++) {
if (x[i] - input_domain <= 0)
continue;
cm_analog_set_perm_bkpt(x[i] - input_domain);
cm_analog_set_perm_bkpt(x[i]);
cm_analog_set_perm_bkpt(x[i] + input_domain);
}
}
else {
double dtlo, dthi, dt;
for (i=1; i<size-1; i++) {
dtlo = x[i] - x[i-1];
dthi = x[i+1] - x[i];
dt = ((dtlo >= dthi) ? dthi : dtlo) * input_domain;
if (x[i] - dt <= 0)
continue;
cm_analog_set_perm_bkpt(x[i] - dt);
cm_analog_set_perm_bkpt(x[i]);
cm_analog_set_perm_bkpt(x[i] + dt);
}
}
/* debug printout
for (i = 0; i < size; i++)
fprintf(stderr, "%e ", y[i]);
@ -229,22 +270,6 @@ void cm_pwlts(ARGS) /* structure holding parms,
}
/* See if input_domain is absolute...if so, test against */
/* breakpoint segments for violation of 50% rule... */
if (PARAM(fraction) == MIF_FALSE) {
if ( 3 < size ) {
for (i=1; i<(size-2); i++) {
/* Test for overlap...0.999999999 factor is to */
/* prevent floating point problems with comparison. */
if ( (test1 = x[i+1] - x[i]) <
(test2 = 0.999999999 * (2.0 * input_domain)) ) {
cm_message_send(limit_error);
}
}
}
}
/* Retrieve x_input value as current simulation time. */
x_input = TIME;