ivector-transform.cc File Reference

#include "base/kaldi-common.h"
#include "util/common-utils.h"
#include "gmm/am-diag-gmm.h"
#include "ivector/ivector-extractor.h"
#include "util/kaldi-thread.h"

Include dependency graph for ivector-transform.cc:

Go to the source code of this file.

Functions
int	main (int argc, char *argv[])

Function Documentation

main()

int main	(	int	argc,
		char *	argv[]
	)

Definition at line 28 of file ivector-transform.cc.

References VectorBase< Real >::AddMatVec(), VectorBase< Real >::CopyColFromMat(), VectorBase< Real >::Dim(), SequentialTableReader< Holder >::Done(), ParseOptions::GetArg(), KALDI_ASSERT, KALDI_LOG, SequentialTableReader< Holder >::Key(), kaldi::kNoTrans, SequentialTableReader< Holder >::Next(), ParseOptions::NumArgs(), MatrixBase< Real >::NumCols(), MatrixBase< Real >::NumRows(), ParseOptions::PrintUsage(), ParseOptions::Read(), kaldi::ReadKaldiObject(), SequentialTableReader< Holder >::Value(), kaldi::VecVec(), and TableWriter< Holder >::Write().

                                 {
  using namespace kaldi;
  typedef kaldi::int32 int32;
  try {
    const char *usage =
        "Multiplies iVectors (on the left) by a supplied transformation matrix\n"
        "\n"
        "Usage:  ivector-transform [options] <matrix-in> <ivector-rspecifier>"
        "<ivector-wspecifier>\n"
        "e.g.: \n"
        " ivector-transform transform.mat ark:ivectors.ark ark:transformed_ivectors.ark\n";

    ParseOptions po(usage);

    po.Read(argc, argv);

    if (po.NumArgs() != 3) {
      po.PrintUsage();
      exit(1);
    }

    std::string matrix_rxfilename = po.GetArg(1),
        ivector_rspecifier = po.GetArg(2),
        ivector_wspecifier = po.GetArg(3);


    Matrix<BaseFloat> transform;
    ReadKaldiObject(matrix_rxfilename, &transform);

    int32 num_done = 0;

    // The following quantities will be needed if we're doing
    // an affine transform (i.e. linear plus an offset)
    SubMatrix<BaseFloat> linear_term(transform,
                                     0, transform.NumRows(),
                                     0, transform.NumCols() - 1);
    Vector<BaseFloat> constant_term(transform.NumRows());
    constant_term.CopyColFromMat(transform, transform.NumCols() - 1);

    Vector<double> sum(transform.NumRows());
    double sumsq = 0.0;

    SequentialBaseFloatVectorReader ivector_reader(ivector_rspecifier);
    BaseFloatVectorWriter ivector_writer(ivector_wspecifier);

    for (; !ivector_reader.Done(); ivector_reader.Next()) {
      std::string key = ivector_reader.Key();
      const Vector<BaseFloat> &ivector = ivector_reader.Value();

      Vector<BaseFloat> transformed_ivector(transform.NumRows());
      if (ivector.Dim() == transform.NumCols()) {
        transformed_ivector.AddMatVec(1.0, transform, kNoTrans, ivector, 0.0);
      } else {
        KALDI_ASSERT(ivector.Dim() == transform.NumCols() - 1);
        transformed_ivector.CopyFromVec(constant_term);
        transformed_ivector.AddMatVec(1.0, linear_term, kNoTrans, ivector, 1.0);
      }
      sum.AddVec(1.0, transformed_ivector);
      sumsq += VecVec(transformed_ivector, transformed_ivector);
      ivector_writer.Write(key, transformed_ivector);
      num_done++;
    }

    KALDI_LOG << "Processed " << num_done << " iVectors.";
    if (num_done != 0) {
      sum.Scale(1.0 / num_done);
      sumsq /= num_done;
      BaseFloat mean_length = sum.Norm(2.0),
          variance = sumsq - VecVec(sum, sum),
          avg_len = sqrt(variance),
          norm_length = avg_len / sqrt(transform.NumRows());
      KALDI_LOG << "Norm of mean was " << mean_length
                << " (should be close to zero), length divided by sqrt(dim) was "
                << norm_length << " (should probably be close to one)";
    }
    return (num_done != 0 ? 0 : 1);
  } catch(const std::exception &e) {
    std::cerr << e.what();
    return -1;
  }
}