Extension_RandomGaussian.cpp

/* -----------------------------------------------------------------------------
Copyright 2021 Kevin P. Barry

Licensed under the Apache License, Version 2.0 (the "License");
you may not use this file except in compliance with the License.
You may obtain a copy of the License at

    http://www.apache.org/licenses/LICENSE-2.0

Unless required by applicable law or agreed to in writing, software
distributed under the License is distributed on an "AS IS" BASIS,
WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
See the License for the specific language governing permissions and
limitations under the License.
----------------------------------------------------------------------------- */

// Author: Kevin P. Barry [ta0kira@gmail.com]

#include <random>

#include "category-source.hpp"
#include "Streamlined_RandomGaussian.hpp"
#include "Category_Float.hpp"
#include "Category_Generator.hpp"
#include "Category_RandomGaussian.hpp"
#include "Category_RandomUniform.hpp"

#ifdef ZEOLITE_PUBLIC_NAMESPACE
namespace ZEOLITE_PUBLIC_NAMESPACE {
#endif  // ZEOLITE_PUBLIC_NAMESPACE

BoxedValue CreateValue_RandomGaussian(S<const Type_RandomGaussian> parent, const ParamsArgs& params_args);

struct ExtCategory_RandomGaussian : public Category_RandomGaussian {
};

struct ExtType_RandomGaussian : public Type_RandomGaussian {
  inline ExtType_RandomGaussian(Category_RandomGaussian& p, Params<0>::Type params) : Type_RandomGaussian(p, params) {}

  ReturnTuple Call_new(const ParamsArgs& params_args) const final {
    TRACE_FUNCTION("RandomGaussian.new")
    const PrimFloat Var_arg1 = (params_args.GetArg(0)).AsFloat();
    const PrimFloat Var_arg2 = (params_args.GetArg(1)).AsFloat();
    if (Var_arg2 <= 0) {
      FAIL() << "Invalid standard deviation " << Var_arg2;
    }
    return ReturnTuple(CreateValue_RandomGaussian(PARAM_SELF, params_args));
  }
};

struct ExtValue_RandomGaussian : public Value_RandomGaussian {
  inline ExtValue_RandomGaussian(S<const Type_RandomGaussian> p, const ParamsArgs& params_args)
    : Value_RandomGaussian(std::move(p)),
      generator_((int) (long) this),
      distribution_(params_args.GetArg(0).AsFloat(), params_args.GetArg(1).AsFloat()) {}

  ReturnTuple Call_generate(const ParamsArgs& params_args) final {
    TRACE_FUNCTION("RandomGaussian.generate")
    return ReturnTuple(Box_Float(distribution_(generator_)));
  }

  ReturnTuple Call_setSeed(const ParamsArgs& params_args) final {
    TRACE_FUNCTION("RandomGaussian.setSeed")
    distribution_.reset();
    generator_.seed(params_args.GetArg(0).AsInt());
    return ReturnTuple(VAR_SELF);
  }

  std::default_random_engine generator_;
  std::normal_distribution<double> distribution_;
};

Category_RandomGaussian& CreateCategory_RandomGaussian() {
  static auto& category = *new ExtCategory_RandomGaussian();
  return category;
}

S<const Type_RandomGaussian> CreateType_RandomGaussian(const Params<0>::Type& params) {
  static const auto cached = S_get(new ExtType_RandomGaussian(CreateCategory_RandomGaussian(), Params<0>::Type()));
  return cached;
}

void RemoveType_RandomGaussian(const Params<0>::Type& params) {}

BoxedValue CreateValue_RandomGaussian(S<const Type_RandomGaussian> parent, const ParamsArgs& params_args) {
  return BoxedValue::New<ExtValue_RandomGaussian>(std::move(parent), params_args);
}

#ifdef ZEOLITE_PUBLIC_NAMESPACE
}  // namespace ZEOLITE_PUBLIC_NAMESPACE
using namespace ZEOLITE_PUBLIC_NAMESPACE;
#endif  // ZEOLITE_PUBLIC_NAMESPACE