/usr/include/octave-4.2.2/octave/svd.h is in liboctave-dev 4.2.2-1ubuntu1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
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 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 | /*
Copyright (C) 2016-2017 Carnë Draug
Copyright (C) 1994-2016 John W. Eaton
This file is part of Octave.
Octave is free software; you can redistribute it and/or modify it
under the terms of the GNU General Public License as published by the
Free Software Foundation; either version 3 of the License, or (at your
option) any later version.
Octave is distributed in the hope that it will be useful, but WITHOUT
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
for more details.
You should have received a copy of the GNU General Public License
along with Octave; see the file COPYING. If not, see
<http://www.gnu.org/licenses/>.
*/
#if ! defined (octave_svd_h)
#define octave_svd_h 1
#include "octave-config.h"
#include <vector>
namespace octave
{
namespace math
{
template <typename T>
class
svd
{
public:
typedef typename T::real_diag_matrix_type DM_T;
enum class Type
{
std,
economy,
sigma_only
};
enum class Driver
{
GESVD,
GESDD
};
svd (void)
: m_type (), m_driver (), left_sm (), sigma (), right_sm ()
{ }
svd (const T& a, svd::Type type = svd::Type::std,
svd::Driver driver = svd::Driver::GESVD);
svd (const svd& a)
: m_type (a.m_type), m_driver (a.m_driver), left_sm (a.left_sm),
sigma (a.sigma), right_sm (a.right_sm)
{ }
svd& operator = (const svd& a)
{
if (this != &a)
{
m_type = a.m_type;
left_sm = a.left_sm;
sigma = a.sigma;
right_sm = a.right_sm;
m_driver = a.m_driver;
}
return *this;
}
~svd (void) { }
T left_singular_matrix (void) const;
DM_T singular_values (void) const { return sigma; }
T right_singular_matrix (void) const;
private:
typedef typename T::element_type P;
typedef typename DM_T::element_type DM_P;
svd::Type m_type;
svd::Driver m_driver;
T left_sm;
DM_T sigma;
T right_sm;
void gesvd (char& jobu, char& jobv, octave_idx_type m, octave_idx_type n,
P* tmp_data, octave_idx_type m1, DM_P* s_vec, P* u, P* vt,
octave_idx_type nrow_vt1, std::vector<P>& work,
octave_idx_type& lwork, octave_idx_type& info);
void gesdd (char& jobz, octave_idx_type m, octave_idx_type n,
P* tmp_data, octave_idx_type m1, DM_P* s_vec, P* u, P* vt,
octave_idx_type nrow_vt1, std::vector<P>& work,
octave_idx_type& lwork,
octave_idx_type* iwork, octave_idx_type& info);
};
}
}
#endif
|