src/Plate/Plate_SampledCurveConstraint.cxx

   1 // Created on: 1998-05-18
   2 // Created by: Andre LIEUTIER
   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 <Plate_SampledCurveConstraint.hxx>
  19 #include <Plate_SequenceOfPinpointConstraint.hxx>
  20 #include <Standard_DimensionMismatch.hxx>
  21
  22 static inline Standard_Real B0( Standard_Real t)
  23 {
  24   Standard_Real s = t;
  25   if(s<0.) s = -s;
  26   s = 1. - s;
  27   if(s<0.) s = 0;
  28   return s;
  29 }
  30
  31 Plate_SampledCurveConstraint::Plate_SampledCurveConstraint(const Plate_SequenceOfPinpointConstraint &SOPPC,
  32                                                            const Standard_Integer n)
  33 :myLXYZC(n,SOPPC.Length())
  34 {
  35   Standard_Integer m = SOPPC.Length();
  36
  37   if (n > m)  throw Standard_DimensionMismatch();
  38   for(Standard_Integer index =1; index <= m;index++)
  39     myLXYZC.SetPPC(index,SOPPC(index));
  40
  41   Standard_Real ratio = Standard_Real(n+1) /Standard_Real(m+1);
  42   for (Standard_Integer i=1;i<=n;i++)
  43     for (Standard_Integer j=1;j<=m;j++)
  44       {
  45         myLXYZC.SetCoeff(i,j,B0(ratio*j - i));
  46       }
  47 }