blob: 17b4b293411de89142dd2219d108d8d2daa4b082 (
plain)
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
|
/*
Copyright 2018 Google Inc. All Rights Reserved.
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
Unless required by applicable law or agreed to in writing, software
distributed under the License is distributed on an "AS-IS" BASIS,
WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
See the License for the specific language governing permissions and
limitations under the License.
*/
#ifndef RESONANCE_AUDIO_UTILS_PSEUDOINVERSE_H_
#define RESONANCE_AUDIO_UTILS_PSEUDOINVERSE_H_
#include "Eigen/Dense"
namespace vraudio {
// Computes the Moore-Penrose pseudoinverse of |matrix|.
//
// @tparam MatrixType The type of the input matrix (an Eigen::Matrix).
// @param matrix The input matrix to compute the pseudoinverse of.
// @return The Moore-Penrose pseudoinverse of |matrix|.
template <typename MatrixType>
Eigen::Matrix<typename MatrixType::Scalar, MatrixType::ColsAtCompileTime,
MatrixType::RowsAtCompileTime>
Pseudoinverse(const MatrixType& matrix) {
Eigen::JacobiSVD<Eigen::Matrix<typename MatrixType::Scalar, Eigen::Dynamic,
Eigen::Dynamic>> svd(matrix,
Eigen::ComputeThinU |
Eigen::ComputeThinV);
return svd.solve(
Eigen::Matrix<typename MatrixType::Scalar, MatrixType::RowsAtCompileTime,
MatrixType::RowsAtCompileTime>::Identity(matrix.rows(),
matrix.rows()));
}
} // namespace vraudio
#endif // RESONANCE_AUDIO_UTILS_PSEUDOINVERSE_H_
|