ref: ea8f175fa4fda3ae8f6e1c0d07d4a831b288d16b
dir: /dnn/osce.c/
/* Copyright (c) 2023 Amazon
Written by Jan Buethe */
/*
Redistribution and use in source and binary forms, with or without
modification, are permitted provided that the following conditions
are met:
- Redistributions of source code must retain the above copyright
notice, this list of conditions and the following disclaimer.
- Redistributions in binary form must reproduce the above copyright
notice, this list of conditions and the following disclaimer in the
documentation and/or other materials provided with the distribution.
THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER
OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
*/
#ifdef HAVE_CONFIG_H
#include "config.h"
#endif
#include <math.h>
#include "osce.h"
#include "osce_features.h"
#include "os_support.h"
#include "nndsp.h"
#include "float_cast.h"
#include "arch.h"
/*#define OSCE_DEBUG*/
#ifdef OSCE_DEBUG
#include <stdio.h>
/*#define WRITE_FEATURES*/
/*#define DEBUG_LACE*/
/*#define DEBUG_NOLACE*/
#define DEBUG_BBWENET
#define FINIT(fid, name, mode) do{if (fid == NULL) {fid = fopen(name, mode);}} while(0)
#endif
#if 0
#include <stdio.h>
static void print_float_array(FILE *fid, const char *name, const float *array, int n)
{
int i;
for (i = 0; i < n; i++)
{
fprintf(fid, "%s[%d]: %f\n", name, i, array[i]);
}
}
static void print_int_array(FILE *fid, const char *name, const int *array, int n)
{
int i;
for (i = 0; i < n; i++)
{
fprintf(fid, "%s[%d]: %d\n", name, i, array[i]);
}
}
static void print_int8_array(FILE *fid, const char *name, const opus_int8 *array, int n)
{
int i;
for (i = 0; i < n; i++)
{
fprintf(fid, "%s[%d]: %d\n", name, i, array[i]);
}
}
static void print_linear_layer(FILE *fid, const char *name, LinearLayer *layer)
{
int i, n_in, n_out, n_total;
char tmp[256];
n_in = layer->nb_inputs;
n_out = layer->nb_outputs;
n_total = n_in * n_out;
fprintf(fid, "\nprinting layer %s...\n", name);
fprintf(fid, "%s.nb_inputs: %d\n%s.nb_outputs: %d\n", name, n_in, name, n_out);
if (layer->bias !=NULL){}
if (layer->subias !=NULL){}
if (layer->weights !=NULL){}
if (layer->float_weights !=NULL){}
if (layer->bias != NULL) {sprintf(tmp, "%s.bias", name); print_float_array(fid, tmp, layer->bias, n_out);}
if (layer->subias != NULL) {sprintf(tmp, "%s.subias", name); print_float_array(fid, tmp, layer->subias, n_out);}
if (layer->weights != NULL) {sprintf(tmp, "%s.weights", name); print_int8_array(fid, tmp, layer->weights, n_total);}
if (layer->float_weights != NULL) {sprintf(tmp, "%s.float_weights", name); print_float_array(fid, tmp, layer->float_weights, n_total);}
/*if (layer->weights_idx != NULL) {sprintf(tmp, "%s.weights_idx", name); print_float_array(fid, tmp, layer->weights_idx, n_total);}*/
if (layer->diag != NULL) {sprintf(tmp, "%s.diag", name); print_float_array(fid, tmp, layer->diag, n_in);}
if (layer->scale != NULL) {sprintf(tmp, "%s.scale", name); print_float_array(fid, tmp, layer->scale, n_out);}
}
#endif
#ifdef ENABLE_OSCE_TRAINING_DATA
#include <stdio.h>
#endif
#define CLIP(a, min, max) (((a) < (min) ? (min) : (a)) > (max) ? (max) : (a))
extern const WeightArray lacelayers_arrays[];
extern const WeightArray nolacelayers_arrays[];
extern const WeightArray bbwenetlayers_arrays[];
/* LACE */
#ifndef DISABLE_LACE
static void compute_lace_numbits_embedding(float *emb, float numbits, int dim, float min_val, float max_val, int logscale)
{
float x;
(void) dim;
numbits = logscale ? log(numbits) : numbits;
x = CLIP(numbits, min_val, max_val) - (max_val + min_val) / 2;
emb[0] = sin(x * LACE_NUMBITS_SCALE_0 - 0.5f);
emb[1] = sin(x * LACE_NUMBITS_SCALE_1 - 0.5f);
emb[2] = sin(x * LACE_NUMBITS_SCALE_2 - 0.5f);
emb[3] = sin(x * LACE_NUMBITS_SCALE_3 - 0.5f);
emb[4] = sin(x * LACE_NUMBITS_SCALE_4 - 0.5f);
emb[5] = sin(x * LACE_NUMBITS_SCALE_5 - 0.5f);
emb[6] = sin(x * LACE_NUMBITS_SCALE_6 - 0.5f);
emb[7] = sin(x * LACE_NUMBITS_SCALE_7 - 0.5f);
}
static int init_lace(LACE *hLACE, const WeightArray *weights)
{
int ret = 0;
OPUS_CLEAR(hLACE, 1);
celt_assert(weights != NULL);
ret = init_lacelayers(&hLACE->layers, weights);
compute_overlap_window(hLACE->window, LACE_OVERLAP_SIZE);
return ret;
}
static void reset_lace_state(LACEState *state)
{
OPUS_CLEAR(state, 1);
init_adacomb_state(&state->cf1_state);
init_adacomb_state(&state->cf2_state);
init_adaconv_state(&state->af1_state);
}
static void lace_feature_net(
LACE *hLACE,
LACEState *state,
float *output,
const float *features,
const float *numbits,
const int *periods,
int arch
)
{
float input_buffer[IMAX(4 * IMAX(LACE_COND_DIM, LACE_HIDDEN_FEATURE_DIM), LACE_NUM_FEATURES + LACE_PITCH_EMBEDDING_DIM + 2*LACE_NUMBITS_EMBEDDING_DIM)];
float output_buffer[4 * IMAX(LACE_COND_DIM, LACE_HIDDEN_FEATURE_DIM)];
float numbits_embedded[2 * LACE_NUMBITS_EMBEDDING_DIM];
int i_subframe;
compute_lace_numbits_embedding(numbits_embedded, numbits[0], LACE_NUMBITS_EMBEDDING_DIM,
log(LACE_NUMBITS_RANGE_LOW), log(LACE_NUMBITS_RANGE_HIGH), 1);
compute_lace_numbits_embedding(numbits_embedded + LACE_NUMBITS_EMBEDDING_DIM, numbits[1], LACE_NUMBITS_EMBEDDING_DIM,
log(LACE_NUMBITS_RANGE_LOW), log(LACE_NUMBITS_RANGE_HIGH), 1);
/* scaling and dimensionality reduction */
for (i_subframe = 0; i_subframe < 4; i_subframe ++)
{
OPUS_COPY(input_buffer, features + i_subframe * LACE_NUM_FEATURES, LACE_NUM_FEATURES);
OPUS_COPY(input_buffer + LACE_NUM_FEATURES, hLACE->layers.lace_pitch_embedding.float_weights + periods[i_subframe] * LACE_PITCH_EMBEDDING_DIM, LACE_PITCH_EMBEDDING_DIM);
OPUS_COPY(input_buffer + LACE_NUM_FEATURES + LACE_PITCH_EMBEDDING_DIM, numbits_embedded, 2 * LACE_NUMBITS_EMBEDDING_DIM);
compute_generic_conv1d(
&hLACE->layers.lace_fnet_conv1,
output_buffer + i_subframe * LACE_HIDDEN_FEATURE_DIM,
NULL,
input_buffer,
LACE_NUM_FEATURES + LACE_PITCH_EMBEDDING_DIM + 2 * LACE_NUMBITS_EMBEDDING_DIM,
ACTIVATION_TANH,
arch);
}
/* subframe accumulation */
OPUS_COPY(input_buffer, output_buffer, 4 * LACE_HIDDEN_FEATURE_DIM);
compute_generic_conv1d(
&hLACE->layers.lace_fnet_conv2,
output_buffer,
state->feature_net_conv2_state,
input_buffer,
4 * LACE_HIDDEN_FEATURE_DIM,
ACTIVATION_TANH,
arch
);
/* tconv upsampling */
OPUS_COPY(input_buffer, output_buffer, 4 * LACE_COND_DIM);
compute_generic_dense(
&hLACE->layers.lace_fnet_tconv,
output_buffer,
input_buffer,
ACTIVATION_TANH,
arch
);
/* GRU */
OPUS_COPY(input_buffer, output_buffer, 4 * LACE_COND_DIM);
for (i_subframe = 0; i_subframe < 4; i_subframe++)
{
compute_generic_gru(
&hLACE->layers.lace_fnet_gru_input,
&hLACE->layers.lace_fnet_gru_recurrent,
state->feature_net_gru_state,
input_buffer + i_subframe * LACE_COND_DIM,
arch
);
OPUS_COPY(output + i_subframe * LACE_COND_DIM, state->feature_net_gru_state, LACE_COND_DIM);
}
}
static void lace_process_20ms_frame(
LACE* hLACE,
LACEState *state,
float *x_out,
const float *x_in,
const float *features,
const float *numbits,
const int *periods,
int arch
)
{
float feature_buffer[4 * LACE_COND_DIM];
float output_buffer[4 * LACE_FRAME_SIZE];
int i_subframe, i_sample;
#ifdef DEBUG_LACE
static FILE *f_features=NULL, *f_encfeatures=NULL, *f_xin=NULL, *f_xpreemph=NULL, *f_postcf1=NULL;
static FILE *f_postcf2=NULL, *f_postaf1=NULL, *f_xdeemph, *f_numbits, *f_periods;
FINIT(f_features, "debug/c_features.f32", "wb");
FINIT(f_encfeatures, "debug/c_encoded_features.f32", "wb");
FINIT(f_xin, "debug/c_x_in.f32", "wb");
FINIT(f_xpreemph, "debug/c_xpreemph.f32", "wb");
FINIT(f_xdeemph, "debug/c_xdeemph.f32", "wb");
FINIT(f_postcf1, "debug/c_post_cf1.f32", "wb");
FINIT(f_postcf2, "debug/c_post_cf2.f32", "wb");
FINIT(f_postaf1, "debug/c_post_af1.f32", "wb");
FINIT(f_numbits, "debug/c_numbits.f32", "wb");
FINIT(f_periods, "debug/c_periods.s32", "wb");
fwrite(x_in, sizeof(*x_in), 4 * LACE_FRAME_SIZE, f_xin);
fwrite(numbits, sizeof(*numbits), 2, f_numbits);
fwrite(periods, sizeof(*periods), 4, f_periods);
#endif
/* pre-emphasis */
for (i_sample = 0; i_sample < 4 * LACE_FRAME_SIZE; i_sample ++)
{
output_buffer[i_sample] = x_in[i_sample] - LACE_PREEMPH * state->preemph_mem;
state->preemph_mem = x_in[i_sample];
}
/* run feature encoder */
lace_feature_net(hLACE, state, feature_buffer, features, numbits, periods, arch);
#ifdef DEBUG_LACE
fwrite(features, sizeof(*features), 4 * LACE_NUM_FEATURES, f_features);
fwrite(feature_buffer, sizeof(*feature_buffer), 4 * LACE_COND_DIM, f_encfeatures);
fwrite(output_buffer, sizeof(float), 4 * LACE_FRAME_SIZE, f_xpreemph);
#endif
/* 1st comb filtering stage */
for (i_subframe = 0; i_subframe < 4; i_subframe++)
{
adacomb_process_frame(
&state->cf1_state,
output_buffer + i_subframe * LACE_FRAME_SIZE,
output_buffer + i_subframe * LACE_FRAME_SIZE,
feature_buffer + i_subframe * LACE_COND_DIM,
&hLACE->layers.lace_cf1_kernel,
&hLACE->layers.lace_cf1_gain,
&hLACE->layers.lace_cf1_global_gain,
periods[i_subframe],
LACE_COND_DIM,
LACE_FRAME_SIZE,
LACE_OVERLAP_SIZE,
LACE_CF1_KERNEL_SIZE,
LACE_CF1_LEFT_PADDING,
LACE_CF1_FILTER_GAIN_A,
LACE_CF1_FILTER_GAIN_B,
LACE_CF1_LOG_GAIN_LIMIT,
hLACE->window,
arch);
}
#ifdef DEBUG_LACE
fwrite(output_buffer, sizeof(float), 4 * LACE_FRAME_SIZE, f_postcf1);
#endif
/* 2nd comb filtering stage */
for (i_subframe = 0; i_subframe < 4; i_subframe++)
{
adacomb_process_frame(
&state->cf2_state,
output_buffer + i_subframe * LACE_FRAME_SIZE,
output_buffer + i_subframe * LACE_FRAME_SIZE,
feature_buffer + i_subframe * LACE_COND_DIM,
&hLACE->layers.lace_cf2_kernel,
&hLACE->layers.lace_cf2_gain,
&hLACE->layers.lace_cf2_global_gain,
periods[i_subframe],
LACE_COND_DIM,
LACE_FRAME_SIZE,
LACE_OVERLAP_SIZE,
LACE_CF2_KERNEL_SIZE,
LACE_CF2_LEFT_PADDING,
LACE_CF2_FILTER_GAIN_A,
LACE_CF2_FILTER_GAIN_B,
LACE_CF2_LOG_GAIN_LIMIT,
hLACE->window,
arch);
}
#ifdef DEBUG_LACE
fwrite(output_buffer, sizeof(float), 4 * LACE_FRAME_SIZE, f_postcf2);
#endif
/* final adaptive filtering stage */
for (i_subframe = 0; i_subframe < 4; i_subframe++)
{
adaconv_process_frame(
&state->af1_state,
output_buffer + i_subframe * LACE_FRAME_SIZE,
output_buffer + i_subframe * LACE_FRAME_SIZE,
feature_buffer + i_subframe * LACE_COND_DIM,
&hLACE->layers.lace_af1_kernel,
&hLACE->layers.lace_af1_gain,
LACE_COND_DIM,
LACE_FRAME_SIZE,
LACE_OVERLAP_SIZE,
LACE_AF1_IN_CHANNELS,
LACE_AF1_OUT_CHANNELS,
LACE_AF1_KERNEL_SIZE,
LACE_AF1_LEFT_PADDING,
LACE_AF1_FILTER_GAIN_A,
LACE_AF1_FILTER_GAIN_B,
LACE_AF1_SHAPE_GAIN,
hLACE->window,
arch);
}
#ifdef DEBUG_LACE
fwrite(output_buffer, sizeof(float), 4 * LACE_FRAME_SIZE, f_postaf1);
#endif
/* de-emphasis */
for (i_sample = 0; i_sample < 4 * LACE_FRAME_SIZE; i_sample ++)
{
x_out[i_sample] = output_buffer[i_sample] + LACE_PREEMPH * state->deemph_mem;
state->deemph_mem = x_out[i_sample];
}
#ifdef DEBUG_LACE
fwrite(x_out, sizeof(float), 4 * LACE_FRAME_SIZE, f_xdeemph);
#endif
}
#endif /* #ifndef DISABLE_LACE */
/* NoLACE */
#ifndef DISABLE_NOLACE
static void compute_nolace_numbits_embedding(float *emb, float numbits, int dim, float min_val, float max_val, int logscale)
{
float x;
(void) dim;
numbits = logscale ? log(numbits) : numbits;
x = CLIP(numbits, min_val, max_val) - (max_val + min_val) / 2;
emb[0] = sin(x * NOLACE_NUMBITS_SCALE_0 - 0.5f);
emb[1] = sin(x * NOLACE_NUMBITS_SCALE_1 - 0.5f);
emb[2] = sin(x * NOLACE_NUMBITS_SCALE_2 - 0.5f);
emb[3] = sin(x * NOLACE_NUMBITS_SCALE_3 - 0.5f);
emb[4] = sin(x * NOLACE_NUMBITS_SCALE_4 - 0.5f);
emb[5] = sin(x * NOLACE_NUMBITS_SCALE_5 - 0.5f);
emb[6] = sin(x * NOLACE_NUMBITS_SCALE_6 - 0.5f);
emb[7] = sin(x * NOLACE_NUMBITS_SCALE_7 - 0.5f);
}
static int init_nolace(NoLACE *hNoLACE, const WeightArray *weights)
{
int ret = 0;
OPUS_CLEAR(hNoLACE, 1);
celt_assert(weights != NULL);
ret = init_nolacelayers(&hNoLACE->layers, weights);
compute_overlap_window(hNoLACE->window, NOLACE_OVERLAP_SIZE);
return ret;
}
static void reset_nolace_state(NoLACEState *state)
{
OPUS_CLEAR(state, 1);
init_adacomb_state(&state->cf1_state);
init_adacomb_state(&state->cf2_state);
init_adaconv_state(&state->af1_state);
init_adaconv_state(&state->af2_state);
init_adaconv_state(&state->af3_state);
init_adaconv_state(&state->af4_state);
init_adashape_state(&state->tdshape1_state);
init_adashape_state(&state->tdshape2_state);
init_adashape_state(&state->tdshape3_state);
}
static void nolace_feature_net(
NoLACE *hNoLACE,
NoLACEState *state,
float *output,
const float *features,
const float *numbits,
const int *periods,
int arch
)
{
float input_buffer[4 * IMAX(NOLACE_COND_DIM, NOLACE_HIDDEN_FEATURE_DIM)];
float output_buffer[4 * IMAX(NOLACE_COND_DIM, NOLACE_HIDDEN_FEATURE_DIM)];
float numbits_embedded[2 * NOLACE_NUMBITS_EMBEDDING_DIM];
int i_subframe;
compute_nolace_numbits_embedding(numbits_embedded, numbits[0], NOLACE_NUMBITS_EMBEDDING_DIM,
log(NOLACE_NUMBITS_RANGE_LOW), log(NOLACE_NUMBITS_RANGE_HIGH), 1);
compute_nolace_numbits_embedding(numbits_embedded + NOLACE_NUMBITS_EMBEDDING_DIM, numbits[1], NOLACE_NUMBITS_EMBEDDING_DIM,
log(NOLACE_NUMBITS_RANGE_LOW), log(NOLACE_NUMBITS_RANGE_HIGH), 1);
/* scaling and dimensionality reduction */
for (i_subframe = 0; i_subframe < 4; i_subframe ++)
{
OPUS_COPY(input_buffer, features + i_subframe * NOLACE_NUM_FEATURES, NOLACE_NUM_FEATURES);
OPUS_COPY(input_buffer + NOLACE_NUM_FEATURES, hNoLACE->layers.nolace_pitch_embedding.float_weights + periods[i_subframe] * NOLACE_PITCH_EMBEDDING_DIM, NOLACE_PITCH_EMBEDDING_DIM);
OPUS_COPY(input_buffer + NOLACE_NUM_FEATURES + NOLACE_PITCH_EMBEDDING_DIM, numbits_embedded, 2 * NOLACE_NUMBITS_EMBEDDING_DIM);
compute_generic_conv1d(
&hNoLACE->layers.nolace_fnet_conv1,
output_buffer + i_subframe * NOLACE_HIDDEN_FEATURE_DIM,
NULL,
input_buffer,
NOLACE_NUM_FEATURES + NOLACE_PITCH_EMBEDDING_DIM + 2 * NOLACE_NUMBITS_EMBEDDING_DIM,
ACTIVATION_TANH,
arch);
}
/* subframe accumulation */
OPUS_COPY(input_buffer, output_buffer, 4 * NOLACE_HIDDEN_FEATURE_DIM);
compute_generic_conv1d(
&hNoLACE->layers.nolace_fnet_conv2,
output_buffer,
state->feature_net_conv2_state,
input_buffer,
4 * NOLACE_HIDDEN_FEATURE_DIM,
ACTIVATION_TANH,
arch
);
/* tconv upsampling */
OPUS_COPY(input_buffer, output_buffer, 4 * NOLACE_COND_DIM);
compute_generic_dense(
&hNoLACE->layers.nolace_fnet_tconv,
output_buffer,
input_buffer,
ACTIVATION_TANH,
arch
);
/* GRU */
OPUS_COPY(input_buffer, output_buffer, 4 * NOLACE_COND_DIM);
for (i_subframe = 0; i_subframe < 4; i_subframe++)
{
compute_generic_gru(
&hNoLACE->layers.nolace_fnet_gru_input,
&hNoLACE->layers.nolace_fnet_gru_recurrent,
state->feature_net_gru_state,
input_buffer + i_subframe * NOLACE_COND_DIM,
arch
);
OPUS_COPY(output + i_subframe * NOLACE_COND_DIM, state->feature_net_gru_state, NOLACE_COND_DIM);
}
}
static void nolace_process_20ms_frame(
NoLACE* hNoLACE,
NoLACEState *state,
float *x_out,
const float *x_in,
const float *features,
const float *numbits,
const int *periods,
int arch
)
{
float feature_buffer[4 * NOLACE_COND_DIM];
float feature_transform_buffer[4 * NOLACE_COND_DIM];
float x_buffer1[8 * NOLACE_FRAME_SIZE];
float x_buffer2[8 * NOLACE_FRAME_SIZE];
int i_subframe, i_sample;
NOLACELayers *layers = &hNoLACE->layers;
#ifdef DEBUG_NOLACE
static FILE *f_features=NULL, *f_encfeatures=NULL, *f_xin=NULL, *f_xpreemph=NULL, *f_postcf1=NULL;
static FILE *f_postcf2=NULL, *f_postaf1=NULL, *f_xdeemph, *f_numbits, *f_periods;
static FILE *f_ffpostcf1, *f_fpostcf2, *f_fpostaf1;
FINIT(f_features, "debug/c_features.f32", "wb");
FINIT(f_encfeatures, "debug/c_encoded_features.f32", "wb");
FINIT(f_xin, "debug/c_x_in.f32", "wb");
FINIT(f_xpreemph, "debug/c_xpreemph.f32", "wb");
FINIT(f_xdeemph, "debug/c_xdeemph.f32", "wb");
FINIT(f_postcf1, "debug/c_post_cf1.f32", "wb");
FINIT(f_postcf2, "debug/c_post_cf2.f32", "wb");
FINIT(f_postaf1, "debug/c_post_af1.f32", "wb");
FINIT(f_numbits, "debug/c_numbits.f32", "wb");
FINIT(f_periods, "debug/c_periods.s32", "wb");
fwrite(x_in, sizeof(*x_in), 4 * NOLACE_FRAME_SIZE, f_xin);
fwrite(numbits, sizeof(*numbits), 2, f_numbits);
fwrite(periods, sizeof(*periods), 4, f_periods);
#endif
/* pre-emphasis */
for (i_sample = 0; i_sample < 4 * NOLACE_FRAME_SIZE; i_sample ++)
{
x_buffer1[i_sample] = x_in[i_sample] - NOLACE_PREEMPH * state->preemph_mem;
state->preemph_mem = x_in[i_sample];
}
/* run feature encoder */
nolace_feature_net(hNoLACE, state, feature_buffer, features, numbits, periods, arch);
#ifdef DEBUG_NOLACE
fwrite(features, sizeof(*features), 4 * NOLACE_NUM_FEATURES, f_features);
fwrite(feature_buffer, sizeof(*feature_buffer), 4 * NOLACE_COND_DIM, f_encfeatures);
fwrite(output_buffer, sizeof(float), 4 * NOLACE_FRAME_SIZE, f_xpreemph);
#endif
/* 1st comb filtering stage */
for (i_subframe = 0; i_subframe < 4; i_subframe++)
{
/* modifies signal in place */
adacomb_process_frame(
&state->cf1_state,
x_buffer1 + i_subframe * NOLACE_FRAME_SIZE,
x_buffer1 + i_subframe * NOLACE_FRAME_SIZE,
feature_buffer + i_subframe * NOLACE_COND_DIM,
&hNoLACE->layers.nolace_cf1_kernel,
&hNoLACE->layers.nolace_cf1_gain,
&hNoLACE->layers.nolace_cf1_global_gain,
periods[i_subframe],
NOLACE_COND_DIM,
NOLACE_FRAME_SIZE,
NOLACE_OVERLAP_SIZE,
NOLACE_CF1_KERNEL_SIZE,
NOLACE_CF1_LEFT_PADDING,
NOLACE_CF1_FILTER_GAIN_A,
NOLACE_CF1_FILTER_GAIN_B,
NOLACE_CF1_LOG_GAIN_LIMIT,
hNoLACE->window,
arch);
compute_generic_conv1d(
&layers->nolace_post_cf1,
feature_transform_buffer + i_subframe * NOLACE_COND_DIM,
state->post_cf1_state,
feature_buffer + i_subframe * NOLACE_COND_DIM,
NOLACE_COND_DIM,
ACTIVATION_TANH,
arch);
}
/* update feature buffer */
OPUS_COPY(feature_buffer, feature_transform_buffer, 4 * NOLACE_COND_DIM);
#ifdef DEBUG_NOLACE
fwrite(x_buffer1, sizeof(float), 4 * NOLACE_FRAME_SIZE, f_postcf1);
#endif
/* 2nd comb filtering stage */
for (i_subframe = 0; i_subframe < 4; i_subframe++)
{
/* modifies signal in place */
adacomb_process_frame(
&state->cf2_state,
x_buffer1 + i_subframe * NOLACE_FRAME_SIZE,
x_buffer1 + i_subframe * NOLACE_FRAME_SIZE,
feature_buffer + i_subframe * NOLACE_COND_DIM,
&hNoLACE->layers.nolace_cf2_kernel,
&hNoLACE->layers.nolace_cf2_gain,
&hNoLACE->layers.nolace_cf2_global_gain,
periods[i_subframe],
NOLACE_COND_DIM,
NOLACE_FRAME_SIZE,
NOLACE_OVERLAP_SIZE,
NOLACE_CF2_KERNEL_SIZE,
NOLACE_CF2_LEFT_PADDING,
NOLACE_CF2_FILTER_GAIN_A,
NOLACE_CF2_FILTER_GAIN_B,
NOLACE_CF2_LOG_GAIN_LIMIT,
hNoLACE->window,
arch);
compute_generic_conv1d(
&layers->nolace_post_cf2,
feature_transform_buffer + i_subframe * NOLACE_COND_DIM,
state->post_cf2_state,
feature_buffer + i_subframe * NOLACE_COND_DIM,
NOLACE_COND_DIM,
ACTIVATION_TANH,
arch);
}
/* update feature buffer */
OPUS_COPY(feature_buffer, feature_transform_buffer, 4 * NOLACE_COND_DIM);
#ifdef DEBUG_NOLACE
fwrite(x_buffer1, sizeof(float), 4 * NOLACE_FRAME_SIZE, f_postcf2);
#endif
/* final adaptive filtering stage */
for (i_subframe = 0; i_subframe < 4; i_subframe++)
{
adaconv_process_frame(
&state->af1_state,
x_buffer2 + i_subframe * NOLACE_FRAME_SIZE * NOLACE_AF1_OUT_CHANNELS,
x_buffer1 + i_subframe * NOLACE_FRAME_SIZE,
feature_buffer + i_subframe * NOLACE_COND_DIM,
&hNoLACE->layers.nolace_af1_kernel,
&hNoLACE->layers.nolace_af1_gain,
NOLACE_COND_DIM,
NOLACE_FRAME_SIZE,
NOLACE_OVERLAP_SIZE,
NOLACE_AF1_IN_CHANNELS,
NOLACE_AF1_OUT_CHANNELS,
NOLACE_AF1_KERNEL_SIZE,
NOLACE_AF1_LEFT_PADDING,
NOLACE_AF1_FILTER_GAIN_A,
NOLACE_AF1_FILTER_GAIN_B,
NOLACE_AF1_SHAPE_GAIN,
hNoLACE->window,
arch);
compute_generic_conv1d(
&layers->nolace_post_af1,
feature_transform_buffer + i_subframe * NOLACE_COND_DIM,
state->post_af1_state,
feature_buffer + i_subframe * NOLACE_COND_DIM,
NOLACE_COND_DIM,
ACTIVATION_TANH,
arch);
}
/* update feature buffer */
OPUS_COPY(feature_buffer, feature_transform_buffer, 4 * NOLACE_COND_DIM);
#ifdef DEBUG_NOLACE
fwrite(x_buffer2, sizeof(float), 4 * NOLACE_FRAME_SIZE * NOLACE_AF1_OUT_CHANNELS, f_postaf1);
#endif
/* first shape-mix round */
for (i_subframe = 0; i_subframe < 4; i_subframe++)
{
celt_assert(NOLACE_AF1_OUT_CHANNELS == 2);
/* modifies second channel in place */
adashape_process_frame(
&state->tdshape1_state,
x_buffer2 + i_subframe * NOLACE_AF1_OUT_CHANNELS * NOLACE_FRAME_SIZE + NOLACE_FRAME_SIZE,
x_buffer2 + i_subframe * NOLACE_AF1_OUT_CHANNELS * NOLACE_FRAME_SIZE + NOLACE_FRAME_SIZE,
feature_buffer + i_subframe * NOLACE_COND_DIM,
&layers->nolace_tdshape1_alpha1_f,
&layers->nolace_tdshape1_alpha1_t,
&layers->nolace_tdshape1_alpha2,
NOLACE_TDSHAPE1_FEATURE_DIM,
NOLACE_TDSHAPE1_FRAME_SIZE,
NOLACE_TDSHAPE1_AVG_POOL_K,
1,
arch
);
adaconv_process_frame(
&state->af2_state,
x_buffer1 + i_subframe * NOLACE_FRAME_SIZE * NOLACE_AF2_OUT_CHANNELS,
x_buffer2 + i_subframe * NOLACE_FRAME_SIZE * NOLACE_AF2_IN_CHANNELS,
feature_buffer + i_subframe * NOLACE_COND_DIM,
&hNoLACE->layers.nolace_af2_kernel,
&hNoLACE->layers.nolace_af2_gain,
NOLACE_COND_DIM,
NOLACE_FRAME_SIZE,
NOLACE_OVERLAP_SIZE,
NOLACE_AF2_IN_CHANNELS,
NOLACE_AF2_OUT_CHANNELS,
NOLACE_AF2_KERNEL_SIZE,
NOLACE_AF2_LEFT_PADDING,
NOLACE_AF2_FILTER_GAIN_A,
NOLACE_AF2_FILTER_GAIN_B,
NOLACE_AF2_SHAPE_GAIN,
hNoLACE->window,
arch);
compute_generic_conv1d(
&layers->nolace_post_af2,
feature_transform_buffer + i_subframe * NOLACE_COND_DIM,
state->post_af2_state,
feature_buffer + i_subframe * NOLACE_COND_DIM,
NOLACE_COND_DIM,
ACTIVATION_TANH,
arch);
}
/* update feature buffer */
OPUS_COPY(feature_buffer, feature_transform_buffer, 4 * NOLACE_COND_DIM);
#ifdef DEBUG_NOLACE
fwrite(x_buffer1, sizeof(float), 4 * NOLACE_FRAME_SIZE * NOLACE_AF2_OUT_CHANNELS, f_postaf2);
#endif
/* second shape-mix round */
for (i_subframe = 0; i_subframe < 4; i_subframe++)
{
celt_assert(NOLACE_AF2_OUT_CHANNELS == 2);
/* modifies second channel in place */
adashape_process_frame(
&state->tdshape2_state,
x_buffer1 + i_subframe * NOLACE_AF2_OUT_CHANNELS * NOLACE_FRAME_SIZE + NOLACE_FRAME_SIZE,
x_buffer1 + i_subframe * NOLACE_AF2_OUT_CHANNELS * NOLACE_FRAME_SIZE + NOLACE_FRAME_SIZE,
feature_buffer + i_subframe * NOLACE_COND_DIM,
&layers->nolace_tdshape2_alpha1_f,
&layers->nolace_tdshape2_alpha1_t,
&layers->nolace_tdshape2_alpha2,
NOLACE_TDSHAPE2_FEATURE_DIM,
NOLACE_TDSHAPE2_FRAME_SIZE,
NOLACE_TDSHAPE2_AVG_POOL_K,
1,
arch
);
adaconv_process_frame(
&state->af3_state,
x_buffer2 + i_subframe * NOLACE_FRAME_SIZE * NOLACE_AF3_OUT_CHANNELS,
x_buffer1 + i_subframe * NOLACE_FRAME_SIZE * NOLACE_AF3_IN_CHANNELS,
feature_buffer + i_subframe * NOLACE_COND_DIM,
&hNoLACE->layers.nolace_af3_kernel,
&hNoLACE->layers.nolace_af3_gain,
NOLACE_COND_DIM,
NOLACE_FRAME_SIZE,
NOLACE_OVERLAP_SIZE,
NOLACE_AF3_IN_CHANNELS,
NOLACE_AF3_OUT_CHANNELS,
NOLACE_AF3_KERNEL_SIZE,
NOLACE_AF3_LEFT_PADDING,
NOLACE_AF3_FILTER_GAIN_A,
NOLACE_AF3_FILTER_GAIN_B,
NOLACE_AF3_SHAPE_GAIN,
hNoLACE->window,
arch);
compute_generic_conv1d(
&layers->nolace_post_af3,
feature_transform_buffer + i_subframe * NOLACE_COND_DIM,
state->post_af3_state,
feature_buffer + i_subframe * NOLACE_COND_DIM,
NOLACE_COND_DIM,
ACTIVATION_TANH,
arch);
}
/* update feature buffer */
OPUS_COPY(feature_buffer, feature_transform_buffer, 4 * NOLACE_COND_DIM);
/* third shape-mix round */
for (i_subframe = 0; i_subframe < 4; i_subframe++)
{
celt_assert(NOLACE_AF3_OUT_CHANNELS == 2);
/* modifies second channel in place */
adashape_process_frame(
&state->tdshape3_state,
x_buffer2 + i_subframe * NOLACE_AF3_OUT_CHANNELS * NOLACE_FRAME_SIZE + NOLACE_FRAME_SIZE,
x_buffer2 + i_subframe * NOLACE_AF3_OUT_CHANNELS * NOLACE_FRAME_SIZE + NOLACE_FRAME_SIZE,
feature_buffer + i_subframe * NOLACE_COND_DIM,
&layers->nolace_tdshape3_alpha1_f,
&layers->nolace_tdshape3_alpha1_t,
&layers->nolace_tdshape3_alpha2,
NOLACE_TDSHAPE3_FEATURE_DIM,
NOLACE_TDSHAPE3_FRAME_SIZE,
NOLACE_TDSHAPE3_AVG_POOL_K,
1,
arch
);
adaconv_process_frame(
&state->af4_state,
x_buffer1 + i_subframe * NOLACE_FRAME_SIZE * NOLACE_AF4_OUT_CHANNELS,
x_buffer2 + i_subframe * NOLACE_FRAME_SIZE * NOLACE_AF4_IN_CHANNELS,
feature_buffer + i_subframe * NOLACE_COND_DIM,
&hNoLACE->layers.nolace_af4_kernel,
&hNoLACE->layers.nolace_af4_gain,
NOLACE_COND_DIM,
NOLACE_FRAME_SIZE,
NOLACE_OVERLAP_SIZE,
NOLACE_AF4_IN_CHANNELS,
NOLACE_AF4_OUT_CHANNELS,
NOLACE_AF4_KERNEL_SIZE,
NOLACE_AF4_LEFT_PADDING,
NOLACE_AF4_FILTER_GAIN_A,
NOLACE_AF4_FILTER_GAIN_B,
NOLACE_AF4_SHAPE_GAIN,
hNoLACE->window,
arch);
}
/* de-emphasis */
for (i_sample = 0; i_sample < 4 * NOLACE_FRAME_SIZE; i_sample ++)
{
x_out[i_sample] = x_buffer1[i_sample] + NOLACE_PREEMPH * state->deemph_mem;
state->deemph_mem = x_out[i_sample];
}
#ifdef DEBUG_NOLACE
fwrite(x_out, sizeof(float), 4 * NOLACE_FRAME_SIZE, f_xdeemph);
#endif
}
#endif /* #ifndef DISABLE_NOLACE */
#ifdef ENABLE_OSCE_BWE
#ifndef DISABLE_BBWENET
static void bbwe_feature_net(
BBWENet *hBBWENET,
BBWENetState *state,
float *output,
const float *features,
int num_frames,
int arch
)
{
float input_buffer[4 * BBWENET_FNET_GRU_STATE_SIZE];
float output_buffer[4 * BBWENET_FNET_GRU_STATE_SIZE];
int i_subframe;
int i_frame;
#ifdef DEBUG_BBWENET
static FILE *f_features=NULL, *f_conv1=NULL, *f_conv2=NULL, *f_tconv=NULL, *f_gru=NULL;
FINIT(f_features, "debug/bbwenet_features.f32", "wb");
FINIT(f_conv1, "debug/bbwenet_conv1.f32", "wb");
FINIT(f_conv2, "debug/bbwenet_conv2.f32", "wb");
FINIT(f_tconv, "debug/bbwenet_tconv.f32", "wb");
FINIT(f_gru, "debug/bbwenet_gru.f32", "wb");
fwrite(features, sizeof(*features), num_frames * BBWENET_FEATURE_DIM, f_features);
#endif
/* adjust buffer sizes if any of this breaks */
celt_assert(BBWENET_FNET_GRU_STATE_SIZE == BBWENET_FNET_TCONV_OUT_CHANNELS);
celt_assert(BBWENET_FNET_TCONV_OUT_CHANNELS == BBWENET_FNET_CONV2_OUT_SIZE);
celt_assert(BBWENET_FNET_CONV2_OUT_SIZE == BBWENET_FNET_CONV1_OUT_SIZE);
/* first conv layer */
for (i_frame = 0; i_frame < num_frames; i_frame++)
{
compute_generic_conv1d(
&hBBWENET->layers.bbwenet_fnet_conv1,
output_buffer + i_frame * BBWENET_FNET_CONV1_OUT_SIZE,
state->feature_net_conv1_state,
features + i_frame * BBWENET_FEATURE_DIM,
BBWENET_FEATURE_DIM,
ACTIVATION_TANH,
arch
);
#ifdef DEBUG_BBWENET
fwrite(output_buffer + i_frame * BBWENET_FNET_CONV1_OUT_SIZE, sizeof(float), BBWENET_FNET_CONV1_OUT_SIZE, f_conv1);
#endif
}
OPUS_COPY(input_buffer, output_buffer, num_frames * BBWENET_FNET_CONV1_OUT_SIZE);
/* second conv layer */
for (i_frame = 0; i_frame < num_frames; i_frame++)
{
compute_generic_conv1d(
&hBBWENET->layers.bbwenet_fnet_conv2,
output_buffer + i_frame * BBWENET_FNET_CONV2_OUT_SIZE,
state->feature_net_conv2_state,
input_buffer + i_frame * BBWENET_FNET_CONV1_OUT_SIZE,
BBWENET_FNET_CONV1_OUT_SIZE,
ACTIVATION_TANH,
arch
);
#ifdef DEBUG_BBWENET
fwrite(output_buffer + i_frame * BBWENET_FNET_CONV2_OUT_SIZE, sizeof(float), BBWENET_FNET_CONV2_OUT_SIZE, f_conv2);
#endif
}
OPUS_COPY(input_buffer, output_buffer, num_frames * BBWENET_FNET_CONV2_OUT_SIZE);
/* tconv upsampling*/
for (i_frame = 0; i_frame < num_frames; i_frame++)
{
compute_generic_dense(
&hBBWENET->layers.bbwenet_fnet_tconv,
output_buffer + i_frame * BBWENET_FNET_TCONV_OUT_CHANNELS * BBWENET_FNET_TCONV_STRIDE,
input_buffer + i_frame * BBWENET_FNET_CONV2_OUT_SIZE,
ACTIVATION_TANH,
arch
);
#ifdef DEBUG_BBWENET
fwrite(output_buffer + i_frame * BBWENET_FNET_TCONV_OUT_CHANNELS * BBWENET_FNET_TCONV_STRIDE, sizeof(float), BBWENET_FNET_TCONV_OUT_CHANNELS * BBWENET_FNET_TCONV_STRIDE, f_tconv);
#endif
}
OPUS_COPY(input_buffer, output_buffer, num_frames * BBWENET_FNET_TCONV_OUT_CHANNELS * BBWENET_FNET_TCONV_STRIDE);
/* GRU */
celt_assert(BBWENET_FNET_TCONV_STRIDE == 2)
for (i_subframe = 0; i_subframe < BBWENET_FNET_TCONV_STRIDE * num_frames; i_subframe ++)
{
compute_generic_gru(
&hBBWENET->layers.bbwenet_fnet_gru_input,
&hBBWENET->layers.bbwenet_fnet_gru_recurrent,
state->feature_net_gru_state,
input_buffer + i_subframe * BBWENET_FNET_TCONV_OUT_CHANNELS,
arch
);
#ifdef DEBUG_BBWENET
fwrite(state->feature_net_gru_state, sizeof(float), BBWENET_FNET_GRU_STATE_SIZE, f_gru);
#endif
OPUS_COPY(output + i_subframe * BBWENET_FNET_GRU_STATE_SIZE, state->feature_net_gru_state, BBWENET_FNET_GRU_STATE_SIZE);
}
}
static float hq_2x_even[3] = {0.026641845703125, 0.228668212890625, -0.4036407470703125};
static float hq_2x_odd[3] = {0.104583740234375, 0.3932037353515625, -0.152496337890625};
static float frac_01_24[8] = {
0.00576782, -0.01831055, 0.01882935, 0.9328308,
0.09143066, -0.04196167, 0.01296997, -0.00140381
};
static float frac_17_24[8] = {
-3.14331055e-03, 2.73437500e-02, -1.06414795e-01, 3.64685059e-01,
8.03863525e-01, -1.02233887e-01, 1.61437988e-02, -1.22070312e-04
};
static float frac_09_24[8] = {
-0.00146484, 0.02313232, -0.12072754, 0.7315979,
0.4621277, -0.12075806, 0.0295105 , -0.00326538
};
static void apply_valin_activation(float *x, int len)
{
int i;
for (i = 0; i < len; i++)
{
x[i] *= sin(log(fabs(x[i]) + 1e-6f));
}
}
#define DELAY_SAMPLES 8 /* ToDo: this probably should be 7, bug in python code? */
static void interpol_3_2(resamp_state *state, float *x_out, const float *x_in, int num_samples)
{
int i_sample, i_out = 0;
float buffer[8 * BBWENET_FRAME_SIZE16 + DELAY_SAMPLES];
celt_assert(num_samples > 1);
celt_assert(num_samples < 8 * BBWENET_FRAME_SIZE16);
celt_assert(num_samples % 2 == 0);
OPUS_COPY(buffer, state->interpol_buffer, DELAY_SAMPLES);
OPUS_COPY(buffer + DELAY_SAMPLES, x_in, num_samples);
for (i_sample = 0; i_sample < num_samples; i_sample+=2)
{
x_out[i_out++] = buffer[i_sample + 0] * frac_01_24[0] +
buffer[i_sample + 1] * frac_01_24[1] +
buffer[i_sample + 2] * frac_01_24[2] +
buffer[i_sample + 3] * frac_01_24[3] +
buffer[i_sample + 4] * frac_01_24[4] +
buffer[i_sample + 5] * frac_01_24[5] +
buffer[i_sample + 6] * frac_01_24[6] +
buffer[i_sample + 7] * frac_01_24[7];
x_out[i_out++] = buffer[i_sample + 0] * frac_17_24[0] +
buffer[i_sample + 1] * frac_17_24[1] +
buffer[i_sample + 2] * frac_17_24[2] +
buffer[i_sample + 3] * frac_17_24[3] +
buffer[i_sample + 4] * frac_17_24[4] +
buffer[i_sample + 5] * frac_17_24[5] +
buffer[i_sample + 6] * frac_17_24[6] +
buffer[i_sample + 7] * frac_17_24[7];
x_out[i_out++] = buffer[i_sample + 1] * frac_09_24[0] +
buffer[i_sample + 2] * frac_09_24[1] +
buffer[i_sample + 3] * frac_09_24[2] +
buffer[i_sample + 4] * frac_09_24[3] +
buffer[i_sample + 5] * frac_09_24[4] +
buffer[i_sample + 6] * frac_09_24[5] +
buffer[i_sample + 7] * frac_09_24[6] +
buffer[i_sample + 8] * frac_09_24[7];
}
/* copy last samples to buffer */
OPUS_COPY(state->interpol_buffer, buffer + num_samples, DELAY_SAMPLES);
}
static void upsamp_2x(resamp_state *state, float *x_out, const float *x_in, int num_samples)
{
float buffer [4 * BBWENET_FRAME_SIZE16];
float *S_even = state->upsamp_buffer[0];
float *S_odd = state->upsamp_buffer[1];
int k;
float x, X, Y, tmp1, tmp2, tmp3;
celt_assert(num_samples > 1);
celt_assert(num_samples < 4 * BBWENET_FRAME_SIZE16);
OPUS_COPY(buffer, x_in, num_samples);
for (k = 0; k < num_samples; k++)
{
x = buffer[k];
/* even sample, first pass, */
Y = x - S_even[0];
X = Y * hq_2x_even[0];
tmp1 = S_even[0] + X;
S_even[0] = x + X;
/* ...second pass, */
Y = tmp1 - S_even[1];
X = Y * hq_2x_even[1];
tmp2 = S_even[1] + X;
S_even[1] = tmp1 + X;
/* ...third pass */
Y = tmp2 - S_even[2];
X = Y * (1 + hq_2x_even[2]);
tmp3 = S_even[2] + X;
S_even[2] = tmp2 + X;
x_out[2 * k] = tmp3;
/* odd sample, first pass, */
Y = x - S_odd[0];
X = Y * hq_2x_odd[0];
tmp1 = S_odd[0] + X;
S_odd[0] = x + X;
/* ...second pass, */
Y = tmp1 - S_odd[1];
X = Y * hq_2x_odd[1];
tmp2 = S_odd[1] + X;
S_odd[1] = tmp1 + X;
/* ...third pass */
Y = tmp2 - S_odd[2];
X = Y * (1 + hq_2x_odd[2]);
tmp3 = S_odd[2] + X;
S_odd[2] = tmp2 + X;
x_out[2 * k + 1] = tmp3;
}
}
static void bbwenet_process_frames(
BBWENet *hBBWENET,
BBWENetState *state,
float *x_out,
const float *x_in,
const float *features,
int num_frames,
int arch
)
{
float latent_features[4 * BBWENET_COND_DIM];
int i_subframe, num_subframes = 2 * num_frames, i_channel;
float x_buffer1[3 * 3 * 4 * 3*BBWENET_FRAME_SIZE16] = {0}; /* 3x3 channels, 4 subframes, 48 kHz */
float x_buffer2[3 * 3 * 4 * 3*BBWENET_FRAME_SIZE16] = {0};
BBWENETLayers *layers = &hBBWENET->layers;
#ifdef DEBUG_BBWENET
static FILE *f_latent=NULL, *f_xin=NULL, *f_af1_1=NULL, *f_af1_2=NULL, *f_af1_3=NULL;
static FILE *f_up2_1=NULL, *f_up2_2=NULL, *f_up2_3=NULL, *f2_up_shape=NULL, *f2_up_func=NULL;
static FILE *f_af2_1=NULL, *f_af2_2=NULL, *f_af2_3=NULL;
static FILE *f_up15_1=NULL, *f_up15_2=NULL, *f_up15_3=NULL;
static FILE *f_up15_shape=NULL, *f_up15_func=NULL;
static FILE *f_af3_1=NULL;
FINIT(f_latent, "dnn/torch/osce/debugdump/feature_net_gru.f32", "rb");
FINIT(f_xin, "debug/bbwenet_x_in.f32", "wb");
FINIT(f_af1_1, "debug/bbwenet_af1_1.f32", "wb");
FINIT(f_af1_2, "debug/bbwenet_af1_2.f32", "wb");
FINIT(f_af1_3, "debug/bbwenet_af1_3.f32", "wb");
FINIT(f_up2_1, "debug/bbwenet_up2_1.f32", "wb");
FINIT(f_up2_2, "debug/bbwenet_up2_2.f32", "wb");
FINIT(f_up2_3, "debug/bbwenet_up2_3.f32", "wb");
FINIT(f2_up_func, "debug/bbwenet_up2_func.f32", "wb");
FINIT(f2_up_shape, "debug/bbwenet_up2_shape.f32", "wb");
FINIT(f_af2_1, "debug/bbwenet_af2_1.f32", "wb");
FINIT(f_af2_2, "debug/bbwenet_af2_2.f32", "wb");
FINIT(f_af2_3, "debug/bbwenet_af2_3.f32", "wb");
FINIT(f_up15_1, "debug/bbwenet_up15_1.f32", "wb");
FINIT(f_up15_2, "debug/bbwenet_up15_2.f32", "wb");
FINIT(f_up15_3, "debug/bbwenet_up15_3.f32", "wb");
FINIT(f_up15_shape, "debug/bbwenet_up15_shape.f32", "wb");
FINIT(f_up15_func, "debug/bbwenet_up15_func.f32", "wb");
FINIT(f_af3_1, "debug/bbwenet_af3_1.f32", "wb");
fwrite(x_in, sizeof(*x_in), num_subframes * BBWENET_AF1_FRAME_SIZE, f_xin);
#endif
/* feature net */
bbwe_feature_net(hBBWENET, state, latent_features, features, num_frames, arch);
#ifdef DEBUG_BBWENET
if (f_latent != NULL){
fread(latent_features, sizeof(*latent_features), num_subframes * BBWENET_COND_DIM, f_latent);
}
#endif
/* signal net
* first adaptive filtering stage, three output channels */
for (i_subframe = 0; i_subframe < num_subframes; i_subframe++)
{
adaconv_process_frame(
&state->af1_state,
x_buffer1 + i_subframe * BBWENET_AF1_FRAME_SIZE * BBWENET_AF1_OUT_CHANNELS,
x_in + i_subframe * BBWENET_AF1_FRAME_SIZE,
latent_features + i_subframe * BBWENET_COND_DIM,
&layers->bbwenet_af1_kernel,
&layers->bbwenet_af1_gain,
BBWENET_COND_DIM,
BBWENET_AF1_FRAME_SIZE,
BBWENET_AF1_OVERLAP_SIZE,
BBWENET_AF1_IN_CHANNELS,
BBWENET_AF1_OUT_CHANNELS,
BBWENET_AF1_KERNEL_SIZE,
BBWENET_AF1_LEFT_PADDING,
BBWENET_AF1_FILTER_GAIN_A,
BBWENET_AF1_FILTER_GAIN_B,
BBWENET_AF1_SHAPE_GAIN,
hBBWENET->window16,
arch);
#ifdef DEBUG_BBWENET
fwrite(x_buffer1 + i_subframe * BBWENET_AF1_FRAME_SIZE * BBWENET_AF1_OUT_CHANNELS, sizeof(float), BBWENET_AF1_FRAME_SIZE, f_af1_1);
fwrite(x_buffer1 + i_subframe * BBWENET_AF1_FRAME_SIZE * BBWENET_AF1_OUT_CHANNELS + BBWENET_AF1_FRAME_SIZE, sizeof(float), BBWENET_AF1_FRAME_SIZE, f_af1_2);
fwrite(x_buffer1 + i_subframe * BBWENET_AF1_FRAME_SIZE * BBWENET_AF1_OUT_CHANNELS + 2 * BBWENET_AF1_FRAME_SIZE, sizeof(float), BBWENET_AF1_FRAME_SIZE, f_af1_3);
#endif
}
/* 1st round of non-linear extension */
for (i_subframe = 0; i_subframe < num_subframes; i_subframe++)
{
/* 2x upsampling on individual channels */
celt_assert(BBWENET_AF1_OUT_CHANNELS == 3);
celt_assert(2 * BBWENET_AF1_FRAME_SIZE == BBWENET_TDSHAPE1_FRAME_SIZE);
for (i_channel = 0; i_channel < 3; i_channel ++)
{
upsamp_2x(
&state->resampler_state[i_channel],
x_buffer2 + i_subframe * BBWENET_TDSHAPE1_FRAME_SIZE * BBWENET_AF1_OUT_CHANNELS + i_channel * BBWENET_TDSHAPE1_FRAME_SIZE,
x_buffer1 + i_subframe * BBWENET_AF1_FRAME_SIZE * BBWENET_AF1_OUT_CHANNELS + i_channel * BBWENET_AF1_FRAME_SIZE,
BBWENET_AF1_FRAME_SIZE
);
}
#ifdef DEBUG_BBWENET
fwrite(x_buffer2 + i_subframe * BBWENET_AF1_OUT_CHANNELS * BBWENET_TDSHAPE1_FRAME_SIZE, sizeof(float), BBWENET_TDSHAPE1_FRAME_SIZE, f_up2_1);
fwrite(x_buffer2 + i_subframe * BBWENET_AF1_OUT_CHANNELS * BBWENET_TDSHAPE1_FRAME_SIZE + BBWENET_TDSHAPE1_FRAME_SIZE, sizeof(float), BBWENET_TDSHAPE1_FRAME_SIZE, f_up2_2);
fwrite(x_buffer2 + i_subframe * BBWENET_AF1_OUT_CHANNELS * BBWENET_TDSHAPE1_FRAME_SIZE + 2 * BBWENET_TDSHAPE1_FRAME_SIZE, sizeof(float), BBWENET_TDSHAPE1_FRAME_SIZE, f_up2_3);
#endif
/* tdshape on second channel (in place) */
adashape_process_frame(
&state->tdshape1_state,
x_buffer2 + i_subframe * BBWENET_AF1_OUT_CHANNELS * BBWENET_TDSHAPE1_FRAME_SIZE + BBWENET_TDSHAPE1_FRAME_SIZE,
x_buffer2 + i_subframe * BBWENET_AF1_OUT_CHANNELS * BBWENET_TDSHAPE1_FRAME_SIZE + BBWENET_TDSHAPE1_FRAME_SIZE,
latent_features + i_subframe * BBWENET_COND_DIM,
&layers->bbwenet_tdshape1_alpha1_f,
&layers->bbwenet_tdshape1_alpha1_t,
&layers->bbwenet_tdshape1_alpha2,
BBWENET_TDSHAPE1_FEATURE_DIM,
BBWENET_TDSHAPE1_FRAME_SIZE,
BBWENET_TDSHAPE1_AVG_POOL_K,
BBWENET_TDSHAPE1_INTERPOLATE_K,
arch
);
#ifdef DEBUG_BBWENET
fwrite(x_buffer2 + i_subframe * BBWENET_AF1_OUT_CHANNELS * BBWENET_TDSHAPE1_FRAME_SIZE + BBWENET_TDSHAPE1_FRAME_SIZE, sizeof(float), BBWENET_TDSHAPE1_FRAME_SIZE, f2_up_shape);
#endif
/* non-linear activation of third channel (in place)*/
apply_valin_activation(
x_buffer2 + i_subframe * BBWENET_AF1_OUT_CHANNELS * BBWENET_TDSHAPE1_FRAME_SIZE + 2 * BBWENET_TDSHAPE1_FRAME_SIZE,
BBWENET_TDSHAPE1_FRAME_SIZE
);
#ifdef DEBUG_BBWENET
fwrite(x_buffer2 + i_subframe * BBWENET_AF1_OUT_CHANNELS * BBWENET_TDSHAPE1_FRAME_SIZE + 2 * BBWENET_TDSHAPE1_FRAME_SIZE, sizeof(float), BBWENET_TDSHAPE1_FRAME_SIZE, f2_up_func);
#endif
}
/* mixing */
for (i_subframe = 0; i_subframe < num_subframes; i_subframe++)
{
adaconv_process_frame(
&state->af2_state,
x_buffer1 + i_subframe * BBWENET_AF2_FRAME_SIZE * BBWENET_AF2_OUT_CHANNELS,
x_buffer2 + i_subframe * BBWENET_AF2_FRAME_SIZE * BBWENET_AF1_OUT_CHANNELS,
latent_features + i_subframe * BBWENET_COND_DIM,
&layers->bbwenet_af2_kernel,
&layers->bbwenet_af2_gain,
BBWENET_COND_DIM,
BBWENET_AF2_FRAME_SIZE,
BBWENET_AF2_OVERLAP_SIZE,
BBWENET_AF2_IN_CHANNELS,
BBWENET_AF2_OUT_CHANNELS,
BBWENET_AF2_KERNEL_SIZE,
BBWENET_AF2_LEFT_PADDING,
BBWENET_AF2_FILTER_GAIN_A,
BBWENET_AF2_FILTER_GAIN_B,
BBWENET_AF2_SHAPE_GAIN,
hBBWENET->window32,
arch);
#ifdef DEBUG_BBWENET
fwrite(x_buffer1 + i_subframe * BBWENET_AF2_FRAME_SIZE * BBWENET_AF2_OUT_CHANNELS, sizeof(float), BBWENET_AF2_FRAME_SIZE, f_af2_1);
fwrite(x_buffer1 + i_subframe * BBWENET_AF2_FRAME_SIZE * BBWENET_AF2_OUT_CHANNELS + BBWENET_AF2_FRAME_SIZE, sizeof(float), BBWENET_AF2_FRAME_SIZE, f_af2_2);
fwrite(x_buffer1 + i_subframe * BBWENET_AF2_FRAME_SIZE * BBWENET_AF2_OUT_CHANNELS + 2 * BBWENET_AF2_FRAME_SIZE, sizeof(float), BBWENET_AF2_FRAME_SIZE, f_af2_3);
#endif
}
/* second round of extension */
for (i_subframe = 0; i_subframe < num_subframes; i_subframe++)
{
/* 1.5x interpolation on individual channels */
celt_assert(BBWENET_AF2_OUT_CHANNELS == 3);
celt_assert(3 * BBWENET_AF2_FRAME_SIZE == 2 * BBWENET_TDSHAPE2_FRAME_SIZE);
for (i_channel = 0; i_channel < 3; i_channel ++)
{
interpol_3_2(
&state->resampler_state[i_channel],
x_buffer2 + i_subframe * BBWENET_AF3_FRAME_SIZE * BBWENET_AF2_OUT_CHANNELS + i_channel * BBWENET_TDSHAPE2_FRAME_SIZE,
x_buffer1 + i_subframe * BBWENET_TDSHAPE1_FRAME_SIZE * BBWENET_AF2_OUT_CHANNELS + i_channel * BBWENET_TDSHAPE1_FRAME_SIZE,
BBWENET_TDSHAPE1_FRAME_SIZE
);
}
#ifdef DEBUG_BBWENET
fwrite(x_buffer2 + i_subframe * BBWENET_AF2_OUT_CHANNELS * BBWENET_TDSHAPE2_FRAME_SIZE, sizeof(float), BBWENET_TDSHAPE2_FRAME_SIZE, f_up15_1);
fwrite(x_buffer2 + i_subframe * BBWENET_AF2_OUT_CHANNELS * BBWENET_TDSHAPE2_FRAME_SIZE + BBWENET_TDSHAPE2_FRAME_SIZE, sizeof(float), BBWENET_TDSHAPE2_FRAME_SIZE, f_up15_2);
fwrite(x_buffer2 + i_subframe * BBWENET_AF2_OUT_CHANNELS * BBWENET_TDSHAPE2_FRAME_SIZE + 2 * BBWENET_TDSHAPE2_FRAME_SIZE, sizeof(float), BBWENET_TDSHAPE2_FRAME_SIZE, f_up15_3);
#endif
/* tdshape on second channel (in place) */
adashape_process_frame(
&state->tdshape2_state,
x_buffer2 + i_subframe * BBWENET_AF2_OUT_CHANNELS * BBWENET_TDSHAPE2_FRAME_SIZE + BBWENET_TDSHAPE2_FRAME_SIZE,
x_buffer2 + i_subframe * BBWENET_AF2_OUT_CHANNELS * BBWENET_TDSHAPE2_FRAME_SIZE + BBWENET_TDSHAPE2_FRAME_SIZE,
latent_features + i_subframe * BBWENET_COND_DIM,
&layers->bbwenet_tdshape2_alpha1_f,
&layers->bbwenet_tdshape2_alpha1_t,
&layers->bbwenet_tdshape2_alpha2,
BBWENET_TDSHAPE2_FEATURE_DIM,
BBWENET_TDSHAPE2_FRAME_SIZE,
BBWENET_TDSHAPE2_AVG_POOL_K,
BBWENET_TDSHAPE2_INTERPOLATE_K,
arch
);
#ifdef DEBUG_BBWENET
fwrite(x_buffer2 + i_subframe * BBWENET_AF2_OUT_CHANNELS * BBWENET_TDSHAPE2_FRAME_SIZE + BBWENET_TDSHAPE2_FRAME_SIZE, sizeof(float), BBWENET_TDSHAPE2_FRAME_SIZE, f_up15_shape);
#endif
/* non-linear activation of third channel (in place)*/
apply_valin_activation(
x_buffer2 + i_subframe * BBWENET_AF2_OUT_CHANNELS * BBWENET_TDSHAPE2_FRAME_SIZE + 2 * BBWENET_TDSHAPE2_FRAME_SIZE,
BBWENET_TDSHAPE2_FRAME_SIZE
);
#ifdef DEBUG_BBWENET
fwrite(x_buffer2 + i_subframe * BBWENET_AF2_OUT_CHANNELS * BBWENET_TDSHAPE2_FRAME_SIZE + 2 * BBWENET_TDSHAPE2_FRAME_SIZE, sizeof(float), BBWENET_TDSHAPE2_FRAME_SIZE, f_up15_func);
#endif
}
/* final mixing */
celt_assert(BBWENET_AF3_OUT_CHANNELS == 1);
for (i_subframe = 0; i_subframe < num_subframes; i_subframe++)
{
adaconv_process_frame(
&state->af3_state,
x_out + i_subframe * BBWENET_AF3_FRAME_SIZE,
x_buffer2 + i_subframe * BBWENET_TDSHAPE2_FRAME_SIZE * BBWENET_AF2_OUT_CHANNELS,
latent_features + i_subframe * BBWENET_COND_DIM,
&layers->bbwenet_af3_kernel,
&layers->bbwenet_af3_gain,
BBWENET_COND_DIM,
BBWENET_AF3_FRAME_SIZE,
BBWENET_AF3_OVERLAP_SIZE,
BBWENET_AF3_IN_CHANNELS,
BBWENET_AF3_OUT_CHANNELS,
BBWENET_AF3_KERNEL_SIZE,
BBWENET_AF3_LEFT_PADDING,
BBWENET_AF3_FILTER_GAIN_A,
BBWENET_AF3_FILTER_GAIN_B,
BBWENET_AF3_SHAPE_GAIN,
hBBWENET->window48,
arch);
}
#ifdef DEBUG_BBWENET
fwrite(x_out, sizeof(float), num_subframes * BBWENET_AF3_FRAME_SIZE, f_af3_1);
#endif
}
static void reset_bbwenet_state(BBWENetState *state)
{
OPUS_CLEAR(state, 1);
init_adaconv_state(&state->af1_state);
init_adaconv_state(&state->af2_state);
init_adaconv_state(&state->af3_state);
init_adashape_state(&state->tdshape1_state);
init_adashape_state(&state->tdshape2_state);
}
static int init_bbwenet(BBWENet *hBBWENET, const WeightArray *weights)
{
int ret = 0;
OPUS_CLEAR(hBBWENET, 1);
celt_assert(weights != NULL);
ret = init_bbwenetlayers(&hBBWENET->layers, weights);
compute_overlap_window(hBBWENET->window16, BBWENET_AF1_OVERLAP_SIZE);
compute_overlap_window(hBBWENET->window32, BBWENET_AF2_OVERLAP_SIZE);
compute_overlap_window(hBBWENET->window48, BBWENET_AF3_OVERLAP_SIZE);
return ret;
}
#endif
#endif /* ENABLE_OSCE_BWE */
/* API */
void osce_reset(silk_OSCE_struct *hOSCE, int method)
{
OSCEState *state = &hOSCE->state;
OPUS_CLEAR(&hOSCE->features, 1);
switch(method)
{
case OSCE_METHOD_NONE:
break;
#ifndef DISABLE_LACE
case OSCE_METHOD_LACE:
reset_lace_state(&state->lace);
break;
#endif
#ifndef DISABLE_NOLACE
case OSCE_METHOD_NOLACE:
reset_nolace_state(&state->nolace);
break;
#endif
default:
celt_assert(0 && "method not defined"); /* Question: return error code? */
}
hOSCE->method = method;
hOSCE->features.reset = 2;
}
#ifdef ENABLE_OSCE_BWE
void osce_bwe_reset(silk_OSCE_BWE_struct *hOSCEBWE)
{
int k;
OPUS_CLEAR(&hOSCEBWE->features, 1);
#if 1
/* weird python initialization: Fix eventually! */
for (k = 0; k <= OSCE_BWE_MAX_INSTAFREQ_BIN; k ++)
{
hOSCEBWE->features.last_spec[2*k] = 1e-9;
}
#endif
reset_bbwenet_state(&hOSCEBWE->state.bbwenet);
}
#endif /* ENABLE_OSCE_BWE */
int osce_load_models(OSCEModel *model, const void *data, int len)
{
int ret = 0;
WeightArray *list;
if (data != NULL && len)
{
/* init from buffer */
parse_weights(&list, data, len);
#ifndef DISABLE_LACE
if (ret == 0) {ret = init_lace(&model->lace, list);}
#endif
#ifndef DISABLE_NOLACE
if (ret == 0) {ret = init_nolace(&model->nolace, list);}
#endif
#ifdef ENABLE_OSCE_BWE
#ifndef DISABLE_BBWENET
if (ret == 0) {ret = init_bbwenet(&model->bbwenet, list);}
#endif
#endif /* ENABLE_OSCE_BWE */
free(list);
} else
{
#ifdef USE_WEIGHTS_FILE
return -1;
#else
#ifndef DISABLE_LACE
if (ret == 0) {ret = init_lace(&model->lace, lacelayers_arrays);}
#endif
#ifndef DISABLE_NOLACE
if (ret == 0) {ret = init_nolace(&model->nolace, nolacelayers_arrays);}
#endif
#ifdef ENABLE_OSCE_BWE
#ifndef DISABLE_BBWENET
if (ret == 0) {ret = init_bbwenet(&model->bbwenet, bbwenetlayers_arrays);}
#endif
#endif /* ENABLE_OSCE_BWE */
#endif /* USE_WEIGHTS_FILE */
}
ret = ret ? -1 : 0;
return ret;
}
#ifdef ENABLE_OSCE_BWE
void osce_bwe(
OSCEModel *model, /* I OSCE model struct */
silk_OSCE_BWE_struct *psOSCEBWE, /* I/O OSCE BWE state */
opus_int16 xq48[], /* O bandwidth-extended speech */
opus_int16 xq16[], /* I Decoded speech */
opus_int32 xq16_len, /* I Length of xq16 in samples */
int arch /* I Run-time architecture */
)
{
float in_buffer[320];
float out_buffer[3*320];
float features[2 * OSCE_BWE_FEATURE_DIM];
int num_frames, i;
/* currently restricting to 10 or 20-ms frames */
celt_assert(xq16_len == 160 || xq16_len == 320);
num_frames = xq16_len / 160;
/* scale input */
for (i = 0; i < xq16_len; i++)
{
in_buffer[i] = ((float) xq16[i]) * (1.f/32768.f);
}
osce_bwe_calculate_features(&psOSCEBWE->features, features, xq16, xq16_len);
#if 0
/* just upsampling for now */
upsamp_2x(&psOSCEBWE->state.bbwenet.resampler_state[0], out_buffer, in_buffer, xq16_len);
interpol_3_2(&psOSCEBWE->state.bbwenet.resampler_state[0], out_buffer, out_buffer, 2 * xq16_len);
#else
/* process frames */
bbwenet_process_frames(
&model->bbwenet,
&psOSCEBWE->state.bbwenet,
out_buffer,
in_buffer,
features,
num_frames,
arch
);
#endif
/* scale and delay output */
OPUS_COPY(xq48, psOSCEBWE->state.bbwenet.outbut_buffer, OSCE_BWE_OUTPUT_DELAY);
for (i = 0; i < 3 * xq16_len - OSCE_BWE_OUTPUT_DELAY; i++)
{
float tmp = 32768.f * out_buffer[i];
if (tmp > 32767.f) tmp = 32767.f;
if (tmp < -32767.f) tmp = -32767.f;
xq48[i + OSCE_BWE_OUTPUT_DELAY] = float2int(tmp);
}
for (i = 0; i < OSCE_BWE_OUTPUT_DELAY; i++)
{
float tmp = 32768.f * out_buffer[3 * xq16_len - OSCE_BWE_OUTPUT_DELAY + i];
if (tmp > 32767.f) tmp = 32767.f;
if (tmp < -32767.f) tmp = -32767.f;
psOSCEBWE->state.bbwenet.outbut_buffer[i] = float2int(tmp);
}
}
#endif
void osce_enhance_frame(
OSCEModel *model, /* I OSCE model struct */
silk_decoder_state *psDec, /* I/O Decoder state */
silk_decoder_control *psDecCtrl, /* I Decoder control */
opus_int16 xq[], /* I/O Decoded speech */
opus_int32 num_bits, /* I Size of SILK payload in bits */
int arch /* I Run-time architecture */
)
{
float in_buffer[320];
float out_buffer[320];
float features[4 * OSCE_FEATURE_DIM];
float numbits[2];
int periods[4];
int i;
int method;
/* enhancement only implemented for 20 ms frame at 16kHz */
if (psDec->fs_kHz != 16 || psDec->nb_subfr != 4)
{
osce_reset(&psDec->osce, psDec->osce.method);
return;
}
osce_calculate_features(psDec, psDecCtrl, features, numbits, periods, xq, num_bits);
/* scale input */
for (i = 0; i < 320; i++)
{
in_buffer[i] = ((float) xq[i]) * (1.f/32768.f);
}
if (model->loaded)
method = psDec->osce.method;
else
method = OSCE_METHOD_NONE;
switch(method)
{
case OSCE_METHOD_NONE:
OPUS_COPY(out_buffer, in_buffer, 320);
break;
#ifndef DISABLE_LACE
case OSCE_METHOD_LACE:
lace_process_20ms_frame(&model->lace, &psDec->osce.state.lace, out_buffer, in_buffer, features, numbits, periods, arch);
break;
#endif
#ifndef DISABLE_NOLACE
case OSCE_METHOD_NOLACE:
nolace_process_20ms_frame(&model->nolace, &psDec->osce.state.nolace, out_buffer, in_buffer, features, numbits, periods, arch);
break;
#endif
default:
celt_assert(0 && "method not defined");
}
#ifdef ENABLE_OSCE_TRAINING_DATA
int k;
static FILE *flpc = NULL;
static FILE *fgain = NULL;
static FILE *fltp = NULL;
static FILE *fperiod = NULL;
static FILE *fnoisy16k = NULL;
static FILE* f_numbits = NULL;
static FILE* f_numbits_smooth = NULL;
if (flpc == NULL) {flpc = fopen("features_lpc.f32", "wb");}
if (fgain == NULL) {fgain = fopen("features_gain.f32", "wb");}
if (fltp == NULL) {fltp = fopen("features_ltp.f32", "wb");}
if (fperiod == NULL) {fperiod = fopen("features_period.s16", "wb");}
if (fnoisy16k == NULL) {fnoisy16k = fopen("noisy_16k.s16", "wb");}
if(f_numbits == NULL) {f_numbits = fopen("features_num_bits.s32", "wb");}
if (f_numbits_smooth == NULL) {f_numbits_smooth = fopen("features_num_bits_smooth.f32", "wb");}
fwrite(&num_bits, sizeof(num_bits), 1, f_numbits);
fwrite(&(psDec->osce.features.numbits_smooth), sizeof(psDec->osce.features.numbits_smooth), 1, f_numbits_smooth);
for (k = 0; k < psDec->nb_subfr; k++)
{
float tmp;
int16_t itmp;
float lpc_buffer[16] = {0};
opus_int16 *A_Q12, *B_Q14;
(void) num_bits;
(void) arch;
/* gain */
tmp = (float) psDecCtrl->Gains_Q16[k] / (1UL << 16);
fwrite(&tmp, sizeof(tmp), 1, fgain);
/* LPC */
A_Q12 = psDecCtrl->PredCoef_Q12[ k >> 1 ];
for (i = 0; i < psDec->LPC_order; i++)
{
lpc_buffer[i] = (float) A_Q12[i] / (1U << 12);
}
fwrite(lpc_buffer, sizeof(lpc_buffer[0]), 16, flpc);
/* LTP */
B_Q14 = &psDecCtrl->LTPCoef_Q14[ k * LTP_ORDER ];
for (i = 0; i < 5; i++)
{
tmp = (float) B_Q14[i] / (1U << 14);
fwrite(&tmp, sizeof(tmp), 1, fltp);
}
/* periods */
itmp = psDec->indices.signalType == TYPE_VOICED ? psDecCtrl->pitchL[ k ] : 0;
fwrite(&itmp, sizeof(itmp), 1, fperiod);
}
fwrite(xq, psDec->nb_subfr * psDec->subfr_length, sizeof(xq[0]), fnoisy16k);
#endif
if (psDec->osce.features.reset > 1)
{
OPUS_COPY(out_buffer, in_buffer, 320);
psDec->osce.features.reset --;
}
else if (psDec->osce.features.reset)
{
osce_cross_fade_10ms(out_buffer, in_buffer, 320);
psDec->osce.features.reset = 0;
}
/* scale output */
for (i = 0; i < 320; i++)
{
float tmp = 32768.f * out_buffer[i];
if (tmp > 32767.f) tmp = 32767.f;
if (tmp < -32767.f) tmp = -32767.f;
xq[i] = float2int(tmp);
}
}
#if 0
#include <stdio.h>
void lace_feature_net_compare(
const char * prefix,
int num_frames,
LACE* hLACE
)
{
char in_feature_file[256];
char out_feature_file[256];
char numbits_file[256];
char periods_file[256];
char message[512];
int i_frame, i_feature;
float mse;
float in_features[4 * LACE_NUM_FEATURES];
float out_features[4 * LACE_COND_DIM];
float out_features2[4 * LACE_COND_DIM];
float numbits[2];
int periods[4];
init_lace(hLACE);
FILE *f_in_features, *f_out_features, *f_numbits, *f_periods;
strcpy(in_feature_file, prefix);
strcat(in_feature_file, "_in_features.f32");
f_in_features = fopen(in_feature_file, "rb");
if (f_in_features == NULL)
{
sprintf(message, "could not open file %s", in_feature_file);
perror(message);
exit(1);
}
strcpy(out_feature_file, prefix);
strcat(out_feature_file, "_out_features.f32");
f_out_features = fopen(out_feature_file, "rb");
if (f_out_features == NULL)
{
sprintf(message, "could not open file %s", out_feature_file);
perror(message);
exit(1);
}
strcpy(periods_file, prefix);
strcat(periods_file, "_periods.s32");
f_periods = fopen(periods_file, "rb");
if (f_periods == NULL)
{
sprintf(message, "could not open file %s", periods_file);
perror(message);
exit(1);
}
strcpy(numbits_file, prefix);
strcat(numbits_file, "_numbits.f32");
f_numbits = fopen(numbits_file, "rb");
if (f_numbits == NULL)
{
sprintf(message, "could not open file %s", numbits_file);
perror(message);
exit(1);
}
for (i_frame = 0; i_frame < num_frames; i_frame ++)
{
if(fread(in_features, sizeof(float), 4 * LACE_NUM_FEATURES, f_in_features) != 4 * LACE_NUM_FEATURES)
{
fprintf(stderr, "could not read frame %d from in_features\n", i_frame);
exit(1);
}
if(fread(out_features, sizeof(float), 4 * LACE_COND_DIM, f_out_features) != 4 * LACE_COND_DIM)
{
fprintf(stderr, "could not read frame %d from out_features\n", i_frame);
exit(1);
}
if(fread(periods, sizeof(int), 4, f_periods) != 4)
{
fprintf(stderr, "could not read frame %d from periods\n", i_frame);
exit(1);
}
if(fread(numbits, sizeof(float), 2, f_numbits) != 2)
{
fprintf(stderr, "could not read frame %d from numbits\n", i_frame);
exit(1);
}
lace_feature_net(hLACE, out_features2, in_features, numbits, periods);
float mse = 0;
for (int i = 0; i < 4 * LACE_COND_DIM; i ++)
{
mse += pow(out_features[i] - out_features2[i], 2);
}
mse /= (4 * LACE_COND_DIM);
printf("rmse: %f\n", sqrt(mse));
}
fclose(f_in_features);
fclose(f_out_features);
fclose(f_numbits);
fclose(f_periods);
}
void lace_demo(
char *prefix,
char *output
)
{
char feature_file[256];
char numbits_file[256];
char periods_file[256];
char x_in_file[256];
char message[512];
int i_frame;
float mse;
float features[4 * LACE_NUM_FEATURES];
float numbits[2];
int periods[4];
float x_in[4 * LACE_FRAME_SIZE];
int16_t x_out[4 * LACE_FRAME_SIZE];
float buffer[4 * LACE_FRAME_SIZE];
LACE hLACE;
int frame_counter = 0;
FILE *f_features, *f_numbits, *f_periods, *f_x_in, *f_x_out;
init_lace(&hLACE);
strcpy(feature_file, prefix);
strcat(feature_file, "_features.f32");
f_features = fopen(feature_file, "rb");
if (f_features == NULL)
{
sprintf(message, "could not open file %s", feature_file);
perror(message);
exit(1);
}
strcpy(x_in_file, prefix);
strcat(x_in_file, "_x_in.f32");
f_x_in = fopen(x_in_file, "rb");
if (f_x_in == NULL)
{
sprintf(message, "could not open file %s", x_in_file);
perror(message);
exit(1);
}
f_x_out = fopen(output, "wb");
if (f_x_out == NULL)
{
sprintf(message, "could not open file %s", output);
perror(message);
exit(1);
}
strcpy(periods_file, prefix);
strcat(periods_file, "_periods.s32");
f_periods = fopen(periods_file, "rb");
if (f_periods == NULL)
{
sprintf(message, "could not open file %s", periods_file);
perror(message);
exit(1);
}
strcpy(numbits_file, prefix);
strcat(numbits_file, "_numbits.f32");
f_numbits = fopen(numbits_file, "rb");
if (f_numbits == NULL)
{
sprintf(message, "could not open file %s", numbits_file);
perror(message);
exit(1);
}
printf("processing %s\n", prefix);
while (fread(x_in, sizeof(float), 4 * LACE_FRAME_SIZE, f_x_in) == 4 * LACE_FRAME_SIZE)
{
printf("\rframe: %d", frame_counter++);
if(fread(features, sizeof(float), 4 * LACE_NUM_FEATURES, f_features) != 4 * LACE_NUM_FEATURES)
{
fprintf(stderr, "could not read frame %d from features\n", i_frame);
exit(1);
}
if(fread(periods, sizeof(int), 4, f_periods) != 4)
{
fprintf(stderr, "could not read frame %d from periods\n", i_frame);
exit(1);
}
if(fread(numbits, sizeof(float), 2, f_numbits) != 2)
{
fprintf(stderr, "could not read frame %d from numbits\n", i_frame);
exit(1);
}
lace_process_20ms_frame(
&hLACE,
buffer,
x_in,
features,
numbits,
periods
);
for (int n=0; n < 4 * LACE_FRAME_SIZE; n ++)
{
float tmp = (1UL<<15) * buffer[n];
tmp = CLIP(tmp, -32768, 32767);
x_out[n] = (int16_t) round(tmp);
}
fwrite(x_out, sizeof(int16_t), 4 * LACE_FRAME_SIZE, f_x_out);
}
printf("\ndone!\n");
fclose(f_features);
fclose(f_numbits);
fclose(f_periods);
fclose(f_x_in);
fclose(f_x_out);
}
void nolace_demo(
char *prefix,
char *output
)
{
char feature_file[256];
char numbits_file[256];
char periods_file[256];
char x_in_file[256];
char message[512];
int i_frame;
float mse;
float features[4 * LACE_NUM_FEATURES];
float numbits[2];
int periods[4];
float x_in[4 * LACE_FRAME_SIZE];
int16_t x_out[4 * LACE_FRAME_SIZE];
float buffer[4 * LACE_FRAME_SIZE];
NoLACE hNoLACE;
int frame_counter = 0;
FILE *f_features, *f_numbits, *f_periods, *f_x_in, *f_x_out;
init_nolace(&hNoLACE);
strcpy(feature_file, prefix);
strcat(feature_file, "_features.f32");
f_features = fopen(feature_file, "rb");
if (f_features == NULL)
{
sprintf(message, "could not open file %s", feature_file);
perror(message);
exit(1);
}
strcpy(x_in_file, prefix);
strcat(x_in_file, "_x_in.f32");
f_x_in = fopen(x_in_file, "rb");
if (f_x_in == NULL)
{
sprintf(message, "could not open file %s", x_in_file);
perror(message);
exit(1);
}
f_x_out = fopen(output, "wb");
if (f_x_out == NULL)
{
sprintf(message, "could not open file %s", output);
perror(message);
exit(1);
}
strcpy(periods_file, prefix);
strcat(periods_file, "_periods.s32");
f_periods = fopen(periods_file, "rb");
if (f_periods == NULL)
{
sprintf(message, "could not open file %s", periods_file);
perror(message);
exit(1);
}
strcpy(numbits_file, prefix);
strcat(numbits_file, "_numbits.f32");
f_numbits = fopen(numbits_file, "rb");
if (f_numbits == NULL)
{
sprintf(message, "could not open file %s", numbits_file);
perror(message);
exit(1);
}
printf("processing %s\n", prefix);
while (fread(x_in, sizeof(float), 4 * LACE_FRAME_SIZE, f_x_in) == 4 * LACE_FRAME_SIZE)
{
printf("\rframe: %d", frame_counter++);
if(fread(features, sizeof(float), 4 * LACE_NUM_FEATURES, f_features) != 4 * LACE_NUM_FEATURES)
{
fprintf(stderr, "could not read frame %d from features\n", i_frame);
exit(1);
}
if(fread(periods, sizeof(int), 4, f_periods) != 4)
{
fprintf(stderr, "could not read frame %d from periods\n", i_frame);
exit(1);
}
if(fread(numbits, sizeof(float), 2, f_numbits) != 2)
{
fprintf(stderr, "could not read frame %d from numbits\n", i_frame);
exit(1);
}
nolace_process_20ms_frame(
&hNoLACE,
buffer,
x_in,
features,
numbits,
periods
);
for (int n=0; n < 4 * LACE_FRAME_SIZE; n ++)
{
float tmp = (1UL<<15) * buffer[n];
tmp = CLIP(tmp, -32768, 32767);
x_out[n] = (int16_t) round(tmp);
}
fwrite(x_out, sizeof(int16_t), 4 * LACE_FRAME_SIZE, f_x_out);
}
printf("\ndone!\n");
fclose(f_features);
fclose(f_numbits);
fclose(f_periods);
fclose(f_x_in);
fclose(f_x_out);
}
int main()
{
#if 0
LACE hLACE;
lace_feature_net_compare("testvec2/lace", 5, &hLACE);
lace_demo("testdata/test9", "out_lace_c_9kbps.pcm");
lace_demo("testdata/test6", "out_lace_c_6kbps.pcm");
#endif
nolace_demo("testdata/test9", "out_nolace_c_9kbps.pcm");
}
#endif
/*gcc -I ../include -I . -I ../silk -I ../celt osce.c nndsp.c lace_data.c nolace_data.c nnet.c parse_lpcnet_weights.c -lm -o lacetest*/