DEFINE_STANDARD_ALLOC
+ //! given a MultiLine, this algorithm computes the least
+ //! square resolution using the Householder-QR method.
+ //! If the first and/or the last point is a constraint
+ //! point, the value of the tangency or curvature is
+ //! computed in the resolution.
+ //! NbPol is the number of control points wanted
+ //! for the approximating curves.
+ //! The system to solve is the following:
+ //! A X = B.
+ //! Where A is the Bernstein matrix computed with the
+ //! parameters, B the points coordinates and X the poles
+ //! solutions.
+ //! The matrix A is the same for each coordinate x, y and z
+ //! and is also the same for each MultiLine point because
+ //! they are approximated in parallel(so with the same
+ //! parameter, only the vector B changes).
Standard_EXPORT GeomInt_ParLeastSquareOfMyGradientbisOfTheComputeLineOfWLApprox(const GeomInt_TheMultiLineOfWLApprox& SSP, const Standard_Integer FirstPoint, const Standard_Integer LastPoint, const AppParCurves_Constraint FirstCons, const AppParCurves_Constraint LastCons, const math_Vector& Parameters, const Standard_Integer NbPol);
+ //! Initializes the fields of the object.
Standard_EXPORT GeomInt_ParLeastSquareOfMyGradientbisOfTheComputeLineOfWLApprox(const GeomInt_TheMultiLineOfWLApprox& SSP, const Standard_Integer FirstPoint, const Standard_Integer LastPoint, const AppParCurves_Constraint FirstCons, const AppParCurves_Constraint LastCons, const Standard_Integer NbPol);
+ //! given a MultiLine, this algorithm computes the least
+ //! square resolution using the Householder-QR method.
+ //! If the first and/or the last point is a constraint
+ //! point, the value of the tangency or curvature is
+ //! computed in the resolution.
+ //! Deg is the degree wanted for the approximating curves.
+ //! The system to solve is the following:
+ //! A X = B.
+ //! Where A is the BSpline functions matrix computed with
+ //! <parameters>, B the points coordinates and X the poles
+ //! solutions.
+ //! The matrix A is the same for each coordinate x, y and z
+ //! and is also the same for each MultiLine point because
+ //! they are approximated in parallel(so with the same
+ //! parameter, only the vector B changes).
Standard_EXPORT GeomInt_ParLeastSquareOfMyGradientbisOfTheComputeLineOfWLApprox(const GeomInt_TheMultiLineOfWLApprox& SSP, const TColStd_Array1OfReal& Knots, const TColStd_Array1OfInteger& Mults, const Standard_Integer FirstPoint, const Standard_Integer LastPoint, const AppParCurves_Constraint FirstCons, const AppParCurves_Constraint LastCons, const math_Vector& Parameters, const Standard_Integer NbPol);
+ //! Initializes the fields of the object.
Standard_EXPORT GeomInt_ParLeastSquareOfMyGradientbisOfTheComputeLineOfWLApprox(const GeomInt_TheMultiLineOfWLApprox& SSP, const TColStd_Array1OfReal& Knots, const TColStd_Array1OfInteger& Mults, const Standard_Integer FirstPoint, const Standard_Integer LastPoint, const AppParCurves_Constraint FirstCons, const AppParCurves_Constraint LastCons, const Standard_Integer NbPol);
+ //! Is used after having initialized the fields.
+ //! The case "CurvaturePoint" is not treated in this method.
Standard_EXPORT void Perform (const math_Vector& Parameters);
+ //! Is used after having initialized the fields.
Standard_EXPORT void Perform (const math_Vector& Parameters, const Standard_Real l1, const Standard_Real l2);
+ //! Is used after having initialized the fields.
+ //! <V1t> is the tangent vector at the first point.
+ //! <V2t> is the tangent vector at the last point.
Standard_EXPORT void Perform (const math_Vector& Parameters, const math_Vector& V1t, const math_Vector& V2t, const Standard_Real l1, const Standard_Real l2);
+ //! Is used after having initialized the fields.
+ //! <V1t> is the tangent vector at the first point.
+ //! <V2t> is the tangent vector at the last point.
+ //! <V1c> is the tangent vector at the first point.
+ //! <V2c> is the tangent vector at the last point.
Standard_EXPORT void Perform (const math_Vector& Parameters, const math_Vector& V1t, const math_Vector& V2t, const math_Vector& V1c, const math_Vector& V2c, const Standard_Real l1, const Standard_Real l2);
+ //! returns True if all has been correctly done.
Standard_EXPORT Standard_Boolean IsDone() const;
+ //! returns the result of the approximation, i.e. all the
+ //! Curves.
+ //! An exception is raised if NotDone.
Standard_EXPORT AppParCurves_MultiCurve BezierValue();
+ //! returns the result of the approximation, i.e. all the
+ //! Curves.
+ //! An exception is raised if NotDone.
Standard_EXPORT const AppParCurves_MultiBSpCurve& BSplineValue();
+ //! returns the function matrix used to approximate the
+ //! set.
Standard_EXPORT const math_Matrix& FunctionMatrix() const;
+ //! returns the derivative function matrix used
+ //! to approximate the set.
Standard_EXPORT const math_Matrix& DerivativeFunctionMatrix() const;
+ //! returns the maximum errors between the MultiLine
+ //! and the approximation curves. F is the sum of the square
+ //! distances. Grad is the derivative vector of the
+ //! function F.
Standard_EXPORT void ErrorGradient (math_Vector& Grad, Standard_Real& F, Standard_Real& MaxE3d, Standard_Real& MaxE2d);
+ //! returns the distances between the points of the
+ //! multiline and the approximation curves.
Standard_EXPORT const math_Matrix& Distance();
+ //! returns the maximum errors between the MultiLine
+ //! and the approximation curves. F is the sum of the square
+ //! distances.
Standard_EXPORT void Error (Standard_Real& F, Standard_Real& MaxE3d, Standard_Real& MaxE2d);
+ //! returns the value (P2 - P1)/ V1 if the first point
+ //! was a tangency point.
Standard_EXPORT Standard_Real FirstLambda() const;
+ //! returns the value (PN - PN-1)/ VN if the last point
+ //! was a tangency point.
Standard_EXPORT Standard_Real LastLambda() const;
+ //! returns the matrix of points value.
Standard_EXPORT const math_Matrix& Points() const;
+ //! returns the matrix of resulting control points value.
Standard_EXPORT const math_Matrix& Poles() const;
+ //! Returns the indexes of the first non null values of
+ //! A and DA.
+ //! The values are non null from Index(ieme point) +1
+ //! to Index(ieme point) + degree +1.
Standard_EXPORT const math_IntegerVector& KIndex() const;
protected:
+ //! is used by the constuctors above.
Standard_EXPORT void Init (const GeomInt_TheMultiLineOfWLApprox& SSP, const Standard_Integer FirstPoint, const Standard_Integer LastPoint);
+ //! returns the number of second member columns.
+ //! Is used internally to initialize the fields.
Standard_EXPORT Standard_Integer NbBColumns (const GeomInt_TheMultiLineOfWLApprox& SSP) const;
+ //! returns the first point beeing fitted.
Standard_EXPORT Standard_Integer TheFirstPoint (const AppParCurves_Constraint FirstCons, const Standard_Integer FirstPoint) const;
+ //! returns the last point beeing fitted.
Standard_EXPORT Standard_Integer TheLastPoint (const AppParCurves_Constraint LastCons, const Standard_Integer LastPoint) const;
+ //! Affects the fields in the case of a constraint point.
Standard_EXPORT void Affect (const GeomInt_TheMultiLineOfWLApprox& SSP, const Standard_Integer Index, AppParCurves_Constraint& Cons, math_Vector& Vt, math_Vector& Vc);
Standard_EXPORT void ComputeFunction (const math_Vector& Parameters);
Standard_EXPORT void SearchIndex (math_IntegerVector& Index);
+ //! computes internal matrixes for the resolution
Standard_EXPORT void MakeTAA (math_Vector& TheA, math_Vector& TheB);
+ //! computes internal matrixes for the resolution
Standard_EXPORT void MakeTAA (math_Vector& TheA);
+ //! computes internal matrixes for the resolution
Standard_EXPORT void MakeTAA (math_Vector& TheA, math_Matrix& TheB);