Bump.h

#pragma once
 
#include "HashGrid.h"
 
#include <array>
 
namespace environment::field::aggregation
{
 
constexpr int kBumpExponent = 2;
 
constexpr int kPolarRadial3D = 3;
constexpr int kPolarRadial2D = 4;
constexpr int kPolarAngularInside = 16;
constexpr int kPolarAngularOutside = 12;
constexpr int kRectGaussPoints = 6;
 
struct Bump
{
   double cx     = 0.0; 
   double cy     = 0.0; 
   double cz     = 0.0; 
   double a      = 1.0; 
   double b      = 1.0; 
   double c      = 1.0; 
   double cosPhi = 1.0; 
   double sinPhi = 0.0; 
   double amp    = 0.0; 
};
 
constexpr double kMassFactor = 32.0 * 3.14159265358979323846 / 105.0;
 
[[nodiscard]] double AmplitudeForMass(double biomass, double a, double b, double c);
 
[[nodiscard]] double Mass(const Bump& bump);
 
[[nodiscard]] Box SupportBox(const Bump& bump);
 
[[nodiscard]] double Value(const Bump& bump, double x, double y, double z);
 
[[nodiscard]] double Projection(const Bump& bump, double x, double y);
 
[[nodiscard]] double SliceMass(const Bump& bump, double z);
 
[[nodiscard]] double DiskIntegral(const Bump& bump, double x, double y, double z, double radius);
 
[[nodiscard]] double CircleIntegral(const Bump& bump, double x, double y, double radius);
 
[[nodiscard]] double RectIntegral(const Bump& bump, double xMin, double xMax, double yMin, double yMax);
 
} // namespace environment::field::aggregation