0024428: Implementation of LGPL license
[occt.git] / src / CSLib / CSLib_NormalPolyDef.cxx
1 // Created on: 1996-08-23
2 // Created by: Benoit TANNIOU
3 // Copyright (c) 1996-1999 Matra Datavision
4 // Copyright (c) 1999-2014 OPEN CASCADE SAS
5 //
6 // This file is part of Open CASCADE Technology software library.
7 //
8 // This library is free software; you can redistribute it and / or modify it
9 // under the terms of the GNU Lesser General Public version 2.1 as published
10 // by the Free Software Foundation, with special exception defined in the file
11 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
12 // distribution for complete text of the license and disclaimer of any warranty.
13 //
14 // Alternatively, this file may be used under the terms of Open CASCADE
15 // commercial license or contractual agreement.
16
17 #include <CSLib_NormalPolyDef.ixx>
18 #include <PLib.hxx>
19
20 //=============================================================================
21 CSLib_NormalPolyDef::CSLib_NormalPolyDef(const Standard_Integer k0,
22                                                 const TColStd_Array1OfReal& li)
23 //=============================================================================
24  :myTABli(0,k0)
25 {
26   myK0=k0;
27   for(Standard_Integer i=0;i<=k0;i++)
28     myTABli(i)=li(i);
29 }
30
31 //=============================================================================
32 Standard_Boolean CSLib_NormalPolyDef::Value(const Standard_Real X,
33                                                  Standard_Real& F)
34 //=============================================================================
35 {
36   F=0.0;
37   Standard_Real co,si;
38   co=cos(X);
39   si=sin(X);
40
41   for(Standard_Integer i=0;i<=myK0;i++){
42     F=F+PLib::Bin(myK0,i)*pow(co,i)*pow(si,(myK0-i))*myTABli(i);
43   }
44   return Standard_True;
45 }
46
47 //=============================================================================
48 Standard_Boolean CSLib_NormalPolyDef::Derivative(const Standard_Real X,
49                                                       Standard_Real& D)
50 //=============================================================================
51 {
52   D=0.0;
53   Standard_Real co,si;
54   co=cos(X);
55   si=sin(X);
56   for(Standard_Integer i=0;i<=myK0;i++){
57     D=D+PLib::Bin(myK0,i)*pow(co,(i-1))*pow(si,(myK0-i-1))*(myK0*co*co-i);
58   }
59   return Standard_True;
60 }
61
62 //=============================================================================
63 Standard_Boolean CSLib_NormalPolyDef::Values(const Standard_Real X,
64                                                   Standard_Real& F,
65                                                   Standard_Real& D)
66 //=============================================================================
67 {
68   F=0;
69   D=0;
70   Standard_Real co,si;
71   co=cos(X);
72   si=sin(X);
73   for(Standard_Integer i=0;i<=myK0;i++){
74     F=F+PLib::Bin(myK0,i)*pow(co,i)*pow(si,(myK0-i))*myTABli(i);
75     D=D+PLib::Bin(myK0,i)*pow(co,(i-1))*
76         pow(si,(myK0-i-1))*(myK0*co*co-i)*myTABli(i);
77    }
78   return Standard_True;
79 }
80
81
82
83
84