epilectrik
/
pyNgSpice
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity
Browse Source

XSPICE core: remove memory leak

h_vogt 14 years ago
parent
8d374157fe
commit
cb220b2d6f
2 changed files with 201 additions and 204 deletions
Whitespace
Show all changes Ignore whitespace when comparing lines Ignore changes in amount of whitespace Ignore changes in whitespace at EOL
Split View
Diff Options
Show Stats Download Patch File Download Diff File
  1. 393
      src/xspice/icm/xtradev/core/cfunc.mod
  2. 12
      src/xspice/icm/xtradev/core/ifspec.ifs

393
src/xspice/icm/xtradev/core/cfunc.mod
View File

@ -8,18 +8,18 @@ Georgia Tech Research Corporation, Atlanta, Ga. 30332
All Rights Reserved
PROJECT A-8503-405
AUTHORS
AUTHORS
24 Apr 1991 Jeffrey P. Murray
MODIFICATIONS
MODIFICATIONS
24 Apr 1991 Jeffrey P. Murray
26 Sep 1991 Jeffrey P. Murray
SUMMARY
@ -27,19 +27,19 @@ SUMMARY
code model.
INTERFACES
INTERFACES
FILE ROUTINE CALLED
FILE ROUTINE CALLED
CMmacros.h cm_message_send();
CMmacros.h cm_message_send();
CMutil.c void cm_smooth_corner();
CMutil.c void cm_smooth_corner();
REFERENCED FILES
Inputs from and outputs to ARGS structure.
NON-STANDARD FEATURES
@ -51,7 +51,7 @@ NON-STANDARD FEATURES
#include "core.h"
/*=== CONSTANTS ========================*/
@ -62,27 +62,38 @@ NON-STANDARD FEATURES
/*=== LOCAL VARIABLES & TYPEDEFS =======*/
/*=== LOCAL VARIABLES & TYPEDEFS =======*/
typedef struct {
double *H_array; /* the storage array for the
control vector (cntl_array) */
double *B_array; /* the storage array for the
pulse width array (pw_array) */
Boolean_t tran_init; /* for initialization of phase1) */
} Local_Data_t;
/*=== FUNCTION PROTOTYPE DEFINITIONS ===*/
/*==============================================================================
FUNCTION cm_core()
AUTHORS
AUTHORS
24 Apr 1991 Jeffrey P. Murray
MODIFICATIONS
MODIFICATIONS
24 Apr 1991 Jeffrey P. Murray
26 Sep 1991 Jeffrey P. Murray
@ -91,21 +102,21 @@ SUMMARY
This function implements the core code model.
INTERFACES
INTERFACES
FILE ROUTINE CALLED
FILE ROUTINE CALLED
CMmacros.h cm_message_send();
CMmacros.h cm_message_send();
CMutil.c void cm_smooth_corner();
CMutil.c void cm_smooth_corner();
RETURNED VALUE
Returns inputs and outputs via ARGS structure.
GLOBAL VARIABLES
NONE
NON-STANDARD FEATURES
@ -117,51 +128,51 @@ NON-STANDARD FEATURES
#include <stdlib.h>
/*=== CM_CORE ROUTINE ===*/
/*******************************************************************/
/* */
/* CORE Model: */
/* */
/* */
/* CORE Model: */
/* */
/* The core model is designed to operate in one of two modes. */
/* The first of these, and the one most likely to be used by */
/* the engineer, is a modified version of the pwl model. This */
/* behavior occurs when the model is in pwl mode (the default). */
/* If the model is set to hyst mode, its behavior mimics that of */
/* the hysteresis block. The following provides additional */
/* detail: */
/* */
/* detail: */
/* */
/* PWL Mode */
/* */
/* */
/* In pwl mode, the core model is a modified version of the */
/* PWL model... */
/* it has a single two-terminal input/output, and accepts as */
/* input the mmf value, represented by a voltage. Its output is */
/* a flux value, which is represented as a current. Additional */
/* inputs include the cross-sectional area of the physical */
/* core, and the median length of the core, seen from the */
/* perspective of the flux that traverses it. */
/* */
/* it has a single two-terminal input/output, and accepts as */
/* input the mmf value, represented by a voltage. Its output is */
/* a flux value, which is represented as a current. Additional */
/* inputs include the cross-sectional area of the physical */
/* core, and the median length of the core, seen from the */
/* perspective of the flux that traverses it. */
/* */
/* The core model in pwl mode DOES NOT include hysteresis... */
/* current thinking is that such provides */
/* little benefit to the designer, aside from the ability to */
/* calculate eddy losses in a modeled device...the nonlinear */
/* B vs. H behavior, however, is of great importance. */
/* */
/* Note that the user must input a piece-wise-linear */
/* description, in the form of a series of coordinate B vs. H */
/* values, in order to model a particular core material type. */
/* Such curves may be found in textbooks, or from manufacturer's */
/* databooks. In this model, the "x" values are assumed to */
/* represent the magnetic field (H), and the "y" values are */
/* assumed to represent the flux density (B). */
/* */
/* little benefit to the designer, aside from the ability to */
/* calculate eddy losses in a modeled device...the nonlinear */
/* B vs. H behavior, however, is of great importance. */
/* */
/* Note that the user must input a piece-wise-linear */
/* description, in the form of a series of coordinate B vs. H */
/* values, in order to model a particular core material type. */
/* Such curves may be found in textbooks, or from manufacturer's */
/* databooks. In this model, the "x" values are assumed to */
/* represent the magnetic field (H), and the "y" values are */
/* assumed to represent the flux density (B). */
/* */
/* Hyst Mode */
/* */
/* */
/* In hyst mode, the core model is a modified version of the */
/* HYST code model... */
/* it has a single two-terminal input/output, and accepts as */
/* input the mmf value, represented by a voltage. Its output is */
/* a flux value, which is represented as a current. Additional */
/* it has a single two-terminal input/output, and accepts as */
/* input the mmf value, represented by a voltage. Its output is */
/* a flux value, which is represented as a current. Additional */
/* inputs include the input high and low values for the */
/* hysteretic behavior, and the output high and low values. */
/* Also, a value of hysteresis must be included, as must an */
@ -176,12 +187,12 @@ NON-STANDARD FEATURES
/* the hysteresis capability will be of only nominal benefit to */
/* the engineer, as it will not typically allow for as accurate */
/* tailoring of the response as is possible in the pwl mode. */
/* */
/* 4/24/91 J.P.Murray */
/* Last modified: 10/24/91 */
/* */
/* 4/24/91 J.P.Murray */
/* Last modified: 10/24/91 */
/*******************************************************************/
void cm_core(ARGS) /* structure holding parms,
void cm_core(ARGS) /* structure holding parms,
inputs, outputs, etc. */
{
@ -213,7 +224,7 @@ void cm_core(ARGS) /* structure holding parms,
double length; /* length of core (in meters) */
Mif_Complex_t ac_gain;
char *allocation_error="\n***ERROR***\nCORE: Allocation calloc failed!\n";
char *limit_error="\n***ERROR***\nCORE: Violation of 50% rule in breakpoints!\n";
@ -222,34 +233,33 @@ void cm_core(ARGS) /* structure holding parms,
/*** The following declarations pertain to HYSTERESIS mode... ***/
double in, /* input to hysteresis block */
out, /* output from hysteresis block */
in_low, /* lower input value for hyst=0 at which
the transfer curve changes from constant
out, /* output from hysteresis block */
in_low, /* lower input value for hyst=0 at which
the transfer curve changes from constant
to linear */
in_high, /* upper input value for hyst=0 at which
the transfer curve changes from constant
in_high, /* upper input value for hyst=0 at which
the transfer curve changes from constant
to linear */
hyst, /* the hysteresis value (see above diagram) */
out_lower_limit, /* the minimum output value from the block */
out_upper_limit, /* the maximum output value from the block */
slope, /* calculated rise and fall slope for the block */
pout_pin, /* partial derivative of output w.r.t. input */
x_rise_linear, /* = in_low + hyst */
x_rise_zero, /* = in_high + hyst */
x_fall_linear, /* = in_high - hyst */
x_fall_zero; /* = in_low - hyst */
Boolean_t *hyst_state, /* TRUE => input is on lower leg of
hyst, /* the hysteresis value (see above diagram) */
out_lower_limit, /* the minimum output value from the block */
out_upper_limit, /* the maximum output value from the block */
slope, /* calculated rise and fall slope for the block */
pout_pin, /* partial derivative of output w.r.t. input */
x_rise_linear, /* = in_low + hyst */
x_rise_zero, /* = in_high + hyst */
x_fall_linear, /* = in_high - hyst */
x_fall_zero; /* = in_low - hyst */
Boolean_t *hyst_state, /* TRUE => input is on lower leg of
hysteresis curve, between -infinity
and in_high + hyst.
FALSE => input is on upper leg
of hysteresis curve, between
of hysteresis curve, between
in_low - hyst and +infinity */
*old_hyst_state; /* previous value of *hyst_state */
*old_hyst_state; /* previous value of *hyst_state */
Local_Data_t *loc; /* Pointer to local static data, not to be included
in the state vector */
/* Retrieve mode parameter... */
@ -261,7 +271,7 @@ void cm_core(ARGS) /* structure holding parms,
/******** pwl mode *****************/
if ( HYSTERESIS != mode ) {
if ( HYSTERESIS != mode ) {
/* Retrieve frequently used parameters... */
@ -269,43 +279,53 @@ void cm_core(ARGS) /* structure holding parms,
area = PARAM(area);
length = PARAM(length);
size = PARAM_SIZE(H_array);
size = PARAM_SIZE(H_array);
/* Allocate storage for breakpoint domain & range values */
H = (double *) calloc((size_t) size, sizeof(double));
if (!H) {
cm_message_send(allocation_error);
return;
}
B = (double *) calloc((size_t) size, sizeof(double));
if (!B) {
cm_message_send(allocation_error);
return;
if(INIT==1) {
/*** allocate static storage for *loc ***/
STATIC_VAR (locdata) = calloc (1 , sizeof ( Local_Data_t ));
loc = STATIC_VAR (locdata);
/* Allocate storage for breakpoint domain & range values */
H = loc->H_array = (double *) calloc((size_t) size, sizeof(double));
if (!H) {
cm_message_send(allocation_error);
return;
}
B = loc->B_array = (double *) calloc((size_t) size, sizeof(double));
if (!B) {
cm_message_send(allocation_error);
return;
}
}
loc = STATIC_VAR (locdata);
H = loc->H_array;
B = loc->B_array;
/* Retrieve H and B values. */
/* Retrieve H and B values. */
for (i=0; i<size; i++) {
*(H+i) = PARAM(H_array[i]);
*(B+i) = PARAM(B_array[i]);
}
}
/* See if input_domain is absolute...if so, test against */
/* breakpoint segments for violation of 50% rule... */
if (PARAM(fraction) == MIF_FALSE) {
for (i=0; i<(size-1); i++) {
for (i=0; i<(size-1); i++) {
if ( (*(H+i+1) - *(H+i)) < (2.0*input_domain) ) {
cm_message_send(limit_error);
return;
}
cm_message_send(limit_error);
return;
}
}
}
/* Retrieve mmf_input value. */
/* Retrieve mmf_input value. */
mmf_input = INPUT(mc);
/* Calculate H_input value from mmf_input... */
H_input = mmf_input / length;
@ -316,23 +336,21 @@ void cm_core(ARGS) /* structure holding parms,
if (H_input <= (*(H+1) + *H)/2.0) {/*** H_input below lowest midpoint ***/
dout_din = (*(B+1) - *B)/(*(H+1) - *H);
B_out = *B + (H_input - *H) * dout_din;
}
else {
if (H_input >= (*(H+size-2) + *(H+size-1))/2.0) {
/*** H_input above highest midpoint ***/
} else {
if (H_input >= (*(H+size-2) + *(H+size-1))/2.0) {
/*** H_input above highest midpoint ***/
dout_din = (*(B+size-1) - *(B+size-2)) /
(*(H+size-1) - *(H+size-2));
(*(H+size-1) - *(H+size-2));
B_out = *(B+size-1) + (H_input - *(H+size-1)) * dout_din;
}
else { /*** H_input within bounds of end midpoints... ***/
/*** must determine position progressively & then ***/
/*** calculate required output. ***/
} else { /*** H_input within bounds of end midpoints... ***/
/*** must determine position progressively & then ***/
/*** calculate required output. ***/
for (i=1; i<size; i++) {
if (H_input < (*(H+i) + *(H+i+1))/2.0) {
/* approximate position known... */
if (H_input < (*(H+i) + *(H+i+1))/2.0) {
/* approximate position known... */
lower_seg = (*(H+i) - *(H+i-1));
upper_seg = (*(H+i+1) - *(H+i));
@ -340,17 +358,17 @@ void cm_core(ARGS) /* structure holding parms,
/* Calculate input_domain about this region's breakpoint.*/
if (PARAM(fraction) == MIF_TRUE) { /* Translate input_domain */
/* into an absolute.... */
/* into an absolute.... */
if ( lower_seg <= upper_seg ) /* Use lower */
/* segment */
/* for % calc.*/
/* segment */
/* for % calc.*/
input_domain = input_domain * lower_seg;
else /* Use upper */
/* segment */
/* for % calc.*/
/* segment */
/* for % calc.*/
input_domain = input_domain * upper_seg;
}
}
/* Set up threshold values about breakpoint... */
threshold_lower = *(H+i) - input_domain;
threshold_upper = *(H+i) + input_domain;
@ -360,29 +378,27 @@ void cm_core(ARGS) /* structure holding parms,
if (H_input < threshold_lower) { /* Lower linear region */
dout_din = (*(B+i) - *(B+i-1))/lower_seg;
B_out = *(B+i) + (H_input - *(H+i)) * dout_din;
}
else {
} else {
if (H_input < threshold_upper) { /* Parabolic region */
lower_slope = (*(B+i) - *(B+i-1))/lower_seg;
upper_slope = (*(B+i+1) - *(B+i))/upper_seg;
cm_smooth_corner(H_input,*(H+i),*(B+i),input_domain,
lower_slope,upper_slope,&B_out,&dout_din);
}
else { /* Upper linear region */
lower_slope,upper_slope,&B_out,&dout_din);
} else { /* Upper linear region */
dout_din = (*(B+i+1) - *(B+i))/upper_seg;
B_out = *(B+i) + (H_input - *(H+i)) * dout_din;
}
}
break; /* Break search loop...H_input has been found, */
/* and B_out and dout_din have been assigned. */
/* and B_out and dout_din have been assigned. */
}
}
}
}
/* Calculate value of flux_out... */
flux_out = B_out * area;
/* Adjust dout_din value to reflect area and length multipliers... */
dout_din = dout_din * area / length;
@ -391,54 +407,51 @@ void cm_core(ARGS) /* structure holding parms,
if(ANALYSIS != MIF_AC) { /* Output DC & Transient Values */
OUTPUT(mc) = flux_out;
PARTIAL(mc,mc) = dout_din;
}
else { /* Output AC Gain */
} else { /* Output AC Gain */
ac_gain.real = dout_din;
ac_gain.imag= 0.0;
AC_GAIN(mc,mc) = ac_gain;
}
}
}
/******** hysteresis mode ******************/
else {
else {
/** Retrieve frequently used parameters... **/
in_low = PARAM(in_low);
in_high = PARAM(in_high);
hyst = PARAM(hyst);
out_lower_limit = PARAM(out_lower_limit);
out_upper_limit = PARAM(out_upper_limit);
out_upper_limit = PARAM(out_upper_limit);
input_domain = PARAM(input_domain);
/** Calculate Hysteresis Linear Region Slopes & Derived Values **/
/* Define slope of rise and fall lines when not being smoothed */
slope = (out_upper_limit - out_lower_limit)/(in_high - in_low);
x_rise_linear = in_low + hyst; /* Breakpoint - x rising to
slope = (out_upper_limit - out_lower_limit)/(in_high - in_low);
x_rise_linear = in_low + hyst; /* Breakpoint - x rising to
linear region */
x_rise_zero = in_high + hyst; /* Breakpoint - x rising to
x_rise_zero = in_high + hyst; /* Breakpoint - x rising to
zero-slope (out_upper_limit) */
x_fall_linear = in_high - hyst; /* Breakpoint - x falling to
x_fall_linear = in_high - hyst; /* Breakpoint - x falling to
linear region */
x_fall_zero = in_low - hyst; /* Breakpoint - x falling to
x_fall_zero = in_low - hyst; /* Breakpoint - x falling to
zero-slope (out_lower_limit) */
/* Set range to absolute value */
if (PARAM(fraction) == MIF_TRUE)
if (PARAM(fraction) == MIF_TRUE)
input_domain = input_domain * (in_high - in_low);
/** Retrieve frequently used inputs... **/
in = INPUT(mc);
@ -449,31 +462,29 @@ void cm_core(ARGS) /* structure holding parms,
/* First pass...allocate storage for previous state. */
/* Also, calculate roughly where the current output */
/* will be and use this value to define current state. */
if (INIT==1) {
cm_analog_alloc(TRUE,sizeof(Boolean_t));
if (INIT==1) {
cm_analog_alloc(TRUE,sizeof(Boolean_t));
hyst_state = (Boolean_t *) cm_analog_get_ptr(TRUE,0);
old_hyst_state = (Boolean_t *) cm_analog_get_ptr(TRUE,1);
if (in < x_rise_zero + input_domain) { /* Set state to X_RISING */
*old_hyst_state = X_RISING;
*old_hyst_state = X_RISING;
} else {
*old_hyst_state = X_FALLING;
}
else {
*old_hyst_state = X_FALLING;
}
}
else { /* Allocation not necessary...retrieve previous values */
} else { /* Allocation not necessary...retrieve previous values */
hyst_state = (Boolean_t *) cm_analog_get_ptr(TRUE,0); /* Set out pointer to current
time storage */
old_hyst_state = (Boolean_t *) cm_analog_get_ptr(TRUE,1); /* Set old-output-state pointer
to previous time storage */
hyst_state = (Boolean_t *) cm_analog_get_ptr(TRUE,0); /* Set out pointer to current
time storage */
old_hyst_state = (Boolean_t *) cm_analog_get_ptr(TRUE,1); /* Set old-output-state pointer
to previous time storage */
}
/** Set *hyst_out = *old_hyst_out, unless changed below...
we don't need the last iteration value of *hyst_state. **/
@ -486,29 +497,25 @@ void cm_core(ARGS) /* structure holding parms,
/*** Calculate value of hyst_state, pout_pin.... ***/
if (*old_hyst_state == X_RISING) { /* Assume calculations on lower */
/* hysteresis section (x rising) */
/* hysteresis section (x rising) */
if ( in <= x_rise_linear - input_domain ) { /* Output @ lower limit */
out = out_lower_limit;
pout_pin = 0.0;
}
else {
} else {
if ( in <= x_rise_linear + input_domain ) { /* lower smoothing region */
cm_smooth_corner(in,x_rise_linear,out_lower_limit,input_domain,
0.0,slope,&out,&pout_pin);
}
else {
if (in <= x_rise_zero - input_domain) { /* Rising linear region */
0.0,slope,&out,&pout_pin);
} else {
if (in <= x_rise_zero - input_domain) { /* Rising linear region */
out = (in - x_rise_linear)*slope + out_lower_limit;
pout_pin = slope;
}
else {
if (in <= x_rise_zero + input_domain) { /* Upper smoothing region */
} else {
if (in <= x_rise_zero + input_domain) { /* Upper smoothing region */
cm_smooth_corner(in,x_rise_zero,out_upper_limit,input_domain,
slope,0.0,&out,&pout_pin);
}
else { /* input has transitioned to X_FALLING region... */
slope,0.0,&out,&pout_pin);
} else { /* input has transitioned to X_FALLING region... */
out = out_upper_limit;
pout_pin = 0.0;
*hyst_state = X_FALLING;
@ -516,30 +523,25 @@ void cm_core(ARGS) /* structure holding parms,
}
}
}
}
else { /* Assume calculations on upper hysteresis section (x falling) */
} else { /* Assume calculations on upper hysteresis section (x falling) */
if ( in >= x_fall_linear + input_domain ) { /* Output @ upper limit */
out = out_upper_limit;
pout_pin = 0.0;
}
else {
} else {
if ( in >= x_fall_linear - input_domain ) { /* Upper smoothing region */
cm_smooth_corner(in,x_fall_linear,out_upper_limit,input_domain,
slope,0.0,&out,&pout_pin);
}
else {
if (in >= x_fall_zero + input_domain) { /* Falling linear region */
slope,0.0,&out,&pout_pin);
} else {
if (in >= x_fall_zero + input_domain) { /* Falling linear region */
out = (in - x_fall_zero)*slope + out_lower_limit;
pout_pin = slope;
}
else {
if (in >= x_fall_zero - input_domain) { /* Lower smoothing region */
} else {
if (in >= x_fall_zero - input_domain) { /* Lower smoothing region */
cm_smooth_corner(in,x_fall_zero,out_lower_limit,input_domain,
0.0,slope,&out,&pout_pin);
}
else { /* input has transitioned to X_RISING region... */
0.0,slope,&out,&pout_pin);
} else { /* input has transitioned to X_RISING region... */
out = out_lower_limit;
pout_pin = 0.0;
*hyst_state = X_RISING;
@ -552,12 +554,11 @@ void cm_core(ARGS) /* structure holding parms,
if (ANALYSIS != MIF_AC) { /* DC & Transient Analyses */
OUTPUT(mc) = out;
PARTIAL(mc,mc) = pout_pin;
}
else { /* AC Analysis */
OUTPUT(mc) = out;
PARTIAL(mc,mc) = pout_pin;
} else { /* AC Analysis */
ac_gain.real = pout_pin;
ac_gain.imag= 0.0;
AC_GAIN(mc,mc) = ac_gain;

12
src/xspice/icm/xtradev/core/ifspec.ifs
View File

@ -131,13 +131,9 @@ Vector: no
Vector_Bounds: -
Null_Allowed: yes
STATIC_VAR_TABLE:
Static_Var_Name: locdata
Description: "local static data"
Data_Type: pointer
Write Preview
Loading…
Cancel
Save