src/Standard/Standard_math.cxx

   1 // Copyright (c) 1998-1999 Matra Datavision
   2 // Copyright (c) 1999-2012 OPEN CASCADE SAS
   3 //
   4 // The content of this file is subject to the Open CASCADE Technology Public
   5 // License Version 6.5 (the "License"). You may not use the content of this file
   6 // except in compliance with the License. Please obtain a copy of the License
   7 // at http://www.opencascade.org and read it completely before using this file.
   8 //
   9 // The Initial Developer of the Original Code is Open CASCADE S.A.S., having its
  10 // main offices at: 1, place des Freres Montgolfier, 78280 Guyancourt, France.
  11 //
  12 // The Original Code and all software distributed under the License is
  13 // distributed on an "AS IS" basis, without warranty of any kind, and the
  14 // Initial Developer hereby disclaims all such warranties, including without
  15 // limitation, any warranties of merchantability, fitness for a particular
  16 // purpose or non-infringement. Please see the License for the specific terms
  17 // and conditions governing the rights and limitations under the License.
  18
  19 // -------------------------------------------------------------------------
  20 // Definition des fonctions mathematiques hyperboliques: acosh, asinh, atanh
  21 // pour Windows/NT uniquement.
  22 // -------------------------------------------------------------------------
  23 #ifdef WNT
  24 #include <math.h>
  25 #include <Standard_math.hxx>
  26
  27 #ifndef _Standard_MacroHeaderFile
  28 # include <Standard_Macro.hxx>
  29 #endif  // _Standard_MacroHeaderFile
  30
  31 Standard_EXPORT double  __cdecl acosh( double X)
  32 {
  33         double  res;
  34         res = log(X + sqrt(X * X - 1));
  35         return res;
  36 };
  37 Standard_EXPORT double  __cdecl  asinh( double X)
  38 {
  39         double  res;
  40 //  Modified by Sergey KHROMOV - Mon Nov 11 16:27:11 2002 Begin
  41 //  Correction of the formula to avoid numerical problems.
  42 //      res = log(X + sqrt(X * X + 1));
  43         if (X > 0.)
  44         res = log(X + sqrt(X * X + 1));
  45         else
  46           res = -log(sqrt(X * X + 1) - X);
  47 //  Modified by Sergey KHROMOV - Mon Nov 11 16:27:13 2002 End
  48         return res;
  49 };
  50 Standard_EXPORT double __cdecl  atanh( double X)
  51 {
  52         double res;
  53         res = log((1 + X) / (1 - X)) / 2;
  54         return res;
  55 };
  56 #endif