OCCT

// Copyright (c) 1995-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 <BRepGProp_Domain.hxx>
#include <BRepGProp_Face.hxx>
#include <BRepGProp_Gauss.hxx>
#include <BRepGProp_Vinert.hxx>
#include <gp_Pln.hxx>
#include <gp_Pnt.hxx>

//=================================================================================================

BRepGProp_Vinert::BRepGProp_Vinert()
    : myEpsilon(0.0)
{
}

//=================================================================================================

BRepGProp_Vinert::BRepGProp_Vinert(BRepGProp_Face& theSurface,
                                   const gp_Pnt&   theLocation,
                                   const double    theEps)
{
  SetLocation(theLocation);
  Perform(theSurface, theEps);
}

//=================================================================================================

BRepGProp_Vinert::BRepGProp_Vinert(BRepGProp_Face&   theSurface,
                                   BRepGProp_Domain& theDomain,
                                   const gp_Pnt&     theLocation,
                                   const double      theEps)
{
  SetLocation(theLocation);
  Perform(theSurface, theDomain, theEps);
}

//=================================================================================================

BRepGProp_Vinert::BRepGProp_Vinert(BRepGProp_Face&   theSurface,
                                   BRepGProp_Domain& theDomain,
                                   const gp_Pnt&     theLocation)
{
  SetLocation(theLocation);
  Perform(theSurface, theDomain);
}

//=================================================================================================

BRepGProp_Vinert::BRepGProp_Vinert(const BRepGProp_Face& theSurface, const gp_Pnt& theLocation)
{
  SetLocation(theLocation);
  Perform(theSurface);
}

//=================================================================================================

BRepGProp_Vinert::BRepGProp_Vinert(BRepGProp_Face& theSurface,
                                   const gp_Pnt&   theOrigin,
                                   const gp_Pnt&   theLocation,
                                   const double    theEps)
{
  SetLocation(theLocation);
  Perform(theSurface, theOrigin, theEps);
}

//=================================================================================================

BRepGProp_Vinert::BRepGProp_Vinert(BRepGProp_Face&   theSurface,
                                   BRepGProp_Domain& theDomain,
                                   const gp_Pnt&     theOrigin,
                                   const gp_Pnt&     theLocation,
                                   const double      theEps)
{
  SetLocation(theLocation);
  Perform(theSurface, theDomain, theOrigin, theEps);
}

//=================================================================================================

BRepGProp_Vinert::BRepGProp_Vinert(const BRepGProp_Face& theSurface,
                                   const gp_Pnt&         theOrigin,
                                   const gp_Pnt&         theLocation)
{
  SetLocation(theLocation);
  Perform(theSurface, theOrigin);
}

//=================================================================================================

BRepGProp_Vinert::BRepGProp_Vinert(BRepGProp_Face&   theSurface,
                                   BRepGProp_Domain& theDomain,
                                   const gp_Pnt&     theOrigin,
                                   const gp_Pnt&     theLocation)
{
  SetLocation(theLocation);
  Perform(theSurface, theDomain, theOrigin);
}

//=================================================================================================

BRepGProp_Vinert::BRepGProp_Vinert(BRepGProp_Face& theSurface,
                                   const gp_Pln&   thePlane,
                                   const gp_Pnt&   theLocation,
                                   const double    theEps)
{
  SetLocation(theLocation);
  Perform(theSurface, thePlane, theEps);
}

//=================================================================================================

BRepGProp_Vinert::BRepGProp_Vinert(BRepGProp_Face&   theSurface,
                                   BRepGProp_Domain& theDomain,
                                   const gp_Pln&     thePlane,
                                   const gp_Pnt&     theLocation,
                                   const double      theEps)
{
  SetLocation(theLocation);
  Perform(theSurface, theDomain, thePlane, theEps);
}

//=================================================================================================

BRepGProp_Vinert::BRepGProp_Vinert(const BRepGProp_Face& theSurface,
                                   const gp_Pln&         thePlane,
                                   const gp_Pnt&         theLocation)
{
  SetLocation(theLocation);
  Perform(theSurface, thePlane);
}

//=================================================================================================

BRepGProp_Vinert::BRepGProp_Vinert(BRepGProp_Face&   theSurface,
                                   BRepGProp_Domain& theDomain,
                                   const gp_Pln&     thePlane,
                                   const gp_Pnt&     theLocation)
{
  SetLocation(theLocation);
  Perform(theSurface, theDomain, thePlane);
}

//=================================================================================================

void BRepGProp_Vinert::SetLocation(const gp_Pnt& theLocation)
{
  loc = theLocation;
}

//=================================================================================================

double BRepGProp_Vinert::Perform(BRepGProp_Face& theSurface, const double theEps)
{
  BRepGProp_Domain anEmptyDomain;
  return Perform(theSurface, anEmptyDomain, theEps);
}

//=================================================================================================

double BRepGProp_Vinert::Perform(BRepGProp_Face&   theSurface,
                                 BRepGProp_Domain& theDomain,
                                 const double      theEps)
{
  const double    aCoeff[] = {0.0, 0.0, 0.0};
  BRepGProp_Gauss aGauss(BRepGProp_Gauss::GaussType::Vinert);

  return myEpsilon =
           aGauss.Compute(theSurface, theDomain, loc, theEps, aCoeff, true, dim, g, inertia);
}

//=================================================================================================

void BRepGProp_Vinert::Perform(const BRepGProp_Face& theSurface)
{
  const double    aCoeff[] = {0.0, 0.0, 0.0};
  BRepGProp_Gauss aGauss(BRepGProp_Gauss::GaussType::Vinert);

  myEpsilon = 1.0;
  aGauss.Compute(theSurface, loc, aCoeff, true, dim, g, inertia);
}

//=================================================================================================

void BRepGProp_Vinert::Perform(BRepGProp_Face& theSurface, BRepGProp_Domain& theDomain)
{
  const double    aCoeff[] = {0.0, 0.0, 0.0};
  BRepGProp_Gauss aGauss(BRepGProp_Gauss::GaussType::Vinert);

  myEpsilon = 1.0;
  aGauss.Compute(theSurface, theDomain, loc, aCoeff, true, dim, g, inertia);
}

//=================================================================================================

double BRepGProp_Vinert::Perform(BRepGProp_Face& theSurface,
                                 const gp_Pnt&   theOrigin,
                                 const double    theEps)
{
  BRepGProp_Domain anEmptyDomain;
  return Perform(theSurface, anEmptyDomain, theOrigin, theEps);
}

//=================================================================================================

double BRepGProp_Vinert::Perform(BRepGProp_Face&   theSurface,
                                 BRepGProp_Domain& theDomain,
                                 const gp_Pnt&     theOrigin,
                                 const double      theEps)
{
  const double aCoeff[] = {theOrigin.X() - loc.X(),
                           theOrigin.Y() - loc.Y(),
                           theOrigin.Z() - loc.Z()};

  BRepGProp_Gauss aGauss(BRepGProp_Gauss::GaussType::Vinert);

  return myEpsilon =
           aGauss.Compute(theSurface, theDomain, loc, theEps, aCoeff, true, dim, g, inertia);
}

//=================================================================================================

void BRepGProp_Vinert::Perform(const BRepGProp_Face& theSurface, const gp_Pnt& theOrigin)
{
  BRepGProp_Gauss aGauss(BRepGProp_Gauss::GaussType::Vinert);
  const double    aCoeff[] = {theOrigin.X() - loc.X(),
                              theOrigin.Y() - loc.Y(),
                              theOrigin.Z() - loc.Z()};

  myEpsilon = 1.0;
  aGauss.Compute(theSurface, loc, aCoeff, true, dim, g, inertia);
}

//=================================================================================================

void BRepGProp_Vinert::Perform(BRepGProp_Face&   theSurface,
                               BRepGProp_Domain& theDomain,
                               const gp_Pnt&     theOrigin)
{
  BRepGProp_Gauss aGauss(BRepGProp_Gauss::GaussType::Vinert);
  const double    aCoeff[] = {theOrigin.X() - loc.X(),
                              theOrigin.Y() - loc.Y(),
                              theOrigin.Z() - loc.Z()};

  myEpsilon = 1.0;
  aGauss.Compute(theSurface, theDomain, loc, aCoeff, true, dim, g, inertia);
}

//=================================================================================================

double BRepGProp_Vinert::Perform(BRepGProp_Face& theSurface,
                                 const gp_Pln&   thePlane,
                                 const double    theEps)
{
  BRepGProp_Domain anEmptyDomain;
  return Perform(theSurface, anEmptyDomain, thePlane, theEps);
}

//=================================================================================================

double BRepGProp_Vinert::Perform(BRepGProp_Face&   theSurface,
                                 BRepGProp_Domain& theDomain,
                                 const gp_Pln&     thePlane,
                                 const double      theEps)
{
  double aCoeff[4];
  thePlane.Coefficients(aCoeff[0], aCoeff[1], aCoeff[2], aCoeff[3]);
  aCoeff[3] = aCoeff[3] - aCoeff[0] * loc.X() - aCoeff[1] * loc.Y() - aCoeff[2] * loc.Z();

  BRepGProp_Gauss aGauss(BRepGProp_Gauss::GaussType::Vinert);

  return myEpsilon =
           aGauss.Compute(theSurface, theDomain, loc, theEps, aCoeff, false, dim, g, inertia);
}

//=================================================================================================

void BRepGProp_Vinert::Perform(const BRepGProp_Face& theSurface, const gp_Pln& thePlane)
{
  BRepGProp_Gauss aGauss(BRepGProp_Gauss::GaussType::Vinert);
  double          aCoeff[4];

  thePlane.Coefficients(aCoeff[0], aCoeff[1], aCoeff[2], aCoeff[3]);

  aCoeff[3] = aCoeff[3] - aCoeff[0] * loc.X() - aCoeff[1] * loc.Y() - aCoeff[2] * loc.Z();

  myEpsilon = 1.0;
  aGauss.Compute(theSurface, loc, aCoeff, false, dim, g, inertia);
}

//=================================================================================================

void BRepGProp_Vinert::Perform(BRepGProp_Face&   theSurface,
                               BRepGProp_Domain& theDomain,
                               const gp_Pln&     thePlane)
{
  BRepGProp_Gauss aGauss(BRepGProp_Gauss::GaussType::Vinert);
  double          aCoeff[4];

  thePlane.Coefficients(aCoeff[0], aCoeff[1], aCoeff[2], aCoeff[3]);

  aCoeff[3] = aCoeff[3] - aCoeff[0] * loc.X() - aCoeff[1] * loc.Y() - aCoeff[2] * loc.Z();

  myEpsilon = 1.0;
  aGauss.Compute(theSurface, theDomain, loc, aCoeff, false, dim, g, inertia);
}

//=================================================================================================

double BRepGProp_Vinert::GetEpsilon()
{
  return myEpsilon;
}