# $!/bin/bash

# source the settings
. ./path.sh

. utils/parse_options.sh || exit 1

# temporary files
mfccdir=$decode_dir/mfcc
feat_scp=$mfccdir/raw_mfcc.scp

mkdir -p $mfccdir

compute-mfcc-feats  --verbose=0 --config=$MFCC scp:$wav_scp \
  ark,scp:$mfccdir/raw_mfcc.ark,$feat_scp || exit 1;

# # For debugging
# # cgdb -q -x .gdbinit --args \
# compute-mfcc-feats  --verbose=0 --config=$MFCC scp:$wav_scp \
#   ark,t,scp:$mfccdir/raw_mfcc.ark.txt,${feat_scp}.txt || exit 1;
# # For debugging
# add-deltas "scp,s,cs:$feat_scp" "ark,t:$mfccdir/dd_mfcc.ark.txt"

# check the splice-feats use the default left=right=4 context
if [ -z $MAT ] ; then
  # no LDA matrix -> use delta-delta
  feats="ark,s,cs:copy-feats scp:$feat_scp ark:- | add-deltas ark:- ark:- |"
else
  # LDA matrix specified -> using it
  feats="ark,s,cs:copy-feats scp:$feat_scp ark:- | splice-feats ark:- ark:- | transform-feats $MAT ark:- ark:- |"
fi

gmm-latgen-faster --verbose=0 --max-mem=500000000 \
    --beam=$beam --lattice-beam=$latbeam --max-active=$max_active \
    --allow-partial=true --word-symbol-table=$WST \
    $AM $HCLG "$feats" \
    "ark:|gzip -c > $lattice"

lattice-best-path --verbose=0 --lm-scale=15 --word-symbol-table=$WST \
    "ark:gunzip -c $lattice|" ark,t:$gmm_latgen_faster_tra || exit 1;

cat $gmm_latgen_faster_tra | utils/int2sym.pl -f 2- $WST \
    > $gmm_latgen_faster_tra_txt || exit 1

# reference is named based on wav_scp
./build_reference.py $wav_scp $decode_dir
reference=$decode_dir/`basename $wav_scp`.tra

echo; echo "Reference"; echo
cat $reference
echo; echo "Decoded"; echo
cat $gmm_latgen_faster_tra_txt

compute-wer --text --mode=present ark:$reference ark,p:$gmm_latgen_faster_tra_txt