/* | | 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. |*//* | Project:     WCDMA simulation environment | Module:      C-level Test bed for downlink discrete channel routines | Author:      Tommi Makelainen, Nokia/NIST | Date:        February 17, 1999 | | History: |              February 17, 1999 Tommi Makelainen |                      Initial version. | |              May 12, 1999 Tommi Makelainen |                      Created separate tx-rx tester. |                      Tests transmission path from Tx modulation |                      to Rx demodulation (data modulation, spreading, |                      channel, despreading and demodulation). | */#include <stdio.h>#include <math.h>#include "chcoding.h"#include "chdecoding.h"#include "discmod.h"#include "discdemod.h"#include "channel.h"#include "rake.h"#include "blockcollect.h"#define INPUT_SIZE 10#define PROB_VECTOR_LEN 2int main(int argc, char *argv[]){    int I_out[3*256*INPUT_SIZE];    double I_out_2[3*256*INPUT_SIZE];    int I_out_3[3*256*INPUT_SIZE];    double I_out_4[3*256*INPUT_SIZE];    int Q_out[3*256*INPUT_SIZE];    double Q_out_2[3*256*INPUT_SIZE];    int Q_out_3[3*256*INPUT_SIZE];    double Q_out_4[3*256*INPUT_SIZE];    int I_symbs_len, Q_symbs_len;#if 0    int pn_code[] = {1, 1, -1, -1};    int pn_code_len = 4;    int sf = 4;#endif#if 1    int pn_code[] = {1, 1, 1, 1, -1, -1, -1, -1, 1, 1, 1, 1, -1, -1, -1, -1,                     1, 1, 1, 1, -1, -1, -1, -1, 1, 1, 1, 1, -1, -1, -1, -1,                     1, 1, 1, 1, -1, -1, -1, -1, 1, 1, 1, 1, -1, -1, -1, -1,                     1, 1, 1, 1, -1, -1, -1, -1, 1, 1, 1, 1, -1, -1, -1, -1,                     1, 1, 1, 1, -1, -1, -1, -1, 1, 1, 1, 1, -1, -1, -1, -1,                     1, 1, 1, 1, -1, -1, -1, -1, 1, 1, 1, 1, -1, -1, -1, -1,                     1, 1, 1, 1, -1, -1, -1, -1, 1, 1, 1, 1, -1, -1, -1, -1,                     1, 1, 1, 1, -1, -1, -1, -1, 1, 1, 1, 1, -1, -1, -1, -1,                     1, 1, 1, 1, -1, -1, -1, -1, 1, 1, 1, 1, -1, -1, -1, -1,                     1, 1, 1, 1, -1, -1, -1, -1, 1, 1, 1, 1, -1, -1, -1, -1,                     1, 1, 1, 1, -1, -1, -1, -1, 1, 1, 1, 1, -1, -1, -1, -1,                     1, 1, 1, 1, -1, -1, -1, -1, 1, 1, 1, 1, -1, -1, -1, -1,                     1, 1, 1, 1, -1, -1, -1, -1, 1, 1, 1, 1, -1, -1, -1, -1,                     1, 1, 1, 1, -1, -1, -1, -1, 1, 1, 1, 1, -1, -1, -1, -1,                     1, 1, 1, 1, -1, -1, -1, -1, 1, 1, 1, 1, -1, -1, -1, -1,                     1, 1, 1, 1, -1, -1, -1, -1, 1, 1, 1, 1, -1, -1, -1, -1};    int pn_code_len = 256;    int sf = 256;#endif#if 0    int pn_code[] = {1, 1, 1, 1, -1, -1, -1, -1, 1, 1, 1, 1, -1, -1, -1, -1,                     1, 1, 1, 1, -1, -1, -1, -1, 1, 1, 1, 1, -1, -1, -1, -1,                     1, 1, 1, 1, -1, -1, -1, -1, 1, 1, 1, 1, -1, -1, -1, -1,                     1, 1, 1, 1, -1, -1, -1, -1, 1, 1, 1, 1, -1, -1, -1, -1,                     1, 1, 1, 1, -1, -1, -1, -1, 1, 1, 1, 1, -1, -1, -1, -1,                     1, 1, 1, 1, -1, -1, -1, -1, 1, 1, 1, 1, -1, -1, -1, -1,                     1, 1, 1, 1, -1, -1, -1, -1, 1, 1, 1, 1, -1, -1, -1, -1,                     1, 1, 1, 1, -1, -1, -1, -1, 1, 1, 1, 1, -1, -1, -1, -1};    int pn_code_len = 128;    int sf = 128;#endif     int tx_coded_symbol_bits[3*INPUT_SIZE];    int rx_coded_symbol_bits[3*INPUT_SIZE];    int ch_delays[8] = {0, 2, 5,7, 8, 8, 12, 16};    int ch_delays2[8] = {5, 8, 16,3};    int *ch_de_ptr;    int nTaps = 3;    int nImpulse = 1;    int nFrames = 0;    double ch_amplitudes[] = {0.7, 0.1, 0.2, 0.2, 0.15, 0.1, 0.1, 0.1};    double ch_amplitudes2[] = {0.1, 0.1, 0.8, 0.2};    int channel_delays[16];    double esno = 16; /* in dB */    double power = 0; /* milliWatts */    double r_temppi[3*INPUT_SIZE];    int in_temppi[10000];    int ch_id0 = 0, ch_id1 = 1;    int coder_instance, decoder_instance;    int rake_instance1, rake_instance2;    int channel_instance1, channel_instance2;    int symbuf_instance;    int ready;    double out_data[20];    double path_select_threshold;    int pilot_len;    int round;    double impulse_prob[1] = { 1.0 };    double impulse_cumul[1];    int channel_no = 0;    int bit_delay=0;    int delayI, delayQ;    int nFingers;    int diff;    int i, j, k, m, tmp_value;    ch_id0 = wcdma_channel_init(ch_amplitudes, ch_delays, nTaps,                                nImpulse, impulse_prob, nFrames,                                esno, power, pn_code_len);    round = 0;    m = 0;    ch_de_ptr = ch_delays;    path_select_threshold = 0.18;    nFingers = 2;    pilot_len = 20;    rake_instance1= wcdma_rake_init(pilot_len, path_select_threshold,                                    nFingers);    rake_instance2= wcdma_rake_init(pilot_len, path_select_threshold,                                    nFingers);    symbuf_instance = wcdma_symbolbuffer_init(20);    for (k=0; k < 8; k++) {    printf("round %d\n", round);/*    if ( (round % 8) == 1 ) ch_de_ptr = ch_delays;    if (round == 2) ch_de_ptr = ch_delays2;*/    if ((round % 2) == 0) {        for (i=0; i < INPUT_SIZE; i++) {            tmp_value = i % 2;            tx_coded_symbol_bits[i] = tmp_value > 0 ? 0 : 1;        }    } else {        for (i=0; i < INPUT_SIZE; i++) {            tmp_value = i % 2;            tx_coded_symbol_bits[i] = tmp_value > 0 ? 1 : 0;        }    }    for (i=0; i < INPUT_SIZE; i++)         in_temppi[k*INPUT_SIZE+i] = tx_coded_symbol_bits[i];    round++;    wcdma_dl_mod(tx_coded_symbol_bits, INPUT_SIZE, pn_code, pn_code_len,                 sf, I_out+pilot_len, Q_out+pilot_len);    wcdma_channel(I_out, Q_out, sf*(INPUT_SIZE/2)+pilot_len, ch_id0,                  ch_amplitudes, channel_delays, I_out_2, Q_out_2);    delayI = delayQ = 0;    wcdma_rake_receiver(rake_instance1, I_out_2, sf*(INPUT_SIZE/2)+pilot_len,                         pn_code, pn_code_len, sf,                        ch_de_ptr, ch_amplitudes, nTaps,                        &I_symbs_len, &delayI, I_out_4);    wcdma_rake_receiver(rake_instance2, Q_out_2, sf*(INPUT_SIZE/2)+pilot_len,                        pn_code, pn_code_len, sf,                        ch_de_ptr, ch_amplitudes, nTaps,                        &Q_symbs_len, &delayQ, Q_out_4);    for (i=0, j=0; i < INPUT_SIZE/2; i++, j += 2) {        rx_coded_symbol_bits[j] = I_out_3[i] = (I_out_4[i] > 0) ? 1 : 0;         rx_coded_symbol_bits[j+1] = Q_out_3[i] = (Q_out_4[i] > 0) ? 1 : 0;     }    ready = 0;    for (i=0; i < I_symbs_len+Q_symbs_len; i++)          r_temppi[i] = (double)rx_coded_symbol_bits[i];    wcdma_symbolbuffer(symbuf_instance, r_temppi,             I_symbs_len+Q_symbs_len, &ready, out_data);    if (ready == 1) {        printf(" round %d output: ", k);        for (j=0; j < 20; j++) {            printf(" %d", (out_data[j] > 0) ? 1 : 0 );            if ( ((m*20+j) % INPUT_SIZE) == 15 ) printf(" x ");        }        printf("\n");        printf(" round %d input : ", k);        for (j=0; j < 20; j++) {            printf(" %d", in_temppi[m*20+j]);            if ( ((m*20+j) % INPUT_SIZE) == 15 ) printf(" x ");        }        m++;        printf("\n");    }#if 0    for (i=0; i < (INPUT_SIZE); i++) {      printf("input %d output %d\n",             tx_coded_symbol_bits[i], rx_coded_symbol_bits[i]);    }#endif    printf(" ---- \n");    } /* for k *//*    printf("\n\tChecking difference of bits before and after convolution.\n");    diff=0;    bit_delay = (delayI == 0) ? 0 : (delayI/pn_code_len + 1);    for (i=0; i < (INPUT_SIZE-bit_delay); i++) {      printf("input %d output %d\n",             tx_coded_symbol_bits[i], rx_coded_symbol_bits[i+bit_delay]);      diff += (tx_coded_symbol_bits[i] == rx_coded_symbol_bits[i+bit_delay])               ? 0 : 1;    }    printf("\t\tDifference is %d bit places.\n", diff);*/    printf("\n");    wcdma_channel_free(ch_id0);    wcdma_channel_free(ch_id1);    wcdma_symbolbuffer_free(0);    return(0);}