[occt.git] / src / StdPrs / StdPrs_ToolTriangulatedShape.cxx

// Created on: 1993-10-27
// Created by: Jean-LOuis FRENKEL
// Copyright (c) 1993-1999 Matra Datavision
// Copyright (c) 1999-2014 OPEN CASCADE SAS
//
// This file is part of Open CASCADE Technology software library.
//
// This library is free software; you can redistribute it and/or modify it under
// the terms of the GNU Lesser General Public License version 2.1 as published
// by the Free Software Foundation, with special exception defined in the file
// OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
// distribution for complete text of the license and disclaimer of any warranty.
//
// Alternatively, this file may be used under the terms of Open CASCADE
// commercial license or contractual agreement.

#include <StdPrs_ToolTriangulatedShape.hxx>

#include <BRepMesh_DiscretFactory.hxx>
#include <BRepMesh_DiscretRoot.hxx>
#include <BRepTools.hxx>
#include <BRep_Tool.hxx>
#include <GeomAbs_SurfaceType.hxx>
#include <GeomLib.hxx>
#include <gp_XYZ.hxx>
#include <Poly.hxx>
#include <Poly_Connect.hxx>
#include <Poly_Triangulation.hxx>
#include <Precision.hxx>
#include <Prs3d.hxx>
#include <Prs3d_Drawer.hxx>
#include <TColgp_Array1OfPnt.hxx>
#include <TColgp_Array1OfPnt2d.hxx>
#include <TopAbs_Orientation.hxx>
#include <TopLoc_Location.hxx>
#include <TShort_HArray1OfShortReal.hxx>
#include <TShort_Array1OfShortReal.hxx>
#include <TopExp_Explorer.hxx>
#include <TopoDS.hxx>
#include <TopoDS_Face.hxx>

//=======================================================================
//function : IsTriangulated
//purpose  :
//=======================================================================
Standard_Boolean StdPrs_ToolTriangulatedShape::IsTriangulated (const TopoDS_Shape& theShape)
{
  TopLoc_Location aLocDummy;
  for (TopExp_Explorer aFaceIter (theShape, TopAbs_FACE); aFaceIter.More(); aFaceIter.Next())
  {
    const TopoDS_Face&                aFace = TopoDS::Face (aFaceIter.Current());
    const Handle(Poly_Triangulation)& aTri  = BRep_Tool::Triangulation (aFace, aLocDummy);
    if (aTri.IsNull())
    {
      return Standard_False;
    }
  }
  return Standard_True;
}

//=======================================================================
//function : IsClosed
//purpose  :
//=======================================================================
Standard_Boolean StdPrs_ToolTriangulatedShape::IsClosed (const TopoDS_Shape& theShape)
{
  if (theShape.IsNull())
  {
    return Standard_True;
  }

  switch (theShape.ShapeType())
  {
    case TopAbs_COMPOUND:
    case TopAbs_COMPSOLID:
    default:
    {
      // check that compound consists of closed solids
      for (TopoDS_Iterator anIter (theShape); anIter.More(); anIter.Next())
      {
        const TopoDS_Shape& aShape = anIter.Value();
        if (!IsClosed (aShape))
        {
          return Standard_False;
        }
      }
      return Standard_True;
    }
    case TopAbs_SOLID:
    {
      // Check for non-manifold topology first of all:
      // have to use BRep_Tool::IsClosed() because it checks the face connectivity
      // inside the shape
      if (!BRep_Tool::IsClosed (theShape))
        return Standard_False;

      for (TopoDS_Iterator anIter (theShape); anIter.More(); anIter.Next())
      {
        const TopoDS_Shape& aShape = anIter.Value();
        if (aShape.IsNull())
        {
          continue;
        }

        if (aShape.ShapeType() == TopAbs_FACE)
        {
          // invalid solid
          return Standard_False;
        }
        else if (!IsTriangulated (aShape))
        {
          // mesh contains holes
          return Standard_False;
        }
      }
      return Standard_True;
    }
    case TopAbs_SHELL:
    case TopAbs_FACE:
    {
      // free faces / shell are not allowed
      return Standard_False;
    }
    case TopAbs_WIRE:
    case TopAbs_EDGE:
    case TopAbs_VERTEX:
    {
      // ignore
      return Standard_True;
    }
  }
}

//=======================================================================
//function : ComputeNormals
//purpose  :
//=======================================================================
void StdPrs_ToolTriangulatedShape::ComputeNormals (const TopoDS_Face& theFace,
                                                   const Handle(Poly_Triangulation)& theTris,
                                                   Poly_Connect& thePolyConnect)
{
  if (theTris.IsNull()
   || theTris->HasNormals())
  {
    return;
  }

  // take in face the surface location
  const TopoDS_Face    aZeroFace = TopoDS::Face (theFace.Located (TopLoc_Location()));
  Handle(Geom_Surface) aSurf     = BRep_Tool::Surface (aZeroFace);
  const Poly_Array1OfTriangle& aTriangles = theTris->Triangles();
  if (!theTris->HasUVNodes() || aSurf.IsNull())
  {
    // compute normals by averaging triangulation normals sharing the same vertex
    Poly::ComputeNormals (theTris);
    return;
  }

  const Standard_Real aTol = Precision::Confusion();
  Handle(TShort_HArray1OfShortReal) aNormals = new TShort_HArray1OfShortReal (1, theTris->NbNodes() * 3);
  const TColgp_Array1OfPnt2d& aNodesUV = theTris->UVNodes();
  Standard_Integer aTri[3];
  const TColgp_Array1OfPnt& aNodes = theTris->Nodes();
  gp_Dir aNorm;
  for (Standard_Integer aNodeIter = aNodes.Lower(); aNodeIter <= aNodes.Upper(); ++aNodeIter)
  {
    // try to retrieve normal from real surface first, when UV coordinates are available
    if (GeomLib::NormEstim (aSurf, aNodesUV.Value (aNodeIter), aTol, aNorm) > 1)
    {
      if (thePolyConnect.Triangulation() != theTris)
      {
        thePolyConnect.Load (theTris);
      }

      // compute flat normals
      gp_XYZ eqPlan (0.0, 0.0, 0.0);
      for (thePolyConnect.Initialize (aNodeIter); thePolyConnect.More(); thePolyConnect.Next())
      {
        aTriangles (thePolyConnect.Value()).Get (aTri[0], aTri[1], aTri[2]);
        const gp_XYZ v1 (aNodes (aTri[1]).Coord() - aNodes (aTri[0]).Coord());
        const gp_XYZ v2 (aNodes (aTri[2]).Coord() - aNodes (aTri[1]).Coord());
        const gp_XYZ vv = v1 ^ v2;
        const Standard_Real aMod = vv.Modulus();
        if (aMod >= aTol)
        {
          eqPlan += vv / aMod;
        }
      }
      const Standard_Real aModMax = eqPlan.Modulus();
      aNorm = (aModMax > aTol) ? gp_Dir (eqPlan) : gp::DZ();
    }

    const Standard_Integer anId = (aNodeIter - aNodes.Lower()) * 3;
    aNormals->SetValue (anId + 1, (Standard_ShortReal )aNorm.X());
    aNormals->SetValue (anId + 2, (Standard_ShortReal )aNorm.Y());
    aNormals->SetValue (anId + 3, (Standard_ShortReal )aNorm.Z());
  }
  theTris->SetNormals (aNormals);
}

//=======================================================================
//function : Normal
//purpose  :
//=======================================================================
void StdPrs_ToolTriangulatedShape::Normal (const TopoDS_Face&  theFace,
                                           Poly_Connect&       thePolyConnect,
                                           TColgp_Array1OfDir& theNormals)
{
  const Handle(Poly_Triangulation)& aPolyTri = thePolyConnect.Triangulation();
  if (!aPolyTri->HasNormals())
  {
    ComputeNormals (theFace, aPolyTri, thePolyConnect);
  }

  const TColgp_Array1OfPnt&       aNodes   = aPolyTri->Nodes();
  const TShort_Array1OfShortReal& aNormals = aPolyTri->Normals();
  const Standard_ShortReal*       aNormArr = &aNormals.First();
  for (Standard_Integer aNodeIter = aNodes.Lower(); aNodeIter <= aNodes.Upper(); ++aNodeIter)
  {
    const Standard_Integer anId = 3 * (aNodeIter - aNodes.Lower());
    const gp_Dir aNorm (aNormArr[anId + 0],
                        aNormArr[anId + 1],
                        aNormArr[anId + 2]);
    theNormals (aNodeIter) = aNorm;
  }

  if (theFace.Orientation() == TopAbs_REVERSED)
  {
    for (Standard_Integer aNodeIter = aNodes.Lower(); aNodeIter <= aNodes.Upper(); ++aNodeIter)
    {
      theNormals.ChangeValue (aNodeIter).Reverse();
    }
  }
}

//=======================================================================
//function : IsTessellated
//purpose  :
//=======================================================================
Standard_Boolean StdPrs_ToolTriangulatedShape::IsTessellated (const TopoDS_Shape&         theShape,
                                                              const Handle(Prs3d_Drawer)& theDrawer)
{
  return BRepTools::Triangulation (theShape, Prs3d::GetDeflection (theShape, theDrawer));
}

// =======================================================================
// function : Tessellate
// purpose  :
// =======================================================================
Standard_Boolean StdPrs_ToolTriangulatedShape::Tessellate (const TopoDS_Shape&         theShape,
                                                           const Handle(Prs3d_Drawer)& theDrawer)
{
  Standard_Boolean wasRecomputed = Standard_False;
  // Check if it is possible to avoid unnecessary recomputation of shape triangulation
  if (IsTessellated (theShape, theDrawer))
  {
    return wasRecomputed;
  }

  Standard_Real aDeflection = Prs3d::GetDeflection (theShape, theDrawer);

  // retrieve meshing tool from Factory
  Handle(BRepMesh_DiscretRoot) aMeshAlgo = BRepMesh_DiscretFactory::Get().Discret (theShape,
                                                                                   aDeflection,
                                                                                   theDrawer->HLRAngle());
  if (!aMeshAlgo.IsNull())
  {
    aMeshAlgo->Perform();
    wasRecomputed = Standard_True;
  }

  return wasRecomputed;
}

// =======================================================================
// function : ClearOnOwnDeflectionChange
// purpose  :
// =======================================================================
void StdPrs_ToolTriangulatedShape::ClearOnOwnDeflectionChange (const TopoDS_Shape&         theShape,
                                                               const Handle(Prs3d_Drawer)& theDrawer,
                                                               const Standard_Boolean      theToResetCoeff)
{
  if (!theDrawer->IsAutoTriangulation()
    || theShape.IsNull())
  {
    return;
  }

  const Standard_Boolean isOwnDeviationAngle       = theDrawer->HasOwnDeviationAngle();
  const Standard_Boolean isOwnDeviationCoefficient = theDrawer->HasOwnDeviationCoefficient();
  const Standard_Real anAngleNew  = theDrawer->DeviationAngle();
  const Standard_Real anAnglePrev = theDrawer->PreviousDeviationAngle();
  const Standard_Real aCoeffNew   = theDrawer->DeviationCoefficient();
  const Standard_Real aCoeffPrev  = theDrawer->PreviousDeviationCoefficient();
  if ((!isOwnDeviationAngle       || Abs (anAngleNew - anAnglePrev) <= Precision::Angular())
   && (!isOwnDeviationCoefficient || Abs (aCoeffNew  - aCoeffPrev)  <= Precision::Confusion()))
  {
    return;
  }

  BRepTools::Clean (theShape);
  if (theToResetCoeff)
  {
    theDrawer->UpdatePreviousDeviationAngle();
    theDrawer->UpdatePreviousDeviationCoefficient();
  }
}
Commit	Line	Data
b311480e	1	// Created on: 1993-10-27
	2	// Created by: Jean-LOuis FRENKEL
	3	// Copyright (c) 1993-1999 Matra Datavision
973c2be1	4	// Copyright (c) 1999-2014 OPEN CASCADE SAS
b311480e	5	//
973c2be1	6	// This file is part of Open CASCADE Technology software library.
b311480e	7	//
d5f74e42	8	// This library is free software; you can redistribute it and/or modify it under
d5f74e42	9	// the terms of the GNU Lesser General Public License version 2.1 as published
973c2be1	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.
b311480e	13	//
973c2be1	14	// Alternatively, this file may be used under the terms of Open CASCADE
973c2be1	15	// commercial license or contractual agreement.
b311480e	16
5ad8c033	17	#include <StdPrs_ToolTriangulatedShape.hxx>
9447f912	18
5ad8c033	19	#include <BRepMesh_DiscretFactory.hxx>
	20	#include <BRepMesh_DiscretRoot.hxx>
	21	#include <BRepTools.hxx>
9447f912	22	#include <BRep_Tool.hxx>
9447f912	23	#include <GeomAbs_SurfaceType.hxx>
9447f912	24	#include <GeomLib.hxx>
5ad8c033	25	#include <gp_XYZ.hxx>
dc2cc135	26	#include <Poly.hxx>
9447f912	27	#include <Poly_Connect.hxx>
7fd59977	28	#include <Poly_Triangulation.hxx>
9447f912	29	#include <Precision.hxx>
5ad8c033	30	#include <Prs3d.hxx>
5ad8c033	31	#include <Prs3d_Drawer.hxx>
7fd59977	32	#include <TColgp_Array1OfPnt.hxx>
7fd59977	33	#include <TColgp_Array1OfPnt2d.hxx>
7fd59977	34	#include <TopAbs_Orientation.hxx>
7fd59977	35	#include <TopLoc_Location.hxx>
	36	#include <TShort_HArray1OfShortReal.hxx>
	37	#include <TShort_Array1OfShortReal.hxx>
fc9b36d6	38	#include <TopExp_Explorer.hxx>
fc9b36d6	39	#include <TopoDS.hxx>
5ad8c033	40	#include <TopoDS_Face.hxx>
7fd59977	41
4769a395	42	//=======================================================================
5ad8c033	43	//function : IsTriangulated
4769a395	44	//purpose :
4769a395	45	//=======================================================================
5ad8c033	46	Standard_Boolean StdPrs_ToolTriangulatedShape::IsTriangulated (const TopoDS_Shape& theShape)
3b1817a9	47	{
4769a395	48	TopLoc_Location aLocDummy;
4769a395	49	for (TopExp_Explorer aFaceIter (theShape, TopAbs_FACE); aFaceIter.More(); aFaceIter.Next())
3b1817a9	50	{
4769a395	51	const TopoDS_Face& aFace = TopoDS::Face (aFaceIter.Current());
	52	const Handle(Poly_Triangulation)& aTri = BRep_Tool::Triangulation (aFace, aLocDummy);
	53	if (aTri.IsNull())
3b1817a9	54	{
4769a395	55	return Standard_False;
3b1817a9	56	}
3b1817a9	57	}
4769a395	58	return Standard_True;
3b1817a9	59	}
3b1817a9	60
7fd59977	61	//=======================================================================
7fd59977	62	//function : IsClosed
9447f912	63	//purpose :
7fd59977	64	//=======================================================================
5ad8c033	65	Standard_Boolean StdPrs_ToolTriangulatedShape::IsClosed (const TopoDS_Shape& theShape)
7fd59977	66	{
3b1817a9	67	if (theShape.IsNull())
	68	{
	69	return Standard_True;
	70	}
	71
	72	switch (theShape.ShapeType())
	73	{
	74	case TopAbs_COMPOUND:
	75	case TopAbs_COMPSOLID:
	76	default:
	77	{
	78	// check that compound consists of closed solids
	79	for (TopoDS_Iterator anIter (theShape); anIter.More(); anIter.Next())
	80	{
	81	const TopoDS_Shape& aShape = anIter.Value();
	82	if (!IsClosed (aShape))
	83	{
	84	return Standard_False;
	85	}
	86	}
	87	return Standard_True;
	88	}
	89	case TopAbs_SOLID:
	90	{
4769a395	91	// Check for non-manifold topology first of all:
	92	// have to use BRep_Tool::IsClosed() because it checks the face connectivity
	93	// inside the shape
	94	if (!BRep_Tool::IsClosed (theShape))
	95	return Standard_False;
	96
3b1817a9	97	for (TopoDS_Iterator anIter (theShape); anIter.More(); anIter.Next())
	98	{
	99	const TopoDS_Shape& aShape = anIter.Value();
	100	if (aShape.IsNull())
	101	{
	102	continue;
	103	}
	104
4769a395	105	if (aShape.ShapeType() == TopAbs_FACE)
3b1817a9	106	{
	107	// invalid solid
	108	return Standard_False;
	109	}
4769a395	110	else if (!IsTriangulated (aShape))
3b1817a9	111	{
	112	// mesh contains holes
	113	return Standard_False;
	114	}
	115	}
	116	return Standard_True;
	117	}
	118	case TopAbs_SHELL:
	119	case TopAbs_FACE:
	120	{
	121	// free faces / shell are not allowed
	122	return Standard_False;
	123	}
	124	case TopAbs_WIRE:
	125	case TopAbs_EDGE:
	126	case TopAbs_VERTEX:
	127	{
	128	// ignore
	129	return Standard_True;
	130	}
	131	}
7fd59977	132	}
7fd59977	133
7fd59977	134	//=======================================================================
450c83ad	135	//function : ComputeNormals
9447f912	136	//purpose :
7fd59977	137	//=======================================================================
450c83ad	138	void StdPrs_ToolTriangulatedShape::ComputeNormals (const TopoDS_Face& theFace,
	139	const Handle(Poly_Triangulation)& theTris,
	140	Poly_Connect& thePolyConnect)
7fd59977	141	{
450c83ad	142	if (theTris.IsNull()
450c83ad	143	\|\| theTris->HasNormals())
9447f912	144	{
9447f912	145	return;
7fd59977	146	}
7fd59977	147
9447f912	148	// take in face the surface location
dc2cc135	149	const TopoDS_Face aZeroFace = TopoDS::Face (theFace.Located (TopLoc_Location()));
dc2cc135	150	Handle(Geom_Surface) aSurf = BRep_Tool::Surface (aZeroFace);
450c83ad	151	const Poly_Array1OfTriangle& aTriangles = theTris->Triangles();
dc2cc135	152	if (!theTris->HasUVNodes() \|\| aSurf.IsNull())
	153	{
	154	// compute normals by averaging triangulation normals sharing the same vertex
	155	Poly::ComputeNormals (theTris);
	156	return;
	157	}
	158
	159	const Standard_Real aTol = Precision::Confusion();
	160	Handle(TShort_HArray1OfShortReal) aNormals = new TShort_HArray1OfShortReal (1, theTris->NbNodes() * 3);
	161	const TColgp_Array1OfPnt2d& aNodesUV = theTris->UVNodes();
9447f912	162	Standard_Integer aTri[3];
450c83ad	163	const TColgp_Array1OfPnt& aNodes = theTris->Nodes();
f4064435	164	gp_Dir aNorm;
9447f912	165	for (Standard_Integer aNodeIter = aNodes.Lower(); aNodeIter <= aNodes.Upper(); ++aNodeIter)
	166	{
	167	// try to retrieve normal from real surface first, when UV coordinates are available
dc2cc135	168	if (GeomLib::NormEstim (aSurf, aNodesUV.Value (aNodeIter), aTol, aNorm) > 1)
9447f912	169	{
450c83ad	170	if (thePolyConnect.Triangulation() != theTris)
	171	{
	172	thePolyConnect.Load (theTris);
	173	}
	174
9447f912	175	// compute flat normals
	176	gp_XYZ eqPlan (0.0, 0.0, 0.0);
	177	for (thePolyConnect.Initialize (aNodeIter); thePolyConnect.More(); thePolyConnect.Next())
	178	{
	179	aTriangles (thePolyConnect.Value()).Get (aTri[0], aTri[1], aTri[2]);
	180	const gp_XYZ v1 (aNodes (aTri[1]).Coord() - aNodes (aTri[0]).Coord());
	181	const gp_XYZ v2 (aNodes (aTri[2]).Coord() - aNodes (aTri[1]).Coord());
	182	const gp_XYZ vv = v1 ^ v2;
	183	const Standard_Real aMod = vv.Modulus();
	184	if (aMod >= aTol)
	185	{
	186	eqPlan += vv / aMod;
	187	}
7fd59977	188	}
9447f912	189	const Standard_Real aModMax = eqPlan.Modulus();
f4064435	190	aNorm = (aModMax > aTol) ? gp_Dir (eqPlan) : gp::DZ();
7fd59977	191	}
7fd59977	192
9447f912	193	const Standard_Integer anId = (aNodeIter - aNodes.Lower()) * 3;
f4064435	194	aNormals->SetValue (anId + 1, (Standard_ShortReal )aNorm.X());
	195	aNormals->SetValue (anId + 2, (Standard_ShortReal )aNorm.Y());
	196	aNormals->SetValue (anId + 3, (Standard_ShortReal )aNorm.Z());
7fd59977	197	}
450c83ad	198	theTris->SetNormals (aNormals);
f4064435	199	}
	200
	201	//=======================================================================
	202	//function : Normal
	203	//purpose :
	204	//=======================================================================
	205	void StdPrs_ToolTriangulatedShape::Normal (const TopoDS_Face& theFace,
	206	Poly_Connect& thePolyConnect,
	207	TColgp_Array1OfDir& theNormals)
	208	{
	209	const Handle(Poly_Triangulation)& aPolyTri = thePolyConnect.Triangulation();
	210	if (!aPolyTri->HasNormals())
	211	{
450c83ad	212	ComputeNormals (theFace, aPolyTri, thePolyConnect);
f4064435	213	}
	214
	215	const TColgp_Array1OfPnt& aNodes = aPolyTri->Nodes();
	216	const TShort_Array1OfShortReal& aNormals = aPolyTri->Normals();
	217	const Standard_ShortReal* aNormArr = &aNormals.First();
	218	for (Standard_Integer aNodeIter = aNodes.Lower(); aNodeIter <= aNodes.Upper(); ++aNodeIter)
	219	{
	220	const Standard_Integer anId = 3 * (aNodeIter - aNodes.Lower());
	221	const gp_Dir aNorm (aNormArr[anId + 0],
	222	aNormArr[anId + 1],
	223	aNormArr[anId + 2]);
	224	theNormals (aNodeIter) = aNorm;
	225	}
7fd59977	226
9447f912	227	if (theFace.Orientation() == TopAbs_REVERSED)
	228	{
	229	for (Standard_Integer aNodeIter = aNodes.Lower(); aNodeIter <= aNodes.Upper(); ++aNodeIter)
	230	{
	231	theNormals.ChangeValue (aNodeIter).Reverse();
7fd59977	232	}
7fd59977	233	}
7fd59977	234	}
5ad8c033	235
	236	//=======================================================================
	237	//function : IsTessellated
	238	//purpose :
	239	//=======================================================================
	240	Standard_Boolean StdPrs_ToolTriangulatedShape::IsTessellated (const TopoDS_Shape& theShape,
	241	const Handle(Prs3d_Drawer)& theDrawer)
	242	{
	243	return BRepTools::Triangulation (theShape, Prs3d::GetDeflection (theShape, theDrawer));
	244	}
	245
	246	// =======================================================================
	247	// function : Tessellate
	248	// purpose :
	249	// =======================================================================
	250	Standard_Boolean StdPrs_ToolTriangulatedShape::Tessellate (const TopoDS_Shape& theShape,
	251	const Handle(Prs3d_Drawer)& theDrawer)
	252	{
	253	Standard_Boolean wasRecomputed = Standard_False;
	254	// Check if it is possible to avoid unnecessary recomputation of shape triangulation
	255	if (IsTessellated (theShape, theDrawer))
	256	{
	257	return wasRecomputed;
	258	}
	259
	260	Standard_Real aDeflection = Prs3d::GetDeflection (theShape, theDrawer);
	261
	262	// retrieve meshing tool from Factory
	263	Handle(BRepMesh_DiscretRoot) aMeshAlgo = BRepMesh_DiscretFactory::Get().Discret (theShape,
	264	aDeflection,
	265	theDrawer->HLRAngle());
	266	if (!aMeshAlgo.IsNull())
	267	{
	268	aMeshAlgo->Perform();
	269	wasRecomputed = Standard_True;
	270	}
	271
	272	return wasRecomputed;
	273	}
83b0f13a	274
	275	// =======================================================================
	276	// function : ClearOnOwnDeflectionChange
	277	// purpose :
	278	// =======================================================================
	279	void StdPrs_ToolTriangulatedShape::ClearOnOwnDeflectionChange (const TopoDS_Shape& theShape,
	280	const Handle(Prs3d_Drawer)& theDrawer,
	281	const Standard_Boolean theToResetCoeff)
	282	{
	283	if (!theDrawer->IsAutoTriangulation()
	284	\|\| theShape.IsNull())
	285	{
	286	return;
	287	}
	288
	289	const Standard_Boolean isOwnDeviationAngle = theDrawer->HasOwnDeviationAngle();
	290	const Standard_Boolean isOwnDeviationCoefficient = theDrawer->HasOwnDeviationCoefficient();
	291	const Standard_Real anAngleNew = theDrawer->DeviationAngle();
	292	const Standard_Real anAnglePrev = theDrawer->PreviousDeviationAngle();
	293	const Standard_Real aCoeffNew = theDrawer->DeviationCoefficient();
	294	const Standard_Real aCoeffPrev = theDrawer->PreviousDeviationCoefficient();
	295	if ((!isOwnDeviationAngle \|\| Abs (anAngleNew - anAnglePrev) <= Precision::Angular())
	296	&& (!isOwnDeviationCoefficient \|\| Abs (aCoeffNew - aCoeffPrev) <= Precision::Confusion()))
	297	{
	298	return;
	299	}
	300
	301	BRepTools::Clean (theShape);
	302	if (theToResetCoeff)
	303	{
	304	theDrawer->UpdatePreviousDeviationAngle();
	305	theDrawer->UpdatePreviousDeviationCoefficient();
	306	}
	307	}