sphere_n.hpp

Go to the documentation of this file.

#pragma once
 
#include <array>
#include <memory>
#include <mutex>
#include <span>
#include <vector>
 
#include "lds.hpp"
 
#ifndef M_PI
#    define M_PI 3.14159265358979323846264338327950288
#endif
 
namespace ldsgen {
 
    constexpr double PI = M_PI;
    constexpr double HALF_PI = M_PI / 2.0;
 
    constexpr std::size_t TABLE_SIZE = 300;
 
    extern const std::array<double, TABLE_SIZE> X;
    extern const std::array<double, TABLE_SIZE> NEG_COSINE;
    extern const std::array<double, TABLE_SIZE> SINE;
    extern const std::array<double, TABLE_SIZE> F2;
 
    std::vector<double> linspace(double start, double stop, unsigned int num);
 
    double simple_interp(double x_value, std::span<const double> x_points,
                         std::span<const double> y_points);
 
    std::vector<double> get_tp(unsigned int n);
 
    class SphereGen {
      public:
        SphereGen() = default;
        SphereGen(const SphereGen&) = default;
        SphereGen(SphereGen&&) noexcept = default;
        SphereGen& operator=(const SphereGen&) = default;
        SphereGen& operator=(SphereGen&&) noexcept = default;
        virtual ~SphereGen() = default;
        virtual std::vector<double> pop() = 0;
        virtual void reseed(unsigned long seed) = 0;
    };
 
    class Sphere3 {
      public:
        explicit Sphere3(std::span<const unsigned long> base);
        std::array<double, 4> pop();
        void reseed(unsigned long seed);
 
      private:
        VdCorput vdc_;
        Sphere sphere2_;
        mutable std::mutex mutex_;
    };
 
    class SphereWrapper : public SphereGen {
      public:
        explicit SphereWrapper(std::span<const unsigned long> base);
        std::vector<double> pop() override;
        void reseed(unsigned long seed) override;
 
      private:
        Sphere sphere_;
        mutable std::mutex mutex_;
    };
 
    class SphereN : public SphereGen {
      public:
        explicit SphereN(std::span<const unsigned long> base);
        std::vector<double> pop() override;
        void reseed(unsigned long seed) override;
 
      private:
        VdCorput vdc_;
        std::unique_ptr<SphereGen> s_gen_;
        unsigned int n_;
        std::vector<double> tp_;
        double range_;
        mutable std::mutex mutex_;
    };
}  // namespace ldsgen