X-Git-Url: http://git.salome-platform.org/gitweb/?a=blobdiff_plain;f=src%2FGeomAPI%2FGeomAPI_Pln.cpp;h=c52dda9d659e869e48c55bc1c6e19a1aeebe1977;hb=ea593bc59e7e9461f6c4e2afd3f24d621edb1011;hp=5780b3d819df8a23ed80928a37b03453e752cc06;hpb=c1bc28d34120d2206f9db487718630c3da7acde7;p=modules%2Fshaper.git diff --git a/src/GeomAPI/GeomAPI_Pln.cpp b/src/GeomAPI/GeomAPI_Pln.cpp index 5780b3d81..c52dda9d6 100644 --- a/src/GeomAPI/GeomAPI_Pln.cpp +++ b/src/GeomAPI/GeomAPI_Pln.cpp @@ -1,41 +1,169 @@ -// File: GeomAPI_Pln.cpp -// Created: 23 Apr 2014 -// Author: Mikhail PONIKAROV +// Copyright (C) 2014-2019 CEA/DEN, EDF R&D +// +// This library is free software; you can redistribute it and/or +// modify it under the terms of the GNU Lesser General Public +// License as published by the Free Software Foundation; either +// version 2.1 of the License, or (at your option) any later version. +// +// This library is distributed in the hope that it will be useful, +// but WITHOUT ANY WARRANTY; without even the implied warranty of +// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +// Lesser General Public License for more details. +// +// You should have received a copy of the GNU Lesser General Public +// License along with this library; if not, write to the Free Software +// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA +// +// See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com +// -#include +#include +#include #include #include +#include +#include -#include +#include -using namespace std; +#include -GeomAPI_Pln::GeomAPI_Pln(const boost::shared_ptr& thePoint, - const boost::shared_ptr& theNormal) -: GeomAPI_Interface(new gp_Pln(thePoint->impl(), - theNormal->impl())) +GeomAPI_Pln::GeomAPI_Pln(const std::shared_ptr& theAxis) +: GeomAPI_Interface(new gp_Ax3(theAxis->impl())) { } -GeomAPI_Pln::GeomAPI_Pln( - const double theA, const double theB, const double theC, const double theD) -: GeomAPI_Interface(new gp_Pln(theA, theB, theC, theD)) +GeomAPI_Pln::GeomAPI_Pln(const std::shared_ptr& thePoint, + const std::shared_ptr& theNormal) + : GeomAPI_Interface(new gp_Pln(thePoint->impl(), theNormal->impl())) { } -boost::shared_ptr GeomAPI_Pln::location() +GeomAPI_Pln::GeomAPI_Pln(const double theA, const double theB, const double theC, const double theD) + : GeomAPI_Interface(new gp_Pln(theA, theB, theC, theD)) +{ +} + +std::shared_ptr GeomAPI_Pln::location() const { gp_Pnt aLoc = impl().Location(); - return boost::shared_ptr(new GeomAPI_Pnt(aLoc.X(), aLoc.Y(), aLoc.Z())); + return std::shared_ptr(new GeomAPI_Pnt(aLoc.X(), aLoc.Y(), aLoc.Z())); } -boost::shared_ptr GeomAPI_Pln::direction() +std::shared_ptr GeomAPI_Pln::direction() const { const gp_Dir& aDir = impl().Axis().Direction(); - return boost::shared_ptr(new GeomAPI_Dir(aDir.X(), aDir.Y(), aDir.Z())); + return std::shared_ptr(new GeomAPI_Dir(aDir.X(), aDir.Y(), aDir.Z())); +} + +std::shared_ptr GeomAPI_Pln::xDirection() const +{ + const gp_Dir& aDir = impl().XAxis().Direction(); + return std::shared_ptr(new GeomAPI_Dir(aDir.X(), aDir.Y(), aDir.Z())); } void GeomAPI_Pln::coefficients(double& theA, double& theB, double& theC, double& theD) { impl().Coefficients(theA, theB, theC, theD); } + +bool GeomAPI_Pln::isCoincident(const std::shared_ptr thePlane, + const double theTolerance) +{ + if(!thePlane.get()) { + return false; + } + + const gp_Pln& aMyPln = impl(); + const gp_Pln& anOtherPln = thePlane->impl(); + return (aMyPln.Contains(anOtherPln.Location(), theTolerance) && + aMyPln.Axis().IsParallel(anOtherPln.Axis(), theTolerance)); +} + +std::shared_ptr + GeomAPI_Pln::intersect(const std::shared_ptr& theLine) const +{ + std::shared_ptr aLineDir = theLine->direction()->xyz(); + std::shared_ptr aLineLoc = theLine->location()->xyz(); + + std::shared_ptr aNormal = direction()->xyz(); + std::shared_ptr aLocation = location()->xyz(); + + double aDot = aNormal->dot(aLineDir); + if (Abs(aDot) < Precision::SquareConfusion()) + return std::shared_ptr(); + + double aParam = aNormal->dot(aLocation->decreased(aLineLoc)) / aDot; + return std::shared_ptr( + new GeomAPI_Pnt(aLineLoc->added(aLineDir->multiplied(aParam)))); +} + +std::shared_ptr + GeomAPI_Pln::project(const std::shared_ptr& thePoint) const +{ + std::shared_ptr aNormal = direction()->xyz(); + std::shared_ptr aLocation = location()->xyz(); + + std::shared_ptr aVec = thePoint->xyz()->decreased(aLocation); + double aDot = aNormal->dot(aVec); + std::shared_ptr aProjection = + aLocation->added(aVec->decreased(aNormal->multiplied(aDot))); + return std::shared_ptr(new GeomAPI_Pnt(aProjection)); +} + +double GeomAPI_Pln::distance(const std::shared_ptr thePlane) const +{ + const gp_Pln& aMyPln = impl(); + const gp_Pln& anOtherPln = thePlane->impl(); + + return aMyPln.Distance(anOtherPln); +} + +double GeomAPI_Pln::distance(const std::shared_ptr thePoint) const +{ + const gp_Pln& aMyPln = impl(); + const gp_Pnt& aPnt = thePoint->impl(); + return aMyPln.Distance(aPnt); +} + +void GeomAPI_Pln::translate(const std::shared_ptr theDir, double theDist) +{ + gp_Vec aVec(theDir->impl()); + aVec.Normalize(); + aVec.Multiply(theDist); + implPtr()->Translate(aVec); +} + +std::shared_ptr + GeomAPI_Pln::intersect(const std::shared_ptr thePlane) const +{ + std::shared_ptr aRes; + + if(!thePlane.get()) { + return aRes; + } + + const gp_Pln& aMyPln = impl(); + const gp_Pln& anOtherPln = thePlane->impl(); + + IntAna_QuadQuadGeo aQuad(aMyPln, anOtherPln, Precision::Confusion(), Precision::Confusion()); + + if(aQuad.IsDone() != Standard_True) { + return aRes; + } + + if(aQuad.NbSolutions() != 1) { + return aRes; + } + + gp_Lin aLin = aQuad.Line(1); + gp_Pnt aLoc = aLin.Location(); + gp_Dir aDir = aLin.Direction(); + + std::shared_ptr aGeomLoc(new GeomAPI_Pnt(aLoc.X(), aLoc.Y(), aLoc.Z())); + std::shared_ptr aGeomDir(new GeomAPI_Dir(aDir.X(), aDir.Y(), aDir.Z())); + + aRes.reset(new GeomAPI_Lin(aGeomLoc, aGeomDir)); + + return aRes; +}