Home | History | Annotate | Download | only in cfront 
      1 /*---------------------------------------------------------------------------*
      2  *  chelfep.c  *
      3  *                                                                           *
      4  *  Copyright 2007, 2008 Nuance Communciations, Inc.                               *
      5  *                                                                           *
      6  *  Licensed under the Apache License, Version 2.0 (the 'License');          *
      7  *  you may not use this file except in compliance with the License.         *
      8  *                                                                           *
      9  *  You may obtain a copy of the License at                                  *
     10  *      http://www.apache.org/licenses/LICENSE-2.0                           *
     11  *                                                                           *
     12  *  Unless required by applicable law or agreed to in writing, software      *
     13  *  distributed under the License is distributed on an 'AS IS' BASIS,        *
     14  *  WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. *
     15  *  See the License for the specific language governing permissions and      *
     16  *  limitations under the License.                                           *
     17  *                                                                           *
     18  *---------------------------------------------------------------------------*/
     19 
     20 
     21 #ifndef _RTT
     22 #include <stdio.h>
     23 #endif
     24 #include <stdlib.h>
     25 #include <string.h>
     26 
     27 
     28 #ifdef unix
     29 #include <unistd.h>
     30 #endif
     31 #include <assert.h>
     32 
     33 #ifndef _RTT
     34 #include "duk_io.h"
     35 #endif
     36 #include "voicing.h"
     37 #include "front.h"
     38 #include "portable.h"
     39 
     40 
     41 #include "../clib/memmove.h"
     42 #include "sh_down.h"
     43 
     44 #define SMOOTH_C0_FOR_VOICING 1
     45 
     46 #if SMOOTH_C0_FOR_VOICING
     47 static featdata smoothed_c0(front_cep *cepobj, front_channel *channel);
     48 #endif
     49 
     50 int make_frame(front_channel *channel, front_wave *waveobj,
     51                front_freq *freqobj, front_cep *cepobj,
     52                voicing_info *voice,
     53                samdata *inFramesWorth, samdata *refFramesWorth,
     54                int num_samples,
     55                featdata *framdata, featdata *voicedata)
     56 {
     57 #if SMOOTH_C0_FOR_VOICING
     58   featdata smooth_c0;
     59 #endif
     60   if (freqobj->do_filterbank_input)
     61     /* This memmove is in filterbank_emulation  */
     62     MEMMOVE(channel->cep + (channel->mel_dim + 1), channel->cep,
     63             (Q2 - 1) *(channel->mel_dim + 1), sizeof(float));
     64   else
     65   {
     66 		/* not fb input */
     67 
     68     /*  JFR data processing is no longer supported. BP 21-Jul98. */
     69 
     70     /* 2. CEP data processing */
     71     filterbank_emulation(channel, waveobj, freqobj, cepobj,
     72                          inFramesWorth, refFramesWorth, num_samples);
     73     /* if doing fb dump then skip frame making. A top level will
     74         read channel->fbo and dump it. */
     75     if (freqobj->do_filterbank_dump)
     76       return True;
     77   }
     78   cepstrum_params(channel, waveobj, freqobj, cepobj);
     79 
     80   /* 4. Delta CEP data processing */
     81   (void) make_std_frame(channel,  cepobj, framdata);
     82   if (!channel->frame_valid)
     83     return (channel->frame_valid);
     84 
     85 
     86   /* 5. Voicing analysis */
     87   if (channel->frame_valid)
     88   {
     89     if (voice != NULL)
     90     {
     91 #if SMOOTH_C0_FOR_VOICING
     92       if (cepobj->do_smooth_c0)
     93       {
     94         smooth_c0 = smoothed_c0(cepobj, channel);
     95         *voicedata = (featdata) voicing_analysis(voice, smooth_c0, NULL);
     96       }
     97       else
     98         *voicedata = (featdata) voicing_analysis(voice, framdata[0], NULL);
     99 
    100 #else
    101       *voicedata = (featdata) voicing_analysis(voice, framdata[0], NULL);
    102 #endif
    103     }
    104     if (cepobj->do_skip_even_frames)
    105       channel->frame_valid = (channel->frame_count) % 2;
    106   }
    107   return (channel->frame_valid);
    108 }
    109 
    110 static featdata smoothed_c0(front_cep *cepobj, front_channel *channel)
    111 {
    112   cepdata  val;
    113   featdata fval;
    114   bigdata  scval;
    115 
    116   if (channel->frame_count <= 1)
    117     val = channel->cep[DELTA * (channel->mel_dim+1)];
    118   else
    119   {
    120     val = (channel->cep[(DELTA-1) * (channel->mel_dim+1)] >> 2)
    121           + (channel->cep[DELTA * (channel->mel_dim+1)] >> 1)
    122           + (channel->cep[(DELTA+1) * (channel->mel_dim+1)] >> 2);
    123   }
    124 
    125   /*  Now scaling and byteranging
    126   */
    127   /* Now take the costable scaling off the ceps. */
    128   ASSERT((cepobj->melA_scale[0] *(float)SHIFT_DOWN(val, COSINE_TABLE_SHIFT))
    129          < LONG_MAX);
    130   ASSERT((cepobj->melA_scale[0] *(float)SHIFT_DOWN(val, COSINE_TABLE_SHIFT))
    131          > -LONG_MAX);
    132 
    133   scval = (bigdata)(SHIFT_DOWN((bigdata)cepobj->melA_scale[0]
    134                                * (bigdata) SHIFT_DOWN(val, COSINE_TABLE_SHIFT)
    135                                + (bigdata)cepobj->melB_scale[0], BYTERANGE_SHIFT + LOG_SCALE_SHIFT));
    136   fval = (featdata) MAKEBYTE(scval);
    137   return (fval);
    138 }
    139 
    140