src/Expr/Expr_ArgSinh.cxx

   1 //static const char* sccsid = "@(#)Expr_ArgSinh.cxx     3.2 95/01/10"; // Do not delete this line. Used by sccs.
   2 // Copyright:   Matra-Datavision 1991
   3 // File:        Expr_ArgSinh.cxx
   4 // Created:     Mon May 27 12:10:02 1991
   5 // Author:      Arnaud BOUZY
   6 //              <adn>
   7
   8 #include <Expr_ArgSinh.ixx>
   9 #include <Expr_NumericValue.hxx>
  10 #include <Expr_Sinh.hxx>
  11 #include <Expr_Division.hxx>
  12 #include <Expr_Square.hxx>
  13 #include <Expr_Sum.hxx>
  14 #include <Expr_SquareRoot.hxx>
  15 #include <Expr_Operators.hxx>
  16 #include <Expr.hxx>
  17
  18 Expr_ArgSinh::Expr_ArgSinh (const Handle(Expr_GeneralExpression)& exp)
  19 {
  20   CreateOperand(exp);
  21 }
  22
  23 Handle(Expr_GeneralExpression) Expr_ArgSinh::ShallowSimplified () const
  24 {
  25   Handle(Expr_GeneralExpression) op = Operand();
  26   if (op->IsKind(STANDARD_TYPE(Expr_NumericValue))) {
  27     Handle(Expr_NumericValue) valop = Handle(Expr_NumericValue)::DownCast(op);
  28     return new Expr_NumericValue(ASinh(valop->GetValue()));
  29   }
  30   if (op->IsKind(STANDARD_TYPE(Expr_Sinh))) {
  31     return op->SubExpression(1);
  32   }
  33   Handle(Expr_ArgSinh) me = this;
  34   return me;
  35 }
  36
  37 Handle(Expr_GeneralExpression) Expr_ArgSinh::Copy () const
  38 {
  39   return  new Expr_ArgSinh(Expr::CopyShare(Operand()));
  40 }
  41
  42 Standard_Boolean Expr_ArgSinh::IsIdentical (const Handle(Expr_GeneralExpression)& Other) const
  43 {
  44   if (!Other->IsKind(STANDARD_TYPE(Expr_ArgSinh))) {
  45     return Standard_False;
  46   }
  47   Handle(Expr_GeneralExpression) op = Operand();
  48   return op->IsIdentical(Other->SubExpression(1));
  49 }
  50
  51 Standard_Boolean Expr_ArgSinh::IsLinear () const
  52 {
  53   if (ContainsUnknowns()) {
  54     return Standard_False;
  55   }
  56   return Standard_True;
  57 }
  58
  59 Handle(Expr_GeneralExpression) Expr_ArgSinh::Derivative (const Handle(Expr_NamedUnknown)& X) const
  60 {
  61   if (!Contains(X)) {
  62     return new Expr_NumericValue(0.0);
  63   }
  64   Handle(Expr_GeneralExpression) op = Operand();
  65   Handle(Expr_GeneralExpression) derop = op->Derivative(X);
  66
  67   Handle(Expr_Square) sq = new Expr_Square(Expr::CopyShare(op));
  68   // X2 + 1
  69   Handle(Expr_Sum) thesum = sq->ShallowSimplified() + 1.0;
  70
  71   // sqrt(X2 + 1)
  72   Handle(Expr_SquareRoot) theroot = new Expr_SquareRoot(thesum->ShallowSimplified());
  73
  74   // ArgSinh'(F(X)) = F'(X)/sqrt(F(X)2+1)
  75   Handle(Expr_Division) thediv = derop / theroot->ShallowSimplified();
  76
  77   return thediv->ShallowSimplified();
  78 }
  79
  80 Standard_Real Expr_ArgSinh::Evaluate(const Expr_Array1OfNamedUnknown& vars, const TColStd_Array1OfReal& vals) const
  81 {
  82   Standard_Real val = Operand()->Evaluate(vars,vals);
  83   return ::Log(val + ::Sqrt(::Square(val)+1.0));
  84 }
  85
  86 TCollection_AsciiString Expr_ArgSinh::String() const
  87 {
  88   TCollection_AsciiString str("ASinh(");
  89   str += Operand()->String();
  90   str += ")";
  91   return str;
  92 }