matrix_algo.h

/********************************************************************
 * Copyright (C) 2015 Liangliang Nan <liangliang.nan@gmail.com>
 * https://3d.bk.tudelft.nl/liangliang/
 *
 * This file is part of Easy3D. If it is useful in your research/work,
 * I would be grateful if you show your appreciation by citing it:
 * ------------------------------------------------------------------
 *      Liangliang Nan.
 *      Easy3D: a lightweight, easy-to-use, and efficient C++ library
 *      for processing and rendering 3D data.
 *      Journal of Open Source Software, 6(64), 3255, 2021.
 * ------------------------------------------------------------------
 *
 * Easy3D is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License Version 3
 * as published by the Free Software Foundation.
 *
 * Easy3D 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 this program. If not, see <http://www.gnu.org/licenses/>.
 ********************************************************************/
 
 
#ifndef EASY3D_CORE_MATRIX_ALGO_H
#define EASY3D_CORE_MATRIX_ALGO_H
 
#include <easy3d/core/matrix.h>
 
 
namespace easy3d {
 
    using MATRIX = easy3d::Matrix<double>;
 
    double determinant(const MATRIX &A);
 
 
    bool inverse(const MATRIX&A, MATRIX&invA);
 
 
    MATRIX inverse(const MATRIX&A);
 
 
    void svd_decompose(const MATRIX&A, MATRIX&U, MATRIX&S, MATRIX&V);
 
 
    bool solve_least_squares(const MATRIX&A, const std::vector<double> &b, std::vector<double> &x);
}
 
#endif // EASY3D_MATRIX_ALGOTHMS_H