29     const f32 EPS = 2e-4f;
44     NN_MATH_INLINE f32 
45     spow(f32 x)  in spow()
51     NN_MATH_INLINE f32 
52     spow<1>(f32 x)  in spow()
72     SolveEquation2(f32* root, /*f32 a==1,*/ f32 b, f32 c)  in SolveEquation2()
78                 f32 r_c = FSqrt(-c);  in SolveEquation2()
95             f32 A = b / 2;  in SolveEquation2()
96             f32 B = c / spow<2>(A);  in SolveEquation2()
97             f32 D = 1 - B;  in SolveEquation2()
101                 f32 C = -1 - FSqrt(D);  in SolveEquation2()
133     f32 
134     SolveEquation3(/*f32* root, f32 a==1,*/ f32 b, f32 c, f32 d)  in SolveEquation3()
136         f32 q   = (spow<2>(b) - 3 * c) / 9;  in SolveEquation3()
137         f32 r   = (2 * spow<3>(b) - 9 * b * c + 27 * d) / 54;  in SolveEquation3()
138         f32 D   = spow<3>(q) - spow<2>(r);  in SolveEquation3()
139         f32 b_3 = b / 3;  in SolveEquation3()
143             f32 theta   = AcosRad( r / FSqrt(spow<3>(q)) );  in SolveEquation3()
144             f32 theta_3 = theta/3;  in SolveEquation3()
145             f32 r_q2    = -2 * FSqrt(q);  in SolveEquation3()
151             f32 r3_Dr = FCbrt( FSqrt(- D) + FAbs(r) );  in SolveEquation3()
152             f32 xp    = r3_Dr + q / r3_Dr;  in SolveEquation3()
158             f32 xp  = FSqrt(q);  in SolveEquation3()
186 SolveEquation2(f32* root, f32 a, f32 b, f32 c)  in SolveEquation2()
206         f32 A = b / (2 * a);  in SolveEquation2()
207         f32 B = c / (a * A * A);  in SolveEquation2()
208         f32 D = 1 - B;  in SolveEquation2()
212             f32 C = -1 - FSqrt(D);  in SolveEquation2()
234         f32 c_a = - c / a;  in SolveEquation2()
238             f32 r_c_a = FSqrt(c_a);  in SolveEquation2()
272 SolveEquation3(f32* root, f32 a, f32 b, f32 c, f32 d)  in SolveEquation3()
283     f32 q = (spow<2>(b) - 3 * c) / 9;  in SolveEquation3()
284     f32 r = (2 * spow<3>(b) - 9 * b * c + 27 * d) / 54;  in SolveEquation3()
285     f32 D = spow<3>(q) - spow<2>(r);  in SolveEquation3()
287     f32 b_3 = b / 3;  in SolveEquation3()
291         f32 theta   = AcosRad( r / FSqrt(spow<3>(q)) );  in SolveEquation3()
292         f32 theta_3 = theta/3;  in SolveEquation3()
293         f32 r_q2    = -2 * FSqrt(q);  in SolveEquation3()
303         f32 r3_Dr = FCbrt( FSqrt(- D) + FAbs(r) );  in SolveEquation3()
304         f32 xp    = r3_Dr + q / r3_Dr;  in SolveEquation3()
312         f32 xp  = FSqrt(q);  in SolveEquation3()
313         f32 sxp = FCopySign(xp, r);  in SolveEquation3()
340 SolveEquation4(f32* root, f32 a, f32 b, f32 c, f32 d, f32 e)  in SolveEquation4()
345     f32 m, n, y;  in SolveEquation4()
360         f32 p = - 3.f * spow<2>(b) / 8 + c;  in SolveEquation4()
361         f32 q = spow<3>(b) / 8 - b * c / 2 + d;  in SolveEquation4()
362         f32 r = - 3.f * spow<4>(b) / 256 + spow<2>(b) * c / 16 - b * d / 4 + e;  in SolveEquation4()
414     f32 b4    = b / 4;  in SolveEquation4()
418         f32 root01[2];  in SolveEquation4()
429         f32 root23[2];  in SolveEquation4()