NumCpp  2.9.0
A Templatized Header Only C++ Implementation of the Python NumPy Library
pinv.hpp
Go to the documentation of this file.
1 #pragma once
29 
30 #include <string>
31 
33 #include "NumCpp/Core/Types.hpp"
35 #include "NumCpp/Linalg/svd.hpp"
36 #include "NumCpp/NdArray.hpp"
37 
38 namespace nc
39 {
40  namespace linalg
41  {
42  //============================================================================
43  // Method Description:
51  template<typename dtype>
53  {
55 
59  svd(inArray, u, d, v);
60 
61  const auto inShape = inArray.shape();
62  auto dPlus = nc::zeros<double>(inShape.cols, inShape.rows); // transpose
63 
64  for (uint32 i = 0; i < d.shape().rows; ++i)
65  {
66  dPlus(i, i) = 1. / d(i, i);
67  }
68 
69  return v.transpose().dot(dPlus).dot(u.transpose());
70  }
71  } // namespace linalg
72 } // namespace nc
#define STATIC_ASSERT_ARITHMETIC_OR_COMPLEX(dtype)
Definition: StaticAsserts.hpp:49
Shape shape() const noexcept
Definition: NdArrayCore.hpp:4092
NdArray< dtype > transpose() const
Definition: NdArrayCore.hpp:4457
NdArray< dtype > dot(const NdArray< dtype > &inOtherArray) const
Definition: NdArrayCore.hpp:2520
void svd(const NdArray< dtype > &inArray, NdArray< double > &outU, NdArray< double > &outS, NdArray< double > &outVt)
Definition: svd.hpp:53
NdArray< double > pinv(const NdArray< dtype > &inArray)
Definition: pinv.hpp:52
Definition: Coordinate.hpp:45
std::uint32_t uint32
Definition: Types.hpp:40