include/opencascade/GeomInt_ResConstraintOfMyGradientOfTheComputeLineBezierOfWLApprox.hxx

0001 // Created on: 1995-01-27
0002 // Created by: Jacques GOUSSARD
0003 // Copyright (c) 1995-1999 Matra Datavision
0004 // Copyright (c) 1999-2014 OPEN CASCADE SAS
0005 //
0006 // This file is part of Open CASCADE Technology software library.
0007 //
0008 // This library is free software; you can redistribute it and/or modify it under
0009 // the terms of the GNU Lesser General Public License version 2.1 as published
0010 // by the Free Software Foundation, with special exception defined in the file
0011 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
0012 // distribution for complete text of the license and disclaimer of any warranty.
0013 //
0014 // Alternatively, this file may be used under the terms of Open CASCADE
0015 // commercial license or contractual agreement.
0016
0017 #ifndef _GeomInt_ResConstraintOfMyGradientOfTheComputeLineBezierOfWLApprox_HeaderFile
0018 #define _GeomInt_ResConstraintOfMyGradientOfTheComputeLineBezierOfWLApprox_HeaderFile
0019
0020 #include <Standard.hxx>
0021 #include <Standard_DefineAlloc.hxx>
0022 #include <Standard_Handle.hxx>
0023
0024 #include <Standard_Real.hxx>
0025 #include <math_Matrix.hxx>
0026 #include <math_Vector.hxx>
0027 #include <Standard_Integer.hxx>
0028 #include <TColStd_Array1OfInteger.hxx>
0029 #include <AppParCurves_HArray1OfConstraintCouple.hxx>
0030 class Standard_OutOfRange;
0031 class GeomInt_TheMultiLineOfWLApprox;
0032 class GeomInt_TheMultiLineToolOfWLApprox;
0033 class AppParCurves_MultiCurve;
0034 class math_Matrix;
0035
0036 class GeomInt_ResConstraintOfMyGradientOfTheComputeLineBezierOfWLApprox
0037 {
0038 public:
0039   DEFINE_STANDARD_ALLOC
0040
0041   //! Given a MultiLine SSP with constraints points, this
0042   //! algorithm finds the best curve solution to approximate it.
0043   //! The poles from SCurv issued for example from the least
0044   //! squares are used as a guess solution for the uzawa
0045   //! algorithm. The tolerance used in the Uzawa algorithms
0046   //! is Tolerance.
0047   //! A is the Bernstein matrix associated to the MultiLine
0048   //! and DA is the derivative bernstein matrix.(They can come
0049   //! from an approximation with ParLeastSquare.)
0050   //! The MultiCurve is modified. New MultiPoles are given.
0051   Standard_EXPORT GeomInt_ResConstraintOfMyGradientOfTheComputeLineBezierOfWLApprox(
0052     const GeomInt_TheMultiLineOfWLApprox&                 SSP,
0053     AppParCurves_MultiCurve&                              SCurv,
0054     const Standard_Integer                                FirstPoint,
0055     const Standard_Integer                                LastPoint,
0056     const Handle(AppParCurves_HArray1OfConstraintCouple)& Constraints,
0057     const math_Matrix&                                    Bern,
0058     const math_Matrix&                                    DerivativeBern,
0059     const Standard_Real                                   Tolerance = 1.0e-10);
0060
0061   //! returns True if all has been correctly done.
0062   Standard_EXPORT Standard_Boolean IsDone() const;
0063
0064   //! returns the maximum difference value between the curve
0065   //! and the given points.
0066   Standard_EXPORT Standard_Real Error() const;
0067
0068   Standard_EXPORT const math_Matrix& ConstraintMatrix() const;
0069
0070   //! returns the duale variables of the system.
0071   Standard_EXPORT const math_Vector& Duale() const;
0072
0073   //! Returns the derivative of the constraint matrix.
0074   Standard_EXPORT const math_Matrix& ConstraintDerivative(const GeomInt_TheMultiLineOfWLApprox& SSP,
0075                                                           const math_Vector&     Parameters,
0076                                                           const Standard_Integer Deg,
0077                                                           const math_Matrix&     DA);
0078
0079   //! returns the Inverse of Cont*Transposed(Cont), where
0080   //! Cont is the constraint matrix for the algorithm.
0081   Standard_EXPORT const math_Matrix& InverseMatrix() const;
0082
0083 protected:
0084   //! is used internally to create the fields.
0085   Standard_EXPORT Standard_Integer
0086     NbConstraints(const GeomInt_TheMultiLineOfWLApprox&                 SSP,
0087                   const Standard_Integer                                FirstPoint,
0088                   const Standard_Integer                                LastPoint,
0089                   const Handle(AppParCurves_HArray1OfConstraintCouple)& TheConstraints) const;
0090
0091   //! is internally used for the fields creation.
0092   Standard_EXPORT Standard_Integer NbColumns(const GeomInt_TheMultiLineOfWLApprox& SSP,
0093                                              const Standard_Integer                Deg) const;
0094
0095 private:
0096   Standard_Boolean        Done;
0097   Standard_Real           Err;
0098   math_Matrix             Cont;
0099   math_Matrix             DeCont;
0100   math_Vector             Secont;
0101   math_Matrix             CTCinv;
0102   math_Vector             Vardua;
0103   Standard_Integer        IncPass;
0104   Standard_Integer        IncTan;
0105   Standard_Integer        IncCurv;
0106   TColStd_Array1OfInteger IPas;
0107   TColStd_Array1OfInteger ITan;
0108   TColStd_Array1OfInteger ICurv;
0109 };
0110
0111 #endif // _GeomInt_ResConstraintOfMyGradientOfTheComputeLineBezierOfWLApprox_HeaderFile