Home | History | Annotate | Download | only in lib
      1 /*
      2  * Copyright (C) 2008-2009 SVOX AG, Baslerstr. 30, 8048 Zuerich, Switzerland
      3  *
      4  * Licensed under the Apache License, Version 2.0 (the "License");
      5  * you may not use this file except in compliance with the License.
      6  * You may obtain a copy of the License at
      7  *
      8  *     http://www.apache.org/licenses/LICENSE-2.0
      9  *
     10  * Unless required by applicable law or agreed to in writing, software
     11  * distributed under the License is distributed on an "AS IS" BASIS,
     12  * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
     13  * See the License for the specific language governing permissions and
     14  * limitations under the License.
     15  */
     16 /**
     17  * @file picosig2.c
     18  *
     19  * Signal Generation PU - Internal functions - Implementation
     20  *
     21  * Copyright (C) 2008-2009 SVOX AG, Baslerstr. 30, 8048 Zuerich, Switzerland
     22  * All rights reserved.
     23  *
     24  * History:
     25  * - 2009-04-20 -- initial version
     26  *
     27  */
     28 #include "picoos.h"
     29 #include "picodsp.h"
     30 #include "picosig2.h"
     31 #include "picofftsg.h"
     32 
     33 #ifdef __cplusplus
     34 extern "C" {
     35 #endif
     36 #if 0
     37 }
     38 #endif
     39 /*---------------------------------------------------------------------------
     40  * INTERNAL FUNCTIONS DECLARATION
     41  *---------------------------------------------------------------------------*/
     42 static void gen_hann2(sig_innerobj_t *sig_inObj);
     43 static void get_simple_excitation(sig_innerobj_t *sig_inObj,
     44         picoos_int16 *nextPeak);
     45 static void enh_wind_init(sig_innerobj_t *sig_inObj);
     46 static void init_rand(sig_innerobj_t *sig_inObj);
     47 static void get_trig(picoos_int32 ang, picoos_int32 *table, picoos_int32 *cs,
     48         picoos_int32 *sn);
     49 
     50 /*---------------------------------------------------------------------------
     51  * PICO SYSTEM FUNCTIONS
     52  *---------------------------------------------------------------------------*/
     53 /**
     54  * allocation of DSP memory for SIG PU
     55  * @param   mm : memory manager
     56  * @param   sig_inObj : sig PU internal object of the sub-object
     57  * @return  PICO_OK : allocation successful
     58  * @return  PICO_ERR_OTHER : allocation NOT successful
     59  * @callgraph
     60  * @callergraph
     61  */
     62 pico_status_t sigAllocate(picoos_MemoryManager mm, sig_innerobj_t *sig_inObj)
     63 {
     64     picoos_int16 *data_i;
     65     picoos_int32 *d32;
     66     picoos_int32 nCount;
     67 
     68     sig_inObj->int_vec22 =
     69     sig_inObj->int_vec23 =
     70     sig_inObj->int_vec24 =
     71     sig_inObj->int_vec25 =
     72     sig_inObj->int_vec26 =
     73     sig_inObj->int_vec28 =
     74     sig_inObj->int_vec29 =
     75     sig_inObj->int_vec38 =
     76     sig_inObj->int_vec30 =
     77     sig_inObj->int_vec31 =
     78     sig_inObj->int_vec32 =
     79     sig_inObj->int_vec33 =
     80     sig_inObj->int_vec34 =
     81     sig_inObj->int_vec35 =
     82     sig_inObj->int_vec36 =
     83     sig_inObj->int_vec37 =
     84     sig_inObj->int_vec38 =
     85     sig_inObj->int_vec39 =
     86     sig_inObj->int_vec40 = NULL;
     87 
     88     sig_inObj->sig_vec1 = NULL;
     89 
     90     sig_inObj->idx_vect1 = sig_inObj->idx_vect2 = sig_inObj->idx_vect4 = NULL;
     91     sig_inObj->idx_vect5 = sig_inObj->idx_vect6 = sig_inObj->idx_vect7 =
     92     sig_inObj->idx_vect8 = sig_inObj->idx_vect9 = NULL;
     93     sig_inObj->ivalue17 = sig_inObj->ivalue18 = 0;
     94 
     95     /*-----------------------------------------------------------------
     96      * Memory allocations
     97      * NOTE : it would be far better to do a single allocation
     98      *          and to do pointer initialization inside this routine
     99      * ------------------------------------------------------------------*/
    100     data_i = (picoos_int16 *) picoos_allocate(mm, sizeof(picoos_int16)
    101             * PICODSP_FFTSIZE);
    102     if (NULL == data_i) {
    103         sigDeallocate(mm, sig_inObj);
    104         return PICO_ERR_OTHER;
    105     }
    106     sig_inObj->idx_vect1 = data_i;
    107 
    108     data_i = (picoos_int16 *) picoos_allocate(mm, sizeof(picoos_int16)
    109             * PICODSP_HFFTSIZE_P1);
    110     if (NULL == data_i) {
    111         sigDeallocate(mm, sig_inObj);
    112         return PICO_ERR_OTHER;
    113     }
    114     sig_inObj->idx_vect2 = data_i;
    115 
    116     data_i = (picoos_int16 *) picoos_allocate(mm, sizeof(picoos_int16)
    117             * PICODSP_FFTSIZE);
    118     if (NULL == data_i) {
    119         sigDeallocate(mm, sig_inObj);
    120         return PICO_ERR_OTHER;
    121     }
    122     sig_inObj->idx_vect4 = data_i;
    123 
    124     data_i = (picoos_int16 *) picoos_allocate(mm, sizeof(picoos_int16)
    125             * PICODSP_FFTSIZE);
    126     if (NULL == data_i) {
    127         sigDeallocate(mm, sig_inObj);
    128         return PICO_ERR_OTHER;
    129     }
    130     sig_inObj->idx_vect5 = data_i;
    131 
    132     data_i = (picoos_int16 *) picoos_allocate(mm, sizeof(picoos_int16)
    133             * PICODSP_FFTSIZE);
    134     if (NULL == data_i) {
    135         sigDeallocate(mm, sig_inObj);
    136         return PICO_ERR_OTHER;
    137     }
    138     sig_inObj->idx_vect6 = data_i;
    139 
    140     data_i = (picoos_int16 *) picoos_allocate(mm, sizeof(picoos_int16)
    141             * PICODSP_HFFTSIZE_P1);
    142     if (NULL == data_i) {
    143         sigDeallocate(mm, sig_inObj);
    144         return PICO_ERR_OTHER;
    145     }
    146     sig_inObj->idx_vect7 = data_i;
    147 
    148     data_i = (picoos_int16 *) picoos_allocate(mm, sizeof(picoos_int16)
    149             * PICODSP_MAX_EX);
    150     if (NULL == data_i) {
    151         sigDeallocate(mm, sig_inObj);
    152         return PICO_ERR_OTHER;
    153     }
    154     sig_inObj->idx_vect8 = data_i;
    155 
    156     data_i = (picoos_int16 *) picoos_allocate(mm, sizeof(picoos_int16)
    157             * PICODSP_MAX_EX);
    158     if (data_i == NULL) {
    159         sigDeallocate(mm, sig_inObj);
    160         return PICO_ERR_OTHER;
    161     }
    162     sig_inObj->idx_vect9 = data_i;
    163 
    164     d32 = (picoos_int32 *) picoos_allocate(mm, sizeof(picoos_int32)
    165             * PICODSP_FFTSIZE);
    166     if (NULL == d32) {
    167         sigDeallocate(mm, sig_inObj);
    168         return PICO_ERR_OTHER;
    169     }
    170     sig_inObj->int_vec22 = d32;
    171     d32 = (picoos_int32 *) picoos_allocate(mm, sizeof(picoos_int32)
    172             * PICODSP_FFTSIZE);
    173     if (NULL == d32) {
    174         sigDeallocate(mm, sig_inObj);
    175         return PICO_ERR_OTHER;
    176     }
    177     sig_inObj->int_vec23 = d32;
    178     d32 = (picoos_int32 *) picoos_allocate(mm, sizeof(picoos_int32)
    179             * PICODSP_FFTSIZE);
    180     if (NULL == d32) {
    181         sigDeallocate(mm, sig_inObj);
    182         return PICO_ERR_OTHER;
    183     }
    184     sig_inObj->int_vec24 = d32;
    185     d32 = (picoos_int32 *) picoos_allocate(mm, sizeof(picoos_int32)
    186             * PICODSP_FFTSIZE);
    187     if (NULL == d32) {
    188         sigDeallocate(mm, sig_inObj);
    189         return PICO_ERR_OTHER;
    190     }
    191     sig_inObj->int_vec25 = d32;
    192     d32 = (picoos_int32 *) picoos_allocate(mm, sizeof(picoos_int32)
    193             * PICODSP_FFTSIZE * 2);
    194     if (NULL == d32) {
    195         sigDeallocate(mm, sig_inObj);
    196         return PICO_ERR_OTHER;
    197     }
    198     sig_inObj->int_vec26 = d32;
    199 
    200     d32 = (picoos_int32 *) picoos_allocate(mm, sizeof(picoos_int32)
    201             * PICODSP_FFTSIZE);
    202     if (NULL == d32) {
    203         sigDeallocate(mm, sig_inObj);
    204         return PICO_ERR_OTHER;
    205     }
    206     sig_inObj->int_vec28 = d32;
    207     d32 = (picoos_int32 *) picoos_allocate(mm, sizeof(picoos_int32)
    208             * PICODSP_FFTSIZE);
    209     if (NULL == d32) {
    210         sigDeallocate(mm, sig_inObj);
    211         return PICO_ERR_OTHER;
    212     }
    213     sig_inObj->int_vec29 = d32;
    214     d32 = (picoos_int32 *) picoos_allocate(mm, sizeof(picoos_int32)
    215             * PICODSP_FFTSIZE);
    216     if (NULL == d32) {
    217         sigDeallocate(mm, sig_inObj);
    218         return PICO_ERR_OTHER;
    219     }
    220     sig_inObj->int_vec38 = d32;
    221     d32 = (picoos_int32 *) picoos_allocate(mm, sizeof(picoos_int32)
    222             * PICODSP_FFTSIZE);
    223     if (NULL == d32) {
    224         sigDeallocate(mm, sig_inObj);
    225         return PICO_ERR_OTHER;
    226     }
    227     sig_inObj->int_vec30 = d32;
    228     d32 = (picoos_int32 *) picoos_allocate(mm, sizeof(picoos_int32)
    229             * PICODSP_FFTSIZE);
    230     if (NULL == d32) {
    231         sigDeallocate(mm, sig_inObj);
    232         return PICO_ERR_OTHER;
    233     }
    234     sig_inObj->int_vec31 = d32;
    235 
    236     d32 = (picoos_int32 *) picoos_allocate(mm, sizeof(picoos_int32)
    237             * PICODSP_FFTSIZE);
    238     if (NULL == d32) {
    239         sigDeallocate(mm, sig_inObj);
    240         return PICO_ERR_OTHER;
    241     }
    242     sig_inObj->int_vec32 = d32;
    243     d32 = (picoos_int32 *) picoos_allocate(mm, sizeof(picoos_int32)
    244             * PICODSP_FFTSIZE);
    245     if (NULL == d32) {
    246         sigDeallocate(mm, sig_inObj);
    247         return PICO_ERR_OTHER;
    248     }
    249     sig_inObj->int_vec33 = d32;
    250     d32 = (picoos_int32 *) picoos_allocate(mm, sizeof(picoos_int32)
    251             * PICODSP_N_RAND_TABLE);
    252     if (NULL == d32) {
    253         sigDeallocate(mm, sig_inObj);
    254         return PICO_ERR_OTHER;
    255     }
    256     sig_inObj->int_vec34 = d32;
    257     d32 = (picoos_int32 *) picoos_allocate(mm, sizeof(picoos_int32)
    258             * PICODSP_N_RAND_TABLE);
    259     if (NULL == d32) {
    260         sigDeallocate(mm, sig_inObj);
    261         return PICO_ERR_OTHER;
    262     }
    263     sig_inObj->int_vec35 = d32;
    264     d32 = (picoos_int32 *) picoos_allocate(mm, sizeof(picoos_int32)
    265             * PICODSP_N_RAND_TABLE);
    266     if (NULL == d32) {
    267         sigDeallocate(mm, sig_inObj);
    268         return PICO_ERR_OTHER;
    269     }
    270     sig_inObj->int_vec36 = d32;
    271     d32 = (picoos_int32 *) picoos_allocate(mm, sizeof(picoos_int32)
    272             * PICODSP_N_RAND_TABLE);
    273     if (NULL == d32) {
    274         sigDeallocate(mm, sig_inObj);
    275         return PICO_ERR_OTHER;
    276     }
    277     sig_inObj->int_vec37 = d32;
    278 
    279     d32 = (picoos_int32 *) picoos_allocate(mm, sizeof(picoos_int32)
    280             * PICODSP_HFFTSIZE_P1);
    281     if (NULL == d32) {
    282         sigDeallocate(mm, sig_inObj);
    283         return PICO_ERR_OTHER;
    284     }
    285     sig_inObj->int_vec39 = d32;
    286     d32 = (picoos_int32 *) picoos_allocate(mm, sizeof(picoos_int32) * (1
    287             + PICODSP_COS_TABLE_LEN));
    288     if (NULL == d32) {
    289         sigDeallocate(mm, sig_inObj);
    290         return PICO_ERR_OTHER;
    291     }
    292     sig_inObj->int_vec40 = d32;
    293 
    294     for (nCount = 0; nCount < CEPST_BUFF_SIZE; nCount++) {
    295         d32 = (picoos_int32 *) picoos_allocate(mm, sizeof(picoos_int32) * (PICODSP_CEPORDER));
    296         if (NULL == d32) {
    297             sigDeallocate(mm, sig_inObj);
    298             return PICO_ERR_OTHER;
    299         }
    300         sig_inObj->int_vec41[nCount] = d32;
    301     }
    302 
    303     for (nCount = 0; nCount < PHASE_BUFF_SIZE; nCount++) {
    304         d32 = (picoos_int32 *) picoos_allocate(mm, sizeof(picoos_int32) * (PICODSP_PHASEORDER));
    305         if (NULL == d32) {
    306             sigDeallocate(mm, sig_inObj);
    307             return PICO_ERR_OTHER;
    308         }
    309         sig_inObj->int_vec42[nCount] = d32;
    310     }
    311 
    312     d32 = (picoos_int32 *) picoos_allocate(mm, sizeof(picoos_int32)
    313             * PICODSP_FFTSIZE * 2); /* - fixed point */
    314     if (NULL == d32) {
    315         sigDeallocate(mm, sig_inObj);
    316         return PICO_ERR_OTHER;
    317     }
    318     sig_inObj->sig_vec1 = d32;
    319 
    320     return PICO_OK;
    321 }/*sigAllocate*/
    322 
    323 /**
    324  * frees DSP memory for SIG PU
    325  * @param   mm : memory manager
    326  * @param   sig_inObj : sig PU internal object of the sub-object
    327  * @return  void
    328  * @callgraph
    329  * @callergraph
    330  */
    331 void sigDeallocate(picoos_MemoryManager mm, sig_innerobj_t *sig_inObj)
    332 {
    333     picoos_int32 nCount;
    334     /*-----------------------------------------------------------------
    335      * Memory de-allocations
    336      * ------------------------------------------------------------------*/
    337     if (NULL != sig_inObj->idx_vect1)
    338         picoos_deallocate(mm, (void *) &(sig_inObj->idx_vect1));
    339     if (NULL != sig_inObj->idx_vect2)
    340         picoos_deallocate(mm, (void *) &(sig_inObj->idx_vect2));
    341     if (NULL != sig_inObj->idx_vect4)
    342         picoos_deallocate(mm, (void *) &(sig_inObj->idx_vect4));
    343     if (NULL != sig_inObj->idx_vect5)
    344         picoos_deallocate(mm, (void *) &(sig_inObj->idx_vect5));
    345     if (NULL != sig_inObj->idx_vect6)
    346         picoos_deallocate(mm, (void *) &(sig_inObj->idx_vect6));
    347     if (NULL != sig_inObj->idx_vect7)
    348         picoos_deallocate(mm, (void *) &(sig_inObj->idx_vect7));
    349     if (NULL != sig_inObj->idx_vect8)
    350         picoos_deallocate(mm, (void *) &(sig_inObj->idx_vect8));
    351     if (NULL != sig_inObj->idx_vect9)
    352         picoos_deallocate(mm, (void *) &(sig_inObj->idx_vect9));
    353 
    354     if (NULL != sig_inObj->int_vec22)
    355         picoos_deallocate(mm, (void *) &(sig_inObj->int_vec22));
    356     if (NULL != sig_inObj->int_vec23)
    357         picoos_deallocate(mm, (void *) &(sig_inObj->int_vec23));
    358     if (NULL != sig_inObj->int_vec24)
    359         picoos_deallocate(mm, (void *) &(sig_inObj->int_vec24));
    360     if (NULL != sig_inObj->int_vec25)
    361         picoos_deallocate(mm, (void *) &(sig_inObj->int_vec25));
    362     if (NULL != sig_inObj->int_vec26)
    363         picoos_deallocate(mm, (void *) &(sig_inObj->int_vec26));
    364     if (NULL != sig_inObj->int_vec28)
    365         picoos_deallocate(mm, (void *) &(sig_inObj->int_vec28));
    366     if (NULL != sig_inObj->int_vec29)
    367         picoos_deallocate(mm, (void *) &(sig_inObj->int_vec29));
    368     if (NULL != sig_inObj->int_vec38)
    369         picoos_deallocate(mm, (void *) &(sig_inObj->int_vec38));
    370     if (NULL != sig_inObj->int_vec30)
    371         picoos_deallocate(mm, (void *) &(sig_inObj->int_vec30));
    372     if (NULL != sig_inObj->int_vec31)
    373         picoos_deallocate(mm, (void *) &(sig_inObj->int_vec31));
    374     if (NULL != sig_inObj->int_vec32)
    375         picoos_deallocate(mm, (void *) &(sig_inObj->int_vec32));
    376     if (NULL != sig_inObj->int_vec33)
    377         picoos_deallocate(mm, (void *) &(sig_inObj->int_vec33));
    378     if (NULL != sig_inObj->int_vec34)
    379         picoos_deallocate(mm, (void *) &(sig_inObj->int_vec34));
    380     if (NULL != sig_inObj->int_vec35)
    381         picoos_deallocate(mm, (void *) &(sig_inObj->int_vec35));
    382     if (NULL != sig_inObj->int_vec36)
    383         picoos_deallocate(mm, (void *) &(sig_inObj->int_vec36));
    384     if (NULL != sig_inObj->int_vec37)
    385         picoos_deallocate(mm, (void *) &(sig_inObj->int_vec37));
    386     if (NULL != sig_inObj->int_vec39)
    387         picoos_deallocate(mm, (void *) &(sig_inObj->int_vec39));
    388     if (NULL != sig_inObj->int_vec40)
    389         picoos_deallocate(mm, (void *) &(sig_inObj->int_vec40));
    390 
    391     for (nCount = 0; nCount < CEPST_BUFF_SIZE; nCount++) {
    392         if (NULL != sig_inObj->int_vec41[nCount]) {
    393             picoos_deallocate(mm, (void *) &(sig_inObj->int_vec41[nCount]));
    394         }
    395     }
    396 
    397     for (nCount = 0; nCount < PHASE_BUFF_SIZE; nCount++) {
    398         if (NULL != sig_inObj->int_vec42[nCount]) {
    399             picoos_deallocate(mm, (void *) &(sig_inObj->int_vec42[nCount]));
    400         }
    401     }
    402 
    403     if (NULL != sig_inObj->sig_vec1) {
    404         picoos_deallocate(mm, (void *) &(sig_inObj->sig_vec1));
    405     }
    406 }/*sigDeAllocate*/
    407 
    408 /**
    409  * initializes all memory neededed by DSP at instance creation time
    410  * @param   sig_inObj : sig PU internal object of the sub-object
    411  * @return  void
    412  * @callgraph
    413  * @callergraph
    414  */
    415 void sigDspInitialize(sig_innerobj_t *sig_inObj, picoos_int32 resetMode)
    416 {
    417     picoos_int32 i, j;
    418     picoos_int32 *pnt;
    419 
    420     if (resetMode == PICO_RESET_SOFT) {
    421         /*minimal initialization when receiving a soft reset */
    422         return;
    423     }
    424     /*-----------------------------------------------------------------
    425      * Initialization
    426      * ------------------------------------------------------------------*/
    427     sig_inObj->warp_p = PICODSP_FREQ_WARP_FACT;
    428     sig_inObj->VCutoff_p = PICODSP_V_CUTOFF_FREQ; /*voicing cut off frequency in Hz (will be modeled in the future)*/
    429     sig_inObj->UVCutoff_p = PICODSP_UV_CUTOFF_FREQ;/*unvoiced frames only (periodize lowest components to mask bad voicing transitions)*/
    430     sig_inObj->Fs_p = PICODSP_SAMP_FREQ; /*Sampling freq*/
    431 
    432     sig_inObj->m1_p = PICODSP_CEPORDER;
    433     sig_inObj->m2_p = PICODSP_FFTSIZE; /*also initializes windowLen*/
    434     sig_inObj->framesz_p = PICODSP_DISPLACE; /*1/4th of the frame size = displacement*/
    435     sig_inObj->hfftsize_p = PICODSP_H_FFTSIZE; /*half of the FFT size*/
    436     sig_inObj->voxbnd_p = (picoos_int32) ((picoos_single) sig_inObj->hfftsize_p
    437             / ((picoos_single) sig_inObj->Fs_p / (picoos_single) 2)
    438             * (picoos_single) sig_inObj->VCutoff_p);
    439     sig_inObj->voxbnd2_p
    440             = (picoos_int32) ((picoos_single) sig_inObj->hfftsize_p
    441                     / ((picoos_single) sig_inObj->Fs_p / (picoos_single) 2)
    442                     * (picoos_single) sig_inObj->UVCutoff_p);
    443     sig_inObj->hop_p = sig_inObj->framesz_p;
    444     sig_inObj->nextPeak_p = (((int) (PICODSP_FFTSIZE))
    445             / ((int) PICODSP_DISPLACE) - 1) * sig_inObj->hop_p;
    446     sig_inObj->phId_p = 0; /*phonetic id*/
    447     sig_inObj->E_p = (picoos_single) 0.0f;
    448     sig_inObj->F0_p = (picoos_single) 0.0f;
    449     sig_inObj->voiced_p = 0;
    450     sig_inObj->nV = sig_inObj->nU = 0;
    451     sig_inObj->sMod_p = (picoos_single) 1.0f;
    452 
    453     /*cleanup vectors*/
    454     for (i = 0; i < 2 * PICODSP_FFTSIZE; i++) {
    455         sig_inObj->sig_vec1[i] = 0;
    456         sig_inObj->int_vec26[i] = 0; /*wav buff cleanup */
    457     }
    458 
    459     for (i = 0; i < PICODSP_FFTSIZE; i++) {
    460         sig_inObj->idx_vect1[i] = sig_inObj->idx_vect4[i]
    461                 = sig_inObj->idx_vect5[i] = sig_inObj->idx_vect6[i] = 0;
    462         sig_inObj->int_vec32[i] = sig_inObj->int_vec33[i] = 0;
    463     }
    464 
    465     for (i = 0; i < PICODSP_HFFTSIZE_P1; i++) {
    466         sig_inObj->idx_vect2[i] = (picoos_int16) 0;
    467     }
    468 
    469     for (i = 0; i < CEPST_BUFF_SIZE; i++) {
    470         sig_inObj->F0Buff[i]=0;
    471         sig_inObj->PhIdBuff[i]=0;
    472         sig_inObj->VoicingBuff[i]=0;
    473         sig_inObj->FuVBuff[i]=0;
    474         if (NULL != sig_inObj->CepBuff[i]) {
    475             pnt = sig_inObj->CepBuff[i];
    476             for (j = 0; j < PICODSP_CEPORDER; j++) {
    477                 pnt[j] = 0;
    478             }
    479         }
    480     }
    481 
    482     for (i = 0; i < PHASE_BUFF_SIZE; i++) {
    483         if (NULL != sig_inObj->int_vec42[i]) {
    484             pnt = sig_inObj->int_vec42[i];
    485             for (j = 0; j < PICODSP_PHASEORDER; j++) {
    486                 pnt[j] = 0;
    487             }
    488         }
    489     }
    490     sig_inObj->n_available=0;
    491     /*---------------------------------------------
    492      Init    formant enhancement window
    493      hanning window,
    494      Post Filter Hermite's interpolator Matrix
    495      Mel-2-Lin lookup tables
    496      ---------------------------------------------*/
    497     enh_wind_init(sig_inObj); /*creates the formant enhancement window*/
    498     init_rand(sig_inObj);
    499     gen_hann2(sig_inObj);
    500     mel_2_lin_init(sig_inObj);
    501 
    502 }/*sigDspInitialize*/
    503 
    504 /*-------------------------------------------------------------------------------
    505  PROCESSING FUNCTIONS : CALLED WITHIN sigStep (cfr. picosig.c)
    506  --------------------------------------------------------------------------------*/
    507 /**
    508  * convert from mel scale to linear scale
    509  * @param   sig_inObj : sig PU internal object of the sub-object
    510  * @param   scmeanMGC : mean value of the MGC
    511  * @return  void
    512  * @callgraph
    513  * @callergraph
    514  * @remarks translated from matlab code to c-code
    515  * Input
    516  * - c1 : input mfcc vector (ceporder=m1, real)
    517  * - m1 : input order
    518  * - A,B,D : lookup tables
    519  * - m2 : output order
    520  * - Xr,Xi (m2=FFT size, real) temporary arrays for FFT
    521  * - WNr,WNi (m2=FFT size, real) cos and sin precalculated tables
    522  * Output
    523  * - Xr (m2=FFT size, real) linear cepstral vector
    524  */
    525 void mel_2_lin_lookup(sig_innerobj_t *sig_inObj, picoos_uint32 scmeanMGC)
    526 {
    527     /*Local vars*/
    528     picoos_int16 nI, k;
    529     picoos_int32 delta, term1, term2;
    530 
    531     /*Local vars to be linked with sig data object*/
    532     picoos_int32 *c1, *XXr;
    533     picoos_single K1;
    534     picoos_int32 *D, K2, shift;
    535     picoos_int16 m1, *A, m2, m4, voiced, i;
    536 
    537     /*Link local variables with sig data object*/
    538     c1 = sig_inObj->wcep_pI;
    539     m1 = sig_inObj->m1_p;
    540     m2 = PICODSP_FFTSIZE;
    541     m4 = m2 >> 1;
    542 
    543     A = sig_inObj->A_p;
    544     D = sig_inObj->d_p;
    545 
    546     XXr = sig_inObj->wcep_pI;
    547     voiced = sig_inObj->voiced_p;
    548 
    549     shift = 27 - scmeanMGC;
    550     K2 = 1 << shift;
    551     K1 = (picoos_single) PICODSP_START_FLOAT_NORM * K2;
    552     XXr[0] = (picoos_int32) ((picoos_single) c1[0] * K1);
    553     for (nI = 1; nI < m1; nI++) {
    554         XXr[nI] = c1[nI] << shift;
    555     }
    556     i = sizeof(picoos_int32) * (PICODSP_FFTSIZE + 1 - m1);
    557     picoos_mem_set(XXr + m1, 0, i);
    558     dfct_nmf(m4, XXr); /* DFCT directly in fixed point */
    559 
    560     /* *****************************************************************************************
    561      Linear frequency scale envelope through interpolation.
    562      Two additions and one multiplication per entry.
    563 
    564      Optimization of linear interpolation algorithm
    565      - Start from 1 and stop at PICODSP_H_FFTSIZE-1 because 0 and PICODSP_H_FFTSIZE are invariant points
    566      - B[k]=A[k]+1 except for 0 and PICODSP_H_FFTSIZE
    567      - get rid of extra -1 operation by adapting the table A[]
    568 
    569      *******************************************************************************************/
    570     for (nI = 1; nI < PICODSP_H_FFTSIZE; nI++) {
    571         k = A[nI];
    572         term2 = XXr[k];
    573         term1 = XXr[k + 1];
    574         delta = term1 - term2;
    575         XXr[nI] = term2 + ((D[nI] * delta) >> 5); /* ok because nI<=A[nI] <=B[nI] */
    576     }
    577 }/*mel_2_lin_lookup*/
    578 
    579 /**
    580  * calculate phase
    581  * @remarks voiced phase taken from phase codebook and smoothed,
    582  * @remarks unvoiced phase - random
    583  * @param   sig_inObj : sig PU internal object of the sub-object
    584  * @return  void
    585  * @callgraph
    586  * @callergraph
    587  */
    588 void phase_spec2(sig_innerobj_t *sig_inObj)
    589 {
    590     picoos_int16 nI, iRand, firstUV;
    591     picoos_int32 *tmp1, *tmp2;
    592     picoos_int16 VOXBND_M1;
    593     picoos_int32 *spect, *ang;
    594     picoos_int16 voiced, m2;
    595     picoos_int32 *co, *so, *c, *s, voxbnd, voxbnd2;
    596     picoos_int16 i,j, k, n_comp;
    597     picoos_int16 *Pvoxbnd;
    598     picoos_int32 *phs_p2, *phs_p1, *phs_n1, *phs_n2;
    599     picoos_int32 *phs;
    600 
    601     /*Link local variables to sig data object*/
    602     spect = sig_inObj->wcep_pI; /* spect_p;*/
    603     /* current spect scale : times PICODSP_FIX_SCALE1 */
    604     ang = sig_inObj->ang_p;
    605     voxbnd = (picoos_int32) (sig_inObj->voxbnd_p * sig_inObj->voicing);
    606     voxbnd2 = sig_inObj->voxbnd2_p;
    607     voiced = sig_inObj->voiced_p;
    608     m2 = sig_inObj->m2_p;
    609     VOXBND_M1 = voxbnd - 1;
    610     firstUV = 1;
    611 
    612     /*code starts here*/
    613     if (voiced == 1) {
    614         firstUV = voxbnd;
    615         Pvoxbnd =  sig_inObj->VoxBndBuff;
    616         n_comp   = Pvoxbnd[2];
    617         phs_p2 = sig_inObj->PhsBuff[0];
    618         phs_p1 = sig_inObj->PhsBuff[1];
    619         phs    = sig_inObj->PhsBuff[2];
    620         phs_n1 = sig_inObj->PhsBuff[3];
    621         phs_n2 = sig_inObj->PhsBuff[4];
    622 
    623         /* find and smooth components which have full context */
    624         j = n_comp;
    625         for (i=0; i<5; i++) {
    626             if (Pvoxbnd[i]<j) j = Pvoxbnd[i];
    627         }
    628         for (i=0; i<j; i++) {
    629             ang[i] = -(((phs_p2[i]+phs_p1[i]+phs[i]+phs_n1[i]+phs_n2[i])<<6) / 5);
    630         }
    631 
    632         /* find and smooth components which at least one component on each side */
    633         k = n_comp;
    634         if (Pvoxbnd[2]<k) k = Pvoxbnd[2];
    635         if (Pvoxbnd[4]<k) k = Pvoxbnd[4];
    636         for (i=j; i<k; i++) {  /* smooth using only two surrounding neighbours */
    637                 ang[i] = -(((phs_p1[i]+phs[i]+phs_n1[i])<<6) / 3);
    638         }
    639 
    640         /* handle rest of components - at least one side does not exist */
    641         for (i=k; i<n_comp; i++) {
    642             ang[i] = -(phs[i]<<6); /* - simple copy without smoothing */
    643         }
    644 
    645         /*Phase unwrap - cumsum */
    646         tmp1 = &(ang[1]);
    647         tmp2 = &(ang[0]);
    648         /* current ang scale : PICODSP_M_PI = PICODSP_FIX_SCALE2 */
    649         FAST_DEVICE(VOXBND_M1,*(tmp1++)+=*(tmp2)-PICODSP_FIX_SCALE2;*(tmp2)=(*tmp2>=0)?(*tmp2)>>PICODSP_SHIFT_FACT4:-((-(*tmp2))>>PICODSP_SHIFT_FACT4);tmp2++);
    650         *tmp2 = (*tmp2 >= 0) ? (*tmp2) >> PICODSP_SHIFT_FACT4 : -((-(*tmp2))
    651                 >> PICODSP_SHIFT_FACT4); /*ang[voxbnd-1]/=2;*/
    652     }
    653 
    654     /* now for the unvoiced part */
    655     iRand = sig_inObj->iRand;
    656     c = sig_inObj->randCosTbl + iRand;
    657     s = sig_inObj->randSinTbl + iRand;
    658     co = sig_inObj->outCosTbl + firstUV;
    659     so = sig_inObj->outSinTbl + firstUV;
    660     for (nI = firstUV; nI < PICODSP_HFFTSIZE_P1 - 1; nI++) {
    661         *co++ = *c++;
    662         *so++ = *s++;
    663     }
    664     *co = 1;
    665     *so = 0;
    666     sig_inObj->iRand += (PICODSP_HFFTSIZE_P1 - firstUV);
    667     if (sig_inObj->iRand > PICODSP_N_RAND_TABLE - PICODSP_HFFTSIZE_P1)
    668         sig_inObj->iRand = 1 + sig_inObj->iRand + PICODSP_HFFTSIZE_P1
    669             - PICODSP_N_RAND_TABLE;
    670 }/*phase_spec2*/
    671 
    672 /**
    673  * Prepare Envelope spectrum for inverse FFT
    674  * @param   sig_inObj : sig PU internal object of the sub-object
    675  * @return  void
    676  * @remarks make phase bilateral -->> angh (FFT size, real)
    677  * @remarks combine in complex input vector for IFFT F = e**(spet/2+j*ang)
    678  * @remarks Compute energy -->> E (scalar, real)
    679  * @callgraph
    680  * @callergraph
    681  * Input
    682  * - spect (FFT size, real)
    683  * - ang (half FFT size -1, real)
    684  * - m2    fftsize
    685  *  - WNr,WNi (FFT size, real) the tabulated sine and cosine values
    686  *  - brev (FFT size, real) the tabulated bit reversal indexes
    687  * Output
    688  * - Fr, Fi (FFT size, complex) the envelope spectrum
    689  * - E (scalar, real) the energy
    690  */
    691 void env_spec(sig_innerobj_t *sig_inObj)
    692 {
    693 
    694     picoos_int16 nI;
    695     picoos_int32 fcX, fsX, fExp, voxbnd;
    696     picoos_int32 *spect, *ang, *ctbl;
    697     picoos_int16 voiced, prev_voiced;
    698     picoos_int32 *co, *so;
    699     picoos_int32 *Fr, *Fi;
    700     picoos_single mult;
    701 
    702     /*Link local variables to sig object*/
    703     spect = sig_inObj->wcep_pI; /*spect_p*/
    704     /*  current spect scale : times PICODSP_FIX_SCALE1 */
    705     ang = sig_inObj->ang_p;
    706     /*  current spect scale : PICODSP_M_PI =  PICODSP_FIX_SCALE2 */
    707     Fr = sig_inObj->F2r_p;
    708     Fi = sig_inObj->F2i_p;
    709     voiced = sig_inObj->voiced_p;
    710     prev_voiced = sig_inObj->prevVoiced_p;
    711     voxbnd = (picoos_int32) (sig_inObj->voxbnd_p * sig_inObj->voicing);
    712     ctbl = sig_inObj->cos_table;
    713     /*  ctbl scale : times 4096 */
    714     mult = PICODSP_ENVSPEC_K1 / PICODSP_FIX_SCALE1;
    715 
    716     /*remove dc from real part*/
    717     if (sig_inObj->F0_p > 120) {
    718         spect[0] = spect[1] = 0;
    719         spect[2] /= PICODSP_ENVSPEC_K2;
    720     } else {
    721         spect[0] = 0;
    722     }
    723 
    724     /* if using rand table, use sin and cos tables as well */
    725     if (voiced || (prev_voiced)) {
    726         /*Envelope becomes a complex exponential : F=exp(.5*spect + j*angh);*/
    727         for (nI = 0; nI < voxbnd; nI++) {
    728             get_trig(ang[nI], ctbl, &fcX, &fsX);
    729             fExp = (picoos_int32) EXP((double)spect[nI]*mult);
    730             Fr[nI] = fExp * fcX;
    731             Fi[nI] = fExp * fsX;
    732         }
    733         /*         ao=sig_inObj->ang_p+(picoos_int32)voxbnd; */
    734         co = sig_inObj->outCosTbl + voxbnd;
    735         so = sig_inObj->outSinTbl + voxbnd;
    736 
    737         for (nI = voxbnd; nI < PICODSP_HFFTSIZE_P1; nI++) {
    738             fcX = *co++;
    739             fsX = *so++;
    740             fExp = (picoos_int32) EXP((double)spect[nI]*mult);
    741             Fr[nI] = fExp * fcX;
    742             Fi[nI] = fExp * fsX;
    743         }
    744     } else {
    745         /*ao=sig_inObj->ang_p+1;*/
    746         co = sig_inObj->outCosTbl + 1;
    747         so = sig_inObj->outSinTbl + 1;
    748         for (nI = 1; nI < PICODSP_HFFTSIZE_P1; nI++) {
    749             fcX = *co++;
    750             fsX = *so++;
    751             fExp = (picoos_int32) EXP((double)spect[nI]*mult);
    752 
    753             Fr[nI] = fExp * fcX;
    754             Fi[nI] = fExp * fsX;
    755         }
    756     }
    757 
    758 }/*env_spec*/
    759 
    760 /**
    761  * Calculates the impulse response of the comlpex spectrum through inverse rFFT
    762  * @param   sig_inObj : sig PU internal object of the sub-object
    763  * @return  void
    764  * @remarks Imp corresponds with the real part of the FFT
    765  * @callgraph
    766  * @callergraph
    767  * Input
    768  * - Fr, Fi (FFT size, real & imaginary) the complex envelope spectrum (only first half of spectrum)
    769  * Output
    770  * - Imp: impulse response (length: m2)
    771  * - E (scalar, real) RMS value
    772  */
    773 void impulse_response(sig_innerobj_t *sig_inObj)
    774 {
    775     /*Define local variables*/
    776     picoos_single f;
    777     picoos_int16 nI, nn, m2, m4, voiced;
    778     picoos_single *E;
    779     picoos_int32 *norm_window; /* - fixed point */
    780     picoos_int32 *fr, *Fr, *Fi, *t1, ff; /* - fixed point */
    781 
    782     /*Link local variables with sig object*/
    783     m2 = sig_inObj->m2_p;
    784     m4 = m2 >> 1;
    785     Fr = sig_inObj->F2r_p;
    786     Fi = sig_inObj->F2i_p;
    787     norm_window = sig_inObj->norm_window_p;
    788     E = &(sig_inObj->E_p); /*as pointer: value will be modified*/
    789     voiced = sig_inObj->voiced_p;
    790     fr = sig_inObj->imp_p;
    791 
    792     /*Inverse FFT*/
    793     for (nI = 0, nn = 0; nI < m4; nI++, nn += 2) {
    794         fr[nn] = Fr[nI]; /* - fixed point */
    795     }
    796 
    797     fr[1] = (picoos_int32) (Fr[m4]);
    798     for (nI = 1, nn = 3; nI < m4; nI++, nn += 2) {
    799         fr[nn] = -Fi[nI]; /* - fixed point */
    800     }
    801 
    802     rdft(m2, -1, fr);
    803     /*window, normalize and differentiate*/
    804     *E = norm_result(m2, fr, norm_window);
    805 
    806     if (*E > 0) {
    807         f = *E * PICODSP_FIXRESP_NORM;
    808     } else {
    809         f = 20; /*PICODSP_FIXRESP_NORM*/
    810     }
    811     ff = (picoos_int32) f;
    812     if (ff < 1)
    813         ff = 1;
    814     /*normalize impulse response*/
    815     t1 = fr;FAST_DEVICE(PICODSP_FFTSIZE,*(t1++) /= ff;); /* - fixed point */
    816 
    817 } /* impulse_response */
    818 
    819 /**
    820  * time domain pitch synchronous overlap add over two frames (when no voicing transition)
    821  * @param    sig_inObj : sig PU internal object of the sub-object
    822  * @return  void
    823  * @remarks Special treatment at voicing boundaries
    824  * @remarks Introduced to get rid of time-domain aliasing (and additional speed up)
    825  * @callgraph
    826  * @callergraph
    827  */
    828 void td_psola2(sig_innerobj_t *sig_inObj)
    829 {
    830     picoos_int16 nI;
    831     picoos_int16 hop, m2, *nextPeak, voiced;
    832     picoos_int32 *t1, *t2;
    833     picoos_int16 cnt;
    834     picoos_int32 *fr, *v1, ff, f;
    835     picoos_int16 a, i;
    836     picoos_int32 *window;
    837     picoos_int16 s = (picoos_int16) 1;
    838     window = sig_inObj->window_p;
    839 
    840     /*Link local variables with sig object*/
    841     hop = sig_inObj->hop_p;
    842     m2 = sig_inObj->m2_p;
    843     nextPeak = &(sig_inObj->nextPeak_p);
    844     voiced = sig_inObj->voiced_p;
    845     fr = sig_inObj->imp_p;
    846     /*toggle the pointers and initialize signal vector */
    847     v1 = sig_inObj->sig_vec1;
    848 
    849     t1 = v1;
    850     FAST_DEVICE(PICODSP_FFTSIZE-PICODSP_DISPLACE,*(t1++)=0;);
    851     t1 = &(v1[PICODSP_FFTSIZE - PICODSP_DISPLACE]);
    852     t2 = &(v1[PICODSP_FFTSIZE]);
    853     FAST_DEVICE(PICODSP_FFTSIZE, *(t1++)=*(t2++););
    854     t1 = &(v1[2 * PICODSP_FFTSIZE - PICODSP_DISPLACE]);FAST_DEVICE(PICODSP_DISPLACE,*(t1++)=0;);
    855     /*calculate excitation points*/
    856     get_simple_excitation(sig_inObj, nextPeak);
    857 
    858     /*TD-PSOLA based on excitation vector */
    859     if ((sig_inObj->nU == 0) && (sig_inObj->voiced_p == 1)) {
    860         /* purely voiced */
    861         for (nI = 0; nI < sig_inObj->nV; nI++) {
    862             f = sig_inObj->EnV[nI];
    863             a = 0;
    864             cnt = PICODSP_FFTSIZE;
    865             ff = (f * window[sig_inObj->LocV[nI]]) >> PICODSP_SHIFT_FACT1;
    866             t1 = &(v1[a + sig_inObj->LocV[nI]]);
    867             t2 = &(fr[a]);
    868             if (cnt > 0)FAST_DEVICE(cnt,*(t1++)+=*(t2++)*ff;);
    869         }
    870     } else if ((sig_inObj->nV == 0) && (sig_inObj->voiced_p == 0)) {
    871         /* PURELY UNVOICED*/
    872         for (nI = 0; nI < sig_inObj->nU; nI++) {
    873             f = sig_inObj->EnU[nI];
    874             s = -s; /*reverse order to reduce the periodicity effect*/
    875             if (s == 1) {
    876                 a = 0;
    877                 cnt = PICODSP_FFTSIZE;
    878                 ff = (f * window[sig_inObj->LocU[nI]]) >> PICODSP_SHIFT_FACT1;
    879                 t1 = &(v1[a + sig_inObj->LocU[nI]]);
    880                 t2 = &(fr[a]);
    881                 if (cnt > 0)FAST_DEVICE(cnt,*(t1++)+=*(t2++)*ff; );
    882             } else { /*s==-1*/
    883                 a = 0;
    884                 cnt = PICODSP_FFTSIZE;
    885                 ff = (f * window[sig_inObj->LocU[nI]]) >> PICODSP_SHIFT_FACT1;
    886                 t1 = &(v1[(m2 - 1 - a) + sig_inObj->LocU[nI]]);
    887                 t2 = &(fr[a]);
    888                 if (cnt > 0)FAST_DEVICE(cnt,*(t1--)+=*(t2++)*ff; );
    889             }
    890         }
    891     } else if (sig_inObj->VoicTrans == 0) {
    892         /*voicing transition from unvoiced to voiced*/
    893         for (nI = 0; nI < sig_inObj->nV; nI++) {
    894             f = sig_inObj->EnV[nI];
    895             a = 0;
    896             cnt = PICODSP_FFTSIZE;
    897             ff = (f * window[sig_inObj->LocV[nI]]) >> PICODSP_SHIFT_FACT1;
    898             t1 = &(v1[a + sig_inObj->LocV[nI]]);
    899             t2 = &(fr[a]);
    900             if (cnt > 0)FAST_DEVICE(cnt,*(t1++)+=*(t2++)*ff;);
    901         }
    902         /*add remaining stuff from unvoiced part*/
    903         for (nI = 0; nI < sig_inObj->nU; nI++) {
    904             f = sig_inObj->EnU[nI];
    905             s = -s; /*reverse order to reduce the periodicity effect*/
    906             if (s == 1) {
    907                 a = 0;
    908                 cnt = PICODSP_FFTSIZE;
    909                 ff = (f * window[sig_inObj->LocU[nI]]) >> PICODSP_SHIFT_FACT1;
    910                 t1 = &(v1[a + sig_inObj->LocU[nI]]);
    911                 t2 = &(sig_inObj->ImpResp_p[a]); /*saved impulse response*/
    912                 if (cnt > 0)FAST_DEVICE(cnt,*(t1++)+=*(t2++)*ff; );
    913             } else {
    914                 a = 0;
    915                 cnt = PICODSP_FFTSIZE;
    916                 ff = (f * window[sig_inObj->LocU[nI]]) >> PICODSP_SHIFT_FACT1;
    917                 t1 = &(v1[(m2 - 1 - a) + sig_inObj->LocU[nI]]);
    918                 t2 = &(sig_inObj->ImpResp_p[a]);
    919                 if (cnt > 0)FAST_DEVICE(cnt,*(t1--)+=*(t2++)*ff; );
    920             }
    921         }
    922     } else {
    923         /*voiced to unvoiced*/
    924         for (nI = 0; nI < sig_inObj->nU; nI++) {
    925             f = sig_inObj->EnU[nI];
    926             s = -s; /*reverse order to reduce the periodicity effect*/
    927             if (s > 0) {
    928                 a = 0;
    929                 cnt = PICODSP_FFTSIZE;
    930                 ff = (f * window[sig_inObj->LocU[nI]]) >> PICODSP_SHIFT_FACT1;
    931                 t1 = &(v1[a + sig_inObj->LocU[nI]]);
    932                 t2 = &(fr[a]);
    933                 if (cnt > 0)FAST_DEVICE(cnt,*(t1++)+=*(t2++)*ff; );
    934             } else {
    935                 a = 0;
    936                 cnt = PICODSP_FFTSIZE;
    937                 ff = (f * window[sig_inObj->LocU[nI]]) >> PICODSP_SHIFT_FACT1;
    938                 t1 = &(v1[(m2 - 1 - a) + sig_inObj->LocU[nI]]);
    939                 t2 = &(fr[a]);
    940                 if (cnt > 0)FAST_DEVICE(cnt,*(t1--)+=*(t2++)*ff; );
    941             }
    942         }
    943         /*add remaining stuff from voiced part*/
    944         for (nI = 0; nI < sig_inObj->nV; nI++) {
    945             f = sig_inObj->EnV[nI];
    946             a = 0;
    947             cnt = PICODSP_FFTSIZE;
    948             ff = (f * window[sig_inObj->LocV[nI]]) >> PICODSP_SHIFT_FACT1;
    949             t1 = &(v1[a + sig_inObj->LocV[nI]]);
    950             t2 = &(sig_inObj->ImpResp_p[a]);
    951             if (cnt > 0)FAST_DEVICE(cnt,*(t1++)+=*(t2++)*ff;);
    952         }
    953     }
    954 
    955     t1 = sig_inObj->sig_vec1;
    956     for (i = 0; i < PICODSP_FFTSIZE; i++, t1++) {
    957         if (*t1 >= 0)
    958             *t1 >>= PICODSP_SHIFT_FACT5;
    959         else
    960             *t1 = -((-*t1) >> PICODSP_SHIFT_FACT5);
    961     }
    962 
    963 }/*td_psola2*/
    964 
    965 /**
    966  * overlap + add summing of impulse responses on the final destination sample buffer
    967  * @param    sig_inObj : sig PU internal object of the sub-object
    968  * @return  void
    969  * @remarks Special treatment at voicing boundaries
    970  * @remarks Introduced to get rid of time-domain aliasing (and additional speed up)
    971  * Input
    972  * - wlet : the generic impulse response (FFT size, real)
    973  * - window : the windowing funcion (FFT size, real) fixed
    974  * - WavBuff : the destination buffer with past samples (FFT size*2, short)
    975  * - m2 : fftsize
    976  * Output
    977  * - WavBuff : the destination buffer with updated samples (FFT size*2, short)
    978  * @callgraph
    979  * @callergraph
    980  */
    981 void overlap_add(sig_innerobj_t *sig_inObj)
    982 {
    983     /*Local variables*/
    984     picoos_int32 *w, *v;
    985 
    986     /*Link local variables with sig object*/
    987     w = sig_inObj->WavBuff_p;
    988     v = sig_inObj->sig_vec1;
    989 
    990     FAST_DEVICE(PICODSP_FFTSIZE, *(w++)+=*(v++)<<PICODSP_SHIFT_FACT6;);
    991 
    992 }/*overlap_add*/
    993 
    994 /*-------------------------------------------------------------------------------
    995  INITIALIZATION AND INTERNAL    FUNCTIONS
    996  --------------------------------------------------------------------------------*/
    997 /**
    998  * Hanning window initialization
    999  * @param    sig_inObj : sig PU internal object of the sub-object
   1000  * @return  PICO_OK
   1001  * @callgraph
   1002  * @callergraph
   1003  */
   1004 static void gen_hann2(sig_innerobj_t *sig_inObj)
   1005 {
   1006     picoos_int32 *hann;
   1007     picoos_int32 *norm;
   1008     /*link local variables with sig object*/
   1009     hann = sig_inObj->window_p;
   1010     norm = sig_inObj->norm_window_p;
   1011 
   1012     norm[0] = 80224;
   1013     norm[1] = 320832;
   1014     norm[2] = 721696;
   1015     norm[3] = 1282560;
   1016     norm[4] = 2003104;
   1017     norm[5] = 2882880;
   1018     norm[6] = 3921376;
   1019     norm[7] = 5117984;
   1020     norm[8] = 6471952;
   1021     norm[9] = 7982496;
   1022     norm[10] = 9648720;
   1023     norm[11] = 11469616;
   1024     norm[12] = 13444080;
   1025     norm[13] = 15570960;
   1026     norm[14] = 17848976;
   1027     norm[15] = 20276752;
   1028     norm[16] = 22852864;
   1029     norm[17] = 25575744;
   1030     norm[18] = 28443776;
   1031     norm[19] = 31455264;
   1032     norm[20] = 34608368;
   1033     norm[21] = 37901248;
   1034     norm[22] = 41331904;
   1035     norm[23] = 44898304;
   1036     norm[24] = 48598304;
   1037     norm[25] = 52429696;
   1038     norm[26] = 56390192;
   1039     norm[27] = 60477408;
   1040     norm[28] = 64688944;
   1041     norm[29] = 69022240;
   1042     norm[30] = 73474720;
   1043     norm[31] = 78043744;
   1044     norm[32] = 82726544;
   1045     norm[33] = 87520352;
   1046     norm[34] = 92422272;
   1047     norm[35] = 97429408;
   1048     norm[36] = 102538752;
   1049     norm[37] = 107747248;
   1050     norm[38] = 113051776;
   1051     norm[39] = 118449184;
   1052     norm[40] = 123936224;
   1053     norm[41] = 129509648;
   1054     norm[42] = 135166080;
   1055     norm[43] = 140902192;
   1056     norm[44] = 146714528;
   1057     norm[45] = 152599584;
   1058     norm[46] = 158553904;
   1059     norm[47] = 164573888;
   1060     norm[48] = 170655936;
   1061     norm[49] = 176796448;
   1062     norm[50] = 182991712;
   1063     norm[51] = 189238064;
   1064     norm[52] = 195531744;
   1065     norm[53] = 201868992;
   1066     norm[54] = 208246016;
   1067     norm[55] = 214659040;
   1068     norm[56] = 221104176;
   1069     norm[57] = 227577616;
   1070     norm[58] = 234075488;
   1071     norm[59] = 240593872;
   1072     norm[60] = 247128912;
   1073     norm[61] = 253676688;
   1074     norm[62] = 260233280;
   1075     norm[63] = 266794768;
   1076     norm[64] = 273357248;
   1077     norm[65] = 279916768;
   1078     norm[66] = 286469440;
   1079     norm[67] = 293011360;
   1080     norm[68] = 299538560;
   1081     norm[69] = 306047168;
   1082     norm[70] = 312533312;
   1083     norm[71] = 318993088;
   1084     norm[72] = 325422656;
   1085     norm[73] = 331818144;
   1086     norm[74] = 338175744;
   1087     norm[75] = 344491680;
   1088     norm[76] = 350762176;
   1089     norm[77] = 356983424;
   1090     norm[78] = 363151808;
   1091     norm[79] = 369263520;
   1092     norm[80] = 375315008;
   1093     norm[81] = 381302592;
   1094     norm[82] = 387222720;
   1095     norm[83] = 393071872;
   1096     norm[84] = 398846528;
   1097     norm[85] = 404543232;
   1098     norm[86] = 410158560;
   1099     norm[87] = 415689216;
   1100     norm[88] = 421131840;
   1101     norm[89] = 426483200;
   1102     norm[90] = 431740096;
   1103     norm[91] = 436899392;
   1104     norm[92] = 441958016;
   1105     norm[93] = 446912928;
   1106     norm[94] = 451761152;
   1107     norm[95] = 456499840;
   1108     norm[96] = 461126080;
   1109     norm[97] = 465637152;
   1110     norm[98] = 470030400;
   1111     norm[99] = 474303104;
   1112     norm[100] = 478452800;
   1113     norm[101] = 482476960;
   1114     norm[102] = 486373184;
   1115     norm[103] = 490139200;
   1116     norm[104] = 493772640;
   1117     norm[105] = 497271424;
   1118     norm[106] = 500633440;
   1119     norm[107] = 503856704;
   1120     norm[108] = 506939200;
   1121     norm[109] = 509879168;
   1122     norm[110] = 512674880;
   1123     norm[111] = 515324544;
   1124     norm[112] = 517826688;
   1125     norm[113] = 520179776;
   1126     norm[114] = 522382368;
   1127     norm[115] = 524433184;
   1128     norm[116] = 526331008;
   1129     norm[117] = 528074688;
   1130     norm[118] = 529663200;
   1131     norm[119] = 531095552;
   1132     norm[120] = 532370944;
   1133     norm[121] = 533488576;
   1134     norm[122] = 534447808;
   1135     norm[123] = 535248000;
   1136     norm[124] = 535888768;
   1137     norm[125] = 536369664;
   1138     norm[126] = 536690432;
   1139     norm[127] = 536850880;
   1140     norm[128] = 536850880;
   1141     norm[129] = 536690432;
   1142     norm[130] = 536369664;
   1143     norm[131] = 535888768;
   1144     norm[132] = 535248000;
   1145     norm[133] = 534447808;
   1146     norm[134] = 533488576;
   1147     norm[135] = 532370944;
   1148     norm[136] = 531095552;
   1149     norm[137] = 529663200;
   1150     norm[138] = 528074688;
   1151     norm[139] = 526331008;
   1152     norm[140] = 524433216;
   1153     norm[141] = 522382368;
   1154     norm[142] = 520179776;
   1155     norm[143] = 517826688;
   1156     norm[144] = 515324544;
   1157     norm[145] = 512674880;
   1158     norm[146] = 509879168;
   1159     norm[147] = 506939200;
   1160     norm[148] = 503856704;
   1161     norm[149] = 500633472;
   1162     norm[150] = 497271424;
   1163     norm[151] = 493772672;
   1164     norm[152] = 490139200;
   1165     norm[153] = 486373184;
   1166     norm[154] = 482476992;
   1167     norm[155] = 478452800;
   1168     norm[156] = 474303104;
   1169     norm[157] = 470030400;
   1170     norm[158] = 465637184;
   1171     norm[159] = 461126080;
   1172     norm[160] = 456499840;
   1173     norm[161] = 451761152;
   1174     norm[162] = 446912960;
   1175     norm[163] = 441958016;
   1176     norm[164] = 436899424;
   1177     norm[165] = 431740096;
   1178     norm[166] = 426483200;
   1179     norm[167] = 421131840;
   1180     norm[168] = 415689216;
   1181     norm[169] = 410158560;
   1182     norm[170] = 404543232;
   1183     norm[171] = 398846528;
   1184     norm[172] = 393071872;
   1185     norm[173] = 387222720;
   1186     norm[174] = 381302592;
   1187     norm[175] = 375315008;
   1188     norm[176] = 369263552;
   1189     norm[177] = 363151808;
   1190     norm[178] = 356983456;
   1191     norm[179] = 350762176;
   1192     norm[180] = 344491712;
   1193     norm[181] = 338175776;
   1194     norm[182] = 331818144;
   1195     norm[183] = 325422656;
   1196     norm[184] = 318993088;
   1197     norm[185] = 312533312;
   1198     norm[186] = 306047168;
   1199     norm[187] = 299538560;
   1200     norm[188] = 293011360;
   1201     norm[189] = 286469472;
   1202     norm[190] = 279916800;
   1203     norm[191] = 273357248;
   1204     norm[192] = 266794784;
   1205     norm[193] = 260233280;
   1206     norm[194] = 253676688;
   1207     norm[195] = 247128928;
   1208     norm[196] = 240593888;
   1209     norm[197] = 234075488;
   1210     norm[198] = 227577632;
   1211     norm[199] = 221104192;
   1212     norm[200] = 214659040;
   1213     norm[201] = 208246032;
   1214     norm[202] = 201868992;
   1215     norm[203] = 195531744;
   1216     norm[204] = 189238080;
   1217     norm[205] = 182991728;
   1218     norm[206] = 176796448;
   1219     norm[207] = 170655952;
   1220     norm[208] = 164573888;
   1221     norm[209] = 158553920;
   1222     norm[210] = 152599600;
   1223     norm[211] = 146714528;
   1224     norm[212] = 140902208;
   1225     norm[213] = 135166096;
   1226     norm[214] = 129509648;
   1227     norm[215] = 123936240;
   1228     norm[216] = 118449184;
   1229     norm[217] = 113051776;
   1230     norm[218] = 107747248;
   1231     norm[219] = 102538752;
   1232     norm[220] = 97429424;
   1233     norm[221] = 92422288;
   1234     norm[222] = 87520352;
   1235     norm[223] = 82726544;
   1236     norm[224] = 78043744;
   1237     norm[225] = 73474736;
   1238     norm[226] = 69022240;
   1239     norm[227] = 64688944;
   1240     norm[228] = 60477424;
   1241     norm[229] = 56390192;
   1242     norm[230] = 52429696;
   1243     norm[231] = 48598304;
   1244     norm[232] = 44898304;
   1245     norm[233] = 41331904;
   1246     norm[234] = 37901248;
   1247     norm[235] = 34608384;
   1248     norm[236] = 31455264;
   1249     norm[237] = 28443792;
   1250     norm[238] = 25575744;
   1251     norm[239] = 22852864;
   1252     norm[240] = 20276752;
   1253     norm[241] = 17848976;
   1254     norm[242] = 15570960;
   1255     norm[243] = 13444080;
   1256     norm[244] = 11469616;
   1257     norm[245] = 9648720;
   1258     norm[246] = 7982512;
   1259     norm[247] = 6471952;
   1260     norm[248] = 5117984;
   1261     norm[249] = 3921376;
   1262     norm[250] = 2882880;
   1263     norm[251] = 2003104;
   1264     norm[252] = 1282560;
   1265     norm[253] = 721696;
   1266     norm[254] = 320832;
   1267     norm[255] = 80224;
   1268     hann[0] = 0;
   1269     hann[1] = 0;
   1270     hann[2] = 1;
   1271     hann[3] = 2;
   1272     hann[4] = 3;
   1273     hann[5] = 5;
   1274     hann[6] = 7;
   1275     hann[7] = 9;
   1276     hann[8] = 12;
   1277     hann[9] = 15;
   1278     hann[10] = 18;
   1279     hann[11] = 21;
   1280     hann[12] = 25;
   1281     hann[13] = 29;
   1282     hann[14] = 34;
   1283     hann[15] = 38;
   1284     hann[16] = 43;
   1285     hann[17] = 48;
   1286     hann[18] = 54;
   1287     hann[19] = 59;
   1288     hann[20] = 66;
   1289     hann[21] = 72;
   1290     hann[22] = 78;
   1291     hann[23] = 85;
   1292     hann[24] = 92;
   1293     hann[25] = 100;
   1294     hann[26] = 107;
   1295     hann[27] = 115;
   1296     hann[28] = 123;
   1297     hann[29] = 131;
   1298     hann[30] = 140;
   1299     hann[31] = 148;
   1300     hann[32] = 157;
   1301     hann[33] = 166;
   1302     hann[34] = 176;
   1303     hann[35] = 185;
   1304     hann[36] = 195;
   1305     hann[37] = 205;
   1306     hann[38] = 215;
   1307     hann[39] = 225;
   1308     hann[40] = 236;
   1309     hann[41] = 247;
   1310     hann[42] = 257;
   1311     hann[43] = 268;
   1312     hann[44] = 279;
   1313     hann[45] = 291;
   1314     hann[46] = 302;
   1315     hann[47] = 313;
   1316     hann[48] = 325;
   1317     hann[49] = 337;
   1318     hann[50] = 349;
   1319     hann[51] = 360;
   1320     hann[52] = 372;
   1321     hann[53] = 385;
   1322     hann[54] = 397;
   1323     hann[55] = 409;
   1324     hann[56] = 421;
   1325     hann[57] = 434;
   1326     hann[58] = 446;
   1327     hann[59] = 458;
   1328     hann[60] = 471;
   1329     hann[61] = 483;
   1330     hann[62] = 496;
   1331     hann[63] = 508;
   1332     hann[64] = 521;
   1333     hann[65] = 533;
   1334     hann[66] = 546;
   1335     hann[67] = 558;
   1336     hann[68] = 571;
   1337     hann[69] = 583;
   1338     hann[70] = 596;
   1339     hann[71] = 608;
   1340     hann[72] = 620;
   1341     hann[73] = 632;
   1342     hann[74] = 645;
   1343     hann[75] = 657;
   1344     hann[76] = 669;
   1345     hann[77] = 680;
   1346     hann[78] = 692;
   1347     hann[79] = 704;
   1348     hann[80] = 715;
   1349     hann[81] = 727;
   1350     hann[82] = 738;
   1351     hann[83] = 749;
   1352     hann[84] = 760;
   1353     hann[85] = 771;
   1354     hann[86] = 782;
   1355     hann[87] = 792;
   1356     hann[88] = 803;
   1357     hann[89] = 813;
   1358     hann[90] = 823;
   1359     hann[91] = 833;
   1360     hann[92] = 842;
   1361     hann[93] = 852;
   1362     hann[94] = 861;
   1363     hann[95] = 870;
   1364     hann[96] = 879;
   1365     hann[97] = 888;
   1366     hann[98] = 896;
   1367     hann[99] = 904;
   1368     hann[100] = 912;
   1369     hann[101] = 920;
   1370     hann[102] = 927;
   1371     hann[103] = 934;
   1372     hann[104] = 941;
   1373     hann[105] = 948;
   1374     hann[106] = 954;
   1375     hann[107] = 961;
   1376     hann[108] = 966;
   1377     hann[109] = 972;
   1378     hann[110] = 977;
   1379     hann[111] = 982;
   1380     hann[112] = 987;
   1381     hann[113] = 992;
   1382     hann[114] = 996;
   1383     hann[115] = 1000;
   1384     hann[116] = 1003;
   1385     hann[117] = 1007;
   1386     hann[118] = 1010;
   1387     hann[119] = 1012;
   1388     hann[120] = 1015;
   1389     hann[121] = 1017;
   1390     hann[122] = 1019;
   1391     hann[123] = 1020;
   1392     hann[124] = 1022;
   1393     hann[125] = 1023;
   1394     hann[126] = 1023;
   1395     hann[127] = 1023;
   1396     hann[128] = 1023;
   1397     hann[129] = 1023;
   1398     hann[130] = 1023;
   1399     hann[131] = 1022;
   1400     hann[132] = 1020;
   1401     hann[133] = 1019;
   1402     hann[134] = 1017;
   1403     hann[135] = 1015;
   1404     hann[136] = 1012;
   1405     hann[137] = 1010;
   1406     hann[138] = 1007;
   1407     hann[139] = 1003;
   1408     hann[140] = 1000;
   1409     hann[141] = 996;
   1410     hann[142] = 992;
   1411     hann[143] = 987;
   1412     hann[144] = 982;
   1413     hann[145] = 977;
   1414     hann[146] = 972;
   1415     hann[147] = 966;
   1416     hann[148] = 961;
   1417     hann[149] = 954;
   1418     hann[150] = 948;
   1419     hann[151] = 941;
   1420     hann[152] = 934;
   1421     hann[153] = 927;
   1422     hann[154] = 920;
   1423     hann[155] = 912;
   1424     hann[156] = 904;
   1425     hann[157] = 896;
   1426     hann[158] = 888;
   1427     hann[159] = 879;
   1428     hann[160] = 870;
   1429     hann[161] = 861;
   1430     hann[162] = 852;
   1431     hann[163] = 842;
   1432     hann[164] = 833;
   1433     hann[165] = 823;
   1434     hann[166] = 813;
   1435     hann[167] = 803;
   1436     hann[168] = 792;
   1437     hann[169] = 782;
   1438     hann[170] = 771;
   1439     hann[171] = 760;
   1440     hann[172] = 749;
   1441     hann[173] = 738;
   1442     hann[174] = 727;
   1443     hann[175] = 715;
   1444     hann[176] = 704;
   1445     hann[177] = 692;
   1446     hann[178] = 680;
   1447     hann[179] = 669;
   1448     hann[180] = 657;
   1449     hann[181] = 645;
   1450     hann[182] = 632;
   1451     hann[183] = 620;
   1452     hann[184] = 608;
   1453     hann[185] = 596;
   1454     hann[186] = 583;
   1455     hann[187] = 571;
   1456     hann[188] = 558;
   1457     hann[189] = 546;
   1458     hann[190] = 533;
   1459     hann[191] = 521;
   1460     hann[192] = 508;
   1461     hann[193] = 496;
   1462     hann[194] = 483;
   1463     hann[195] = 471;
   1464     hann[196] = 458;
   1465     hann[197] = 446;
   1466     hann[198] = 434;
   1467     hann[199] = 421;
   1468     hann[200] = 409;
   1469     hann[201] = 397;
   1470     hann[202] = 385;
   1471     hann[203] = 372;
   1472     hann[204] = 360;
   1473     hann[205] = 349;
   1474     hann[206] = 337;
   1475     hann[207] = 325;
   1476     hann[208] = 313;
   1477     hann[209] = 302;
   1478     hann[210] = 291;
   1479     hann[211] = 279;
   1480     hann[212] = 268;
   1481     hann[213] = 257;
   1482     hann[214] = 247;
   1483     hann[215] = 236;
   1484     hann[216] = 225;
   1485     hann[217] = 215;
   1486     hann[218] = 205;
   1487     hann[219] = 195;
   1488     hann[220] = 185;
   1489     hann[221] = 176;
   1490     hann[222] = 166;
   1491     hann[223] = 157;
   1492     hann[224] = 148;
   1493     hann[225] = 140;
   1494     hann[226] = 131;
   1495     hann[227] = 123;
   1496     hann[228] = 115;
   1497     hann[229] = 107;
   1498     hann[230] = 100;
   1499     hann[231] = 92;
   1500     hann[232] = 85;
   1501     hann[233] = 78;
   1502     hann[234] = 72;
   1503     hann[235] = 66;
   1504     hann[236] = 59;
   1505     hann[237] = 54;
   1506     hann[238] = 48;
   1507     hann[239] = 43;
   1508     hann[240] = 38;
   1509     hann[241] = 34;
   1510     hann[242] = 29;
   1511     hann[243] = 25;
   1512     hann[244] = 21;
   1513     hann[245] = 18;
   1514     hann[246] = 15;
   1515     hann[247] = 12;
   1516     hann[248] = 9;
   1517     hann[249] = 7;
   1518     hann[250] = 5;
   1519     hann[251] = 3;
   1520     hann[252] = 2;
   1521     hann[253] = 1;
   1522     hann[254] = 0;
   1523     hann[255] = 0;
   1524 
   1525 } /* gen_hann2 */
   1526 
   1527 /**
   1528  * Creates the non uniform enhancement window (bilinearly mapped hannning window)
   1529  * in order to reduce the formant enhancement in the high spectrum
   1530  * @param    mm : memory manager
   1531  * @param    sig_inObj : sig PU internal object of the sub-object
   1532  * @return   void
   1533  * @callgraph
   1534  * @callergraph
   1535  * @remarks the outgput is based on the matlab script below\n
   1536  matlab script
   1537  -------------
   1538  function makeEnhWind(alpha)
   1539  N=129; % =HFFTSIZEE_P1
   1540  s=(bilinmap(alpha,N));
   1541  h=[hann(N)'];
   1542  x=(1:N)/N;
   1543  W=interp1(x,h,s); W(1)=0;
   1544 
   1545  fid=fopen('enhwind.txt','wt');
   1546  fprintf(fid,'picoos_int16 enh_wind_init(sig_innerobj_t *sig_inObj) {\n');
   1547  for i=1:N,
   1548  fprintf(fid,['   sig_inObj->enhwind[' int2str(i) '] = (picoos_single)' num2str(W(i),'%0.7g') ';\n']);
   1549  end;
   1550  fprintf(fid,'   return PICO_OK;\n');
   1551  fprintf(fid,' }; \n');
   1552  fclose(fid);
   1553  %figure(1); plot(x*8,W);
   1554  */
   1555 static void enh_wind_init(sig_innerobj_t *sig_inObj)
   1556 {
   1557     /*     picoos_int16 i; */
   1558     picoos_int32 *c;
   1559 
   1560     c = sig_inObj->cos_table;
   1561 
   1562     c[0] = 4096;
   1563     c[1] = 4095;
   1564     c[2] = 4095;
   1565     c[3] = 4095;
   1566     c[4] = 4095;
   1567     c[5] = 4095;
   1568     c[6] = 4095;
   1569     c[7] = 4095;
   1570     c[8] = 4094;
   1571     c[9] = 4094;
   1572     c[10] = 4094;
   1573     c[11] = 4093;
   1574     c[12] = 4093;
   1575     c[13] = 4092;
   1576     c[14] = 4092;
   1577     c[15] = 4091;
   1578     c[16] = 4091;
   1579     c[17] = 4090;
   1580     c[18] = 4089;
   1581     c[19] = 4089;
   1582     c[20] = 4088;
   1583     c[21] = 4087;
   1584     c[22] = 4086;
   1585     c[23] = 4085;
   1586     c[24] = 4084;
   1587     c[25] = 4083;
   1588     c[26] = 4082;
   1589     c[27] = 4081;
   1590     c[28] = 4080;
   1591     c[29] = 4079;
   1592     c[30] = 4078;
   1593     c[31] = 4077;
   1594     c[32] = 4076;
   1595     c[33] = 4075;
   1596     c[34] = 4073;
   1597     c[35] = 4072;
   1598     c[36] = 4071;
   1599     c[37] = 4069;
   1600     c[38] = 4068;
   1601     c[39] = 4066;
   1602     c[40] = 4065;
   1603     c[41] = 4063;
   1604     c[42] = 4062;
   1605     c[43] = 4060;
   1606     c[44] = 4058;
   1607     c[45] = 4057;
   1608     c[46] = 4055;
   1609     c[47] = 4053;
   1610     c[48] = 4051;
   1611     c[49] = 4049;
   1612     c[50] = 4047;
   1613     c[51] = 4045;
   1614     c[52] = 4043;
   1615     c[53] = 4041;
   1616     c[54] = 4039;
   1617     c[55] = 4037;
   1618     c[56] = 4035;
   1619     c[57] = 4033;
   1620     c[58] = 4031;
   1621     c[59] = 4029;
   1622     c[60] = 4026;
   1623     c[61] = 4024;
   1624     c[62] = 4022;
   1625     c[63] = 4019;
   1626     c[64] = 4017;
   1627     c[65] = 4014;
   1628     c[66] = 4012;
   1629     c[67] = 4009;
   1630     c[68] = 4007;
   1631     c[69] = 4004;
   1632     c[70] = 4001;
   1633     c[71] = 3999;
   1634     c[72] = 3996;
   1635     c[73] = 3993;
   1636     c[74] = 3990;
   1637     c[75] = 3988;
   1638     c[76] = 3985;
   1639     c[77] = 3982;
   1640     c[78] = 3979;
   1641     c[79] = 3976;
   1642     c[80] = 3973;
   1643     c[81] = 3970;
   1644     c[82] = 3967;
   1645     c[83] = 3963;
   1646     c[84] = 3960;
   1647     c[85] = 3957;
   1648     c[86] = 3954;
   1649     c[87] = 3950;
   1650     c[88] = 3947;
   1651     c[89] = 3944;
   1652     c[90] = 3940;
   1653     c[91] = 3937;
   1654     c[92] = 3933;
   1655     c[93] = 3930;
   1656     c[94] = 3926;
   1657     c[95] = 3923;
   1658     c[96] = 3919;
   1659     c[97] = 3915;
   1660     c[98] = 3912;
   1661     c[99] = 3908;
   1662     c[100] = 3904;
   1663     c[101] = 3900;
   1664     c[102] = 3897;
   1665     c[103] = 3893;
   1666     c[104] = 3889;
   1667     c[105] = 3885;
   1668     c[106] = 3881;
   1669     c[107] = 3877;
   1670     c[108] = 3873;
   1671     c[109] = 3869;
   1672     c[110] = 3864;
   1673     c[111] = 3860;
   1674     c[112] = 3856;
   1675     c[113] = 3852;
   1676     c[114] = 3848;
   1677     c[115] = 3843;
   1678     c[116] = 3839;
   1679     c[117] = 3834;
   1680     c[118] = 3830;
   1681     c[119] = 3826;
   1682     c[120] = 3821;
   1683     c[121] = 3816;
   1684     c[122] = 3812;
   1685     c[123] = 3807;
   1686     c[124] = 3803;
   1687     c[125] = 3798;
   1688     c[126] = 3793;
   1689     c[127] = 3789;
   1690     c[128] = 3784;
   1691     c[129] = 3779;
   1692     c[130] = 3774;
   1693     c[131] = 3769;
   1694     c[132] = 3764;
   1695     c[133] = 3759;
   1696     c[134] = 3754;
   1697     c[135] = 3749;
   1698     c[136] = 3744;
   1699     c[137] = 3739;
   1700     c[138] = 3734;
   1701     c[139] = 3729;
   1702     c[140] = 3723;
   1703     c[141] = 3718;
   1704     c[142] = 3713;
   1705     c[143] = 3708;
   1706     c[144] = 3702;
   1707     c[145] = 3697;
   1708     c[146] = 3691;
   1709     c[147] = 3686;
   1710     c[148] = 3680;
   1711     c[149] = 3675;
   1712     c[150] = 3669;
   1713     c[151] = 3664;
   1714     c[152] = 3658;
   1715     c[153] = 3652;
   1716     c[154] = 3647;
   1717     c[155] = 3641;
   1718     c[156] = 3635;
   1719     c[157] = 3629;
   1720     c[158] = 3624;
   1721     c[159] = 3618;
   1722     c[160] = 3612;
   1723     c[161] = 3606;
   1724     c[162] = 3600;
   1725     c[163] = 3594;
   1726     c[164] = 3588;
   1727     c[165] = 3582;
   1728     c[166] = 3576;
   1729     c[167] = 3570;
   1730     c[168] = 3563;
   1731     c[169] = 3557;
   1732     c[170] = 3551;
   1733     c[171] = 3545;
   1734     c[172] = 3538;
   1735     c[173] = 3532;
   1736     c[174] = 3526;
   1737     c[175] = 3519;
   1738     c[176] = 3513;
   1739     c[177] = 3506;
   1740     c[178] = 3500;
   1741     c[179] = 3493;
   1742     c[180] = 3487;
   1743     c[181] = 3480;
   1744     c[182] = 3473;
   1745     c[183] = 3467;
   1746     c[184] = 3460;
   1747     c[185] = 3453;
   1748     c[186] = 3447;
   1749     c[187] = 3440;
   1750     c[188] = 3433;
   1751     c[189] = 3426;
   1752     c[190] = 3419;
   1753     c[191] = 3412;
   1754     c[192] = 3405;
   1755     c[193] = 3398;
   1756     c[194] = 3391;
   1757     c[195] = 3384;
   1758     c[196] = 3377;
   1759     c[197] = 3370;
   1760     c[198] = 3363;
   1761     c[199] = 3356;
   1762     c[200] = 3348;
   1763     c[201] = 3341;
   1764     c[202] = 3334;
   1765     c[203] = 3326;
   1766     c[204] = 3319;
   1767     c[205] = 3312;
   1768     c[206] = 3304;
   1769     c[207] = 3297;
   1770     c[208] = 3289;
   1771     c[209] = 3282;
   1772     c[210] = 3274;
   1773     c[211] = 3267;
   1774     c[212] = 3259;
   1775     c[213] = 3252;
   1776     c[214] = 3244;
   1777     c[215] = 3236;
   1778     c[216] = 3229;
   1779     c[217] = 3221;
   1780     c[218] = 3213;
   1781     c[219] = 3205;
   1782     c[220] = 3197;
   1783     c[221] = 3190;
   1784     c[222] = 3182;
   1785     c[223] = 3174;
   1786     c[224] = 3166;
   1787     c[225] = 3158;
   1788     c[226] = 3150;
   1789     c[227] = 3142;
   1790     c[228] = 3134;
   1791     c[229] = 3126;
   1792     c[230] = 3117;
   1793     c[231] = 3109;
   1794     c[232] = 3101;
   1795     c[233] = 3093;
   1796     c[234] = 3085;
   1797     c[235] = 3076;
   1798     c[236] = 3068;
   1799     c[237] = 3060;
   1800     c[238] = 3051;
   1801     c[239] = 3043;
   1802     c[240] = 3034;
   1803     c[241] = 3026;
   1804     c[242] = 3018;
   1805     c[243] = 3009;
   1806     c[244] = 3000;
   1807     c[245] = 2992;
   1808     c[246] = 2983;
   1809     c[247] = 2975;
   1810     c[248] = 2966;
   1811     c[249] = 2957;
   1812     c[250] = 2949;
   1813     c[251] = 2940;
   1814     c[252] = 2931;
   1815     c[253] = 2922;
   1816     c[254] = 2914;
   1817     c[255] = 2905;
   1818     c[256] = 2896;
   1819     c[257] = 2887;
   1820     c[258] = 2878;
   1821     c[259] = 2869;
   1822     c[260] = 2860;
   1823     c[261] = 2851;
   1824     c[262] = 2842;
   1825     c[263] = 2833;
   1826     c[264] = 2824;
   1827     c[265] = 2815;
   1828     c[266] = 2806;
   1829     c[267] = 2796;
   1830     c[268] = 2787;
   1831     c[269] = 2778;
   1832     c[270] = 2769;
   1833     c[271] = 2760;
   1834     c[272] = 2750;
   1835     c[273] = 2741;
   1836     c[274] = 2732;
   1837     c[275] = 2722;
   1838     c[276] = 2713;
   1839     c[277] = 2703;
   1840     c[278] = 2694;
   1841     c[279] = 2684;
   1842     c[280] = 2675;
   1843     c[281] = 2665;
   1844     c[282] = 2656;
   1845     c[283] = 2646;
   1846     c[284] = 2637;
   1847     c[285] = 2627;
   1848     c[286] = 2617;
   1849     c[287] = 2608;
   1850     c[288] = 2598;
   1851     c[289] = 2588;
   1852     c[290] = 2578;
   1853     c[291] = 2569;
   1854     c[292] = 2559;
   1855     c[293] = 2549;
   1856     c[294] = 2539;
   1857     c[295] = 2529;
   1858     c[296] = 2519;
   1859     c[297] = 2510;
   1860     c[298] = 2500;
   1861     c[299] = 2490;
   1862     c[300] = 2480;
   1863     c[301] = 2470;
   1864     c[302] = 2460;
   1865     c[303] = 2450;
   1866     c[304] = 2439;
   1867     c[305] = 2429;
   1868     c[306] = 2419;
   1869     c[307] = 2409;
   1870     c[308] = 2399;
   1871     c[309] = 2389;
   1872     c[310] = 2379;
   1873     c[311] = 2368;
   1874     c[312] = 2358;
   1875     c[313] = 2348;
   1876     c[314] = 2337;
   1877     c[315] = 2327;
   1878     c[316] = 2317;
   1879     c[317] = 2306;
   1880     c[318] = 2296;
   1881     c[319] = 2286;
   1882     c[320] = 2275;
   1883     c[321] = 2265;
   1884     c[322] = 2254;
   1885     c[323] = 2244;
   1886     c[324] = 2233;
   1887     c[325] = 2223;
   1888     c[326] = 2212;
   1889     c[327] = 2201;
   1890     c[328] = 2191;
   1891     c[329] = 2180;
   1892     c[330] = 2170;
   1893     c[331] = 2159;
   1894     c[332] = 2148;
   1895     c[333] = 2138;
   1896     c[334] = 2127;
   1897     c[335] = 2116;
   1898     c[336] = 2105;
   1899     c[337] = 2094;
   1900     c[338] = 2084;
   1901     c[339] = 2073;
   1902     c[340] = 2062;
   1903     c[341] = 2051;
   1904     c[342] = 2040;
   1905     c[343] = 2029;
   1906     c[344] = 2018;
   1907     c[345] = 2007;
   1908     c[346] = 1997;
   1909     c[347] = 1986;
   1910     c[348] = 1975;
   1911     c[349] = 1964;
   1912     c[350] = 1952;
   1913     c[351] = 1941;
   1914     c[352] = 1930;
   1915     c[353] = 1919;
   1916     c[354] = 1908;
   1917     c[355] = 1897;
   1918     c[356] = 1886;
   1919     c[357] = 1875;
   1920     c[358] = 1864;
   1921     c[359] = 1852;
   1922     c[360] = 1841;
   1923     c[361] = 1830;
   1924     c[362] = 1819;
   1925     c[363] = 1807;
   1926     c[364] = 1796;
   1927     c[365] = 1785;
   1928     c[366] = 1773;
   1929     c[367] = 1762;
   1930     c[368] = 1751;
   1931     c[369] = 1739;
   1932     c[370] = 1728;
   1933     c[371] = 1717;
   1934     c[372] = 1705;
   1935     c[373] = 1694;
   1936     c[374] = 1682;
   1937     c[375] = 1671;
   1938     c[376] = 1659;
   1939     c[377] = 1648;
   1940     c[378] = 1636;
   1941     c[379] = 1625;
   1942     c[380] = 1613;
   1943     c[381] = 1602;
   1944     c[382] = 1590;
   1945     c[383] = 1579;
   1946     c[384] = 1567;
   1947     c[385] = 1555;
   1948     c[386] = 1544;
   1949     c[387] = 1532;
   1950     c[388] = 1520;
   1951     c[389] = 1509;
   1952     c[390] = 1497;
   1953     c[391] = 1485;
   1954     c[392] = 1474;
   1955     c[393] = 1462;
   1956     c[394] = 1450;
   1957     c[395] = 1438;
   1958     c[396] = 1427;
   1959     c[397] = 1415;
   1960     c[398] = 1403;
   1961     c[399] = 1391;
   1962     c[400] = 1379;
   1963     c[401] = 1368;
   1964     c[402] = 1356;
   1965     c[403] = 1344;
   1966     c[404] = 1332;
   1967     c[405] = 1320;
   1968     c[406] = 1308;
   1969     c[407] = 1296;
   1970     c[408] = 1284;
   1971     c[409] = 1272;
   1972     c[410] = 1260;
   1973     c[411] = 1248;
   1974     c[412] = 1237;
   1975     c[413] = 1225;
   1976     c[414] = 1213;
   1977     c[415] = 1201;
   1978     c[416] = 1189;
   1979     c[417] = 1176;
   1980     c[418] = 1164;
   1981     c[419] = 1152;
   1982     c[420] = 1140;
   1983     c[421] = 1128;
   1984     c[422] = 1116;
   1985     c[423] = 1104;
   1986     c[424] = 1092;
   1987     c[425] = 1080;
   1988     c[426] = 1068;
   1989     c[427] = 1056;
   1990     c[428] = 1043;
   1991     c[429] = 1031;
   1992     c[430] = 1019;
   1993     c[431] = 1007;
   1994     c[432] = 995;
   1995     c[433] = 983;
   1996     c[434] = 970;
   1997     c[435] = 958;
   1998     c[436] = 946;
   1999     c[437] = 934;
   2000     c[438] = 921;
   2001     c[439] = 909;
   2002     c[440] = 897;
   2003     c[441] = 885;
   2004     c[442] = 872;
   2005     c[443] = 860;
   2006     c[444] = 848;
   2007     c[445] = 836;
   2008     c[446] = 823;
   2009     c[447] = 811;
   2010     c[448] = 799;
   2011     c[449] = 786;
   2012     c[450] = 774;
   2013     c[451] = 762;
   2014     c[452] = 749;
   2015     c[453] = 737;
   2016     c[454] = 725;
   2017     c[455] = 712;
   2018     c[456] = 700;
   2019     c[457] = 687;
   2020     c[458] = 675;
   2021     c[459] = 663;
   2022     c[460] = 650;
   2023     c[461] = 638;
   2024     c[462] = 625;
   2025     c[463] = 613;
   2026     c[464] = 601;
   2027     c[465] = 588;
   2028     c[466] = 576;
   2029     c[467] = 563;
   2030     c[468] = 551;
   2031     c[469] = 538;
   2032     c[470] = 526;
   2033     c[471] = 513;
   2034     c[472] = 501;
   2035     c[473] = 488;
   2036     c[474] = 476;
   2037     c[475] = 463;
   2038     c[476] = 451;
   2039     c[477] = 438;
   2040     c[478] = 426;
   2041     c[479] = 413;
   2042     c[480] = 401;
   2043     c[481] = 388;
   2044     c[482] = 376;
   2045     c[483] = 363;
   2046     c[484] = 351;
   2047     c[485] = 338;
   2048     c[486] = 326;
   2049     c[487] = 313;
   2050     c[488] = 301;
   2051     c[489] = 288;
   2052     c[490] = 276;
   2053     c[491] = 263;
   2054     c[492] = 251;
   2055     c[493] = 238;
   2056     c[494] = 226;
   2057     c[495] = 213;
   2058     c[496] = 200;
   2059     c[497] = 188;
   2060     c[498] = 175;
   2061     c[499] = 163;
   2062     c[500] = 150;
   2063     c[501] = 138;
   2064     c[502] = 125;
   2065     c[503] = 113;
   2066     c[504] = 100;
   2067     c[505] = 87;
   2068     c[506] = 75;
   2069     c[507] = 62;
   2070     c[508] = 50;
   2071     c[509] = 37;
   2072     c[510] = 25;
   2073     c[511] = 12;
   2074     c[512] = 0;
   2075 } /*enh_wind_init*/
   2076 
   2077 /**
   2078  * Initializes a useful large array of random numbers
   2079  * @param  sig_inObj : sig PU internal object of the sub-object
   2080  * @return void
   2081  * @callgraph
   2082  * @callergraph
   2083  */
   2084 static void init_rand(sig_innerobj_t *sig_inObj)
   2085 {
   2086     picoos_int32 *q = sig_inObj->int_vec34;
   2087     picoos_int32 *r = sig_inObj->int_vec35;
   2088 
   2089     sig_inObj->iRand = 0;
   2090 
   2091     q[0] = -2198;
   2092     r[0] = 3455;
   2093     q[1] = 3226;
   2094     r[1] = -2522;
   2095     q[2] = -845;
   2096     r[2] = 4007;
   2097     q[3] = -1227;
   2098     r[3] = 3907;
   2099     q[4] = -3480;
   2100     r[4] = 2158;
   2101     q[5] = -1325;
   2102     r[5] = -3875;
   2103     q[6] = 2089;
   2104     r[6] = -3522;
   2105     q[7] = -468;
   2106     r[7] = 4069;
   2107     q[8] = 711;
   2108     r[8] = -4033;
   2109     q[9] = 3862;
   2110     r[9] = 1362;
   2111     q[10] = 4054;
   2112     r[10] = -579;
   2113     q[11] = 2825;
   2114     r[11] = 2965;
   2115     q[12] = 2704;
   2116     r[12] = -3076;
   2117     q[13] = 4081;
   2118     r[13] = 344;
   2119     q[14] = -3912;
   2120     r[14] = 1211;
   2121     q[15] = -3541;
   2122     r[15] = 2058;
   2123     q[16] = 2694;
   2124     r[16] = 3084;
   2125     q[17] = 835;
   2126     r[17] = 4009;
   2127     q[18] = -2578;
   2128     r[18] = -3182;
   2129     q[19] = 3205;
   2130     r[19] = 2550;
   2131     q[20] = -4074;
   2132     r[20] = -418;
   2133     q[21] = -183;
   2134     r[21] = -4091;
   2135     q[22] = -2665;
   2136     r[22] = -3110;
   2137     q[23] = -1367;
   2138     r[23] = 3860;
   2139     q[24] = -2266;
   2140     r[24] = -3411;
   2141     q[25] = 3327;
   2142     r[25] = -2387;
   2143     q[26] = -2807;
   2144     r[26] = -2982;
   2145     q[27] = -3175;
   2146     r[27] = -2587;
   2147     q[28] = -4095;
   2148     r[28] = 27;
   2149     q[29] = -811;
   2150     r[29] = -4014;
   2151     q[30] = 4082;
   2152     r[30] = 332;
   2153     q[31] = -2175;
   2154     r[31] = 3470;
   2155     q[32] = 3112;
   2156     r[32] = 2662;
   2157     q[33] = 1168;
   2158     r[33] = -3925;
   2159     q[34] = 2659;
   2160     r[34] = 3115;
   2161     q[35] = 4048;
   2162     r[35] = 622;
   2163     q[36] = 4092;
   2164     r[36] = -165;
   2165     q[37] = -4036;
   2166     r[37] = 697;
   2167     q[38] = 1081;
   2168     r[38] = -3950;
   2169     q[39] = -548;
   2170     r[39] = 4059;
   2171     q[40] = 4038;
   2172     r[40] = 685;
   2173     q[41] = -511;
   2174     r[41] = 4063;
   2175     q[42] = 3317;
   2176     r[42] = -2402;
   2177     q[43] = -3180;
   2178     r[43] = 2580;
   2179     q[44] = 851;
   2180     r[44] = -4006;
   2181     q[45] = 2458;
   2182     r[45] = -3276;
   2183     q[46] = -1453;
   2184     r[46] = 3829;
   2185     q[47] = -3577;
   2186     r[47] = 1995;
   2187     q[48] = -3708;
   2188     r[48] = -1738;
   2189     q[49] = -3890;
   2190     r[49] = 1282;
   2191     q[50] = 4041;
   2192     r[50] = 666;
   2193     q[51] = -3511;
   2194     r[51] = -2108;
   2195     q[52] = -1454;
   2196     r[52] = -3828;
   2197     q[53] = 2124;
   2198     r[53] = 3502;
   2199     q[54] = -3159;
   2200     r[54] = 2606;
   2201     q[55] = 2384;
   2202     r[55] = -3330;
   2203     q[56] = -3767;
   2204     r[56] = -1607;
   2205     q[57] = -4063;
   2206     r[57] = -513;
   2207     q[58] = 3952;
   2208     r[58] = -1075;
   2209     q[59] = -3778;
   2210     r[59] = -1581;
   2211     q[60] = -301;
   2212     r[60] = -4084;
   2213     q[61] = -4026;
   2214     r[61] = 751;
   2215     q[62] = -3346;
   2216     r[62] = 2361;
   2217     q[63] = -2426;
   2218     r[63] = 3299;
   2219     q[64] = 428;
   2220     r[64] = -4073;
   2221     q[65] = 3968;
   2222     r[65] = 1012;
   2223     q[66] = 2900;
   2224     r[66] = -2892;
   2225     q[67] = -263;
   2226     r[67] = 4087;
   2227     q[68] = 4083;
   2228     r[68] = 322;
   2229     q[69] = 2024;
   2230     r[69] = 3560;
   2231     q[70] = 4015;
   2232     r[70] = 808;
   2233     q[71] = -3971;
   2234     r[71] = -1000;
   2235     q[72] = 3785;
   2236     r[72] = -1564;
   2237     q[73] = -3726;
   2238     r[73] = 1701;
   2239     q[74] = -3714;
   2240     r[74] = 1725;
   2241     q[75] = 743;
   2242     r[75] = 4027;
   2243     q[76] = 875;
   2244     r[76] = -4001;
   2245     q[77] = 294;
   2246     r[77] = 4085;
   2247     q[78] = 2611;
   2248     r[78] = 3155;
   2249     q[79] = 2491;
   2250     r[79] = -3251;
   2251     q[80] = 1558;
   2252     r[80] = 3787;
   2253     q[81] = -2063;
   2254     r[81] = -3538;
   2255     q[82] = 3809;
   2256     r[82] = -1505;
   2257     q[83] = -2987;
   2258     r[83] = 2802;
   2259     q[84] = -1955;
   2260     r[84] = 3599;
   2261     q[85] = 1980;
   2262     r[85] = -3585;
   2263     q[86] = -539;
   2264     r[86] = -4060;
   2265     q[87] = -3210;
   2266     r[87] = 2543;
   2267     q[88] = 2415;
   2268     r[88] = -3308;
   2269     q[89] = 1587;
   2270     r[89] = 3775;
   2271     q[90] = -3943;
   2272     r[90] = 1106;
   2273     q[91] = 3476;
   2274     r[91] = 2165;
   2275     q[92] = 2253;
   2276     r[92] = 3420;
   2277     q[93] = -2584;
   2278     r[93] = 3177;
   2279     q[94] = 3804;
   2280     r[94] = -1518;
   2281     q[95] = -3637;
   2282     r[95] = 1883;
   2283     q[96] = 3289;
   2284     r[96] = -2440;
   2285     q[97] = -1621;
   2286     r[97] = 3761;
   2287     q[98] = 1645;
   2288     r[98] = 3751;
   2289     q[99] = -3471;
   2290     r[99] = 2173;
   2291     q[100] = 4071;
   2292     r[100] = -449;
   2293     q[101] = -872;
   2294     r[101] = -4001;
   2295     q[102] = -3897;
   2296     r[102] = 1259;
   2297     q[103] = -3590;
   2298     r[103] = 1970;
   2299     q[104] = -2456;
   2300     r[104] = -3277;
   2301     q[105] = -3004;
   2302     r[105] = 2783;
   2303     q[106] = 2589;
   2304     r[106] = 3173;
   2305     q[107] = 3727;
   2306     r[107] = -1698;
   2307     q[108] = 2992;
   2308     r[108] = 2796;
   2309     q[109] = 794;
   2310     r[109] = -4018;
   2311     q[110] = -918;
   2312     r[110] = 3991;
   2313     q[111] = 1446;
   2314     r[111] = -3831;
   2315     q[112] = 3871;
   2316     r[112] = -1338;
   2317     q[113] = -612;
   2318     r[113] = -4049;
   2319     q[114] = -1566;
   2320     r[114] = -3784;
   2321     q[115] = 672;
   2322     r[115] = -4040;
   2323     q[116] = 3841;
   2324     r[116] = 1422;
   2325     q[117] = 3545;
   2326     r[117] = -2051;
   2327     q[118] = -1982;
   2328     r[118] = -3584;
   2329     q[119] = -3413;
   2330     r[119] = 2263;
   2331     q[120] = -3265;
   2332     r[120] = -2473;
   2333     q[121] = -2876;
   2334     r[121] = -2915;
   2335     q[122] = 4094;
   2336     r[122] = -117;
   2337     q[123] = -269;
   2338     r[123] = 4087;
   2339     q[124] = -4077;
   2340     r[124] = 391;
   2341     q[125] = -3759;
   2342     r[125] = 1626;
   2343     q[126] = 1639;
   2344     r[126] = 3753;
   2345     q[127] = 3041;
   2346     r[127] = -2743;
   2347     q[128] = 5;
   2348     r[128] = 4095;
   2349     q[129] = 2778;
   2350     r[129] = -3009;
   2351     q[130] = 1121;
   2352     r[130] = -3939;
   2353     q[131] = -455;
   2354     r[131] = -4070;
   2355     q[132] = 3532;
   2356     r[132] = 2073;
   2357     q[133] = -143;
   2358     r[133] = -4093;
   2359     q[134] = -2357;
   2360     r[134] = -3349;
   2361     q[135] = 458;
   2362     r[135] = 4070;
   2363     q[136] = -2887;
   2364     r[136] = -2904;
   2365     q[137] = -1104;
   2366     r[137] = 3944;
   2367     q[138] = -2104;
   2368     r[138] = -3513;
   2369     q[139] = 126;
   2370     r[139] = 4094;
   2371     q[140] = -3655;
   2372     r[140] = -1848;
   2373     q[141] = -3896;
   2374     r[141] = 1263;
   2375     q[142] = -3874;
   2376     r[142] = -1327;
   2377     q[143] = 4058;
   2378     r[143] = 553;
   2379     q[144] = -1831;
   2380     r[144] = -3663;
   2381     q[145] = -255;
   2382     r[145] = -4088;
   2383     q[146] = -1211;
   2384     r[146] = 3912;
   2385     q[147] = 445;
   2386     r[147] = 4071;
   2387     q[148] = 2268;
   2388     r[148] = 3410;
   2389     q[149] = -4010;
   2390     r[149] = 833;
   2391     q[150] = 2621;
   2392     r[150] = 3147;
   2393     q[151] = -250;
   2394     r[151] = 4088;
   2395     q[152] = -3409;
   2396     r[152] = -2269;
   2397     q[153] = -2710;
   2398     r[153] = -3070;
   2399     q[154] = 4063;
   2400     r[154] = 518;
   2401     q[155] = -3611;
   2402     r[155] = -1933;
   2403     q[156] = -3707;
   2404     r[156] = -1741;
   2405     q[157] = -1151;
   2406     r[157] = -3930;
   2407     q[158] = 3976;
   2408     r[158] = -983;
   2409     q[159] = -1736;
   2410     r[159] = 3709;
   2411     q[160] = 3669;
   2412     r[160] = 1820;
   2413     q[161] = -143;
   2414     r[161] = 4093;
   2415     q[162] = -3879;
   2416     r[162] = -1313;
   2417     q[163] = -2242;
   2418     r[163] = -3427;
   2419     q[164] = -4095;
   2420     r[164] = 0;
   2421     q[165] = -1159;
   2422     r[165] = -3928;
   2423     q[166] = -3155;
   2424     r[166] = 2611;
   2425     q[167] = -2887;
   2426     r[167] = -2904;
   2427     q[168] = -4095;
   2428     r[168] = 56;
   2429     q[169] = -3861;
   2430     r[169] = -1364;
   2431     q[170] = -2814;
   2432     r[170] = 2976;
   2433     q[171] = -3680;
   2434     r[171] = -1798;
   2435     q[172] = -4094;
   2436     r[172] = -107;
   2437     q[173] = -3626;
   2438     r[173] = 1903;
   2439     q[174] = 3403;
   2440     r[174] = 2278;
   2441     q[175] = -1735;
   2442     r[175] = -3710;
   2443     q[176] = -2126;
   2444     r[176] = -3500;
   2445     q[177] = 3183;
   2446     r[177] = -2577;
   2447     q[178] = -3499;
   2448     r[178] = 2128;
   2449     q[179] = -1736;
   2450     r[179] = 3709;
   2451     q[180] = 2592;
   2452     r[180] = -3170;
   2453     q[181] = 3875;
   2454     r[181] = 1326;
   2455     q[182] = 3596;
   2456     r[182] = 1960;
   2457     q[183] = 3915;
   2458     r[183] = -1202;
   2459     q[184] = 1570;
   2460     r[184] = 3783;
   2461     q[185] = -3319;
   2462     r[185] = -2400;
   2463     q[186] = 4019;
   2464     r[186] = 787;
   2465     q[187] = -187;
   2466     r[187] = 4091;
   2467     q[188] = 1370;
   2468     r[188] = -3859;
   2469     q[189] = -4091;
   2470     r[189] = 199;
   2471     q[190] = 3626;
   2472     r[190] = 1904;
   2473     q[191] = -2943;
   2474     r[191] = 2848;
   2475     q[192] = 56;
   2476     r[192] = 4095;
   2477     q[193] = 2824;
   2478     r[193] = 2966;
   2479     q[194] = -3994;
   2480     r[194] = -904;
   2481     q[195] = 56;
   2482     r[195] = 4095;
   2483     q[196] = -2045;
   2484     r[196] = 3548;
   2485     q[197] = -3653;
   2486     r[197] = 1850;
   2487     q[198] = -2864;
   2488     r[198] = 2927;
   2489     q[199] = -1996;
   2490     r[199] = 3576;
   2491     q[200] = -4061;
   2492     r[200] = 527;
   2493     q[201] = 159;
   2494     r[201] = 4092;
   2495     q[202] = -3363;
   2496     r[202] = 2336;
   2497     q[203] = -4074;
   2498     r[203] = 421;
   2499     q[204] = 2043;
   2500     r[204] = 3549;
   2501     q[205] = 4095;
   2502     r[205] = -70;
   2503     q[206] = -2107;
   2504     r[206] = -3512;
   2505     q[207] = -1973;
   2506     r[207] = 3589;
   2507     q[208] = -3138;
   2508     r[208] = 2631;
   2509     q[209] = -3625;
   2510     r[209] = -1905;
   2511     q[210] = 2413;
   2512     r[210] = 3309;
   2513     q[211] = -50;
   2514     r[211] = -4095;
   2515     q[212] = 2813;
   2516     r[212] = 2976;
   2517     q[213] = -535;
   2518     r[213] = -4060;
   2519     q[214] = 1250;
   2520     r[214] = 3900;
   2521     q[215] = 1670;
   2522     r[215] = -3739;
   2523     q[216] = 1945;
   2524     r[216] = -3604;
   2525     q[217] = -476;
   2526     r[217] = -4068;
   2527     q[218] = -3659;
   2528     r[218] = -1840;
   2529     q[219] = 2745;
   2530     r[219] = 3039;
   2531     q[220] = -674;
   2532     r[220] = -4040;
   2533     q[221] = 2383;
   2534     r[221] = 3330;
   2535     q[222] = 4086;
   2536     r[222] = 274;
   2537     q[223] = -4030;
   2538     r[223] = 730;
   2539     q[224] = 768;
   2540     r[224] = -4023;
   2541     q[225] = 3925;
   2542     r[225] = 1170;
   2543     q[226] = 785;
   2544     r[226] = 4019;
   2545     q[227] = -3101;
   2546     r[227] = -2675;
   2547     q[228] = 4030;
   2548     r[228] = -729;
   2549     q[229] = 3422;
   2550     r[229] = 2249;
   2551     q[230] = -3847;
   2552     r[230] = 1403;
   2553     q[231] = 3902;
   2554     r[231] = -1243;
   2555     q[232] = 2114;
   2556     r[232] = -3507;
   2557     q[233] = -2359;
   2558     r[233] = 3348;
   2559     q[234] = 3754;
   2560     r[234] = -1638;
   2561     q[235] = -4095;
   2562     r[235] = -83;
   2563     q[236] = 2301;
   2564     r[236] = -3388;
   2565     q[237] = 3336;
   2566     r[237] = 2375;
   2567     q[238] = -2045;
   2568     r[238] = 3548;
   2569     q[239] = -413;
   2570     r[239] = -4075;
   2571     q[240] = 1848;
   2572     r[240] = 3655;
   2573     q[241] = 4072;
   2574     r[241] = -437;
   2575     q[242] = 4069;
   2576     r[242] = -463;
   2577     q[243] = 1386;
   2578     r[243] = -3854;
   2579     q[244] = 966;
   2580     r[244] = 3980;
   2581     q[245] = -1684;
   2582     r[245] = -3733;
   2583     q[246] = 2953;
   2584     r[246] = 2837;
   2585     q[247] = -3961;
   2586     r[247] = -1040;
   2587     q[248] = 3512;
   2588     r[248] = -2107;
   2589     q[249] = 1363;
   2590     r[249] = 3862;
   2591     q[250] = 1883;
   2592     r[250] = 3637;
   2593     q[251] = 2657;
   2594     r[251] = 3116;
   2595     q[252] = 2347;
   2596     r[252] = -3356;
   2597     q[253] = -1635;
   2598     r[253] = -3755;
   2599     q[254] = 3170;
   2600     r[254] = 2593;
   2601     q[255] = 2856;
   2602     r[255] = 2935;
   2603     q[256] = 494;
   2604     r[256] = 4066;
   2605     q[257] = 1936;
   2606     r[257] = -3609;
   2607     q[258] = 245;
   2608     r[258] = 4088;
   2609     q[259] = -1211;
   2610     r[259] = -3912;
   2611     q[260] = -3600;
   2612     r[260] = 1952;
   2613     q[261] = 1632;
   2614     r[261] = 3756;
   2615     q[262] = 2341;
   2616     r[262] = 3360;
   2617     q[263] = 186;
   2618     r[263] = -4091;
   2619     q[264] = 4011;
   2620     r[264] = 829;
   2621     q[265] = -3490;
   2622     r[265] = -2143;
   2623     q[266] = 269;
   2624     r[266] = -4087;
   2625     q[267] = -2939;
   2626     r[267] = 2852;
   2627     q[268] = 1600;
   2628     r[268] = 3770;
   2629     q[269] = -3405;
   2630     r[269] = -2275;
   2631     q[270] = -3134;
   2632     r[270] = -2636;
   2633     q[271] = 2642;
   2634     r[271] = -3129;
   2635     q[272] = 3629;
   2636     r[272] = 1898;
   2637     q[273] = 3413;
   2638     r[273] = 2264;
   2639     q[274] = 2050;
   2640     r[274] = 3545;
   2641     q[275] = 988;
   2642     r[275] = -3975;
   2643     q[276] = -660;
   2644     r[276] = 4042;
   2645     q[277] = 978;
   2646     r[277] = -3977;
   2647     q[278] = 1965;
   2648     r[278] = -3593;
   2649     q[279] = -1513;
   2650     r[279] = -3806;
   2651     q[280] = -4076;
   2652     r[280] = 401;
   2653     q[281] = -4094;
   2654     r[281] = -92;
   2655     q[282] = -1914;
   2656     r[282] = 3621;
   2657     q[283] = 2006;
   2658     r[283] = -3570;
   2659     q[284] = -1550;
   2660     r[284] = -3791;
   2661     q[285] = 3774;
   2662     r[285] = -1591;
   2663     q[286] = -3958;
   2664     r[286] = 1052;
   2665     q[287] = -3576;
   2666     r[287] = 1997;
   2667     q[288] = -382;
   2668     r[288] = 4078;
   2669     q[289] = 1288;
   2670     r[289] = 3888;
   2671     q[290] = -2965;
   2672     r[290] = -2825;
   2673     q[291] = 1608;
   2674     r[291] = 3767;
   2675     q[292] = 3052;
   2676     r[292] = -2731;
   2677     q[293] = -622;
   2678     r[293] = 4048;
   2679     q[294] = -3836;
   2680     r[294] = 1434;
   2681     q[295] = -3542;
   2682     r[295] = 2056;
   2683     q[296] = -2648;
   2684     r[296] = 3124;
   2685     q[297] = -1178;
   2686     r[297] = -3922;
   2687     q[298] = -1109;
   2688     r[298] = 3942;
   2689     q[299] = 3910;
   2690     r[299] = 1217;
   2691     q[300] = 1199;
   2692     r[300] = -3916;
   2693     q[301] = -3386;
   2694     r[301] = 2303;
   2695     q[302] = -3453;
   2696     r[302] = 2202;
   2697     q[303] = -2877;
   2698     r[303] = 2914;
   2699     q[304] = 4095;
   2700     r[304] = -47;
   2701     q[305] = 3635;
   2702     r[305] = 1886;
   2703     q[306] = -2134;
   2704     r[306] = -3495;
   2705     q[307] = 613;
   2706     r[307] = -4049;
   2707     q[308] = -2700;
   2708     r[308] = 3079;
   2709     q[309] = 4091;
   2710     r[309] = -195;
   2711     q[310] = 3989;
   2712     r[310] = -927;
   2713     q[311] = -2385;
   2714     r[311] = 3329;
   2715     q[312] = 4094;
   2716     r[312] = -103;
   2717     q[313] = 1044;
   2718     r[313] = -3960;
   2719     q[314] = -1734;
   2720     r[314] = -3710;
   2721     q[315] = 1646;
   2722     r[315] = 3750;
   2723     q[316] = 575;
   2724     r[316] = 4055;
   2725     q[317] = -2629;
   2726     r[317] = -3140;
   2727     q[318] = 3266;
   2728     r[318] = 2471;
   2729     q[319] = 4091;
   2730     r[319] = -194;
   2731     q[320] = -2154;
   2732     r[320] = 3483;
   2733     q[321] = 659;
   2734     r[321] = 4042;
   2735     q[322] = -1785;
   2736     r[322] = -3686;
   2737     q[323] = -717;
   2738     r[323] = -4032;
   2739     q[324] = 4095;
   2740     r[324] = -37;
   2741     q[325] = -2963;
   2742     r[325] = -2827;
   2743     q[326] = -2645;
   2744     r[326] = -3126;
   2745     q[327] = 2619;
   2746     r[327] = -3148;
   2747     q[328] = 1855;
   2748     r[328] = -3651;
   2749     q[329] = -3726;
   2750     r[329] = 1699;
   2751     q[330] = -3437;
   2752     r[330] = 2227;
   2753     q[331] = 2948;
   2754     r[331] = 2842;
   2755     q[332] = -2125;
   2756     r[332] = 3501;
   2757     q[333] = -1700;
   2758     r[333] = 3726;
   2759     q[334] = 4094;
   2760     r[334] = -101;
   2761     q[335] = 2084;
   2762     r[335] = -3525;
   2763     q[336] = 3225;
   2764     r[336] = -2524;
   2765     q[337] = 2220;
   2766     r[337] = 3442;
   2767     q[338] = 3174;
   2768     r[338] = 2588;
   2769     q[339] = 229;
   2770     r[339] = -4089;
   2771     q[340] = -2381;
   2772     r[340] = -3332;
   2773     q[341] = -3677;
   2774     r[341] = -1803;
   2775     q[342] = -3191;
   2776     r[342] = -2567;
   2777     q[343] = 2465;
   2778     r[343] = 3270;
   2779     q[344] = 2681;
   2780     r[344] = -3096;
   2781     q[345] = 975;
   2782     r[345] = -3978;
   2783     q[346] = 2004;
   2784     r[346] = -3572;
   2785     q[347] = -3442;
   2786     r[347] = -2219;
   2787     q[348] = 3676;
   2788     r[348] = -1805;
   2789     q[349] = -3753;
   2790     r[349] = 1638;
   2791     q[350] = 3544;
   2792     r[350] = 2053;
   2793     q[351] = 397;
   2794     r[351] = -4076;
   2795     q[352] = 2221;
   2796     r[352] = 3440;
   2797     q[353] = -302;
   2798     r[353] = 4084;
   2799     q[354] = 4083;
   2800     r[354] = -323;
   2801     q[355] = -2253;
   2802     r[355] = -3420;
   2803     q[356] = -3038;
   2804     r[356] = 2746;
   2805     q[357] = 2884;
   2806     r[357] = 2908;
   2807     q[358] = 4070;
   2808     r[358] = 454;
   2809     q[359] = -1072;
   2810     r[359] = -3953;
   2811     q[360] = 3831;
   2812     r[360] = 1449;
   2813     q[361] = 3663;
   2814     r[361] = -1831;
   2815     q[362] = -1971;
   2816     r[362] = 3590;
   2817     q[363] = 3226;
   2818     r[363] = -2522;
   2819     q[364] = -145;
   2820     r[364] = -4093;
   2821     q[365] = 1882;
   2822     r[365] = -3637;
   2823     q[366] = 529;
   2824     r[366] = 4061;
   2825     q[367] = 2637;
   2826     r[367] = 3133;
   2827     q[368] = -4077;
   2828     r[368] = 389;
   2829     q[369] = 2156;
   2830     r[369] = -3482;
   2831     q[370] = -3276;
   2832     r[370] = 2458;
   2833     q[371] = -2687;
   2834     r[371] = -3090;
   2835     q[372] = 3469;
   2836     r[372] = -2177;
   2837     q[373] = -4093;
   2838     r[373] = -139;
   2839     q[374] = -850;
   2840     r[374] = 4006;
   2841     q[375] = -625;
   2842     r[375] = 4048;
   2843     q[376] = 1110;
   2844     r[376] = -3942;
   2845     q[377] = -3078;
   2846     r[377] = -2702;
   2847     q[378] = -2719;
   2848     r[378] = 3063;
   2849     q[379] = 742;
   2850     r[379] = 4028;
   2851     q[380] = -3902;
   2852     r[380] = -1245;
   2853     q[381] = 3888;
   2854     r[381] = -1287;
   2855     q[382] = -4081;
   2856     r[382] = 347;
   2857     q[383] = 1070;
   2858     r[383] = 3953;
   2859     q[384] = -996;
   2860     r[384] = -3972;
   2861     q[385] = 4041;
   2862     r[385] = -668;
   2863     q[386] = -2712;
   2864     r[386] = 3069;
   2865     q[387] = -3403;
   2866     r[387] = -2279;
   2867     q[388] = -3320;
   2868     r[388] = -2398;
   2869     q[389] = 3036;
   2870     r[389] = -2749;
   2871     q[390] = 1308;
   2872     r[390] = -3881;
   2873     q[391] = 2256;
   2874     r[391] = 3418;
   2875     q[392] = -1486;
   2876     r[392] = 3816;
   2877     q[393] = -2771;
   2878     r[393] = -3015;
   2879     q[394] = -3883;
   2880     r[394] = -1302;
   2881     q[395] = -3867;
   2882     r[395] = -1349;
   2883     q[396] = 3952;
   2884     r[396] = -1075;
   2885     q[397] = -789;
   2886     r[397] = 4019;
   2887     q[398] = 1458;
   2888     r[398] = 3827;
   2889     q[399] = 3832;
   2890     r[399] = -1446;
   2891     q[400] = -3001;
   2892     r[400] = -2787;
   2893     q[401] = 3463;
   2894     r[401] = 2186;
   2895     q[402] = 3606;
   2896     r[402] = 1942;
   2897     q[403] = 4023;
   2898     r[403] = 764;
   2899     q[404] = 3387;
   2900     r[404] = 2303;
   2901     q[405] = 2648;
   2902     r[405] = -3124;
   2903     q[406] = 1370;
   2904     r[406] = -3860;
   2905     q[407] = -3134;
   2906     r[407] = 2636;
   2907     q[408] = 4051;
   2908     r[408] = -600;
   2909     q[409] = -1977;
   2910     r[409] = -3587;
   2911     q[410] = 3160;
   2912     r[410] = 2605;
   2913     q[411] = 4042;
   2914     r[411] = 659;
   2915     q[412] = 3004;
   2916     r[412] = 2783;
   2917     q[413] = 3370;
   2918     r[413] = 2327;
   2919     q[414] = -419;
   2920     r[414] = -4074;
   2921     q[415] = -1968;
   2922     r[415] = 3591;
   2923     q[416] = -3705;
   2924     r[416] = -1746;
   2925     q[417] = -3331;
   2926     r[417] = -2383;
   2927     q[418] = -3634;
   2928     r[418] = 1888;
   2929     q[419] = -1981;
   2930     r[419] = -3584;
   2931     q[420] = 4069;
   2932     r[420] = -469;
   2933     q[421] = -628;
   2934     r[421] = -4047;
   2935     q[422] = -1900;
   2936     r[422] = 3628;
   2937     q[423] = 1039;
   2938     r[423] = -3961;
   2939     q[424] = 2554;
   2940     r[424] = -3201;
   2941     q[425] = -2955;
   2942     r[425] = 2836;
   2943     q[426] = 2286;
   2944     r[426] = -3398;
   2945     q[427] = -1624;
   2946     r[427] = 3760;
   2947     q[428] = 2213;
   2948     r[428] = 3446;
   2949     q[429] = -3989;
   2950     r[429] = -926;
   2951     q[430] = 192;
   2952     r[430] = -4091;
   2953     q[431] = -723;
   2954     r[431] = 4031;
   2955     q[432] = 2878;
   2956     r[432] = 2913;
   2957     q[433] = -2109;
   2958     r[433] = 3511;
   2959     q[434] = 1463;
   2960     r[434] = -3825;
   2961     q[435] = -741;
   2962     r[435] = -4028;
   2963     q[436] = -1314;
   2964     r[436] = -3879;
   2965     q[437] = 3115;
   2966     r[437] = 2659;
   2967     q[438] = -3160;
   2968     r[438] = -2605;
   2969     q[439] = 1868;
   2970     r[439] = 3644;
   2971     q[440] = -824;
   2972     r[440] = 4012;
   2973     q[441] = 781;
   2974     r[441] = 4020;
   2975     q[442] = -1257;
   2976     r[442] = -3898;
   2977     q[443] = 3331;
   2978     r[443] = -2382;
   2979     q[444] = 1642;
   2980     r[444] = -3752;
   2981     q[445] = 3748;
   2982     r[445] = -1650;
   2983     q[446] = -487;
   2984     r[446] = -4066;
   2985     q[447] = 3085;
   2986     r[447] = -2694;
   2987     q[448] = 4009;
   2988     r[448] = 839;
   2989     q[449] = -2308;
   2990     r[449] = -3383;
   2991     q[450] = 3850;
   2992     r[450] = 1397;
   2993     q[451] = -4078;
   2994     r[451] = -374;
   2995     q[452] = 2989;
   2996     r[452] = -2799;
   2997     q[453] = 3023;
   2998     r[453] = -2762;
   2999     q[454] = 1397;
   3000     r[454] = -3850;
   3001     q[455] = 323;
   3002     r[455] = 4083;
   3003     q[456] = 268;
   3004     r[456] = -4087;
   3005     q[457] = 2414;
   3006     r[457] = 3308;
   3007     q[458] = 3876;
   3008     r[458] = 1322;
   3009     q[459] = -3584;
   3010     r[459] = 1982;
   3011     q[460] = 1603;
   3012     r[460] = 3769;
   3013     q[461] = -1502;
   3014     r[461] = 3810;
   3015     q[462] = 1318;
   3016     r[462] = 3878;
   3017     q[463] = 1554;
   3018     r[463] = -3789;
   3019     q[464] = 2492;
   3020     r[464] = 3250;
   3021     q[465] = -4093;
   3022     r[465] = -154;
   3023     q[466] = 4008;
   3024     r[466] = 842;
   3025     q[467] = -2279;
   3026     r[467] = 3403;
   3027     q[468] = 3013;
   3028     r[468] = 2774;
   3029     q[469] = 2557;
   3030     r[469] = 3199;
   3031     q[470] = 4068;
   3032     r[470] = 475;
   3033     q[471] = 3324;
   3034     r[471] = -2392;
   3035     q[472] = 2653;
   3036     r[472] = -3120;
   3037     q[473] = 796;
   3038     r[473] = 4017;
   3039     q[474] = -1312;
   3040     r[474] = 3880;
   3041     q[475] = 1794;
   3042     r[475] = 3681;
   3043     q[476] = -2347;
   3044     r[476] = -3356;
   3045     q[477] = -4008;
   3046     r[477] = -840;
   3047     q[478] = -3773;
   3048     r[478] = -1592;
   3049     q[479] = 1609;
   3050     r[479] = 3766;
   3051     q[480] = -1564;
   3052     r[480] = -3785;
   3053     q[481] = 3004;
   3054     r[481] = 2784;
   3055     q[482] = 1258;
   3056     r[482] = 3897;
   3057     q[483] = 3729;
   3058     r[483] = 1693;
   3059     q[484] = -4095;
   3060     r[484] = -28;
   3061     q[485] = -4093;
   3062     r[485] = -146;
   3063     q[486] = 1393;
   3064     r[486] = -3851;
   3065     q[487] = 297;
   3066     r[487] = -4085;
   3067     q[488] = 2294;
   3068     r[488] = 3393;
   3069     q[489] = -2562;
   3070     r[489] = 3195;
   3071     q[490] = -1716;
   3072     r[490] = -3718;
   3073     q[491] = 2224;
   3074     r[491] = -3439;
   3075     q[492] = 2032;
   3076     r[492] = 3555;
   3077     q[493] = -2968;
   3078     r[493] = 2822;
   3079     q[494] = 2338;
   3080     r[494] = 3363;
   3081     q[495] = 1584;
   3082     r[495] = -3776;
   3083     q[496] = -3072;
   3084     r[496] = 2708;
   3085     q[497] = -1596;
   3086     r[497] = -3771;
   3087     q[498] = -2256;
   3088     r[498] = -3418;
   3089     q[499] = 4095;
   3090     r[499] = 89;
   3091     q[500] = -1949;
   3092     r[500] = 3602;
   3093     q[501] = 1844;
   3094     r[501] = 3657;
   3095     q[502] = -3375;
   3096     r[502] = 2319;
   3097     q[503] = -1481;
   3098     r[503] = -3818;
   3099     q[504] = 3228;
   3100     r[504] = -2520;
   3101     q[505] = 1116;
   3102     r[505] = 3940;
   3103     q[506] = -2783;
   3104     r[506] = 3004;
   3105     q[507] = 3915;
   3106     r[507] = 1201;
   3107     q[508] = 283;
   3108     r[508] = 4086;
   3109     q[509] = -3732;
   3110     r[509] = 1685;
   3111     q[510] = -433;
   3112     r[510] = -4072;
   3113     q[511] = -3667;
   3114     r[511] = 1823;
   3115     q[512] = 3883;
   3116     r[512] = 1300;
   3117     q[513] = -3742;
   3118     r[513] = 1663;
   3119     q[514] = 4093;
   3120     r[514] = -143;
   3121     q[515] = 3874;
   3122     r[515] = 1328;
   3123     q[516] = -3800;
   3124     r[516] = 1528;
   3125     q[517] = -1257;
   3126     r[517] = 3898;
   3127     q[518] = -1606;
   3128     r[518] = 3767;
   3129     q[519] = 3394;
   3130     r[519] = 2291;
   3131     q[520] = 2255;
   3132     r[520] = 3419;
   3133     q[521] = -4094;
   3134     r[521] = 120;
   3135     q[522] = -3767;
   3136     r[522] = 1606;
   3137     q[523] = 1849;
   3138     r[523] = -3654;
   3139     q[524] = -2883;
   3140     r[524] = 2908;
   3141     q[525] = 3469;
   3142     r[525] = 2176;
   3143     q[526] = 2654;
   3144     r[526] = 3119;
   3145     q[527] = -239;
   3146     r[527] = 4088;
   3147     q[528] = -651;
   3148     r[528] = 4043;
   3149     q[529] = -1140;
   3150     r[529] = 3934;
   3151     q[530] = 328;
   3152     r[530] = -4082;
   3153     q[531] = 3246;
   3154     r[531] = 2497;
   3155     q[532] = 4026;
   3156     r[532] = -753;
   3157     q[533] = -2041;
   3158     r[533] = -3550;
   3159     q[534] = -1154;
   3160     r[534] = 3929;
   3161     q[535] = -2710;
   3162     r[535] = 3070;
   3163     q[536] = -2860;
   3164     r[536] = 2932;
   3165     q[537] = 2097;
   3166     r[537] = 3517;
   3167     q[538] = 3492;
   3168     r[538] = -2140;
   3169     q[539] = 3123;
   3170     r[539] = 2649;
   3171     q[540] = 3360;
   3172     r[540] = 2342;
   3173     q[541] = 2498;
   3174     r[541] = 3245;
   3175     q[542] = 3976;
   3176     r[542] = 982;
   3177     q[543] = -2441;
   3178     r[543] = -3288;
   3179     q[544] = 3601;
   3180     r[544] = 1951;
   3181     q[545] = -4008;
   3182     r[545] = -842;
   3183     q[546] = 1243;
   3184     r[546] = 3902;
   3185     q[547] = 4069;
   3186     r[547] = 466;
   3187     q[548] = -2031;
   3188     r[548] = 3556;
   3189     q[549] = 4077;
   3190     r[549] = 386;
   3191     q[550] = -3112;
   3192     r[550] = -2663;
   3193     q[551] = 4087;
   3194     r[551] = -262;
   3195     q[552] = 4087;
   3196     r[552] = 266;
   3197     q[553] = -3907;
   3198     r[553] = -1228;
   3199     q[554] = -1611;
   3200     r[554] = 3765;
   3201     q[555] = 3066;
   3202     r[555] = -2715;
   3203     q[556] = 2657;
   3204     r[556] = 3117;
   3205     q[557] = 3912;
   3206     r[557] = -1213;
   3207     q[558] = -2531;
   3208     r[558] = -3220;
   3209     q[559] = 3500;
   3210     r[559] = -2127;
   3211     q[560] = -76;
   3212     r[560] = -4095;
   3213     q[561] = 3413;
   3214     r[561] = -2264;
   3215     q[562] = -4071;
   3216     r[562] = -448;
   3217     q[563] = 828;
   3218     r[563] = 4011;
   3219     q[564] = 3664;
   3220     r[564] = 1830;
   3221     q[565] = -1578;
   3222     r[565] = 3779;
   3223     q[566] = 3555;
   3224     r[566] = 2033;
   3225     q[567] = 3868;
   3226     r[567] = -1345;
   3227     q[568] = 4054;
   3228     r[568] = -580;
   3229     q[569] = -4094;
   3230     r[569] = 124;
   3231     q[570] = -3820;
   3232     r[570] = -1477;
   3233     q[571] = -3658;
   3234     r[571] = -1842;
   3235     q[572] = 2595;
   3236     r[572] = 3168;
   3237     q[573] = 3354;
   3238     r[573] = 2350;
   3239     q[574] = -701;
   3240     r[574] = -4035;
   3241     q[575] = -772;
   3242     r[575] = -4022;
   3243     q[576] = 2799;
   3244     r[576] = 2990;
   3245     q[577] = -3632;
   3246     r[577] = 1893;
   3247     q[578] = 310;
   3248     r[578] = 4084;
   3249     q[579] = 3984;
   3250     r[579] = -947;
   3251     q[580] = 3794;
   3252     r[580] = -1542;
   3253     q[581] = -2419;
   3254     r[581] = 3304;
   3255     q[582] = -3916;
   3256     r[582] = 1200;
   3257     q[583] = -3886;
   3258     r[583] = 1292;
   3259     q[584] = -3299;
   3260     r[584] = 2426;
   3261     q[585] = -437;
   3262     r[585] = 4072;
   3263     q[586] = 2053;
   3264     r[586] = -3544;
   3265     q[587] = 3987;
   3266     r[587] = 937;
   3267     q[588] = -789;
   3268     r[588] = -4019;
   3269     q[589] = 4055;
   3270     r[589] = -575;
   3271     q[590] = -3894;
   3272     r[590] = 1270;
   3273     q[591] = 4003;
   3274     r[591] = -864;
   3275     q[592] = -3060;
   3276     r[592] = 2721;
   3277     q[593] = -4009;
   3278     r[593] = 836;
   3279     q[594] = -1655;
   3280     r[594] = -3746;
   3281     q[595] = 3954;
   3282     r[595] = -1067;
   3283     q[596] = -773;
   3284     r[596] = 4022;
   3285     q[597] = -422;
   3286     r[597] = 4074;
   3287     q[598] = -3384;
   3288     r[598] = 2306;
   3289     q[599] = 195;
   3290     r[599] = -4091;
   3291     q[600] = -298;
   3292     r[600] = 4085;
   3293     q[601] = -3988;
   3294     r[601] = 931;
   3295     q[602] = 2014;
   3296     r[602] = -3566;
   3297     q[603] = 3349;
   3298     r[603] = -2357;
   3299     q[604] = 3800;
   3300     r[604] = 1526;
   3301     q[605] = 3858;
   3302     r[605] = 1374;
   3303     q[606] = 2947;
   3304     r[606] = 2844;
   3305     q[607] = -1483;
   3306     r[607] = -3818;
   3307     q[608] = 4056;
   3308     r[608] = -565;
   3309     q[609] = 2612;
   3310     r[609] = -3154;
   3311     q[610] = 2326;
   3312     r[610] = 3371;
   3313     q[611] = -3545;
   3314     r[611] = 2051;
   3315     q[612] = -1001;
   3316     r[612] = -3971;
   3317     q[613] = 3211;
   3318     r[613] = 2541;
   3319     q[614] = -2717;
   3320     r[614] = 3065;
   3321     q[615] = -3159;
   3322     r[615] = -2606;
   3323     q[616] = 2869;
   3324     r[616] = -2922;
   3325     q[617] = -1290;
   3326     r[617] = -3887;
   3327     q[618] = 2479;
   3328     r[618] = 3260;
   3329     q[619] = 3420;
   3330     r[619] = 2252;
   3331     q[620] = 1823;
   3332     r[620] = 3667;
   3333     q[621] = 3368;
   3334     r[621] = 2330;
   3335     q[622] = -3819;
   3336     r[622] = -1480;
   3337     q[623] = 3800;
   3338     r[623] = 1528;
   3339     q[624] = 3773;
   3340     r[624] = 1594;
   3341     q[625] = -189;
   3342     r[625] = -4091;
   3343     q[626] = -4067;
   3344     r[626] = -485;
   3345     q[627] = 2277;
   3346     r[627] = -3404;
   3347     q[628] = -4089;
   3348     r[628] = -233;
   3349     q[629] = -3634;
   3350     r[629] = 1889;
   3351     q[630] = 3292;
   3352     r[630] = 2437;
   3353     q[631] = -530;
   3354     r[631] = 4061;
   3355     q[632] = -3109;
   3356     r[632] = 2666;
   3357     q[633] = -3741;
   3358     r[633] = 1667;
   3359     q[634] = -1903;
   3360     r[634] = -3626;
   3361     q[635] = 3879;
   3362     r[635] = -1315;
   3363     q[636] = 4083;
   3364     r[636] = -315;
   3365     q[637] = -1148;
   3366     r[637] = 3931;
   3367     q[638] = 2630;
   3368     r[638] = 3139;
   3369     q[639] = -4001;
   3370     r[639] = 876;
   3371     q[640] = -2295;
   3372     r[640] = -3392;
   3373     q[641] = 1090;
   3374     r[641] = -3948;
   3375     q[642] = -3024;
   3376     r[642] = 2762;
   3377     q[643] = 2728;
   3378     r[643] = -3054;
   3379     q[644] = -3305;
   3380     r[644] = 2419;
   3381     q[645] = 60;
   3382     r[645] = 4095;
   3383     q[646] = 4048;
   3384     r[646] = -620;
   3385     q[647] = 589;
   3386     r[647] = -4053;
   3387     q[648] = -3867;
   3388     r[648] = 1347;
   3389     q[649] = -2944;
   3390     r[649] = -2847;
   3391     q[650] = -2721;
   3392     r[650] = 3060;
   3393     q[651] = 2928;
   3394     r[651] = 2863;
   3395     q[652] = 801;
   3396     r[652] = 4016;
   3397     q[653] = -3644;
   3398     r[653] = 1870;
   3399     q[654] = -1648;
   3400     r[654] = -3749;
   3401     q[655] = 825;
   3402     r[655] = -4012;
   3403     q[656] = -2036;
   3404     r[656] = -3553;
   3405     q[657] = -1192;
   3406     r[657] = -3918;
   3407     q[658] = 2875;
   3408     r[658] = 2916;
   3409     q[659] = -1831;
   3410     r[659] = -3663;
   3411     q[660] = -2865;
   3412     r[660] = -2926;
   3413     q[661] = -575;
   3414     r[661] = -4055;
   3415     q[662] = -3870;
   3416     r[662] = 1340;
   3417     q[663] = -4080;
   3418     r[663] = -356;
   3419     q[664] = -2176;
   3420     r[664] = -3469;
   3421     q[665] = -2986;
   3422     r[665] = -2803;
   3423     q[666] = 3978;
   3424     r[666] = -972;
   3425     q[667] = 2437;
   3426     r[667] = 3291;
   3427     q[668] = -3528;
   3428     r[668] = 2080;
   3429     q[669] = -3300;
   3430     r[669] = -2425;
   3431     q[670] = 3085;
   3432     r[670] = 2693;
   3433     q[671] = -3700;
   3434     r[671] = -1756;
   3435     q[672] = 3216;
   3436     r[672] = -2535;
   3437     q[673] = 4094;
   3438     r[673] = -91;
   3439     q[674] = 3775;
   3440     r[674] = -1589;
   3441     q[675] = 1097;
   3442     r[675] = -3946;
   3443     q[676] = -152;
   3444     r[676] = -4093;
   3445     q[677] = -3490;
   3446     r[677] = 2142;
   3447     q[678] = 3747;
   3448     r[678] = 1654;
   3449     q[679] = -1490;
   3450     r[679] = -3815;
   3451     q[680] = -3998;
   3452     r[680] = -886;
   3453     q[681] = 3726;
   3454     r[681] = -1700;
   3455     q[682] = -1600;
   3456     r[682] = 3770;
   3457     q[683] = -87;
   3458     r[683] = 4095;
   3459     q[684] = 2538;
   3460     r[684] = -3214;
   3461     q[685] = -4095;
   3462     r[685] = 52;
   3463     q[686] = -3993;
   3464     r[686] = -910;
   3465     q[687] = 4051;
   3466     r[687] = 603;
   3467     q[688] = -1242;
   3468     r[688] = -3902;
   3469     q[689] = 2155;
   3470     r[689] = 3482;
   3471     q[690] = 1270;
   3472     r[690] = 3893;
   3473     q[691] = 1919;
   3474     r[691] = -3618;
   3475     q[692] = -3145;
   3476     r[692] = 2623;
   3477     q[693] = 2475;
   3478     r[693] = 3263;
   3479     q[694] = 2226;
   3480     r[694] = -3437;
   3481     q[695] = -3894;
   3482     r[695] = -1269;
   3483     q[696] = -429;
   3484     r[696] = 4073;
   3485     q[697] = -1346;
   3486     r[697] = 3868;
   3487     q[698] = 1297;
   3488     r[698] = 3885;
   3489     q[699] = 1699;
   3490     r[699] = 3726;
   3491     q[700] = -3375;
   3492     r[700] = 2319;
   3493     q[701] = 1577;
   3494     r[701] = -3779;
   3495     q[702] = -63;
   3496     r[702] = 4095;
   3497     q[703] = 1215;
   3498     r[703] = -3911;
   3499     q[704] = -1492;
   3500     r[704] = 3814;
   3501     q[705] = -1530;
   3502     r[705] = -3799;
   3503     q[706] = 3442;
   3504     r[706] = 2218;
   3505     q[707] = -3867;
   3506     r[707] = -1349;
   3507     q[708] = -3291;
   3508     r[708] = -2437;
   3509     q[709] = -2253;
   3510     r[709] = -3420;
   3511     q[710] = -150;
   3512     r[710] = -4093;
   3513     q[711] = -2686;
   3514     r[711] = -3092;
   3515     q[712] = 3470;
   3516     r[712] = 2175;
   3517     q[713] = -3826;
   3518     r[713] = -1461;
   3519     q[714] = -3148;
   3520     r[714] = 2619;
   3521     q[715] = -3858;
   3522     r[715] = 1375;
   3523     q[716] = -3844;
   3524     r[716] = -1412;
   3525     q[717] = -3652;
   3526     r[717] = 1854;
   3527     q[718] = 4018;
   3528     r[718] = -791;
   3529     q[719] = 179;
   3530     r[719] = -4092;
   3531     q[720] = 3498;
   3532     r[720] = 2129;
   3533     q[721] = -1999;
   3534     r[721] = -3574;
   3535     q[722] = 3531;
   3536     r[722] = 2075;
   3537     q[723] = 4050;
   3538     r[723] = -606;
   3539     q[724] = -1639;
   3540     r[724] = 3753;
   3541     q[725] = -3661;
   3542     r[725] = 1835;
   3543     q[726] = 4039;
   3544     r[726] = 679;
   3545     q[727] = 3561;
   3546     r[727] = 2023;
   3547     q[728] = 528;
   3548     r[728] = 4061;
   3549     q[729] = -634;
   3550     r[729] = -4046;
   3551     q[730] = 364;
   3552     r[730] = -4079;
   3553     q[731] = 2735;
   3554     r[731] = 3048;
   3555     q[732] = 3978;
   3556     r[732] = 973;
   3557     q[733] = -4073;
   3558     r[733] = -427;
   3559     q[734] = -3722;
   3560     r[734] = 1708;
   3561     q[735] = 2356;
   3562     r[735] = -3350;
   3563     q[736] = -1125;
   3564     r[736] = -3938;
   3565     q[737] = 4054;
   3566     r[737] = 580;
   3567     q[738] = 3328;
   3568     r[738] = -2387;
   3569     q[739] = 1439;
   3570     r[739] = -3834;
   3571     q[740] = 1746;
   3572     r[740] = 3705;
   3573     q[741] = 2507;
   3574     r[741] = 3238;
   3575     q[742] = 3839;
   3576     r[742] = -1427;
   3577     q[743] = 488;
   3578     r[743] = -4066;
   3579     q[744] = 1187;
   3580     r[744] = 3920;
   3581     q[745] = 2038;
   3582     r[745] = -3552;
   3583     q[746] = -905;
   3584     r[746] = -3994;
   3585     q[747] = -236;
   3586     r[747] = 4089;
   3587     q[748] = 208;
   3588     r[748] = -4090;
   3589     q[749] = 1660;
   3590     r[749] = 3744;
   3591     q[750] = -4074;
   3592     r[750] = -415;
   3593     q[751] = -2304;
   3594     r[751] = 3385;
   3595     q[752] = -2457;
   3596     r[752] = 3276;
   3597     q[753] = 3302;
   3598     r[753] = 2423;
   3599     q[754] = 1778;
   3600     r[754] = -3689;
   3601     q[755] = 2019;
   3602     r[755] = 3563;
   3603     q[756] = 4037;
   3604     r[756] = 687;
   3605     q[757] = -2365;
   3606     r[757] = 3343;
   3607     q[758] = 5;
   3608     r[758] = -4095;
   3609     q[759] = 160;
   3610     r[759] = -4092;
   3611 
   3612 } /*initRand*/
   3613 
   3614 /**
   3615  * initializes the MEL-2_LINEAR LOOKUP TABLE
   3616  * @param   sig_inObj : sig PU internal object of the sub-object
   3617  * @return  void
   3618  * @remarks translated from matlab code to c-code
   3619  * @callgraph
   3620  * @callergraph
   3621  *
   3622  *  input
   3623  *  - bilinTab : base address of bilinTable destination vector
   3624  *  - alpha : warping factor
   3625  *  - size : size of the vectors to be generated
   3626  *  - A,B : base address of array of indexes for lookup table implementation
   3627  *  - D : base address of delta array for lookup table implementation
   3628  *
   3629  *  output
   3630  *  - bilinTab, A, B, D  : initialized vectors
   3631  */
   3632 void mel_2_lin_init(sig_innerobj_t *sig_inObj)
   3633 {
   3634 
   3635     /*Declare variables tied to  I/O PARAMS formerly passed by value or reference*/
   3636     picoos_single alpha;
   3637     picoos_int32 *D;
   3638     picoos_int32 size;
   3639     picoos_int16 *A;
   3640 
   3641     /*Link local variables with sig data object*/
   3642 
   3643     alpha = sig_inObj->warp_p;
   3644     size = (sig_inObj->hfftsize_p) + 1;
   3645     A = sig_inObj->A_p;
   3646     D = sig_inObj->d_p;
   3647     /*
   3648      fixed point interpolation tables
   3649      scaling factor: 0x20
   3650      corresponding bit shift: 5
   3651      */
   3652 
   3653     A[0] = 0;
   3654     D[0] = 0;
   3655     A[1] = 2;
   3656     D[1] = 14;
   3657     A[2] = 4;
   3658     D[2] = 29;
   3659     A[3] = 7;
   3660     D[3] = 11;
   3661     A[4] = 9;
   3662     D[4] = 24;
   3663     A[5] = 12;
   3664     D[5] = 5;
   3665     A[6] = 14;
   3666     D[6] = 18;
   3667     A[7] = 16;
   3668     D[7] = 30;
   3669     A[8] = 19;
   3670     D[8] = 9;
   3671     A[9] = 21;
   3672     D[9] = 19;
   3673     A[10] = 23;
   3674     D[10] = 29;
   3675     A[11] = 26;
   3676     D[11] = 5;
   3677     A[12] = 28;
   3678     D[12] = 12;
   3679     A[13] = 30;
   3680     D[13] = 19;
   3681     A[14] = 32;
   3682     D[14] = 24;
   3683     A[15] = 34;
   3684     D[15] = 27;
   3685     A[16] = 36;
   3686     D[16] = 30;
   3687     A[17] = 38;
   3688     D[17] = 31;
   3689     A[18] = 40;
   3690     D[18] = 31;
   3691     A[19] = 42;
   3692     D[19] = 29;
   3693     A[20] = 44;
   3694     D[20] = 26;
   3695     A[21] = 46;
   3696     D[21] = 22;
   3697     A[22] = 48;
   3698     D[22] = 17;
   3699     A[23] = 50;
   3700     D[23] = 10;
   3701     A[24] = 52;
   3702     D[24] = 2;
   3703     A[25] = 53;
   3704     D[25] = 24;
   3705     A[26] = 55;
   3706     D[26] = 13;
   3707     A[27] = 57;
   3708     D[27] = 1;
   3709     A[28] = 58;
   3710     D[28] = 20;
   3711     A[29] = 60;
   3712     D[29] = 5;
   3713     A[30] = 61;
   3714     D[30] = 21;
   3715     A[31] = 63;
   3716     D[31] = 4;
   3717     A[32] = 64;
   3718     D[32] = 18;
   3719     A[33] = 65;
   3720     D[33] = 31;
   3721     A[34] = 67;
   3722     D[34] = 11;
   3723     A[35] = 68;
   3724     D[35] = 21;
   3725     A[36] = 69;
   3726     D[36] = 31;
   3727     A[37] = 71;
   3728     D[37] = 7;
   3729     A[38] = 72;
   3730     D[38] = 14;
   3731     A[39] = 73;
   3732     D[39] = 21;
   3733     A[40] = 74;
   3734     D[40] = 27;
   3735     A[41] = 75;
   3736     D[41] = 31;
   3737     A[42] = 77;
   3738     D[42] = 3;
   3739     A[43] = 78;
   3740     D[43] = 6;
   3741     A[44] = 79;
   3742     D[44] = 8;
   3743     A[45] = 80;
   3744     D[45] = 10;
   3745     A[46] = 81;
   3746     D[46] = 10;
   3747     A[47] = 82;
   3748     D[47] = 10;
   3749     A[48] = 83;
   3750     D[48] = 9;
   3751     A[49] = 84;
   3752     D[49] = 8;
   3753     A[50] = 85;
   3754     D[50] = 6;
   3755     A[51] = 86;
   3756     D[51] = 3;
   3757     A[52] = 86;
   3758     D[52] = 31;
   3759     A[53] = 87;
   3760     D[53] = 27;
   3761     A[54] = 88;
   3762     D[54] = 23;
   3763     A[55] = 89;
   3764     D[55] = 18;
   3765     A[56] = 90;
   3766     D[56] = 12;
   3767     A[57] = 91;
   3768     D[57] = 6;
   3769     A[58] = 91;
   3770     D[58] = 31;
   3771     A[59] = 92;
   3772     D[59] = 24;
   3773     A[60] = 93;
   3774     D[60] = 16;
   3775     A[61] = 94;
   3776     D[61] = 8;
   3777     A[62] = 94;
   3778     D[62] = 31;
   3779     A[63] = 95;
   3780     D[63] = 22;
   3781     A[64] = 96;
   3782     D[64] = 13;
   3783     A[65] = 97;
   3784     D[65] = 3;
   3785     A[66] = 97;
   3786     D[66] = 25;
   3787     A[67] = 98;
   3788     D[67] = 14;
   3789     A[68] = 99;
   3790     D[68] = 3;
   3791     A[69] = 99;
   3792     D[69] = 24;
   3793     A[70] = 100;
   3794     D[70] = 13;
   3795     A[71] = 101;
   3796     D[71] = 1;
   3797     A[72] = 101;
   3798     D[72] = 21;
   3799     A[73] = 102;
   3800     D[73] = 8;
   3801     A[74] = 102;
   3802     D[74] = 27;
   3803     A[75] = 103;
   3804     D[75] = 14;
   3805     A[76] = 104;
   3806     D[76] = 1;
   3807     A[77] = 104;
   3808     D[77] = 19;
   3809     A[78] = 105;
   3810     D[78] = 6;
   3811     A[79] = 105;
   3812     D[79] = 24;
   3813     A[80] = 106;
   3814     D[80] = 9;
   3815     A[81] = 106;
   3816     D[81] = 27;
   3817     A[82] = 107;
   3818     D[82] = 12;
   3819     A[83] = 107;
   3820     D[83] = 29;
   3821     A[84] = 108;
   3822     D[84] = 14;
   3823     A[85] = 108;
   3824     D[85] = 31;
   3825     A[86] = 109;
   3826     D[86] = 15;
   3827     A[87] = 109;
   3828     D[87] = 31;
   3829     A[88] = 110;
   3830     D[88] = 16;
   3831     A[89] = 110;
   3832     D[89] = 31;
   3833     A[90] = 111;
   3834     D[90] = 15;
   3835     A[91] = 111;
   3836     D[91] = 31;
   3837     A[92] = 112;
   3838     D[92] = 14;
   3839     A[93] = 112;
   3840     D[93] = 30;
   3841     A[94] = 113;
   3842     D[94] = 13;
   3843     A[95] = 113;
   3844     D[95] = 28;
   3845     A[96] = 114;
   3846     D[96] = 11;
   3847     A[97] = 114;
   3848     D[97] = 26;
   3849     A[98] = 115;
   3850     D[98] = 9;
   3851     A[99] = 115;
   3852     D[99] = 23;
   3853     A[100] = 116;
   3854     D[100] = 6;
   3855     A[101] = 116;
   3856     D[101] = 20;
   3857     A[102] = 117;
   3858     D[102] = 2;
   3859     A[103] = 117;
   3860     D[103] = 16;
   3861     A[104] = 117;
   3862     D[104] = 31;
   3863     A[105] = 118;
   3864     D[105] = 13;
   3865     A[106] = 118;
   3866     D[106] = 27;
   3867     A[107] = 119;
   3868     D[107] = 8;
   3869     A[108] = 119;
   3870     D[108] = 22;
   3871     A[109] = 120;
   3872     D[109] = 4;
   3873     A[110] = 120;
   3874     D[110] = 17;
   3875     A[111] = 120;
   3876     D[111] = 31;
   3877     A[112] = 121;
   3878     D[112] = 13;
   3879     A[113] = 121;
   3880     D[113] = 26;
   3881     A[114] = 122;
   3882     D[114] = 8;
   3883     A[115] = 122;
   3884     D[115] = 21;
   3885     A[116] = 123;
   3886     D[116] = 2;
   3887     A[117] = 123;
   3888     D[117] = 15;
   3889     A[118] = 123;
   3890     D[118] = 29;
   3891     A[119] = 124;
   3892     D[119] = 10;
   3893     A[120] = 124;
   3894     D[120] = 23;
   3895     A[121] = 125;
   3896     D[121] = 4;
   3897     A[122] = 125;
   3898     D[122] = 17;
   3899     A[123] = 125;
   3900     D[123] = 31;
   3901     A[124] = 126;
   3902     D[124] = 12;
   3903     A[125] = 126;
   3904     D[125] = 25;
   3905     A[126] = 127;
   3906     D[126] = 6;
   3907     A[127] = 127;
   3908     D[127] = 19;
   3909     A[128] = 128;
   3910     D[128] = 0;
   3911 
   3912 }/*mel_2_lin_init*/
   3913 
   3914 /**
   3915  * function to be documented
   3916  * @param    ang : ??
   3917  * @param    table : ??
   3918  * @param    cs : ??
   3919  * @param    sn : ??
   3920  * @return  void
   3921  * @callgraph
   3922  * @callergraph
   3923  */
   3924 static void get_trig(picoos_int32 ang, picoos_int32 *table, picoos_int32 *cs,
   3925         picoos_int32 *sn)
   3926 {
   3927     picoos_int32 i, j, k;
   3928 
   3929     i = k = ang >> PICODSP_PI_SHIFT; /*  * PICODSP_COS_TABLE_LEN2/PICODSP_FIX_SCALE2 */
   3930     if (i < 0)
   3931         i = -i;
   3932     j = 1;
   3933     i &= (PICODSP_COS_TABLE_LEN4 - 1);
   3934     if (i > PICODSP_COS_TABLE_LEN2)
   3935         i = PICODSP_COS_TABLE_LEN4 - i;
   3936     if (i > PICODSP_COS_TABLE_LEN) {
   3937         j = -1;
   3938         i = PICODSP_COS_TABLE_LEN2 - i;
   3939     }
   3940     if (j == 1)
   3941         *cs = table[i];
   3942     else
   3943         *cs = -table[i];
   3944 
   3945     i = k - PICODSP_COS_TABLE_LEN;
   3946     if (i < 0)
   3947         i = -i;
   3948     j = 1;
   3949     i &= (PICODSP_COS_TABLE_LEN4 - 1);
   3950     if (i > PICODSP_COS_TABLE_LEN2)
   3951         i = PICODSP_COS_TABLE_LEN4 - i;
   3952     if (i > PICODSP_COS_TABLE_LEN) {
   3953         j = -1;
   3954         i = PICODSP_COS_TABLE_LEN2 - i;
   3955     }
   3956     if (j == 1)
   3957         *sn = table[i];
   3958     else
   3959         *sn = -table[i];
   3960 }/*get_trig*/
   3961 
   3962 /**
   3963  * function to be documented
   3964  * @param    sig_inObj : sig PU internal object of the sub-object
   3965  * @return  void
   3966  * @callgraph
   3967  * @callergraph
   3968  */
   3969 void save_transition_frame(sig_innerobj_t *sig_inObj)
   3970 {
   3971     picoos_int32 *tmp, *tmp2; /*for loop unrolling*/
   3972 
   3973     if (sig_inObj->voiced_p != sig_inObj->prevVoiced_p) {
   3974         sig_inObj->VoicTrans = sig_inObj->prevVoiced_p; /*remember last voicing transition*/
   3975         tmp = sig_inObj->ImpResp_p;
   3976         tmp2 = sig_inObj->imp_p;
   3977         FAST_DEVICE(PICODSP_FFTSIZE,*(tmp++)=*(tmp2++););
   3978         if (sig_inObj->voiced_p == 1)
   3979             sig_inObj->nV = 0;
   3980         else
   3981             sig_inObj->nU = 0; /*to avoid problems in case of very short voiced or unvoiced parts (less than 4 frames long)*/
   3982     }
   3983 }/*save_transition_frame*/
   3984 
   3985 /**
   3986  * calculates an unweighted excitation window
   3987  * @param    sig_inObj : sig PU internal object of the sub-object
   3988  * @param    nextPeak : position of next peak (excitation position)
   3989  * @return  PICO_OK
   3990  * @callgraph
   3991  * @callergraph
   3992  * input
   3993  * - hop : hop size in samples
   3994  * - winlen : excitation window length
   3995  * - E : energy
   3996  * - F0 : pitch
   3997  * - nextPeak : state that remembers next excitation index
   3998  * - Fs - sampling frequency
   3999  * output
   4000  * - LocV, LocU : (MAX_EN size) location of excitation points
   4001  * - EnV,  EnU    : (MAX_EN size) RMS values of excitation (scaled)
   4002  * - nV, nU :    (integers) number of excitation points
   4003  * - nextPeak    new position of lastPeak to calculate next frame
   4004  */
   4005 static void get_simple_excitation(sig_innerobj_t *sig_inObj,
   4006         picoos_int16 *nextPeak)
   4007 {
   4008     /*Define local variables*/
   4009     picoos_int16 nI, nJ, k;
   4010     /* picoos_single    InvSqrt3=(picoos_single)2/(picoos_single)sqrt(3.0); *//*constant*/
   4011     picoos_int32 Ti, sqrtTi;
   4012     picoos_int16 hop, winlen, Fs;
   4013     picoos_single E, F0;
   4014     picoos_int16 voiced;
   4015     picoos_single fact; /*normalization factor*/
   4016     picoos_single rounding = 0.5f;
   4017 
   4018     /*Link local variables to sig object*/
   4019     hop = sig_inObj->hop_p;
   4020     winlen = sig_inObj->m2_p;
   4021     Fs = sig_inObj->Fs_p;
   4022     E = sig_inObj->E_p;
   4023     F0 = sig_inObj->F0_p;
   4024     voiced = sig_inObj->voiced_p;
   4025 
   4026     E = (E > 5) ? 9 : (E > 1) ? 2 * E - 1 : E;
   4027 
   4028 
   4029     /* shift previous excitation window by hop samples*/
   4030     for (nI = 0; nI < sig_inObj->nV; nI++) {
   4031         sig_inObj->LocV[nI] = sig_inObj->LocV[nI] - hop;
   4032     }
   4033     for (nI = 0; nI < sig_inObj->nU; nI++) {
   4034         sig_inObj->LocU[nI] = sig_inObj->LocU[nI] - hop;
   4035     }
   4036 
   4037     /*get rid of the voiced points that fall out of the interval*/
   4038     nI = 0;
   4039     while ((sig_inObj->LocV[nI] < 0) && nI < sig_inObj->nV)
   4040         nI++;
   4041 
   4042     for (nJ = nI; nJ < sig_inObj->nV; nJ++) {
   4043         sig_inObj->LocV[nJ - nI] = sig_inObj->LocV[nJ];
   4044         sig_inObj->EnV[nJ - nI] = sig_inObj->EnV[nJ];
   4045     }
   4046     sig_inObj->nV -= nI;
   4047     /*get rid of the unvoiced points that fall out of the interval */
   4048     nI = 0;
   4049     while ((sig_inObj->LocU[nI] < 0) && nI < sig_inObj->nU)
   4050         nI++;
   4051 
   4052     for (nJ = nI; nJ < sig_inObj->nU; nJ++) {
   4053         sig_inObj->LocU[nJ - nI] = sig_inObj->LocU[nJ];
   4054         sig_inObj->EnU[nJ - nI] = sig_inObj->EnU[nJ];
   4055     }
   4056     sig_inObj->nU -= nI;
   4057 
   4058     *nextPeak -= hop;
   4059     k = *nextPeak;
   4060 
   4061     fact = 3;
   4062     if (voiced == 0) { /*Unvoiced*/
   4063 
   4064         Ti = (picoos_int32) (rounding + (picoos_single) Fs
   4065                 / (picoos_single) sig_inObj->Fuv_p); /* round Period*/
   4066         sqrtTi = (picoos_int32) (E * sqrt((double) Fs
   4067                 / (hop * sig_inObj->Fuv_p)) * fact * PICODSP_GETEXC_K1);
   4068         while (k < winlen) {
   4069             if (k < winlen) {
   4070                 sig_inObj->LocU[sig_inObj->nU] = k;
   4071                 sig_inObj->EnU[sig_inObj->nU] = sqrtTi;
   4072                 sig_inObj->nU++;
   4073                 k += (picoos_int16) Ti;
   4074             }
   4075         }
   4076     } else { /*Voiced*/
   4077         Ti
   4078                 = (picoos_int32) (rounding + (picoos_single) Fs
   4079                         / (picoos_single) F0); /*Period*/
   4080         sqrtTi = (picoos_int32) (E
   4081                 * sqrt((double) Fs / (hop * sig_inObj->F0_p)) * fact
   4082                 * PICODSP_GETEXC_K1);
   4083         while (k < winlen) {
   4084             sig_inObj->LocV[sig_inObj->nV] = k;
   4085             sig_inObj->EnV[sig_inObj->nV] = sqrtTi;
   4086             sig_inObj->nV++;
   4087             k += (picoos_int16) Ti;
   4088         }
   4089     }
   4090     *nextPeak = k;
   4091 
   4092 }/*get_simple_excitation*/
   4093 
   4094 #ifdef __cplusplus
   4095 }
   4096 #endif
   4097 
   4098 /* end picosig2.c */
   4099