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