#!/bin/bash
  
  # Copyright 2014  Guoguo Chen
  # Apache 2.0
  
  # Begin configuration section.
  nj=4
  cmd=run.pl
  tscale=1.0      # transition scale.
  loopscale=0.1   # scale for self-loops.
  cleanup=true
  ngram_order=1
  srilm_options="-wbdiscount"   # By default, use Witten-Bell discounting in SRILM
  # End configuration section.
  
  set -e
  
  echo "$0 $@"
  
  [ -f ./path.sh ] && . ./path.sh
  . parse_options.sh || exit 1;
  
  if [ $# -ne 4 ]; then
    echo "This script builds one decoding graph for each truncated utterance in"
    echo "segmentation. It first calls steps/cleanup/make_utterance_graph.sh to"
    echo "build one decoding graph for each original utterance, which will be"
    echo "shared by the truncated utterances from the same original utterance."
    echo "We assign the decoding graph to each truncated utterance using the scp"
    echo "file so that we can avoid duplicating the graphs on the disk."
    echo ""
    echo "Usage: $0 [options] <data-dir> <lang-dir> <model-dir> <graph-dir>"
    echo " e.g.: $0 data/train_si284_split/ \\"
    echo "                data/lang exp/tri2b/ exp/tri2b/graph_train_si284_split"
    echo ""
    echo "Options:"
    echo "    --ngram-order           # order of n-gram language model"
    echo "    --srilm-options         # options for ngram-count in SRILM tool"
    echo "    --tscale                # transition scale"
    echo "    --loopscale             # scale for self-loops"
    echo "    --cleanup               # if true, removes the intermediate files"
    exit 1;
  fi
  
  data=$1
  lang=$2
  model_dir=$3
  graph_dir=$4
  
  for f in $data/text.orig $data/orig2utt $lang/L_disambig.fst \
    $lang/words.txt $lang/oov.int $model_dir/final.mdl $model_dir/tree; do
    if [ ! -f $f ]; then
      echo "$0: expected $f to exist"
      exit 1;
    fi
  done
  
  utils/lang/check_phones_compatible.sh $lang/phones.txt $model_dir/phones.txt
  
  # If --ngram-order is larger than 1, we will have to use SRILM
  if [ $ngram_order -gt 1 ]; then
    ngram_count=`which ngram-count` || true
    if [ -z $ngram_count ]; then
      if uname -a | grep 64 >/dev/null; then # some kind of 64 bit...
        sdir=$KALDI_ROOT/tools/srilm/bin/i686-m64
      else
        sdir=$KALDI_ROOT/tools/srilm/bin/i686
      fi
      if [ -f $sdir/ngram-count ]; then
        echo Using SRILM tools from $sdir
        export PATH=$PATH:$sdir
      else
        echo You appear to not have SRILM tools installed, either on your path,
        echo or installed in $sdir.  See tools/install_srilm.sh for installation
        echo instructions.
        exit 1
      fi
    fi
  fi
  
  # Creates one graph for each transcript. We parallelize the process a little
  # bit.
  num_lines=`cat $data/text.orig | wc -l`
  if [ $nj -gt $num_lines ]; then
    nj=$num_lines
    echo "$0: Too many number of jobs, using $nj instead"
  fi
  
  mkdir -p $graph_dir/split$nj
  mkdir -p $graph_dir/log
   
  split_texts=""
  for n in $(seq $nj); do
    mkdir -p $graph_dir/split$nj/$n
    split_texts="$split_texts $graph_dir/split$nj/$n/text"
  done
  utils/split_scp.pl $data/text.orig $split_texts
  
  $cmd JOB=1:$nj $graph_dir/log/make_utterance_graph.JOB.log \
    steps/cleanup/make_utterance_graph.sh --cleanup $cleanup \
    --tscale $tscale --loopscale $loopscale \
    --ngram-order $ngram_order --srilm-options "$srilm_options" \
    $graph_dir/split$nj/JOB/text $lang \
    $model_dir $graph_dir/split$nj/JOB || exit 1;
  
  # Copies files from lang directory.
  mkdir -p $graph_dir
  cp -r $lang/* $graph_dir
  
  am-info --print-args=false $model_dir/final.mdl |\
   grep pdfs | awk '{print $NF}' > $graph_dir/num_pdfs
  
  # Creates the graph table.
  cat $graph_dir/split$nj/*/HCLG.fsts.scp > $graph_dir/split$nj/HCLG.fsts.scp
  fstcopy scp:$graph_dir/split$nj/HCLG.fsts.scp \
    "ark,scp:$graph_dir/HCLG.fsts,$graph_dir/tmp.HCLG.fsts.scp"
  
  # The graphs we created above were indexed by the old utterance id. We have to
  # duplicate them for the new utterance id. We do this in the scp file so we do
  # not have to store the duplicated graphs on the disk.
  cat $graph_dir/tmp.HCLG.fsts.scp | perl -e '
    open(O2U, "<$ARGV[0]") || die "Error: fail to open $ARGV[0]
  ";
    while (<STDIN>) {
      chomp;
      @col = split;
      @col == 2 || die "Error: bad line $_
  ";
      $scp{$col[0]} = $col[1];
    }
    while (<O2U>) {
      chomp;
      @col = split;
      @col >= 2 || die "Error: bad line $_
  ";
      defined($scp{$col[0]}) ||
        die "Error: $col[0] not defined in original scp file
  ";
      for ($i = 1; $i < @col; $i += 1) {
        print "$col[$i] $scp{$col[0]}
  "
      }
    }' $data/orig2utt > $graph_dir/HCLG.fsts.scp
  rm $graph_dir/tmp.HCLG.fsts.scp
  
  if $cleanup; then
    rm -r $graph_dir/split$nj
  fi
  
  exit 0;