You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
240 lines
6.6 KiB
240 lines
6.6 KiB
/*.......1.........2.........3.........4.........5.........6.........7.........8
|
|
================================================================================
|
|
|
|
FILE int/cfunc.mod
|
|
|
|
Public Domain
|
|
|
|
Georgia Tech Research Corporation
|
|
Atlanta, Georgia 30332
|
|
PROJECT A-8503-405
|
|
|
|
|
|
AUTHORS
|
|
|
|
6 Nov 1991 Jeffrey P. Murray
|
|
|
|
|
|
MODIFICATIONS
|
|
|
|
2 Oct 1991 Jeffrey P. Murray
|
|
|
|
|
|
SUMMARY
|
|
|
|
This file contains the model-specific routines used to
|
|
functionally describe the int code model.
|
|
|
|
|
|
INTERFACES
|
|
|
|
FILE ROUTINE CALLED
|
|
|
|
CMutil.c void cm_smooth_corner();
|
|
|
|
CM.c void *cm_analog_alloc()
|
|
void *cm_analog_get_ptr()
|
|
int cm_analog_integrate()
|
|
|
|
|
|
REFERENCED FILES
|
|
|
|
Inputs from and outputs to ARGS structure.
|
|
|
|
|
|
NON-STANDARD FEATURES
|
|
|
|
NONE
|
|
|
|
===============================================================================*/
|
|
|
|
/*=== INCLUDE FILES ====================*/
|
|
|
|
#include "int.h"
|
|
|
|
|
|
|
|
/*=== CONSTANTS ========================*/
|
|
|
|
|
|
|
|
|
|
/*=== MACROS ===========================*/
|
|
|
|
|
|
|
|
|
|
/*=== LOCAL VARIABLES & TYPEDEFS =======*/
|
|
|
|
|
|
|
|
|
|
/*=== FUNCTION PROTOTYPE DEFINITIONS ===*/
|
|
|
|
|
|
|
|
|
|
|
|
/*==============================================================================
|
|
|
|
FUNCTION void cm_int()
|
|
|
|
AUTHORS
|
|
|
|
2 Oct 1991 Jeffrey P. Murray
|
|
|
|
MODIFICATIONS
|
|
|
|
NONE
|
|
|
|
SUMMARY
|
|
|
|
This function implements the int code model.
|
|
|
|
INTERFACES
|
|
|
|
FILE ROUTINE CALLED
|
|
|
|
CMutil.c void cm_smooth_corner();
|
|
|
|
CM.c void *cm_analog_alloc()
|
|
void *cm_analog_get_ptr()
|
|
int cm_analog_integrate()
|
|
|
|
RETURNED VALUE
|
|
|
|
Returns inputs and outputs via ARGS structure.
|
|
|
|
GLOBAL VARIABLES
|
|
|
|
NONE
|
|
|
|
NON-STANDARD FEATURES
|
|
|
|
NONE
|
|
|
|
==============================================================================*/
|
|
|
|
/*=== CM_INT ROUTINE ===*/
|
|
|
|
void cm_int(ARGS) /* structure holding parms,
|
|
inputs, outputs, etc. */
|
|
{
|
|
double *out, /* current output */
|
|
*in, /* input */
|
|
in_offset, /* input offset */
|
|
gain, /* gain parameter */
|
|
out_lower_limit, /* output lower limit */
|
|
out_upper_limit, /* output upper limit */
|
|
limit_range, /* range of output below out_upper_limit
|
|
and above out_lower_limit within which
|
|
smoothing will take place */
|
|
out_ic, /* output initial condition - initial output value */
|
|
pout_pin, /* partial derivative of output w.r.t. input */
|
|
pout_gain; /* temporary storage variable for partial
|
|
value returned by smoothing function
|
|
(subsequently multiplied by pout_pin) */
|
|
|
|
Mif_Complex_t ac_gain; /* AC gain */
|
|
|
|
|
|
|
|
/** Retrieve frequently used parameters (used by all analyses)... **/
|
|
|
|
gain = PARAM(gain);
|
|
|
|
|
|
|
|
if (ANALYSIS != MIF_AC) { /**** DC & Transient Analyses ****/
|
|
|
|
/** Retrieve frequently used parameters... **/
|
|
|
|
in_offset = PARAM(in_offset);
|
|
out_lower_limit = PARAM(out_lower_limit);
|
|
out_upper_limit = PARAM(out_upper_limit);
|
|
limit_range = PARAM(limit_range);
|
|
out_ic = PARAM(out_ic);
|
|
|
|
|
|
|
|
/** Test for INIT; if so, allocate storage, otherwise, retrieve
|
|
previous timepoint input value... **/
|
|
|
|
if (INIT==1) { /* First pass...allocate storage for previous value. */
|
|
|
|
cm_analog_alloc(INT1,sizeof(double));
|
|
cm_analog_alloc(INT2,sizeof(double));
|
|
}
|
|
/* retrieve previous value */
|
|
|
|
in = (double *) cm_analog_get_ptr(INT1,0); /* Set out pointer to input storage location */
|
|
out = (double *) cm_analog_get_ptr(INT2,0); /* Set out pointer to output storage location */
|
|
|
|
|
|
/*** Read input value for current time, and calculate pseudo-input ***/
|
|
/*** which includes input offset and gain.... ***/
|
|
|
|
*in = gain*(INPUT(in)+in_offset);
|
|
|
|
/*** Test to see if this is the first timepoint calculation... ***/
|
|
/*** this would imply that TIME equals zero. ***/
|
|
|
|
if ( 0.0 == TIME ) { /*** Test to see if this is the first ***/
|
|
/*** timepoint calculation...if ***/
|
|
*out = out_ic; /*** so, return out_ic. ***/
|
|
pout_pin = 0.0;
|
|
}
|
|
else { /*** Calculate value of integral.... ***/
|
|
cm_analog_integrate(*in,out,&pout_pin);
|
|
}
|
|
|
|
|
|
/*** Smooth output if it is within limit_range of
|
|
out_lower_limit or out_upper_limit. ***/
|
|
|
|
if (*out < (out_lower_limit - limit_range)) { /* At lower limit. */
|
|
*out = out_lower_limit;
|
|
pout_pin = 0.0;
|
|
}
|
|
else {
|
|
if (*out < (out_lower_limit + limit_range)) { /* Lower smoothing range */
|
|
cm_smooth_corner(*out,out_lower_limit,out_lower_limit,limit_range,
|
|
0.0,1.0,out,&pout_gain);
|
|
pout_pin = pout_pin * pout_gain;
|
|
}
|
|
else {
|
|
if (*out > (out_upper_limit + limit_range)) { /* At upper limit */
|
|
*out = out_upper_limit;
|
|
pout_pin = 0.0;
|
|
}
|
|
else {
|
|
if (*out > (out_upper_limit - limit_range)) { /* Upper smoothing region */
|
|
cm_smooth_corner(*out,out_upper_limit,out_upper_limit,limit_range,
|
|
1.0,0.0,out,&pout_gain);
|
|
pout_pin = pout_pin * pout_gain;
|
|
}
|
|
}
|
|
}
|
|
}
|
|
|
|
|
|
|
|
|
|
/** Output values for DC & Transient **/
|
|
|
|
OUTPUT(out) = *out;
|
|
PARTIAL(out,in) = pout_pin;
|
|
|
|
}
|
|
|
|
else { /**** AC Analysis...output (0.0,gain/s) ****/
|
|
ac_gain.real = 0.0;
|
|
ac_gain.imag = -gain / RAD_FREQ;
|
|
AC_GAIN(out,in) = ac_gain;
|
|
}
|
|
}
|
|
|
|
|
|
|
|
|
|
|