// Created on: 1993-10-27 // Created by: Jean-LOuis FRENKEL // Copyright (c) 1993-1999 Matra Datavision // Copyright (c) 1999-2012 OPEN CASCADE SAS // // The content of this file is subject to the Open CASCADE Technology Public // License Version 6.5 (the "License"). You may not use the content of this file // except in compliance with the License. Please obtain a copy of the License // at http://www.opencascade.org and read it completely before using this file. // // The Initial Developer of the Original Code is Open CASCADE S.A.S., having its // main offices at: 1, place des Freres Montgolfier, 78280 Guyancourt, France. // // The Original Code and all software distributed under the License is // distributed on an "AS IS" basis, without warranty of any kind, and the // Initial Developer hereby disclaims all such warranties, including without // limitation, any warranties of merchantability, fitness for a particular // purpose or non-infringement. Please see the License for the specific terms // and conditions governing the rights and limitations under the License. #include #include #include #include #include #include #include #include //#include #include #include #include #include #include #include #include #include //======================================================================= //function : IsClosed //purpose : //======================================================================= Standard_Boolean StdPrs_ToolShadedShape::IsClosed(const TopoDS_Shape& aShape) { return aShape.Closed(); } //======================================================================= //function : Triangulation //purpose : //======================================================================= Handle(Poly_Triangulation) StdPrs_ToolShadedShape::Triangulation (const TopoDS_Face& aFace, TopLoc_Location& loc) { return BRep_Tool::Triangulation(aFace, loc); } //======================================================================= //function : Normal //purpose : //======================================================================= void StdPrs_ToolShadedShape::Normal(const TopoDS_Face& aFace, Poly_Connect& pc, TColgp_Array1OfDir& Nor) { const Handle(Poly_Triangulation)& T = pc.Triangulation(); BRepAdaptor_Surface S; Standard_Boolean hasUV = T->HasUVNodes(); Standard_Integer i; TopLoc_Location l; // move to zero TopoDS_Face zeroFace = TopoDS::Face(aFace.Located(TopLoc_Location())); //take in face the surface location //Handle(Geom_Surface) GS = BRep_Tool::Surface(aFace, l); Handle(Geom_Surface) GS = BRep_Tool::Surface(zeroFace); if(T->HasNormals()) { const TColgp_Array1OfPnt& Nodes = T->Nodes(); const TShort_Array1OfShortReal& Normals = T->Normals(); const Standard_ShortReal * arrN = &(Normals.Value(Normals.Lower())); for( i = Nodes.Lower(); i <= Nodes.Upper(); i++) { Standard_Integer in = 3*(i-Nodes.Lower()); gp_Dir N(arrN[in + 0], arrN[in + 1], arrN[in + 2]); Nor(i) = N; } if (aFace.Orientation() == TopAbs_REVERSED) { for( i = Nodes.Lower(); i <= Nodes.Upper(); i++) { Nor.ChangeValue(i).Reverse(); } } } else if (hasUV && !GS.IsNull()) { Standard_Integer nbNormVal = T->NbNodes() * 3; Handle(TShort_HArray1OfShortReal) Normals = new TShort_HArray1OfShortReal(1, nbNormVal); const TColgp_Array1OfPnt2d& UVNodes = T->UVNodes(); Standard_Real Tol = Precision::Confusion(); for (i = UVNodes.Lower(); i <= UVNodes.Upper(); i++) { if(GeomLib::NormEstim(GS, UVNodes(i), Tol, Nor(i)) > 1) { const TColgp_Array1OfPnt& Nodes = T->Nodes(); Standard_Integer n[3]; const Poly_Array1OfTriangle& triangles = T->Triangles(); gp_XYZ eqPlan(0, 0, 0); Standard_Real modmax = 0.; for (pc.Initialize(i); pc.More(); pc.Next()) { triangles(pc.Value()).Get(n[0], n[1], n[2]); gp_XYZ v1(Nodes(n[1]).Coord()-Nodes(n[0]).Coord()); gp_XYZ v2(Nodes(n[2]).Coord()-Nodes(n[1]).Coord()); gp_XYZ vv = v1^v2; Standard_Real mod = vv.Modulus(); if(mod < Tol) continue; eqPlan += vv/mod; } modmax = eqPlan.Modulus(); if(modmax > Tol) Nor(i) = gp_Dir(eqPlan); else Nor(i) = gp_Dir(0., 0., 1.); } Standard_Integer j = (i - UVNodes.Lower()) * 3; Normals->SetValue(j + 1, (Standard_ShortReal)Nor(i).X()); Normals->SetValue(j + 2, (Standard_ShortReal)Nor(i).Y()); Normals->SetValue(j + 3, (Standard_ShortReal)Nor(i).Z()); if (aFace.Orientation() == TopAbs_REVERSED) (Nor(i)).Reverse(); } T->SetNormals(Normals); } else { Standard_Integer nbNormVal = T->NbNodes() * 3; Handle(TShort_HArray1OfShortReal) Normals = new TShort_HArray1OfShortReal(1, nbNormVal); const TColgp_Array1OfPnt& Nodes = T->Nodes(); Standard_Integer n[3]; const Poly_Array1OfTriangle& triangles = T->Triangles(); Standard_Real Tol = Precision::Confusion(); for (i = Nodes.Lower(); i <= Nodes.Upper(); i++) { gp_XYZ eqPlan(0, 0, 0); for (pc.Initialize(i); pc.More(); pc.Next()) { triangles(pc.Value()).Get(n[0], n[1], n[2]); gp_XYZ v1(Nodes(n[1]).Coord()-Nodes(n[0]).Coord()); gp_XYZ v2(Nodes(n[2]).Coord()-Nodes(n[1]).Coord()); gp_XYZ vv = v1^v2; Standard_Real mod = vv.Modulus(); if(mod < Tol) continue; eqPlan += vv/mod; } Standard_Real modmax = eqPlan.Modulus(); if(modmax > Tol) Nor(i) = gp_Dir(eqPlan); else Nor(i) = gp_Dir(0., 0., 1.); Nor(i) = gp_Dir(eqPlan); Standard_Integer j = (i - Nodes.Lower()) * 3; Normals->SetValue(j + 1, (Standard_ShortReal)Nor(i).X()); Normals->SetValue(j + 2, (Standard_ShortReal)Nor(i).Y()); Normals->SetValue(j + 3, (Standard_ShortReal)Nor(i).Z()); if (aFace.Orientation() == TopAbs_REVERSED) (Nor(i)).Reverse(); } T->SetNormals(Normals); } }