diff --git a/apps/gpu/RotatingGrid/RotatingGrid.cpp b/apps/gpu/RotatingGrid/RotatingGrid.cpp index 628bc2dba508f0e5e30c45a89fd7c3daa4e968aa..8c8dcc06487ad72ab52cfefec5f2ccc305558b4e 100644 --- a/apps/gpu/RotatingGrid/RotatingGrid.cpp +++ b/apps/gpu/RotatingGrid/RotatingGrid.cpp @@ -109,8 +109,8 @@ int main() gridBuilder->addCoarseGrid(-1.5 * L, -0.5 * L, -0.5 * L, 1.5 * L, 0.5 * L, 0.5 * L, dx); - if (rotOrInt == Rot) gridBuilder->addGridRotatingGrid(std::make_shared<Cylinder>(0.0, 0.0, 0.0, 0.3 * L, 1. * L, Cylinder::PrincipalAxis::x)); - if (rotOrInt == Int) gridBuilder->addGrid(std::make_shared<Cylinder>(0.0, 0.0, 0.0, 0.3 * L, 1. * L, Cylinder::PrincipalAxis::x), 1); + if (rotOrInt == Rot) gridBuilder->addGridRotatingGrid(std::make_shared<Cylinder>(0.0, 0.0, 0.0, 0.3 * L, 1. * L, Cylinder::RotationalAxis::x)); + if (rotOrInt == Int) gridBuilder->addGrid(std::make_shared<Cylinder>(0.0, 0.0, 0.0, 0.3 * L, 1. * L, Cylinder::RotationalAxis::x), 1); GridScalingFactory scalingFactory = GridScalingFactory(); scalingFactory.setScalingFactory(GridScalingFactory::GridScaling::ScaleCompressible); diff --git a/src/gpu/GridGenerator/geometries/Cylinder/Cylinder.cpp b/src/gpu/GridGenerator/geometries/Cylinder/Cylinder.cpp index 29dc74bad6b5f3f8d493cf4b09ea4d36f951a7af..40d21f4f9ee545d37641f408134bcf1ba82f8483 100644 --- a/src/gpu/GridGenerator/geometries/Cylinder/Cylinder.cpp +++ b/src/gpu/GridGenerator/geometries/Cylinder/Cylinder.cpp @@ -1,36 +1,36 @@ #include "Cylinder.h" #include <numeric> -Cylinder::Cylinder(double centerX, double centerY, double centerZ, double radius, double height, PrincipalAxis axis) - : center({ centerX, centerY, centerZ }), radius(radius), height(height), principalAxis(axis) +Cylinder::Cylinder(double centerX, double centerY, double centerZ, double radius, double height, RotationalAxis axis) + : center({ centerX, centerY, centerZ }), radius(radius), height(height), rotationalAxis(axis) { } -Cylinder::Cylinder(std::array<double, 3> center, double radius, double height, PrincipalAxis axis) - : center(center), radius(radius), height(height), principalAxis(axis) +Cylinder::Cylinder(std::array<double, 3> center, double radius, double height, RotationalAxis axis) + : center(center), radius(radius), height(height), rotationalAxis(axis) { } SPtr<Object> Cylinder::clone() const { - return std::make_shared<Cylinder>(center, radius, height, principalAxis); + return std::make_shared<Cylinder>(center, radius, height, rotationalAxis); } -double Cylinder::getCentroidCoordinate(PrincipalAxis coordinateDirection) const +double Cylinder::getCentroidCoordinate(RotationalAxis coordinateDirection) const { return center.at(coordinateDirection); } -double Cylinder::getMinimunCoordinate(PrincipalAxis coordinateDirection) const +double Cylinder::getMinimunCoordinate(RotationalAxis coordinateDirection) const { - const auto unitVector = unitVectors.at(principalAxis); + const auto unitVector = unitVectors.at(rotationalAxis); return center.at(coordinateDirection) - 0.5 * height * unitVector.at(coordinateDirection) + radius * (unitVector.at(coordinateDirection) - 1); } -double Cylinder::getMaximumCoordinate(PrincipalAxis coordinateDirection) const +double Cylinder::getMaximumCoordinate(RotationalAxis coordinateDirection) const { - const auto unitVector = unitVectors.at(principalAxis); + const auto unitVector = unitVectors.at(rotationalAxis); return center.at(coordinateDirection) + 0.5 * height * unitVector.at(coordinateDirection) - radius * (unitVector.at(coordinateDirection) - 1); } @@ -90,9 +90,9 @@ double Cylinder::getHeight() const return height; } -Cylinder::PrincipalAxis Cylinder::getPrincipalAxis() const +Cylinder::RotationalAxis Cylinder::getRotationalAxis() const { - return principalAxis; + return rotationalAxis; } bool Cylinder::isInCircle(double delta1, double delta2, double offset) const @@ -110,7 +110,7 @@ bool Cylinder::isPointInObject(const double &x1, const double &x2, const double const double deltaX2 = x2 - center.at(y); const double deltaX3 = x3 - center.at(z); - switch (principalAxis) { + switch (rotationalAxis) { case x: if (deltaX1 > 0.5 * height || deltaX1 < -0.5 * height) return false; return isInCircle(deltaX2, deltaX3, offset); @@ -122,7 +122,7 @@ bool Cylinder::isPointInObject(const double &x1, const double &x2, const double return isInCircle(deltaX1, deltaX2, offset); } - VF_LOG_CRITICAL("Unknown principal axis in Cylinder."); + VF_LOG_CRITICAL("Unknown rotational axis in Cylinder."); return false; } diff --git a/src/gpu/GridGenerator/geometries/Cylinder/Cylinder.h b/src/gpu/GridGenerator/geometries/Cylinder/Cylinder.h index 108d14ef524b9fcab5d238df7ead7e31360703a0..1a3439aae761cfc4858cd2046a7fa18f94ace28b 100644 --- a/src/gpu/GridGenerator/geometries/Cylinder/Cylinder.h +++ b/src/gpu/GridGenerator/geometries/Cylinder/Cylinder.h @@ -42,14 +42,14 @@ class GRIDGENERATOR_EXPORT Cylinder : public Object { public: - enum PrincipalAxis { + enum RotationalAxis { x = 0, y = 1, z = 2, }; - Cylinder(double centerX, double centerY, double centerZ, double radius, double height, PrincipalAxis axis); - Cylinder(std::array<double, 3> center, double radius, double height, PrincipalAxis axis); + Cylinder(double centerX, double centerY, double centerZ, double radius, double height, RotationalAxis axis); + Cylinder(std::array<double, 3> center, double radius, double height, RotationalAxis axis); SPtr<Object> clone() const override; @@ -65,24 +65,24 @@ public: double getRadius() const; double getHeight() const; - PrincipalAxis getPrincipalAxis() const; + RotationalAxis getRotationalAxis() const; bool isPointInObject(const double &x1, const double &x2, const double &x3, const double &minOffset, const double &maxOffset) override; void changeSizeByDelta(double delta) override; private: - double getCentroidCoordinate(PrincipalAxis coordinateDirection) const; - double getMinimunCoordinate(PrincipalAxis coordinateDirection) const; - double getMaximumCoordinate(PrincipalAxis coordinateDirection) const; + double getCentroidCoordinate(RotationalAxis coordinateDirection) const; + double getMinimunCoordinate(RotationalAxis coordinateDirection) const; + double getMaximumCoordinate(RotationalAxis coordinateDirection) const; bool isInCircle(double delta1, double delta2, double offset) const; - const std::map<PrincipalAxis, std::array<double, 3>> unitVectors{ { x, { 1, 0, 0 } }, - { y, { 0, 1, 0 } }, - { z, { 0, 0, 1 } } }; + const std::map<RotationalAxis, std::array<double, 3>> unitVectors{ { x, { 1, 0, 0 } }, + { y, { 0, 1, 0 } }, + { z, { 0, 0, 1 } } }; - PrincipalAxis principalAxis; + RotationalAxis rotationalAxis; const std::array<double, 3> center; double radius; diff --git a/src/gpu/GridGenerator/geometries/Cylinder/CylinderTest.cpp b/src/gpu/GridGenerator/geometries/Cylinder/CylinderTest.cpp index 6c0582f488783adc5fb5678d9b482ba151cbd9a2..cc636a94ac22cb923a8ee29eef0f3358713749a1 100644 --- a/src/gpu/GridGenerator/geometries/Cylinder/CylinderTest.cpp +++ b/src/gpu/GridGenerator/geometries/Cylinder/CylinderTest.cpp @@ -36,7 +36,7 @@ public: std::array<double, 3> center = { 0.1, 0.2, 0.3 }; double radius = 2.0; double height = 8.0; - Cylinder cylinder = Cylinder(center[0], center[1], center[2], radius, height, Cylinder::PrincipalAxis::x); + Cylinder cylinder = Cylinder(center[0], center[1], center[2], radius, height, Cylinder::RotationalAxis::x); }; TEST_F(CylinderTestAxisNormalToX, getCentroid) @@ -97,7 +97,7 @@ protected: std::array<double, 3> center = { 0.1, 0.2, 0.3 }; double radius = 2.0; double height = 8.0; - Cylinder cylinder = Cylinder({ center[0], center[1], center[2] }, radius, height, Cylinder::PrincipalAxis::y); + Cylinder cylinder = Cylinder({ center[0], center[1], center[2] }, radius, height, Cylinder::RotationalAxis::y); }; TEST_F(CylinderTestAxisNormalToY, getCentroid) @@ -158,7 +158,7 @@ protected: std::array<double, 3> center = { 0.1, 0.2, 0.3 }; double radius = 2.0; double height = 8.0; - Cylinder cylinder = Cylinder({ center[0], center[1], center[2] }, radius, height, Cylinder::PrincipalAxis::z); + Cylinder cylinder = Cylinder({ center[0], center[1], center[2] }, radius, height, Cylinder::RotationalAxis::z); }; TEST_F(CylinderTestAxisNormalToZ, getCentroid) diff --git a/src/gpu/GridGenerator/grid/GridBuilder/MultipleGridBuilder.cpp b/src/gpu/GridGenerator/grid/GridBuilder/MultipleGridBuilder.cpp index 58834dabe7431045d39025c234ff67728b8f4817..627f9204effdd5a6c59b88221426503dafda5bd8 100644 --- a/src/gpu/GridGenerator/grid/GridBuilder/MultipleGridBuilder.cpp +++ b/src/gpu/GridGenerator/grid/GridBuilder/MultipleGridBuilder.cpp @@ -248,7 +248,7 @@ SPtr<Grid> MultipleGridBuilder::makeRotatingGrid(SPtr<Cylinder> cylinder, uint l // #TODO: change next line to use staggering cylinder = std::make_shared<Cylinder>(cylinder->getX1Centroid(), cylinder->getX2Centroid(), cylinder->getX3Centroid(), - cylinder->getRadius() + 6.0 * delta, cylinder->getHeight() + 6.0 * delta, cylinder->getPrincipalAxis()); + cylinder->getRadius() + 6.0 * delta, cylinder->getHeight() + 6.0 * delta, cylinder->getRotationalAxis()); SPtr<Grid> newGrid = this->makeGrid(cylinder, staggeredCoordinates[0], staggeredCoordinates[1], diff --git a/src/gpu/GridGenerator/grid/GridBuilder/MultipleGridBuilderTest.cpp b/src/gpu/GridGenerator/grid/GridBuilder/MultipleGridBuilderTest.cpp index 866911b19cb45221bbb3b17ef0536de6144e1989..61e6285be4682ae51df5da51c74ccaf46ab4abd4 100644 --- a/src/gpu/GridGenerator/grid/GridBuilder/MultipleGridBuilderTest.cpp +++ b/src/gpu/GridGenerator/grid/GridBuilder/MultipleGridBuilderTest.cpp @@ -18,7 +18,7 @@ protected: public: void SetUp() override { - cylinder = std::make_shared<Cylinder>(0.0, 0.0, 0.0, 2.0, 8.0, Cylinder::PrincipalAxis::x); + cylinder = std::make_shared<Cylinder>(0.0, 0.0, 0.0, 2.0, 8.0, Cylinder::RotationalAxis::x); gridBuilder.addCoarseGrid(-10. + 0.5 * delta, -10. + 0.5 * delta, -10. + 0.5 * delta, 10. - 0.5 * delta, 10. - 0.5 * delta, 10. - 0.5 * delta, delta); }