Fix T67009, T67013: crash with Intel HD Graphics and some older drivers
[blender.git] / source / blender / blenlib / BLI_math_solvers.h
1 /*
2  * This program is free software; you can redistribute it and/or
3  * modify it under the terms of the GNU General Public License
4  * as published by the Free Software Foundation; either version 2
5  * of the License, or (at your option) any later version.
6  *
7  * This program is distributed in the hope that it will be useful,
8  * but WITHOUT ANY WARRANTY; without even the implied warranty of
9  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
10  * GNU General Public License for more details.
11  *
12  * You should have received a copy of the GNU General Public License
13  * along with this program; if not, write to the Free Software Foundation,
14  * Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
15  *
16  * The Original Code is Copyright (C) 2015 by Blender Foundation
17  * All rights reserved.
18  * */
19
20 #ifndef __BLI_MATH_SOLVERS_H__
21 #define __BLI_MATH_SOLVERS_H__
22
23 /** \file
24  * \ingroup bli
25  */
26
27 #ifdef __cplusplus
28 extern "C" {
29 #endif
30
31 #include "BLI_compiler_attrs.h"
32 #include "BLI_math_inline.h"
33
34 #ifdef BLI_MATH_GCC_WARN_PRAGMA
35 #  pragma GCC diagnostic push
36 #  pragma GCC diagnostic ignored "-Wredundant-decls"
37 #endif
38
39 /********************************** Eigen Solvers *********************************/
40
41 bool BLI_eigen_solve_selfadjoint_m3(const float m3[3][3],
42                                     float r_eigen_values[3],
43                                     float r_eigen_vectors[3][3]);
44
45 void BLI_svd_m3(const float m3[3][3], float r_U[3][3], float r_S[], float r_V[3][3]);
46
47 /***************************** Simple Solvers ************************************/
48
49 bool BLI_tridiagonal_solve(
50     const float *a, const float *b, const float *c, const float *d, float *r_x, const int count);
51 bool BLI_tridiagonal_solve_cyclic(
52     const float *a, const float *b, const float *c, const float *d, float *r_x, const int count);
53
54 /* Generic 3 variable Newton's method solver. */
55 typedef void (*Newton3D_DeltaFunc)(void *userdata, const float x[3], float r_delta[3]);
56 typedef void (*Newton3D_JacobianFunc)(void *userdata, const float x[3], float r_jacobian[3][3]);
57 typedef bool (*Newton3D_CorrectionFunc)(void *userdata,
58                                         const float x[3],
59                                         float step[3],
60                                         float x_next[3]);
61
62 bool BLI_newton3d_solve(Newton3D_DeltaFunc func_delta,
63                         Newton3D_JacobianFunc func_jacobian,
64                         Newton3D_CorrectionFunc func_correction,
65                         void *userdata,
66                         float epsilon,
67                         int max_iterations,
68                         bool trace,
69                         const float x_init[3],
70                         float result[3]);
71
72 #ifdef BLI_MATH_GCC_WARN_PRAGMA
73 #  pragma GCC diagnostic pop
74 #endif
75
76 #ifdef __cplusplus
77 }
78 #endif
79
80 #endif /* __BLI_MATH_SOLVERS_H__ */