a12a8a71bb
The goal is to solve confusion of the "All rights reserved" for licensing
code under an open-source license.
The phrase "All rights reserved" comes from a historical convention that
required this phrase for the copyright protection to apply. This convention
is no longer relevant.
However, even though the phrase has no meaning in establishing the copyright
it has not lost meaning in terms of licensing.
This change makes it so code under the Blender Foundation copyright does
not use "all rights reserved". This is also how the GPL license itself
states how to apply it to the source code:
<one line to give the program's name and a brief idea of what it does.>
Copyright (C) <year> <name of author>
This program is free software ...
This change does not change copyright notice in cases when the copyright
is dual (BF and an author), or just an author of the code. It also does
mot change copyright which is inherited from NaN Holding BV as it needs
some further investigation about what is the proper way to handle it.
60 lines
1.3 KiB
C++
60 lines
1.3 KiB
C++
/* SPDX-License-Identifier: GPL-2.0-or-later
|
|
* Copyright 2015 Blender Foundation */
|
|
|
|
/** \file
|
|
* \ingroup intern_eigen
|
|
*/
|
|
|
|
#ifndef __EIGEN3_SVD_C_API_CC__
|
|
#define __EIGEN3_SVD_C_API_CC__
|
|
|
|
/* Eigen gives annoying huge amount of warnings here, silence them! */
|
|
#if defined(__GNUC__) && !defined(__clang__)
|
|
# pragma GCC diagnostic ignored "-Wlogical-op"
|
|
#endif
|
|
|
|
#ifdef __EIGEN3_SVD_C_API_CC__ /* quiet warning */
|
|
#endif
|
|
|
|
#include <Eigen/Core>
|
|
#include <Eigen/Dense>
|
|
#include <Eigen/SVD>
|
|
|
|
#include "svd.h"
|
|
|
|
using Eigen::JacobiSVD;
|
|
|
|
using Eigen::NoQRPreconditioner;
|
|
|
|
using Eigen::ComputeThinU;
|
|
using Eigen::ComputeThinV;
|
|
|
|
using Eigen::Map;
|
|
using Eigen::MatrixXf;
|
|
using Eigen::VectorXf;
|
|
|
|
using Eigen::Matrix4f;
|
|
|
|
void EIG_svd_square_matrix(const int size, const float *matrix, float *r_U, float *r_S, float *r_V)
|
|
{
|
|
/* Since our matrix is squared, we can use thinU/V. */
|
|
unsigned int flags = (r_U ? ComputeThinU : 0) | (r_V ? ComputeThinV : 0);
|
|
|
|
/* Blender and Eigen matrices are both column-major. */
|
|
JacobiSVD<MatrixXf, NoQRPreconditioner> svd(Map<MatrixXf>((float *)matrix, size, size), flags);
|
|
|
|
if (r_U) {
|
|
Map<MatrixXf>(r_U, size, size) = svd.matrixU();
|
|
}
|
|
|
|
if (r_S) {
|
|
Map<VectorXf>(r_S, size) = svd.singularValues();
|
|
}
|
|
|
|
if (r_V) {
|
|
Map<MatrixXf>(r_V, size, size) = svd.matrixV();
|
|
}
|
|
}
|
|
|
|
#endif /* __EIGEN3_SVD_C_API_CC__ */
|