fedem-solvers/solver/FFaMath_8H_source.html

 // SPDX-FileCopyrightText: 2023 SAP SE

 //

 // SPDX-License-Identifier: Apache-2.0

 //

 // This file is part of FEDEM - https://openfedem.org


 #ifndef FFA_MATH_H

 #define FFA_MATH_H


 #include <cmath>


 #if _MSC_VER > 1300 && _MSC_VER < 1500

 #define hypot _hypot

 #endif


 #ifndef EPS_ZERO

 #define EPS_ZERO 1.0e-15

 #endif


 #ifndef M_PI

 #define M_PI 3.14159265358979323846

 #endif


 #ifndef RAD

 #define RAD(x) (x)*M_PI/180.0

 #endif


 extern double round (double value, int precision);


 #ifdef FFA_DEBUG

 extern double atan3 (double y, double x, const char* func = 0);

 #else

 inline double atan3 (double y, double x, const char* = 0)

 {

   return fabs(y) > EPS_ZERO || fabs(x) > EPS_ZERO ? atan2(y,x) : 0.0;

 }

 #endif


 namespace FFa

 {

   int cubicSolve (double A, double B, double C, double D, double* X);

   int bilinearSolve (const double* A, const double* B, double* X, double* Y);

 }


 #endif

round
double round(double value, int precision)
Rounds off value down to precision significant digits.
Definition: FFaMath.C:31

atan3
double atan3(double y, double x, const char *=0)
Checks for singular rotations before using atan2 library function.
Definition: FFaMath.H:44

EPS_ZERO
#define EPS_ZERO
Division by zero tolerance.
Definition: FFaMath.H:24

FFa
Various math utilities.
Definition: FFaMath.H:52

FFa::bilinearSolve
int bilinearSolve(const double *A, const double *B, double *X, double *Y)
Solves a bi-linear set of equations in X and Y.
Definition: FFaMath.C:158

FFa::cubicSolve
int cubicSolve(double A, double B, double C, double D, double *X)
Solves the cubic equation A*x3 + B*x2 + C*x + D = 0.
Definition: FFaMath.C:61

extctrlsysroutinesmodule::x
integer(ptr), save, private x
Definition: extCtrlSysRoutinesModule.f90:16