EIC code displayed by LXR master/​acts/​Tests/​UnitTests/​Core/​Utilities/​AnyGridViewTests.cpp	Source navigation Diff markup Identifier search General search
	Version: master test Revert   Change

File indexing completed on 2026-10-11 08:09:31

 // This file is part of the ACTS project.
 //
 // Copyright (C) 2016 CERN for the benefit of the ACTS project
 //
 // This Source Code Form is subject to the terms of the Mozilla Public
 // License, v. 2.0. If a copy of the MPL was not distributed with this
 // file, You can obtain one at https://mozilla.org/MPL/2.0/.
 
 #include <boost/test/unit_test.hpp>
 
 #include "Acts/Utilities/AnyGridView.hpp"
 #include "Acts/Utilities/Axis.hpp"
 #include "Acts/Utilities/AxisDefinitions.hpp"
 #include "Acts/Utilities/Grid.hpp"
 #include "Acts/Utilities/TypeTraits.hpp"
 
 #include <stdexcept>
 #include <string>
 #include <utility>
 #include <vector>
 
 using namespace Acts;
 
 namespace ActsTests {
 
 BOOST_AUTO_TEST_SUITE(UtilitiesSuite)
 
 // Helper function to create a 1D grid with int values
 auto createIntGrid1D() {
   Axis a(0.0, 4.0, 4u);
   Grid g(Type<int>, a);
 
   // Fill with values
   for (std::size_t i = 0; i < g.size(); i++) {
     g.at(i) = static_cast<int>(i * 10);
   }
 
   return g;
 }
 
 // Helper function to create a 2D grid with double values
 auto createDoubleGrid2D() {
   Axis a(0.0, 4.0, 4u);
   Axis b(0.0, 6.0, 3u);
   Grid g(Type<double>, std::move(a), std::move(b));
 
   // Fill with values
   for (std::size_t i = 0; i < g.size(); i++) {
     g.at(i) = static_cast<double>(i) * 1.5;
   }
 
   return g;
 }
 
 // Helper function to create a 2D grid with string values
 auto createStringGrid2D() {
   Axis a(0.0, 4.0, 4u);
   Axis b(0.0, 6.0, 3u);
   Grid g(Type<std::string>, std::move(a), std::move(b));
 
   // Fill with values
   for (std::size_t i = 0; i < g.size(); i++) {
     g.at(i) = "Value_" + std::to_string(i);
   }
 
   return g;
 }
 
 // Helper function to create a 3D grid with vector values
 auto createVectorGrid3D() {
   Axis a(0.0, 4.0, 2u);
   Axis b(0.0, 6.0, 2u);
   Axis c(0.0, 8.0, 2u);
   Grid g(Type<std::vector<float>>, std::move(a), std::move(b), std::move(c));
 
   // Fill with values
   for (std::size_t i = 0; i < g.size(); i++) {
     g.at(i) = {static_cast<float>(i), static_cast<float>(i * 2)};
   }
 
   return g;
 }
 
 // Test read and write access to int grid using type safe interface
 BOOST_AUTO_TEST_CASE(IntGridReadWriteAccess) {
   auto grid = createIntGrid1D();
 
   // Create mutable view
   AnyGridView<int> view(grid);
 
   // Test read access
   BOOST_CHECK_EQUAL(view.atLocalBins({1}), 10);
   BOOST_CHECK_EQUAL(view.atLocalBins({2}), 20);
 
   // Test write access
   view.atLocalBins({1}) = 100;
   BOOST_CHECK_EQUAL(view.atLocalBins({1}), 100);
   BOOST_CHECK_EQUAL(grid.atLocalBins({1}), 100);
 }
 
 // Test read-only access to int grid using type safe interface
 BOOST_AUTO_TEST_CASE(IntGridReadOnlyAccess) {
   auto grid = createIntGrid1D();
 
   // Create const view
   AnyGridConstView<int> constView(grid);
 
   // Test read access
   BOOST_CHECK_EQUAL(constView.atLocalBins({1}), 10);
   BOOST_CHECK_EQUAL(constView.atLocalBins({2}), 20);
 
   // Verify that the following would not compile:
   // constView.atLocalBins({1}) = 100; // Should cause compilation error
 }
 
 // Test creation of any grid view from concrete grid type
 BOOST_AUTO_TEST_CASE(CreateFromConcreteGrid) {
   auto doubleGrid = createDoubleGrid2D();
 
   // Create view from concrete grid
   AnyGridView<double> view(doubleGrid);
 
   // Check dimensions
   BOOST_CHECK_EQUAL(view.dimensions(), 2u);
 
   // Check values
   BOOST_CHECK_CLOSE(
       view.atLocalBins({1, 1}),
       1.5 * doubleGrid.multiAxis().getGlobalBinFromLocalBins({1, 1}), 1e-10);
 
   // Modify through view
   view.atLocalBins({1, 1}) = 42.0;
   BOOST_CHECK_CLOSE(doubleGrid.atLocalBins({1, 1}), 42.0, 1e-10);
 }
 
 // Test creation of any grid view from IGrid type
 BOOST_AUTO_TEST_CASE(CreateFromIGrid) {
   auto doubleGrid = createDoubleGrid2D();
   IGrid& iGrid = doubleGrid;
 
   // Create view from IGrid
   AnyGridView<double> view(iGrid);
 
   // Check dimensions
   BOOST_CHECK_EQUAL(view.dimensions(), 2u);
 
   // Check values
   BOOST_CHECK_CLOSE(
       view.atLocalBins({1, 1}),
       1.5 * doubleGrid.multiAxis().getGlobalBinFromLocalBins({1, 1}), 1e-10);
 }
 
 // Test creation of const grid view from const IGrid type
 BOOST_AUTO_TEST_CASE(CreateConstViewFromConstIGrid) {
   auto doubleGrid = createDoubleGrid2D();
   const IGrid& constIGrid = doubleGrid;
 
   // Create const view from const IGrid
   AnyGridConstView<double> constView(constIGrid);
 
   // Check dimensions
   BOOST_CHECK_EQUAL(constView.dimensions(), 2u);
 
   // Check values
   BOOST_CHECK_CLOSE(
       constView.atLocalBins({1, 1}),
       1.5 * doubleGrid.multiAxis().getGlobalBinFromLocalBins({1, 1}), 1e-10);
 }
 
 // Test type mismatch handling when creating from IGrid
 BOOST_AUTO_TEST_CASE(TypeMismatchFromIGrid) {
   auto doubleGrid = createDoubleGrid2D();
   IGrid& iGrid = doubleGrid;
 
   // Try to create int view from double grid
   BOOST_CHECK_THROW((AnyGridView<int>(iGrid)), std::invalid_argument);
 }
 
 // Test type mismatch handling when creating const view from const IGrid
 BOOST_AUTO_TEST_CASE(TypeMismatchFromConstIGrid) {
   auto doubleGrid = createDoubleGrid2D();
   const IGrid& constIGrid = doubleGrid;
 
   // Try to create int view from double grid
   BOOST_CHECK_THROW((AnyGridConstView<int>(constIGrid)), std::invalid_argument);
 }
 
 // Test creation of mutable view from mutable grid
 BOOST_AUTO_TEST_CASE(MutableViewFromMutableGrid) {
   auto stringGrid = createStringGrid2D();
 
   // Create mutable view
   AnyGridView<std::string> view(stringGrid);
 
   // Modify through view
   view.atLocalBins({2, 2}) = "Modified";
 
   // Check modification in original grid
   BOOST_CHECK_EQUAL(stringGrid.atLocalBins({2, 2}), "Modified");
 }
 
 // Test creation of const view from mutable grid
 BOOST_AUTO_TEST_CASE(ConstViewFromMutableGrid) {
   auto stringGrid = createStringGrid2D();
 
   // Create const view from mutable grid
   AnyGridConstView<std::string> constView(stringGrid);
 
   // Check read access
   BOOST_CHECK_EQUAL(constView.atLocalBins({2, 2}),
                     stringGrid.atLocalBins({2, 2}));
 
   // Verify that the following would not compile:
   // constView.atLocalBins({2, 2}) = "Modified"; // Should cause compilation
   // error
 }
 
 // Test creation of const view from const grid
 BOOST_AUTO_TEST_CASE(ConstViewFromConstGrid) {
   const auto stringGrid = createStringGrid2D();
 
   // Create const view from const grid
   AnyGridConstView<std::string> constView(stringGrid);
 
   // Check read access
   BOOST_CHECK_EQUAL(constView.atLocalBins({2, 2}),
                     stringGrid.atLocalBins({2, 2}));
 }
 
 // Test complex type (vector) with both IGrid and concrete grid type
 // construction
 BOOST_AUTO_TEST_CASE(VectorTypeWithBothConstructions) {
   auto vectorGrid = createVectorGrid3D();
 
   // Test with concrete grid type
   AnyGridView<std::vector<float>> concreteView(vectorGrid);
   BOOST_CHECK_EQUAL(concreteView.dimensions(), 3u);
   BOOST_CHECK_EQUAL(concreteView.atLocalBins({1, 1, 1}).size(), 2u);
   BOOST_CHECK_EQUAL(concreteView.atLocalBins({1, 1, 1})[0],
                     vectorGrid.atLocalBins({1, 1, 1})[0]);
 
   // Test with IGrid type
   IGrid& iGrid = vectorGrid;
   AnyGridView<std::vector<float>> iGridView(iGrid);
   BOOST_CHECK_EQUAL(iGridView.dimensions(), 3u);
   BOOST_CHECK_EQUAL(iGridView.atLocalBins({1, 1, 1}).size(), 2u);
   BOOST_CHECK_EQUAL(iGridView.atLocalBins({1, 1, 1})[0],
                     vectorGrid.atLocalBins({1, 1, 1})[0]);
 
   // Test with const IGrid type
   const IGrid& constIGrid = vectorGrid;
   AnyGridConstView<std::vector<float>> constIGridView(constIGrid);
   BOOST_CHECK_EQUAL(constIGridView.dimensions(), 3u);
   BOOST_CHECK_EQUAL(constIGridView.atLocalBins({1, 1, 1}).size(), 2u);
   BOOST_CHECK_EQUAL(constIGridView.atLocalBins({1, 1, 1})[0],
                     vectorGrid.atLocalBins({1, 1, 1})[0]);
 
   // Modify through view
   std::vector<float> newValue = {99.0f, 88.0f};
   concreteView.atLocalBins({1, 1, 1}) = newValue;
   BOOST_CHECK_EQUAL(vectorGrid.atLocalBins({1, 1, 1})[0], 99.0f);
   BOOST_CHECK_EQUAL(vectorGrid.atLocalBins({1, 1, 1})[1], 88.0f);
 }
 
 // Test grid properties access through view
 BOOST_AUTO_TEST_CASE(GridPropertiesAccess) {
   auto doubleGrid = createDoubleGrid2D();
   AnyGridView<double> view(doubleGrid);
 
   // Test dimensions
   BOOST_CHECK_EQUAL(view.dimensions(), 2u);
 
   // Test bin center
   auto center = view.multiAxisAny().getBinCenterAny({1, 1});
   auto expectedCenter = doubleGrid.multiAxis().getBinCenter({1, 1});
   BOOST_CHECK_CLOSE(center[0], expectedCenter[0], 1e-10);
   BOOST_CHECK_CLOSE(center[1], expectedCenter[1], 1e-10);
 
   // Test lower left bin edge
   auto lowerLeft = view.multiAxisAny().getLowerLeftBinEdgeAny({1, 1});
   auto expectedLowerLeft = doubleGrid.multiAxis().getLowerLeftBinEdge({1, 1});
   BOOST_CHECK_CLOSE(lowerLeft[0], expectedLowerLeft[0], 1e-10);
   BOOST_CHECK_CLOSE(lowerLeft[1], expectedLowerLeft[1], 1e-10);
 
   // Test upper right bin edge
   auto upperRight = view.multiAxisAny().getUpperRightBinEdgeAny({1, 1});
   auto expectedUpperRight = doubleGrid.multiAxis().getUpperRightBinEdge({1, 1});
   BOOST_CHECK_CLOSE(upperRight[0], expectedUpperRight[0], 1e-10);
   BOOST_CHECK_CLOSE(upperRight[1], expectedUpperRight[1], 1e-10);
 
   // Test number of local bins
   auto numBins = view.multiAxisAny().getNBinsAny();
   auto expectedNumBins = doubleGrid.multiAxis().getNBins();
   BOOST_CHECK_EQUAL(numBins[0], expectedNumBins[0]);
   BOOST_CHECK_EQUAL(numBins[1], expectedNumBins[1]);
 }
 
 // Test grid properties access through const view from const IGrid
 BOOST_AUTO_TEST_CASE(GridPropertiesAccessConstView) {
   auto doubleGrid = createDoubleGrid2D();
   const IGrid& constIGrid = doubleGrid;
   AnyGridConstView<double> constView(constIGrid);
 
   // Test dimensions
   BOOST_CHECK_EQUAL(constView.dimensions(), 2u);
 
   // Test bin center
   auto center = constView.multiAxisAny().getBinCenterAny({1, 1});
   auto expectedCenter = doubleGrid.multiAxis().getBinCenter({1, 1});
   BOOST_CHECK_CLOSE(center[0], expectedCenter[0], 1e-10);
   BOOST_CHECK_CLOSE(center[1], expectedCenter[1], 1e-10);
 
   // Test lower left bin edge
   auto lowerLeft = constView.multiAxisAny().getLowerLeftBinEdgeAny({1, 1});
   auto expectedLowerLeft = doubleGrid.multiAxis().getLowerLeftBinEdge({1, 1});
   BOOST_CHECK_CLOSE(lowerLeft[0], expectedLowerLeft[0], 1e-10);
   BOOST_CHECK_CLOSE(lowerLeft[1], expectedLowerLeft[1], 1e-10);
 
   // Test upper right bin edge
   auto upperRight = constView.multiAxisAny().getUpperRightBinEdgeAny({1, 1});
   auto expectedUpperRight = doubleGrid.multiAxis().getUpperRightBinEdge({1, 1});
   BOOST_CHECK_CLOSE(upperRight[0], expectedUpperRight[0], 1e-10);
   BOOST_CHECK_CLOSE(upperRight[1], expectedUpperRight[1], 1e-10);
 
   // Test number of local bins
   auto numBins = constView.multiAxisAny().getNBinsAny();
   auto expectedNumBins = doubleGrid.multiAxis().getNBins();
   BOOST_CHECK_EQUAL(numBins[0], expectedNumBins[0]);
   BOOST_CHECK_EQUAL(numBins[1], expectedNumBins[1]);
 }
 
 // Test error cases
 BOOST_AUTO_TEST_CASE(ErrorCases) {
   auto intGrid = createIntGrid1D();
   auto doubleGrid = createDoubleGrid2D();
 
   // Test accessing with wrong number of indices
   AnyGridView<int> intView(intGrid);
   BOOST_CHECK_THROW(intView.atLocalBins({1, 2}), std::invalid_argument);
 
   // Test accessing with out-of-bounds indices
   BOOST_CHECK_THROW(intView.atLocalBins({10}), std::out_of_range);
 
   // Test creating view with wrong value type
   IGrid& iGrid = doubleGrid;
   BOOST_CHECK_THROW(AnyGridView<std::string>{iGrid}, std::invalid_argument);
 
   // Test creating const view with wrong value type
   const IGrid& constIGrid = doubleGrid;
   BOOST_CHECK_THROW(AnyGridConstView<std::string>{constIGrid},
                     std::invalid_argument);
 }
 
 // Test copy operations
 BOOST_AUTO_TEST_CASE(CopyOperations) {
   // Create grids that will live for the duration of the test
   auto doubleGrid1 = createDoubleGrid2D();
   auto doubleGrid2 = createDoubleGrid2D();
   auto doubleGrid3 = createDoubleGrid2D();
 
   // Test copy constructor
   AnyGridView<double> view1(doubleGrid1);
   AnyGridView<double> view2(view1);
   BOOST_CHECK_CLOSE(view1.atLocalBins({1, 1}), view2.atLocalBins({1, 1}),
                     1e-10);
 
   // Test copy assignment
   AnyGridView<double> view3(doubleGrid2);
   view3 = view1;
   BOOST_CHECK_CLOSE(view1.atLocalBins({1, 1}), view3.atLocalBins({1, 1}),
                     1e-10);
 
   // Same for const views
   AnyGridConstView<double> constView1(doubleGrid1);
   AnyGridConstView<double> constView2(constView1);
   BOOST_CHECK_CLOSE(constView1.atLocalBins({1, 1}),
                     constView2.atLocalBins({1, 1}), 1e-10);
 
   auto doubleGrid4 = createDoubleGrid2D();
   AnyGridConstView<double> constView3(doubleGrid4);
   constView3 = constView1;
   BOOST_CHECK_CLOSE(constView1.atLocalBins({1, 1}),
                     constView3.atLocalBins({1, 1}), 1e-10);
 
   // Test const view from const IGrid
   const IGrid& constIGrid = doubleGrid1;
   AnyGridConstView<double> constView4(constIGrid);
   AnyGridConstView<double> constView5(constView4);
   BOOST_CHECK_CLOSE(constView4.atLocalBins({1, 1}),
                     constView5.atLocalBins({1, 1}), 1e-10);
 }
 
 BOOST_AUTO_TEST_CASE(TypeDeduction) {
   auto grid = createIntGrid1D();
   AnyGridView view(grid);
   static_assert(std::is_same_v<decltype(view), AnyGridView<int>>);
 
   const auto constGrid = createIntGrid1D();
   AnyGridConstView<int> constView(constGrid);
   static_assert(std::is_same_v<decltype(constView), AnyGridConstView<int>>);
 }
 
 BOOST_AUTO_TEST_SUITE_END()
 
 }  // namespace ActsTests

This page was automatically generated by the 2.3.7 LXR engine. The LXR team