0027035: General fuse algorithm loses face
[occt.git] / src / GeomLib / GeomLib_LogSample.cxx
1 // Created on: 1998-09-23
2 // Created by: Philippe MANGIN
3 // Copyright (c) 1998-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 under
9 // the terms of the GNU Lesser General Public License 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
18 #include <GeomLib_LogSample.hxx>
19 #include <Standard_OutOfRange.hxx>
20
21 GeomLib_LogSample::GeomLib_LogSample(const Standard_Real A,
22                                       const Standard_Real B,
23                                       const Standard_Integer N)
24                    :math_FunctionSample(A, B, N)
25 {
26   myF = A - 1;
27   myexp = Log(B-A)/N;
28 }
29
30 Standard_Real GeomLib_LogSample::GetParameter(const Standard_Integer Index) const
31 {
32   Standard_Integer aN = NbPoints();
33  
34
35   if ((Index >= aN) || (Index <= 1)) {
36     Standard_Real aA, aB;
37     Bounds(aA, aB);
38     if (Index == 1) return aA;
39     else if (Index == aN) return aB;
40     else Standard_OutOfRange::Raise("GeomLib_LogSample::GetParameter");
41   }
42
43   Standard_Real v = myF + Exp(myexp*Index);
44   return v;
45 }