228a03bff0d8f6f3298ab05818263777138a9c94
[occt.git] / src / Expr / Expr_ArcSine.cxx
1 //static const char* sccsid = "@(#)Expr_ArcSine.cxx     3.2 95/01/10"; // Do not delete this line. Used by sccs.
2 // Copyright:   Matra-Datavision 1991
3 // File:        Expr_ArcSine.cxx
4 // Created:     Mon May 27 11:31:25 1991
5 // Author:      Arnaud BOUZY
6 //              <adn>
7
8 #include <Expr_ArcSine.ixx>
9 #include <Expr_NumericValue.hxx>
10 #include <Expr_Sine.hxx>
11 #include <Expr_Division.hxx>
12 #include <Expr_Square.hxx>
13 #include <Expr_Difference.hxx>
14 #include <Expr_SquareRoot.hxx>
15 #include <Expr_Operators.hxx>
16 #include <Expr.hxx>
17
18 Expr_ArcSine::Expr_ArcSine (const Handle(Expr_GeneralExpression)& exp)
19 {
20   CreateOperand(exp);
21 }
22
23 Handle(Expr_GeneralExpression) Expr_ArcSine::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(ASin(valop->GetValue()));
29   }
30   if (op->IsKind(STANDARD_TYPE(Expr_Sine))) {
31     return op->SubExpression(1);
32   }
33   Handle(Expr_ArcSine) me = this;
34   return me;
35 }
36
37 Handle(Expr_GeneralExpression) Expr_ArcSine::Copy () const 
38 {
39   return  new Expr_ArcSine(Expr::CopyShare(Operand()));
40 }
41
42 Standard_Boolean Expr_ArcSine::IsIdentical (const Handle(Expr_GeneralExpression)& Other) const
43 {
44   if (!Other->IsKind(STANDARD_TYPE(Expr_ArcSine))) {
45     return Standard_False;
46   }
47   Handle(Expr_GeneralExpression) op = Operand();
48   return op->IsIdentical(Other->SubExpression(1));
49 }
50
51 Standard_Boolean Expr_ArcSine::IsLinear () const
52 {
53   if (ContainsUnknowns()) {
54     return Standard_False;
55   }
56   return Standard_True;
57 }
58
59 Handle(Expr_GeneralExpression) Expr_ArcSine::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   // 1 - X2
69   Handle(Expr_Difference) thedif = 1.0 - sq->ShallowSimplified(); 
70
71   // sqrt(1-X2)
72   Handle(Expr_SquareRoot) theroot = new Expr_SquareRoot(thedif->ShallowSimplified());
73
74   // ArcSine'(F(X)) = F'(X)/sqrt(1-F(X)2) 
75   Handle(Expr_Division) thediv = derop / theroot->ShallowSimplified(); 
76
77   return thediv->ShallowSimplified();
78 }
79
80 Standard_Real Expr_ArcSine::Evaluate(const Expr_Array1OfNamedUnknown& vars, const TColStd_Array1OfReal& vals) const
81 {
82   return ::ASin(Operand()->Evaluate(vars,vals));
83 }
84
85 TCollection_AsciiString Expr_ArcSine::String() const
86 {
87   TCollection_AsciiString str("ASin(");
88   str += Operand()->String();
89   str += ")";
90   return str;
91 }