0024624: Lost word in license statement in source files
[occt.git] / src / AppParCurves / AppParCurves_Gradient_BFGS.gxx
1 // Created on: 1999-12-16
2 // Created by: Atelier CAS2000
3 // Copyright (c) 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 // Redefinition de math_BFGS:
18 // ==========================
19
20
21 AppParCurves_Gradient_BFGS::AppParCurves_Gradient_BFGS(math_MultipleVarFunctionWithGradient& F,
22                                                        const math_Vector& StartingPoint,
23                                                        const Standard_Real Tolerance3d,
24                                                        const Standard_Real Tolerance2d,
25                                                        const Standard_Real Eps,
26                                                        const Standard_Integer NbIterations ):
27        math_BFGS(F, Eps, NbIterations, Eps),
28        myTol3d(Tolerance3d),
29        myTol2d(Tolerance2d) 
30 {
31   Perform(F, StartingPoint);
32 }
33
34 Standard_Boolean AppParCurves_Gradient_BFGS::IsSolutionReached(math_MultipleVarFunctionWithGradient& F) const 
35 {
36     AppParCurves_ParFunction *F1 = (AppParCurves_ParFunction*) &F;
37     Standard_Boolean Result, Result2;
38
39     Result = (2.0 * fabs(TheMinimum - PreviousMinimum) <= 
40         1.e-10 * (fabs(TheMinimum) + fabs(PreviousMinimum))+1.e-12);
41     Standard_Real MErr3d = F1->MaxError3d();
42     Standard_Real MErr2d = F1->MaxError2d();
43     
44     Result2 = ((MErr3d <= myTol3d) && (MErr2d <= myTol2d));
45     return (Result || Result2);
46   }
47
48