Blame view
src/matrix/matrix-functions-inl.h
1.82 KB
8dcb6dfcb
first commit
|
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 |
// matrix/matrix-functions-inl.h
// Copyright 2009-2011 Microsoft Corporation
//
// See ../../COPYING for clarification regarding multiple authors
//
// Licensed under the Apache License, Version 2.0 (the "License");
// you may not use this file except in compliance with the License.
// You may obtain a copy of the License at
//
// http://www.apache.org/licenses/LICENSE-2.0
//
// THIS CODE IS PROVIDED *AS IS* BASIS, WITHOUT WARRANTIES OR CONDITIONS OF ANY
// KIND, EITHER EXPRESS OR IMPLIED, INCLUDING WITHOUT LIMITATION ANY IMPLIED
// WARRANTIES OR CONDITIONS OF TITLE, FITNESS FOR A PARTICULAR PURPOSE,
// MERCHANTABLITY OR NON-INFRINGEMENT.
// See the Apache 2 License for the specific language governing permissions and
// limitations under the License.
//
// (*) incorporates, with permission, FFT code from his book
// "Signal Processing with Lapped Transforms", Artech, 1992.
#ifndef KALDI_MATRIX_MATRIX_FUNCTIONS_INL_H_
#define KALDI_MATRIX_MATRIX_FUNCTIONS_INL_H_
namespace kaldi {
//! ComplexMul implements, inline, the complex multiplication b *= a.
template<typename Real> inline void ComplexMul(const Real &a_re, const Real &a_im,
Real *b_re, Real *b_im) {
Real tmp_re = (*b_re * a_re) - (*b_im * a_im);
*b_im = *b_re * a_im + *b_im * a_re;
*b_re = tmp_re;
}
template<typename Real> inline void ComplexAddProduct(const Real &a_re, const Real &a_im,
const Real &b_re, const Real &b_im,
Real *c_re, Real *c_im) {
*c_re += b_re*a_re - b_im*a_im;
*c_im += b_re*a_im + b_im*a_re;
}
template<typename Real> inline void ComplexImExp(Real x, Real *a_re, Real *a_im) {
*a_re = std::cos(x);
*a_im = std::sin(x);
}
} // end namespace kaldi
#endif // KALDI_MATRIX_MATRIX_FUNCTIONS_INL_H_
|