/* | Project:     WCDMA simulation environment | Module:       | Author:       | Date:        June 8, 1999 | | History: |  |              June 8,  1999 Maarit Melvasalo  | | * *  File    : s_crc_remove.c *  Abstract: *  Abstract: * * USER GIVEN PARAMETERS * * 1: N  = the number of input bits * 2: CRC = number of CRC bits * 3: CRC_POL = polynomial for CRC calulation * 4: nFrames = number of frames in a block * * Inputs: *        Input bits  * * Outputs:  *        Input bits + CRC bits          * | | Copyright disclaimer: |   This software was developed at the National Institute of Standards |   and Technology by employees of the Federal Government in the course |   of their official duties. Pursuant to title 17 Section 105 of the |   United States Code this software is not subject to copyright |   protection and is in the public domain. | |   We would appreciate acknowledgement if the software is used. | | */#define S_FUNCTION_NAME s_crc_remove#define S_FUNCTION_LEVEL 2#include <math.h>#include "simstruc.h"#include "tmwtypes.h"#include "config_wcdma.h"#include "wcdma_simulink.h"#include "crc_routines.h"/* USER GIVEN PARAMETERS AND DEFINITIONS /**//*Number of inputs and outputs/**/#define NINPUTS   2#define NOUTPUTS  2/* number of user given parameters /**/#define NPARAMS 4                                      /* Input block size*/#define IN_SIZE(S) ssGetSFcnParam(S,0)  #define nInputs (int_T)(mxGetPr(IN_SIZE(S))[0])    /* Number of CRC bits*/#define CRC_SIZE(S) ssGetSFcnParam(S,1)               #define nCRC (int_T)(mxGetPr(CRC_SIZE(S))[0])         /* CRC Polynomial/**/#define CRC_POL(S) ssGetSFcnParam(S,2)                /* Number of frames in one block/**/#define FRAME_PARAM(S) ssGetSFcnParam(S,3)            #define nFrames (int_T)(mxGetPr(FRAME_PARAM(S))[0]) /* Pointers to Input Ports */#define U(element) (*uPtrs[element])                  #define F(element) (*fPtrs[element])                  /* Sampletime -- defined in config_cdma /**/#define td  nFrames * TD_FRAME               /*====================* * S-function methods * *====================*//* Function: mdlInitializeSizes =============================================== * Abstract: *    The sizes information is used by Simulink to determine the S-function *    block's characteristics (number of inputs, outputs, states, etc.). */static void mdlInitializeSizes(SimStruct *S){  /* Number of expected parameters */  ssSetNumSFcnParams(S, NPARAMS);    if (ssGetNumSFcnParams(S) != ssGetSFcnParamsCount(S))  {      return; }  /* Initialize the input and output port sizes/**/    if (!ssSetNumInputPorts(S, NINPUTS)) return;    ssSetInputPortWidth(S, 0, nInputs);    ssSetInputPortWidth(S, 1, nSide);    ssSetInputPortDirectFeedThrough(S, 0, 1);    ssSetInputPortDirectFeedThrough(S, 1, 1);    if (!ssSetNumOutputPorts(S, NOUTPUTS)) return;    ssSetOutputPortWidth(S, 0, nInputs - nCRC);    ssSetOutputPortWidth(S, 1, 2);      /* Initialize number of sample times and simulink memory vectors/**/    ssSetNumSampleTimes(S, 1);    ssSetNumIWork(S, 1);     /* To Speeds up the simulations/**/    ssSetSFcnParamNotTunable(S,0);     ssSetSFcnParamNotTunable(S,1);      ssSetSFcnParamNotTunable(S,2);     ssSetSFcnParamNotTunable(S,3);     /* Take care when specifying exception free code - see sfuntmpl.doc */    ssSetOptions(S, SS_OPTION_EXCEPTION_FREE_CODE);}/* Function: mdlInitializeSampleTimes ========================================= * Abstract: *    Specifiy the sample time */static void mdlInitializeSampleTimes(SimStruct *S){     ssSetSampleTime(S, 0, td);  ssSetOffsetTime(S, 0, 0.0);}#define MDL_INITIALIZE_CONDITIONS/* Function: mdlInitializeConditions ======================================== * Abstract: *    Initialize the CRC table with given polynomial. *    Set the initial FER rate to zero. */static void mdlInitializeConditions(SimStruct *S){  int_T            crc_poly  = (int_T)mxGetPr(CRC_POL(S))[0];  int_T             *iwork = ssGetIWork(S);    wcdma_crctable_init (crc_poly);  iwork[0] = 0;  /* FER /**/}/* Function: mdlOutputs ======================================================= * Abstract: *         Call the CRC calculation function *         Return the bits with CRC bits *         Return cumulative FER */static void mdlOutputs(SimStruct *S, int_T tid){    real_T            *y    = ssGetOutputPortRealSignal(S,0);    real_T            *crc_out    = ssGetOutputPortRealSignal(S,1);    int_T             *iwork = ssGetIWork(S);    InputRealPtrsType uPtrs = ssGetInputPortRealSignalPtrs(S,0);    InputRealPtrsType fPtrs = ssGetInputPortRealSignalPtrs(S,1);    int_T             input[nInputs];    int_T             crc[nCRC];    uint16_T          crc_value;    int_T             i,tmp;    /* If input is valid data/**/    if (F(0) >0){      tmp = 0 ;      for (i = 0; i < nInputs - nCRC; i++){	input[i] = (int_T)U(i);	y[i] = U(i);      }      for (i = nInputs - nCRC; i < nInputs; i++)	input[i] = (int_T)U(i);    /* If number of CRC bits is > 0/**/      if (nCRC >0 ){	crc_value = (int_T)wcdma_get_crc_value(input, nInputs, crc);	if (((int_T)crc_value) >0) iwork[0]++;/**/	crc_out[0] = iwork[0];	crc_out[1] = crc_value;      }    }    else crc_out[1] = -1; /* ouput is not valid/**/}/* Function: mdlTerminate ===================================================== * Abstract: *    No termination needed, but we are required to have this routine. */static void mdlTerminate(SimStruct *S){}#ifdef  MATLAB_MEX_FILE    /* Is this file being compiled as a MEX-file? */#include "simulink.c"      /* MEX-file interface mechanism */#else#include "cg_sfun.h"       /* Code generation registration function */#endif