mirror of
https://github.com/Open-Cascade-SAS/OCCT.git
synced 2026-05-18 17:06:14 +08:00
b659a6dc19
Data member myIsDivergent was added to class math_FunctionSetRoot. The member is initialized in method Perform start by Standard_False. The member is changed in the method to Standard_True if an approximation point is located outside of search set. Method IsDivergent was added to class math_FunctionSetRoot. The method returns value of myIsDivergent. Parameter theStopOnDivergent with default value Standard_False was added to a constructor and method Perform of class math_FunctionSetRoot. The parameter shows whether return from method Perform if myIsDivergent became Standard_True. Class Extrema_GenExtCS was optimized for dynamic memory consumption: the class surface values are not stored now. Method Perform of class Extrema_GenExtCS was improved for local extreme search in case of iteration algorithm divergence. Size of the algorithm grid cell for initial approximation is reduced by 2 in the case. The count of the reduces is limited. Build fix in *.cdl. Case of extrusion was carried out of the repeat loop. Division method of the 3-dimensional grid (curve param, surface first param, surface second param) by two was replaced by method of moving the grid. The moving method is about two times faster than division method. Now grid moving method is depending on some average linear sizes of the curve and the surface. In case (theStopOnDivergent && myIsDivergent) method math_FunctionSetRoot::Perform sets Done to Standard_False and discard it's last solution. Test case for the bug fix. Test case correction Initial approximation was corrected for search of maximum in Extrema_GenExtCS::Perform. Test case bugs moddata_3 bug23830 was formatted.
270 lines
8.8 KiB
Plaintext
Executable File
270 lines
8.8 KiB
Plaintext
Executable File
-- Created on: 1991-05-14
|
|
-- Created by: Laurent PAINNOT
|
|
-- Copyright (c) 1991-1999 Matra Datavision
|
|
-- Copyright (c) 1999-2012 OPEN CASCADE SAS
|
|
--
|
|
-- The content of this file is subject to the Open CASCADE Technology Public
|
|
-- License Version 6.5 (the "License"). You may not use the content of this file
|
|
-- except in compliance with the License. Please obtain a copy of the License
|
|
-- at http://www.opencascade.org and read it completely before using this file.
|
|
--
|
|
-- The Initial Developer of the Original Code is Open CASCADE S.A.S., having its
|
|
-- main offices at: 1, place des Freres Montgolfier, 78280 Guyancourt, France.
|
|
--
|
|
-- The Original Code and all software distributed under the License is
|
|
-- distributed on an "AS IS" basis, without warranty of any kind, and the
|
|
-- Initial Developer hereby disclaims all such warranties, including without
|
|
-- limitation, any warranties of merchantability, fitness for a particular
|
|
-- purpose or non-infringement. Please see the License for the specific terms
|
|
-- and conditions governing the rights and limitations under the License.
|
|
|
|
|
|
|
|
|
|
class FunctionSetRoot from math
|
|
---Purpose: The math_FunctionSetRoot class calculates the root
|
|
-- of a set of N functions of M variables (N<M, N=M or N>M). Knowing
|
|
-- an initial guess of the solution and using a minimization algorithm, a search
|
|
-- is made in the Newton direction and then in the Gradient direction if there
|
|
-- is no success in the Newton direction. This algorithm can also be
|
|
-- used for functions minimization. Knowledge of all the partial
|
|
-- derivatives (the Jacobian) is required.
|
|
|
|
|
|
uses Vector from math,
|
|
Matrix from math,
|
|
IntegerVector from math,
|
|
FunctionSetWithDerivatives from math,
|
|
OStream from Standard
|
|
|
|
raises NotDone from StdFail,
|
|
DimensionError from Standard
|
|
|
|
is
|
|
|
|
Create(F: in out FunctionSetWithDerivatives;
|
|
Tolerance : Vector;
|
|
NbIterations: Integer = 100)
|
|
---Purpose: is used in a sub-class to initialize correctly all the fields
|
|
-- of this class.
|
|
-- The range (1, F.NbVariables()) must be especially
|
|
-- respected for all vectors and matrix declarations.
|
|
|
|
returns FunctionSetRoot from math;
|
|
|
|
|
|
|
|
Create(F: in out FunctionSetWithDerivatives;
|
|
NbIterations: Integer = 100)
|
|
---Purpose: is used in a sub-class to initialize correctly all the fields
|
|
-- of this class.
|
|
-- The range (1, F.NbVariables()) must be especially
|
|
-- respected for all vectors and matrix declarations.
|
|
-- The method SetTolerance must be called after this
|
|
-- constructor.
|
|
|
|
returns FunctionSetRoot from math;
|
|
|
|
|
|
|
|
Create(F: in out FunctionSetWithDerivatives; StartingPoint: Vector;
|
|
Tolerance: Vector; NbIterations: Integer = 100)
|
|
---Purpose: is used to improve the root of the function F
|
|
-- from the initial guess StartingPoint.
|
|
-- The maximum number of iterations allowed is given by
|
|
-- NbIterations.
|
|
-- In this case, the solution is found when:
|
|
-- abs(Xi - Xi-1)(j) <= Tolerance(j) for all unknowns.
|
|
|
|
returns FunctionSetRoot from math;
|
|
|
|
|
|
|
|
Create(F: in out FunctionSetWithDerivatives; StartingPoint: Vector;
|
|
Tolerance: Vector; infBound, supBound: Vector;
|
|
NbIterations: Integer = 100; theStopOnDivergent : Boolean from Standard = Standard_False)
|
|
---Purpose: is used to improve the root of the function F
|
|
-- from the initial guess StartingPoint.
|
|
-- The maximum number of iterations allowed is given
|
|
-- by NbIterations.
|
|
-- In this case, the solution is found when:
|
|
-- abs(Xi - Xi-1) <= Tolerance for all unknowns.
|
|
|
|
returns FunctionSetRoot from math;
|
|
|
|
|
|
|
|
Delete(me:out) is virtual;
|
|
---C++: alias "Standard_EXPORT virtual ~math_FunctionSetRoot(){Delete();}"
|
|
|
|
SetTolerance(me: in out; Tolerance: Vector)
|
|
---Purpose: Initializes the tolerance values.
|
|
|
|
is static;
|
|
|
|
|
|
|
|
Perform(me: in out; F: in out FunctionSetWithDerivatives;
|
|
StartingPoint: Vector;
|
|
infBound, supBound: Vector; theStopOnDivergent : Boolean from Standard = Standard_False)
|
|
---Purpose: Improves the root of function F from the initial guess
|
|
-- StartingPoint. infBound and supBound may be given to constrain the solution.
|
|
-- Warning
|
|
-- This method is called when computation of the solution is
|
|
-- not performed by the constructors.
|
|
|
|
is static;
|
|
|
|
|
|
IsSolutionReached(me: in out; F: in out FunctionSetWithDerivatives)
|
|
---Purpose: This routine is called at the end of each iteration
|
|
-- to check if the solution was found. It can be redefined
|
|
-- in a sub-class to implement a specific test to stop the
|
|
-- iterations.
|
|
-- In this case, the solution is found when:
|
|
-- abs(Xi - Xi-1) <= Tolerance for all unknowns.
|
|
|
|
returns Boolean is virtual;
|
|
|
|
|
|
IsDone(me)
|
|
---Purpose:
|
|
-- Returns true if the computations are successful, otherwise returns false.
|
|
---C++: inline
|
|
returns Boolean
|
|
is static;
|
|
|
|
|
|
NbIterations(me)
|
|
---Purpose: Returns the number of iterations really done
|
|
-- during the computation of the root.
|
|
-- Exception NotDone is raised if the root was not found.
|
|
---C++: inline
|
|
returns Integer
|
|
raises NotDone
|
|
is static;
|
|
|
|
|
|
StateNumber(me)
|
|
---Purpose: returns the stateNumber (as returned by
|
|
-- F.GetStateNumber()) associated to the root found.
|
|
---C++: inline
|
|
returns Integer
|
|
raises NotDone
|
|
is static;
|
|
|
|
|
|
Root(me)
|
|
---Purpose: Returns the value of the root of function F.
|
|
-- Exception NotDone is raised if the root was not found.
|
|
---C++: inline
|
|
---C++: return const&
|
|
returns Vector
|
|
raises NotDone,
|
|
DimensionError
|
|
is static;
|
|
|
|
|
|
Root(me; Root: out Vector)
|
|
---Purpose: Outputs the root vector in Root.
|
|
-- Exception NotDone is raised if the root was not found.
|
|
-- Exception DimensionError is raised if the range of Root
|
|
-- is not equal to the range of the StartingPoint.
|
|
|
|
raises NotDone from StdFail
|
|
is static;
|
|
|
|
|
|
Derivative(me)
|
|
---Purpose: Returns the matrix value of the derivative at the root.
|
|
-- Exception NotDone is raised if the root was not found.
|
|
---C++: inline
|
|
---C++: return const&
|
|
|
|
returns Matrix
|
|
raises NotDone
|
|
is static;
|
|
|
|
|
|
Derivative(me; Der: out Matrix)
|
|
---Purpose: outputs the matrix value of the derivative
|
|
-- at the root in Der.
|
|
-- Exception NotDone is raised if the root was not found.
|
|
-- Exception DimensionError is raised if the column range
|
|
-- of <Der> is not equal to the range of the startingPoint.
|
|
---C++: inline
|
|
|
|
raises NotDone,
|
|
DimensionError
|
|
is static;
|
|
|
|
|
|
|
|
FunctionSetErrors(me)
|
|
---Purpose: returns the vector value of the error done
|
|
-- on the functions at the root.
|
|
-- Exception NotDone is raised if the root was not found.
|
|
---C++: inline
|
|
---C++: return const&
|
|
|
|
returns Vector
|
|
raises NotDone
|
|
is static;
|
|
|
|
|
|
FunctionSetErrors(me; Err: out Vector)
|
|
---Purpose: outputs the vector value of the error done
|
|
-- on the functions at the root in Err.
|
|
-- Exception NotDone is raised if the root was not found.
|
|
-- Exception DimensionError is raised if the range of Err
|
|
-- is not equal to the range of the StartingPoint.
|
|
|
|
raises NotDone,
|
|
DimensionError
|
|
is static;
|
|
|
|
|
|
Dump(me; o: in out OStream)
|
|
---Purpose: Prints on the stream o information on the current state
|
|
-- of the object.
|
|
-- Is used to redefine the operator <<.
|
|
|
|
is static;
|
|
|
|
IsDivergent(me)
|
|
returns Boolean
|
|
is static;
|
|
|
|
|
|
fields
|
|
|
|
Done : Boolean;
|
|
Delta : Vector is protected;
|
|
Sol : Vector is protected;
|
|
DF : Matrix is protected;
|
|
Kount : Integer;
|
|
State : Integer;
|
|
Tol : Vector is protected;
|
|
Itermax : Integer;
|
|
|
|
-- working tables
|
|
|
|
InfBound : Vector;
|
|
SupBound : Vector;
|
|
|
|
SolSave : Vector;
|
|
GH : Vector;
|
|
DH : Vector;
|
|
DHSave : Vector;
|
|
FF : Vector;
|
|
PreviousSolution : Vector;
|
|
Save : Vector;
|
|
Constraints : IntegerVector;
|
|
Temp1 : Vector;
|
|
Temp2 : Vector;
|
|
Temp3 : Vector;
|
|
Temp4 : Vector;
|
|
myIsDivergent : Boolean from Standard;
|
|
|
|
end FunctionSetRoot;
|