rvoip-codec-core 0.2.2

G.711 and optional G.729A/G.729AB audio codec implementation for RVOIP
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78

//! Provenance: Gain prediction/quantization adapted from ITU G.729 gain predictor and quantizer flow.
//! Q-format: Pitch/code gains and predictors use Q14/Q15 arithmetic with 32-bit intermediates.

use crate::codecs::g729::impls::codec::state::EncoderState;
use crate::codecs::g729::impls::constants::{L_INTER10, L_SUBFR, L_THRESH_ERR};
use crate::codecs::g729::impls::dsp::arith::{add, sub};
use crate::codecs::g729::impls::dsp::arith32::{l_add, l_sub};
use crate::codecs::g729::impls::dsp::oper32::{l_extract, mpy_32_16};
use crate::codecs::g729::impls::dsp::shift::l_shl;
use crate::codecs::g729::impls::dsp::types::{DspContext, Word16, Word32};
use crate::codecs::g729::impls::tables::annexa::TAB_ZONE;

pub(crate) fn test_excitation_error(state: &EncoderState, t0: i16, t0_frac: i16) -> i16 {
    let mut ctx = DspContext::default();
    let t1 = if t0_frac > 0 {
        add(&mut ctx, Word16(t0), Word16(1)).0
    } else {
        t0
    };

    let mut i = sub(&mut ctx, Word16(t1), Word16((L_SUBFR + L_INTER10) as i16)).0;
    if i < 0 {
        i = 0;
    }
    let zone1 = TAB_ZONE[i as usize] as usize;
    i = add(&mut ctx, Word16(t1), Word16((L_INTER10 - 2) as i16)).0;
    let zone2 = TAB_ZONE[i as usize] as usize;

    let mut l_maxloc = Word32(-1);
    for zi in (zone1..=zone2).rev() {
        let l_acc = l_sub(&mut ctx, Word32(state.l_exc_err[zi]), l_maxloc);
        if l_acc.0 > 0 {
            l_maxloc = Word32(state.l_exc_err[zi]);
        }
    }
    if l_sub(&mut ctx, l_maxloc, Word32(L_THRESH_ERR)).0 > 0 {
        1
    } else {
        0
    }
}

pub(crate) fn update_excitation_error(state: &mut EncoderState, gain_pit: i16, t0: i16) {
    let mut ctx = DspContext::default();
    let mut l_worst = Word32(-1);
    let n = sub(&mut ctx, Word16(t0), Word16(L_SUBFR as i16)).0;

    if n < 0 {
        let mut l_temp = Word32(state.l_exc_err[0]);
        for _ in 0..2 {
            let (hi, lo) = l_extract(l_temp);
            l_temp = mpy_32_16(hi, lo, Word16(gain_pit));
            l_temp = l_shl(&mut ctx, l_temp, 1);
            l_temp = l_add(&mut ctx, Word32(0x0000_4000), l_temp);
            if l_sub(&mut ctx, l_temp, l_worst).0 > 0 {
                l_worst = l_temp;
            }
        }
    } else {
        let zone1 = TAB_ZONE[n as usize] as usize;
        let i = sub(&mut ctx, Word16(t0), Word16(1)).0;
        let zone2 = TAB_ZONE[i as usize] as usize;
        for zi in zone1..=zone2 {
            let (hi, lo) = l_extract(Word32(state.l_exc_err[zi]));
            let mut l_temp = mpy_32_16(hi, lo, Word16(gain_pit));
            l_temp = l_shl(&mut ctx, l_temp, 1);
            l_temp = l_add(&mut ctx, Word32(0x0000_4000), l_temp);
            if l_sub(&mut ctx, l_temp, l_worst).0 > 0 {
                l_worst = l_temp;
            }
        }
    }

    for i in (1..4).rev() {
        state.l_exc_err[i] = state.l_exc_err[i - 1];
    }
    state.l_exc_err[0] = l_worst.0;
}