ESS_GUIESS_GUI_DocsESS_GUI_batch_fittingESS_GUI_bumps_abstractionESS_GUI_iss1116ESS_GUI_iss879ESS_GUI_iss959ESS_GUI_openclESS_GUI_orderingESS_GUI_sync_sascalccostrafo411magnetic_scattrelease-4.1.1release-4.1.2release-4.2.2release_4.0.1ticket-1009ticket-1094-headlessticket-1242-2d-resolutionticket-1243ticket-1249ticket885unittest-saveload
Last change
on this file since 45f896f was
f2d6445,
checked in by Mathieu Doucet <doucetm@…>, 17 years ago
|
moving realSpaceModeling to trunk
|
-
Property mode set to
100644
|
File size:
494 bytes
|
Rev | Line | |
---|
[f2d6445] | 1 | #ifndef MATH_UTILS_H |
---|
| 2 | #define MATH_UTILS_H |
---|
| 3 | |
---|
| 4 | /* needed for fabs, sqrt() below */ |
---|
| 5 | #include <cmath> |
---|
| 6 | |
---|
| 7 | |
---|
| 8 | |
---|
| 9 | namespace TNT |
---|
| 10 | { |
---|
| 11 | /** |
---|
| 12 | @returns hypotenuse of real (non-complex) scalars a and b by |
---|
| 13 | avoiding underflow/overflow |
---|
| 14 | using (a * sqrt( 1 + (b/a) * (b/a))), rather than |
---|
| 15 | sqrt(a*a + b*b). |
---|
| 16 | */ |
---|
| 17 | template <class Real> |
---|
| 18 | Real hypot(const Real &a, const Real &b) |
---|
| 19 | { |
---|
| 20 | |
---|
| 21 | if (a== 0) |
---|
| 22 | return abs(b); |
---|
| 23 | else |
---|
| 24 | { |
---|
| 25 | Real c = b/a; |
---|
| 26 | return fabs(a) * sqrt(1 + c*c); |
---|
| 27 | } |
---|
| 28 | } |
---|
| 29 | } /* TNT namespace */ |
---|
| 30 | |
---|
| 31 | |
---|
| 32 | |
---|
| 33 | #endif |
---|
| 34 | /* MATH_UTILS_H */ |
---|
Note: See
TracBrowser
for help on using the repository browser.