Cycles: Fix missing type declaration in OpenCL image
[blender.git] / intern / eigen / intern / svd.cc
1 /*
2  * ***** BEGIN GPL LICENSE BLOCK *****
3  *
4  * This program is free software; you can redistribute it and/or
5  * modify it under the terms of the GNU General Public License
6  * as published by the Free Software Foundation; either version 2
7  * of the License, or (at your option) any later version.
8  *
9  * This program is distributed in the hope that it will be useful,
10  * but WITHOUT ANY WARRANTY; without even the implied warranty of
11  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
12  * GNU General Public License for more details.
13  *
14  * You should have received a copy of the GNU General Public License
15  * along with this program; if not, write to the Free Software Foundation,
16  * Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
17  *
18  * The Original Code is Copyright (C) 2015 Blender Foundation.
19  * All rights reserved.
20  *
21  * Contributor(s): Blender Foundation,
22  *                 Bastien Montagne
23  *
24  * ***** END GPL LICENSE BLOCK *****
25  */
26
27 #ifndef __EIGEN3_SVD_C_API_CC__
28 #define __EIGEN3_SVD_C_API_CC__
29
30 /* Eigen gives annoying huge amount of warnings here, silence them! */
31 #ifdef __GNUC__
32 #  pragma GCC diagnostic ignored "-Wlogical-op"
33 #endif
34
35 #ifdef __EIGEN3_SVD_C_API_CC__  /* quiet warning */
36 #endif
37
38 #include <Eigen/Core>
39 #include <Eigen/SVD>
40
41 #include "svd.h"
42
43 using Eigen::JacobiSVD;
44
45 using Eigen::NoQRPreconditioner;
46
47 using Eigen::ComputeThinU;
48 using Eigen::ComputeThinV;
49
50 using Eigen::MatrixXf;
51 using Eigen::VectorXf;
52 using Eigen::Map;
53
54 void EIG_svd_square_matrix(const int size, const float *matrix, float *r_U, float *r_S, float *r_V)
55 {
56         /* Since our matrix is squared, we can use thinU/V. */
57         unsigned int flags = (r_U ? ComputeThinU : 0) | (r_V ? ComputeThinV : 0);
58
59         /* Blender and Eigen matrices are both column-major. */
60         JacobiSVD<MatrixXf, NoQRPreconditioner> svd(Map<MatrixXf>((float *)matrix, size, size), flags);
61
62         if (r_U) {
63                 Map<MatrixXf>(r_U, size, size) = svd.matrixU();
64         }
65
66         if (r_S) {
67                 Map<VectorXf>(r_S, size) = svd.singularValues();
68         }
69
70         if (r_V) {
71                 Map<MatrixXf>(r_V, size, size) = svd.matrixV();
72         }
73 }
74
75 #endif  /* __EIGEN3_SVD_C_API_CC__ */