const Filter

Author: ProExpertProg

cpp/lib/Filter.cpp

@@ -4,7 +4,7 @@
 
 #include "Filter.hpp"
 namespace ahr {
-void HouLiFilter::operator()(Grid::View::C_XY view) {
+void HouLiFilter::operator()(Grid::View::C_XY view) const {
   grid.for_each_kxky([&](Dim kx, Dim ky) {
     view(kx, ky) *=
         exp(-36.0 * pow(kx_(kx) / grid.KX, 36.0)) * exp(-36.0 * pow(ky_(ky) / grid.KY, 36.0));
@@ -19,7 +19,7 @@ HouLiFilterCached::HouLiFilterCached(Grid const &grid)
   });
 }
 
-void HouLiFilterCached::operator()(Grid::View::C_XY view) {
+void HouLiFilterCached::operator()(Grid::View::C_XY view) const {
   grid.for_each_kxky([&](Dim kx, Dim ky) { view(kx, ky) *= factors(kx, ky); });
 }
 
@@ -33,7 +33,7 @@ HouLiFilterCached1D::HouLiFilterCached1D(Grid const &grid)
   }
 }
 
-void HouLiFilterCached1D::operator()(Grid::View::C_XY view) {
+void HouLiFilterCached1D::operator()(Grid::View::C_XY view) const {
   grid.for_each_kxky([&](Dim kx, Dim ky) {
     // Extra multiplication at runtime for lower memory cost
     view(kx, ky) *= factors_x[kx] * factors_y[ky];

cpp/lib/Filter.hpp

@@ -8,7 +8,7 @@ class HouLiFilter {
 public:
   explicit HouLiFilter(Grid const &grid) : grid(grid) {}
 
-  void operator()(Grid::View::C_XY view);
+  void operator()(Grid::View::C_XY view) const;
 
 protected:
   Grid const &grid;
@@ -23,7 +23,7 @@ class HouLiFilter {
 class HouLiFilterCached : HouLiFilter {
 public:
   explicit HouLiFilterCached(Grid const &grid);
-  void operator()(Grid::View::C_XY view);
+  void operator()(Grid::View::C_XY view) const;
 
 private:
   /// Pre-calculated factors for the Hou-Li filter.
@@ -34,7 +34,7 @@ class HouLiFilterCached : HouLiFilter {
 class HouLiFilterCached1D : protected HouLiFilter {
 public:
   explicit HouLiFilterCached1D(Grid const &grid);
-  void operator()(Grid::View::C_XY view);
+  void operator()(Grid::View::C_XY view) const;
 
 protected:
   std::vector<Real> factors_x, factors_y;