src/SketchPlugin/Test/TestFillet.py

   1 """
   2     TestFillet.py
   3     Unit test of SketchPlugin_Fillet class
   4
   5     SketchPlugin_Fillet
   6         static const std::string MY_CONSTRAINT_FILLET_ID("SketchFillet");
   7         data()->addAttribute(SketchPlugin_Constraint::VALUE(), ModelAPI_AttributeDouble::typeId());
   8         data()->addAttribute(SketchPlugin_Constraint::ENTITY_A(), ModelAPI_AttributeRefAttrList::typeId());
   9         data()->addAttribute(SketchPlugin_Constraint::ENTITY_B(), ModelAPI_AttributeRefList::typeId());
  10         data()->addAttribute(SketchPlugin_Constraint::ENTITY_C(), ModelAPI_AttributeRefAttrList::typeId());
  11
  12 """
  13 from GeomDataAPI import *
  14 from ModelAPI import *
  15 import math
  16 from salome.shaper import model
  17
  18 #=========================================================================
  19 # Auxiliary functions
  20 #=========================================================================
  21 TOLERANCE = 1.e-7
  22
  23 def createSketch1(theSketch):
  24     global aEndPoint1, aEndPoint2
  25     # Initialize sketch by three lines with coincident boundaries
  26     allFeatures = []
  27
  28     aSession.startOperation()
  29     # Line1
  30     aSketchLine1 = theSketch.addFeature("SketchLine")
  31     aStartPoint1 = geomDataAPI_Point2D(aSketchLine1.attribute("StartPoint"))
  32     aEndPoint1   = geomDataAPI_Point2D(aSketchLine1.attribute("EndPoint"))
  33     aStartPoint1.setValue(-10., -10.)
  34     aEndPoint1.setValue(-10., 10.)
  35     allFeatures.append(aSketchLine1)
  36     # Line2
  37     aSketchLine2 = theSketch.addFeature("SketchLine")
  38     aStartPoint2 = geomDataAPI_Point2D(aSketchLine2.attribute("StartPoint"))
  39     aEndPoint2   = geomDataAPI_Point2D(aSketchLine2.attribute("EndPoint"))
  40     aStartPoint2.setValue(-10., 10.)
  41     aEndPoint2.setValue(10., 10.)
  42     allFeatures.append(aSketchLine2)
  43     # Line3
  44     aSketchLine3 = theSketch.addFeature("SketchLine")
  45     aStartPoint3 = geomDataAPI_Point2D(aSketchLine3.attribute("StartPoint"))
  46     aEndPoint3   = geomDataAPI_Point2D(aSketchLine3.attribute("EndPoint"))
  47     aStartPoint3.setValue(10., 10.)
  48     aEndPoint3.setValue(10., -10.)
  49     allFeatures.append(aSketchLine3)
  50     # Coincidence1
  51     aCoincidence1 = theSketch.addFeature("SketchConstraintCoincidence")
  52     aCoincidence1.refattr("ConstraintEntityA").setAttr(aEndPoint1)
  53     aCoincidence1.refattr("ConstraintEntityB").setAttr(aStartPoint2)
  54     # Coincidence2
  55     aCoincidence2 = theSketch.addFeature("SketchConstraintCoincidence")
  56     aCoincidence2.refattr("ConstraintEntityA").setAttr(aEndPoint2)
  57     aCoincidence2.refattr("ConstraintEntityB").setAttr(aStartPoint3)
  58
  59     aSession.finishOperation()
  60     return allFeatures
  61
  62
  63 def createSketch2(theSketch):
  64     global aStartPoint1
  65     # Initialize sketch by line and arc with coincident boundary
  66     allFeatures = []
  67
  68     aSession.startOperation()
  69     # Line
  70     aSketchLine = theSketch.addFeature("SketchLine")
  71     aStartPoint1 = geomDataAPI_Point2D(aSketchLine.attribute("StartPoint"))
  72     aEndPoint1   = geomDataAPI_Point2D(aSketchLine.attribute("EndPoint"))
  73     aStartPoint1.setValue(10., 10.)
  74     aEndPoint1.setValue(30., 15.)
  75     allFeatures.append(aSketchLine)
  76     # Arc
  77     aSketchArc = theSketch.addFeature("SketchArc")
  78     aStartPoint2 = geomDataAPI_Point2D(aSketchArc.attribute("ArcStartPoint"))
  79     aEndPoint2   = geomDataAPI_Point2D(aSketchArc.attribute("ArcEndPoint"))
  80     aCenterPoint = geomDataAPI_Point2D(aSketchArc.attribute("ArcCenter"))
  81     aCenterPoint.setValue(20., 10.)
  82     aStartPoint2.setValue(10., 10.)
  83     aEndPoint2.setValue(20., 0.)
  84     allFeatures.append(aSketchArc)
  85     # Coincidence
  86     aCoincidence = theSketch.addFeature("SketchConstraintCoincidence")
  87     aCoincidence.refattr("ConstraintEntityA").setAttr(aStartPoint1)
  88     aCoincidence.refattr("ConstraintEntityB").setAttr(aStartPoint2)
  89
  90     aSession.finishOperation()
  91     return allFeatures
  92
  93 def checkSmoothness(theSketch):
  94     aPtPtCoincidences = getCoincidences(theSketch)
  95     for coinc in aPtPtCoincidences:
  96         aConnectedFeatures = connectedFeatures(coinc)
  97         assert(len(aConnectedFeatures) == 2)
  98         if aConnectedFeatures[0].getKind() == "SketchArc":
  99             if aConnectedFeatures[1].getKind() == "SketchArc":
 100                 checkArcArcSmoothness(aConnectedFeatures[0], aConnectedFeatures[1])
 101             elif aConnectedFeatures[1].getKind() == "SketchLine":
 102                 checkArcLineSmoothness(aConnectedFeatures[0], aConnectedFeatures[1])
 103         elif aConnectedFeatures[0].getKind() == "SketchLine" and aConnectedFeatures[1].getKind() == "SketchArc":
 104             checkArcLineSmoothness(aConnectedFeatures[1], aConnectedFeatures[0])
 105
 106 def checkArcLineSmoothness(theArc, theLine):
 107     aCenter = geomDataAPI_Point2D(theArc.attribute("ArcCenter"))
 108     aDistance = distancePointLine(aCenter, theLine)
 109     aRadius = arcRadius(theArc)
 110     assert(math.fabs(aRadius - aDistance) < TOLERANCE)
 111
 112 def checkArcArcSmoothness(theArc1, theArc2):
 113     aCenter1 = geomDataAPI_Point2D(theArc1.attribute("ArcCenter"))
 114     aCenter2 = geomDataAPI_Point2D(theArc2.attribute("ArcCenter"))
 115     aDistance = distancePointPoint(aCenter1, aCenter2)
 116     aRadius1 = arcRadius(theArc1)
 117     aRadius2 = arcRadius(theArc2)
 118     aRadSum = aRadius1 + aRadius2
 119     aRadDiff = math.fabs(aRadius1 - aRadius2)
 120     assert(math.fabs(aDistance - aRadSum) < TOLERANCE or math.fabs(aDistance - aRadDiff) < TOLERANCE)
 121
 122 def getCoincidences(theSketch):
 123     aCoincidences = []
 124     for anIndex in range(0, theSketch.numberOfSubs()):
 125         aSubFeature = theSketch.subFeature(anIndex)
 126         if aSubFeature.getKind == "SketchConstraintCoincidence":
 127             anEntityA = aSubFeature.refattr("ConstraintEntityA")
 128             anEntityB = aSubFeature.refattr("ConstraintEntityB")
 129             if not anEntityA.isObject() and not anEntityB.isObject():
 130                 aCoincidences.append(aSubFeature)
 131     return aCoincidences
 132
 133 def connectedFeatures(theCoincidence):
 134     anEntityA = aSubFeature.refattr("ConstraintEntityA")
 135     anEntityB = aSubFeature.refattr("ConstraintEntityB")
 136     return [anEntityA.attr().owner(), anEntityB.attr().owner()]
 137
 138 def arcRadius(theArc):
 139     aCenter = geomDataAPI_Point2D(theArc.attribute("ArcCenter"))
 140     aStart = geomDataAPI_Point2D(theArc.attribute("ArcStartPoint"))
 141     return distancePointPoint(aCenter, aStart)
 142
 143 def distancePointPoint(thePoint1, thePoint2):
 144     return math.hypot(thePoint1.x() - thePoint2.x(), thePoint1.y() - thePoint2.y())
 145
 146 def distancePointLine(thePoint, theLine):
 147     aLineStart = geomDataAPI_Point2D(theLine.attribute("StartPoint"))
 148     aLineEnd = geomDataAPI_Point2D(theLine.attribute("EndPoint"))
 149     aLength = distancePointPoint(aLineStart, aLineEnd)
 150
 151     aDir1x, aDir1y = aLineEnd.x() - aLineStart.x(), aLineEnd.y() - aLineStart.y()
 152     aDir2x, aDir2y = thePoint.x() - aLineStart.x(), thePoint.y() - aLineStart.y()
 153     aCross = aDir1x * aDir2y - aDir1y * aDir2x
 154     return math.fabs(aCross) / aLength
 155
 156
 157 #=========================================================================
 158 # Initialization of the test
 159 #=========================================================================
 160
 161 __updated__ = "2015-09-18"
 162
 163 aSession = ModelAPI_Session.get()
 164 aDocument = aSession.moduleDocument()
 165 #=========================================================================
 166 # Creation of a sketch
 167 #=========================================================================
 168 aSession.startOperation()
 169 aSketchCommonFeature = aDocument.addFeature("Sketch")
 170 aSketchFeature = featureToCompositeFeature(aSketchCommonFeature)
 171 origin = geomDataAPI_Point(aSketchFeature.attribute("Origin"))
 172 origin.setValue(0, 0, 0)
 173 dirx = geomDataAPI_Dir(aSketchFeature.attribute("DirX"))
 174 dirx.setValue(1, 0, 0)
 175 norm = geomDataAPI_Dir(aSketchFeature.attribute("Norm"))
 176 norm.setValue(0, 0, 1)
 177 aSession.finishOperation()
 178 #=========================================================================
 179 # Initialize sketch by three connected lines
 180 #=========================================================================
 181 createSketch1(aSketchFeature)
 182 assert (model.dof(aSketchFeature) == 8)
 183 #=========================================================================
 184 # Create the Fillet
 185 #=========================================================================
 186 aSession.startOperation()
 187 aFillet = aSketchFeature.addFeature("SketchFillet")
 188 aFillet.refattr("fillet_point").setAttr(aEndPoint1);
 189 aSession.finishOperation()
 190 aSession.startOperation()
 191 aFillet = aSketchFeature.addFeature("SketchFillet")
 192 aFillet.refattr("fillet_point").setAttr(aEndPoint2);
 193 aSession.finishOperation()
 194 #=========================================================================
 195 # Verify the objects of fillet are created
 196 #=========================================================================
 197 checkSmoothness(aSketchFeature)
 198 assert model.dof(aSketchFeature) == 14, "PlaneGCS limitation: if you see this message, then maybe PlaneGCS has solved DoF for sketch with fillet correctly (expected DoF = 10, observed = {0}".format(model.dof(aSketchFeature))
 199
 200 #=========================================================================
 201 # Create another sketch
 202 #=========================================================================
 203 aSession.startOperation()
 204 aSketchCommonFeature = aDocument.addFeature("Sketch")
 205 aSketchFeature = featureToCompositeFeature(aSketchCommonFeature)
 206 origin = geomDataAPI_Point(aSketchFeature.attribute("Origin"))
 207 origin.setValue(0, 0, 0)
 208 dirx = geomDataAPI_Dir(aSketchFeature.attribute("DirX"))
 209 dirx.setValue(1, 0, 0)
 210 norm = geomDataAPI_Dir(aSketchFeature.attribute("Norm"))
 211 norm.setValue(0, 1, 0)
 212 aSession.finishOperation()
 213 #=========================================================================
 214 # Initialize sketch by line and arc
 215 #=========================================================================
 216 createSketch2(aSketchFeature)
 217 assert (model.dof(aSketchFeature) == 7)
 218 #=========================================================================
 219 # Create the Fillet
 220 #=========================================================================
 221 aSession.startOperation()
 222 aFillet = aSketchFeature.addFeature("SketchFillet")
 223 aFillet.refattr("fillet_point").setAttr(aStartPoint1)
 224 aSession.finishOperation()
 225 #=========================================================================
 226 # Verify the objects of fillet are created
 227 #=========================================================================
 228 checkSmoothness(aSketchFeature)
 229 assert model.dof(aSketchFeature) == 10, "PlaneGCS limitation: if you see this message, then maybe PlaneGCS has solved DoF for sketch with fillet correctly (expected DoF = 8, observed = {0}".format(model.dof(aSketchFeature))
 230 #=========================================================================
 231 # End of test
 232 #=========================================================================
 233
 234 # TODO: Checking of Python dump has been disabled until the Fillet redesigned.
 235 #assert(model.checkPythonDump())