X-Git-Url: http://git.salome-platform.org/gitweb/?a=blobdiff_plain;ds=sidebyside;f=src%2FGeomAPI%2FGeomAPI_Pln.cpp;h=9ccd1a9508fad214b29e838c3e91efb227fc1a7e;hb=696c11e9e4cb089e1c5496dac79420147d85496a;hp=dfb7e9bd773534d588c2a355cf7188eaf134c743;hpb=cbde248859fb0072f6012907391ea90cfc254574;p=modules%2Fshaper.git diff --git a/src/GeomAPI/GeomAPI_Pln.cpp b/src/GeomAPI/GeomAPI_Pln.cpp index dfb7e9bd7..9ccd1a950 100644 --- a/src/GeomAPI/GeomAPI_Pln.cpp +++ b/src/GeomAPI/GeomAPI_Pln.cpp @@ -1,17 +1,33 @@ -// Copyright (C) 2014-20xx CEA/DEN, EDF R&D +// Copyright (C) 2014-2017 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 +// -// File: GeomAPI_Pln.cpp -// Created: 23 Apr 2014 -// Author: Mikhail PONIKAROV - -#include +#include #include #include #include +#include +#include -#include +#include -using namespace std; +#include GeomAPI_Pln::GeomAPI_Pln(const std::shared_ptr& theAxis) : GeomAPI_Interface(new gp_Ax3(theAxis->impl())) @@ -29,24 +45,31 @@ GeomAPI_Pln::GeomAPI_Pln(const double theA, const double theB, const double theC { } -std::shared_ptr GeomAPI_Pln::location() +std::shared_ptr GeomAPI_Pln::location() const { gp_Pnt aLoc = impl().Location(); return std::shared_ptr(new GeomAPI_Pnt(aLoc.X(), aLoc.Y(), aLoc.Z())); } -std::shared_ptr GeomAPI_Pln::direction() +std::shared_ptr GeomAPI_Pln::direction() const { const gp_Dir& aDir = impl().Axis().Direction(); 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) +bool GeomAPI_Pln::isCoincident(const std::shared_ptr thePlane, + const double theTolerance) { if(!thePlane.get()) { return false; @@ -54,5 +77,94 @@ bool GeomAPI_Pln::isCoincident(const std::shared_ptr thePlane, cons const gp_Pln& aMyPln = impl(); const gp_Pln& anOtherPln = thePlane->impl(); - return (aMyPln.Contains(anOtherPln.Location(), theTolerance) && aMyPln.Axis().IsParallel(anOtherPln.Axis(), theTolerance)); + 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; }