/* | Project:     WCDMA simulation environment | Module:       | Author:       | Date:        MAy 24, 1999 | | History: |              May 23, 1999 Maarit Melvasalo |                      Interleaver separated from channel coding |                      see s_chcoding.c |              May 28, 1999 Maarit Melvasalo |                      sample time changed -- new parameter added  |              June 21, 1999 Maarit Melvasalo |                      Multistage interleaver added |  | |  File    : s_interleaving.c |  Abstract: | | USER GIVEN PARAMETERS | | 1: B   = Number of bits in a frame  | 2: COL = Number of colums for interleaver | 3: nFrames = Number of frames in a block | 4: inter_mode = Inter interleaver mode |   | Inputs: |        input bits  |                | Outputs:  |        Interleaved bits | | DEPENDS ON FILES  |        interleaving.c | | |  | 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_interleaver #define S_FUNCTION_LEVEL 2#include <math.h>#include "simstruc.h"#include "tmwtypes.h"#include "config_wcdma.h"#include "wcdma_simulink.h"/* USER GIVEN PARAMETERS AND DEFINITIONS /**//*Number of input and output ports/**/#define NINPUTS   1#define NOUTPUTS  1/* Number of user given parameters/**/  #define NPARAMS 4      /* Input and Output Size parameters /**/#define B_PARAM(S) ssGetSFcnParam(S,0)           #define bits_frame (int_T)(mxGetPr(B_PARAM(S))[0])   /* Number of columns for block interleaver/**/ #define COL_PARAM(S) ssGetSFcnParam(S,1)   /* Number of frames in one block/**/#define frame(S) ssGetSFcnParam(S,2)#define nFrames (int_T)(mxGetPr(frame(S))[0]) /* Inter frame interleaver mode flag */#define inter_mode(S) ssGetSFcnParam(S,3)  /* Pointer to Input Port */#define U(element) (*uPtrs[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, nFrames*bits_frame);    ssSetInputPortDirectFeedThrough(S, 0, 1);    if (!ssSetNumOutputPorts(S, NOUTPUTS)) return;    ssSetOutputPortWidth(S, 0, nFrames*bits_frame);     /* Initialize number of sample times and simulink work vectors /**/     ssSetNumSampleTimes(S, 1);    ssSetNumIWork(S, 3);     /* 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 sample tiem */static void mdlInitializeSampleTimes(SimStruct *S){  ssSetSampleTime(S, 0, td);    ssSetOffsetTime(S, 0, 0.0);}#define MDL_INITIALIZE_CONDITIONS/* Function: mdlInitializeConditions ======================================== * Abstract: *           For block interleaver calculate number of rows */static void mdlInitializeConditions(SimStruct *S){    int_T             *iwork = ssGetIWork(S);    int_T             cols  = (int_T)(mxGetPr(COL_PARAM(S))[0]);     real_T            mode = mxGetPr(inter_mode(S))[0];     int_T             rows;    /*      If the Block interleaver is selected      /**/    if ( mode <2 ){      iwork[0] = cols;        if(cols > 0) {	rows = (nFrames * bits_frame)/ cols;	if (rows * cols == (nFrames * bits_frame)){	  iwork[1] = rows;	}	/* if there there is something wrong with the cols /**/	else iwork[0] = 0;      }    }    }/* Function: mdlOutputs ======================================================= * Abstract: *           Performs for the given block (nFrames * bits_in_frame *           block interleaving or multistage interleaving *           Depending on the mode parameter * */static void mdlOutputs(SimStruct *S, int_T tid){    real_T            *y    = ssGetOutputPortRealSignal(S,0);    InputRealPtrsType uPtrs = ssGetInputPortRealSignalPtrs(S,0);    int_T             *iwork = ssGetIWork(S);    real_T            mode = mxGetPr(inter_mode(S))[0];     int_T             i;    int_T             inputs[nFrames*bits_frame];    int_T             outputs[nFrames*bits_frame];    /*        Move the input bits to a vector        /**/     for (i=0; i < nFrames*bits_frame; i++)       inputs[i] = (int_T)U(i);	    /*      If the Block interleaver is selected      /**/    if ( mode <2 ){      if(iwork[0] > 0) {		wcdma_block_interleaver(iwork[1], iwork[0], inputs, outputs);      }          }    /*      If the Multistage  interleaver is selected      /**/    else {      wcdma_interframe_interleaver(nFrames*bits_frame,bits_frame,inputs,outputs);     }            /*        Update the return value        /**/    for (i = 0; i < nFrames*bits_frame; i++)       y[i] = outputs[i];}/* 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