cfunc.mod

支持数字元件仿真的SPICE插件

/* $Id: cfunc.tpl,v 1.1 91/03/18 19:01:04 bill Exp $ */
/*.......1.........2.........3.........4.........5.........6.........7.........8
================================================================================

FILE oneshot/cfunc.mod

Copyright 1991
Georgia Tech Research Corporation, Atlanta, Ga. 30332
All Rights Reserved

PROJECT A-8503-405
               

AUTHORS                      

    20 Mar 1991     Harry Li


MODIFICATIONS   

    17 Sep 1991    Jeffrey P. Murray
     2 Oct 1991    Jeffrey P. Murray
                                   

SUMMARY

    This file contains the model-specific routines used to
    functionally describe the oneshot code model.


INTERFACES       

    FILE                 ROUTINE CALLED     

    CMmacros.h           cm_message_send();                   
                             
    CM.c                 void *cm_analog_alloc()
                         void *cm_analog_get_ptr()
                         int cm_analog_set_temp_bkpt()

REFERENCED FILES

    Inputs from and outputs to ARGS structure.
                     

NON-STANDARD FEATURES

    NONE

===============================================================================*/

/*=== INCLUDE FILES ====================*/

#include "oneshot.h"  

                                      

/*=== CONSTANTS ========================*/




/*=== MACROS ===========================*/



  
/*=== LOCAL VARIABLES & TYPEDEFS =======*/                         


    
           
/*=== FUNCTION PROTOTYPE DEFINITIONS ===*/




                   
/*==============================================================================

FUNCTION void cm_oneshot()

AUTHORS                      

    20 Mar 1991     Harry Li


MODIFICATIONS   

    17 Sep 1991    Jeffrey P. Murray
     2 Oct 1991    Jeffrey P. Murray

SUMMARY

    This function implements the oneshot code model.

INTERFACES       

    FILE                 ROUTINE CALLED     

    CMmacros.h           cm_message_send();                   
                             
    CM.c                 void *cm_analog_alloc()
                         void *cm_analog_get_ptr()
                         int cm_analog_set_temp_bkpt()

RETURNED VALUE
    
    Returns inputs and outputs via ARGS structure.

GLOBAL VARIABLES
    
    NONE

NON-STANDARD FEATURES

    NONE

==============================================================================*/

/*=== CM_ONESHOT ROUTINE ===*/

/***************************************************************************************
*
*  This model describes a totally analog oneshot.
*  After a rising edge is detected, the model will
*  output a pulse width specified by the controling 
*  voltage.
*                          HWL 20Mar91
*
*
*
*                              ___________________________________
*                             /<---pulse width --->              :\
*                            / :                  :              : \ 
*         <---rise_delay--> /  :                  :<-fall_delay->:  \
*      ___|________________/   :                  :              :   \____________
*         ^                <-->:                  :              :<-->
*         Trigger         Risetime                               Falltime
*
*
****************************************************************************************/

void cm_oneshot(ARGS)  /* structure holding parms, 
                                   inputs, outputs, etc.     */
{
    int i;             /* generic loop counter index                    */
	int *locked;        /* pointer used to store the locked1 variable */ 
	int locked1;        /* flag which allows the time points to be
						  reset.  value determined by retrig parameter  */
	int cntl_size;     /* size of the control array                     */
	int pw_size;       /* size of the pulse-width array                 */
	int *state;        /* pointer used to store state1 variable         */
	int state1;        /* if state1 = 1, then oneshot has
                          been triggered.  if state1 = 0, no change     */
	int *set;          /* pointer used to store the state of set1       */
	int set1;          /* flag used to set/reset the oneshot            */
    int trig_pos_edge; /* flag used to define positive or negative
                          edge triggering.  1=positive, 0=negative      */

    double *x;         /* pointer used to store the control array       */
	double *y;         /* pointer used to store the pulse-width array   */
	double cntl_input; /* the actual value of the control input         */
	double out;        /* value of the output                           */
	double dout_din;   /* slope of the pw wrt the control voltage       */
	double output_low; /* output low value                              */
	double output_hi;  /* output high value                             */
    double pw;         /* actual value of the pulse-width               */
/*	double del_out;     value of the delay time between triggering 
                          and a change in the output                    */
	double del_rise;    /* value of the delay time between triggering 
                          and a change in the output                    */
	double del_fall;    /* value of the delay time between the end of the
						   pw and a change in the output                */
	double *t1;        /* pointer used to store time1                   */
	double *t2;        /* pointer used to store time2                   */
	double *t3;        /* pointer used to store time3                   */
	double *t4;        /* pointer used to store time4                   */
	double time1;      /* time at which the output first begins to 
                          change (trigger + delay)                      */
	double time2;      /* time2 = time1 + risetime                      */
	double time3;      /* time3 = time2 + pw                            */
	double time4;      /* time4 = time3 + falltime                      */
	double t_rise;     /* risetime                                      */
	double t_fall;     /* falltime                                      */
	double *output_old;/* pointer which stores the previous output      */
    double *clock;     /* pointer which stores the clock                */
	double *old_clock; /* pointer which stores the previous clock       */
	double trig_clk;   /* value at which the clock triggers the oneshot */

    Mif_Complex_t ac_gain;
                                           

    /**** Retrieve frequently used parameters... ****/

    cntl_size = PARAM_SIZE(cntl_array);           
    pw_size = PARAM_SIZE(pw_array);           
    trig_clk = PARAM(clk_trig);
    trig_pos_edge = PARAM(pos_edge_trig);
    output_low = PARAM(out_low);
    output_hi = PARAM(out_high);
    /*del_out = PARAM(delay);*/
	del_rise = PARAM(rise_delay);
	del_fall = PARAM(fall_delay);
    t_rise = PARAM(rise_time);
    t_fall = PARAM(fall_time);

    /* set minimum rise and fall_times */

	if(t_rise < 1e-12){
		t_rise = 1e-12;
    }

	if(t_fall < 1e-12){
		t_fall = 1e-12;
    }

   /* the control array must be the same size as the pulse-width array */ 

	if(cntl_size != pw_size){
		cm_message_send(oneshot_array_error);
		return;
 	}

  if(INIT == 1){  /* first time through, allocate memory */

		t1 = cm_analog_alloc(T1,sizeof(double));
		t2 = cm_analog_alloc(T2,sizeof(double));
		t3 = cm_analog_alloc(T3,sizeof(double));
		t4 = cm_analog_alloc(T4,sizeof(double));
        set = cm_analog_alloc(SET,sizeof(int));
        state = cm_analog_alloc(STATE,sizeof(int));
        clock = cm_analog_alloc(CLOCK,sizeof(double));
		locked = cm_analog_alloc(LOCKED,sizeof(int));
        output_old = cm_analog_alloc(OUTPUT_OLD,sizeof(double));

  } 

  if(ANALYSIS == MIF_DC){

    /* for DC, initialize values and set the output = output_low */

		t1 = cm_analog_get_ptr(T1,0);
		t2 = cm_analog_get_ptr(T2,0);
		t3 = cm_analog_get_ptr(T3,0);
		t4 = cm_analog_get_ptr(T4,0);
		set = cm_analog_get_ptr(SET,0);
		state = cm_analog_get_ptr(STATE,0);
		locked = cm_analog_get_ptr(LOCKED,0);
		output_old = cm_analog_get_ptr(OUTPUT_OLD,0);

		/* initialize time and state values */
		*t1 = -1;
		*t2 = -1;
		*t3 = -1;
		*t4 = -1;
        *set = 0;
		*locked = 0;
        *state = 0;
        *output_old = output_low;

		OUTPUT(out) = output_low;
		if(PORT_NULL(cntl_in) != 1){
		   PARTIAL(out,cntl_in) = 0; 
        }
		if(PORT_NULL(clear) != 1){
		   PARTIAL(out,clear) = 0; 
        }
		PARTIAL(out,clk) = 0; 

  }else

  if(ANALYSIS == MIF_TRAN){

/* retrieve previous values, set them equal to the variables 
   Note that these pointer values are immediately dumped into
   other variables because the previous values can't change- 
   can't rewrite the old values */

	t1 = cm_analog_get_ptr(T1,1);
	t2 = cm_analog_get_ptr(T2,1);
	t3 = cm_analog_get_ptr(T3,1);
	t4 = cm_analog_get_ptr(T4,1);
	set = cm_analog_get_ptr(SET,1);
	state = cm_analog_get_ptr(STATE,1);
	locked = cm_analog_get_ptr(LOCKED,1);
    clock = cm_analog_get_ptr(CLOCK,0);
    old_clock = cm_analog_get_ptr(CLOCK,1);
    output_old = cm_analog_get_ptr(OUTPUT_OLD,1);

	time1 = *t1;
    time2 = *t2;
    time3 = *t3;
    time4 = *t4;
    set1 = *set;
    state1 = *state;
	locked1 = *locked;
   
if((PORT_NULL(clear) != 1) && (INPUT(clear) > trig_clk)){ 

	time1 = -1;
	time2 = -1;
	time3 = -1;
	time4 = -1;
    set1 = 0;
	locked1 = 0;
    state1 = 0;

    OUTPUT(out) = output_low;
	

}else{

    /* Allocate storage for breakpoint domain & freq. range values */

    x = (double *) calloc(cntl_size, sizeof(double));
    if (x == '\0') {
        cm_message_send(oneshot_allocation_error); 
        return;
    }

    y = (double *) calloc(pw_size, sizeof(double));
    if (y == '\0') {
        cm_message_send(oneshot_allocation_error);  
        return;
    }


    /* Retrieve control and pulse-width values. */       
    for (i=0; i<cntl_size; i++) {
        *(x+i) = PARAM(cntl_array[i]);
        *(y+i) = PARAM(pw_array[i]);
    }                       
    

    /* Retrieve cntl_input and clock value. */       

	if(PORT_NULL(cntl_in) != 1){

    	cntl_input = INPUT(cntl_in);

	}else{

		cntl_input = 0;

    }

    *clock = INPUT(clk);




    /* Determine segment boundaries within which cntl_input resides */

    if (cntl_input <= *x) { /* cntl_input below lowest cntl_voltage */
            dout_din = (*(y+1) - *y)/(*(x+1) - *x);
             pw = *y + (cntl_input - *x) * dout_din;

	     if(pw < 0){
	        cm_message_send(oneshot_pw_clamp);
	        pw = 0;
             }
    }
    else 
        /*** cntl_input above highest cntl_voltage ***/
	
        if (cntl_input >= *(x+cntl_size-1)){ 
            dout_din = (*(y+cntl_size-1) - *(y+cntl_size-2)) /
                          (*(x+cntl_size-1) - *(x+cntl_size-2));
            pw = *(y+cntl_size-1) + (cntl_input - *(x+cntl_size-1)) * dout_din;

        } else { /*** cntl_input within bounds of end midpoints...   
                must determine position progressively & then 
                calculate required output.                    ***/

            for (i=0; i<cntl_size-1; i++) {

                if ((cntl_input < *(x+i+1)) && (cntl_input >= *(x+i))){ 
            		
					/* Interpolate to get the correct pulse width value */

					pw = ((cntl_input - *(x+i))/(*(x+i+1) - *(x+i)))* 
							(*(y+i+1)-*(y+i)) + *(y+i); 
                } 
    
            }
        }


if(trig_pos_edge){  /* for a positive edge trigger */

    if(!set1){    /* if set1=0, then look for 
                     1.  a rising edge trigger
                     2.  the clock to be higher than the trigger value */

      if((*clock > *old_clock) && (*clock > trig_clk)){

        state1 = 1;
        set1 = 1;

      }

	}else

            /* look for a neg edge before resetting the trigger */

      if((*clock < *old_clock) && (*clock < trig_clk)){

        set1 = 0;
		
      }

}else{
         /* This stuff belongs to the case where a negative edge
			is needed */

    if(!set1){

      if((*clock < *old_clock) && (*clock < trig_clk)){

        state1 = 1;
        set1 = 1;

      }

	}else
            /* look for a pos edge before resetting the trigger */

    if((*clock > *old_clock) && (*clock > trig_clk)){

        set1 = 0;
    }
}


	/*  I can only set the breakpoints if the state1 is high and
		the output is low, and locked = 0 */
    if((state1) && (*output_old - output_low < 1e-20) && (!locked1)){

/* if state1 is 1, and the output is low, then set the time points
   and the temporary breakpoints */

   	time1 = TIME + del_rise; 
   	time2 = time1 + t_rise;
    time3 = time2 + pw + del_fall;
    time4 = time3 + t_fall;

    if(PARAM(retrig) == MIF_FALSE){

		locked1 = 1;

    }

   if((TIME < time1) || (T(1) == 0)){
    	cm_analog_set_temp_bkpt(time1);
   }

	   cm_analog_set_temp_bkpt(time2);
	   cm_analog_set_temp_bkpt(time3);
	   cm_analog_set_temp_bkpt(time4);

      /* reset the state value */
	   state1 = 0;
	   OUTPUT(out) = output_low;

   }else

	/* state1 = 1, and the output is high,  then just set time3 and time4 
       and their temporary breakpoints This implies that the oneshot was
	   retriggered */	

	 if((state1) && (*output_old - output_hi < 1e-20) && (!locked1)){

        time3 = TIME + pw + del_rise + del_fall + t_rise;
        time4 = time3 + t_fall;

        cm_analog_set_temp_bkpt(time3);
        cm_analog_set_temp_bkpt(time4);

        OUTPUT(out) = output_hi;

        state1 = 0;

    }

/* reset the state if it's 1 and the locked flag is 1.  This 
   means that the clock tried to retrigger the oneshot, but
   the retrig flag prevented it from doing so */

 if((state1) && (locked1)){

        state1 = 0;

 }
   /*  set the value for the output depending on the current time, and
       the values of time1, time2, time3, and time4 */
    if(TIME < time1){

		OUTPUT(out) = output_low;

    }else
	if((time1 <= TIME) && (TIME < time2)){

        OUTPUT(out) = output_low + ((TIME - time1)/(time2 - time1))*
			(output_hi - output_low); 

	}else
	if((time2 <= TIME) && (TIME < time3)){

        OUTPUT(out) = output_hi;

	}else
	if((time3 <= TIME) && (TIME < time4)){
		
	   OUTPUT(out) = output_hi + ((TIME - time3)/(time4 - time3))*
		(output_low - output_hi);

	}else{


  	  OUTPUT(out) = output_low;

	  /* oneshot can now be retriggered, set locked to 0 */
	   if(PARAM(retrig) == MIF_FALSE){

	   locked1 = 0;  
      
	   }
     }

    /* set the variables which need to be stored for the next iteration */
 }
	t1 = cm_analog_get_ptr(T1,0);
	t2 = cm_analog_get_ptr(T2,0);
	t3 = cm_analog_get_ptr(T3,0);
	t4 = cm_analog_get_ptr(T4,0);
	set = cm_analog_get_ptr(SET,0);
	locked = cm_analog_get_ptr(LOCKED,0);
	state = cm_analog_get_ptr(STATE,0);
	output_old = cm_analog_get_ptr(OUTPUT_OLD,0);

	*t1 = time1;
    *t2 = time2;
    *t3 = time3;
    *t4 = time4;
    *set = set1;
    *state = state1;
    *output_old = OUTPUT(out);
	*locked = locked1;

	if(PORT_NULL(cntl_in) != 1){
		PARTIAL(out,cntl_in) = 0;
    }
	if(PORT_NULL(clear) != 1){
		PARTIAL(out,clear) = 0;
    }
    PARTIAL(out,clk) = 0 ;

    } else {                      /* Output AC Gain */

  /* This model has no AC capability */

        ac_gain.real = 0.0; 
        ac_gain.imag= 0.0;
        AC_GAIN(out,clk) = ac_gain;
    }
}