log_range_param.h

Go to the documentation of this file.

// -----------------------------------------------------------------------------
//
// Copyright (C) 2021 CERN & University of Surrey for the benefit of the
// BioDynaMo collaboration. All Rights Reserved.
//
// Licensed under the Apache License, Version 2.0 (the "License");
// you may not use this file except in compliance with the License.
//
// See the LICENSE file distributed with this work for details.
// See the NOTICE file distributed with this work for additional information
// regarding copyright ownership.
//
// -----------------------------------------------------------------------------
 
#ifndef CORE_MULTI_SIMULATION_OPTIMIZATION_PARAM_TYPE_LOG_RANGE_PARAM_H_
#define CORE_MULTI_SIMULATION_OPTIMIZATION_PARAM_TYPE_LOG_RANGE_PARAM_H_
 
#include <cmath>
#include <string>
 
#include "core/multi_simulation/optimization_param_type/optimization_param_type.h"
#include "core/util/log.h"
 
namespace bdm {
 
struct LogRangeParam : public OptimizationParamType {
  LogRangeParam() = default;
  LogRangeParam(const std::string& n, real_t base, real_t min, real_t max,
                real_t stride)
      : OptimizationParamType(n),
        base(base),
        lower_bound(min),
        upper_bound(max),
        stride(stride) {
    Validate();
  };
 
  void Validate() const override {
    if (lower_bound > upper_bound) {
      Log::Fatal("LogRangeParam", "Tried to initialize parameter '", param_name,
                 "' with a lower_bound value higher than upper_bound: ",
                 lower_bound, " > ", upper_bound);
    }
  }
 
  OptimizationParamType* GetCopy() const override {
    return new LogRangeParam(*this);
  }
 
  // Get the nth value
  real_t GetValue(int n) const override {
    real_t exp = lower_bound + n * stride;
    return exp > upper_bound ? std::pow(base, upper_bound)
                             : std::pow(base, exp);
  }
 
  // Returns the number of discrete values that this range contains (including
  // the `lower_bound` and `upper_bound` values)
  uint32_t GetNumElements() const override {
    return std::round(((upper_bound - lower_bound) + stride) / stride);
  }
 
  // The base value
  real_t base = 10;
  // The minimum value
  real_t lower_bound = 0;
  // THe maximum value
  real_t upper_bound = 0;
  // The stride
  real_t stride = 1;
  BDM_CLASS_DEF_OVERRIDE(LogRangeParam, 1);
};
 
}  // namespace bdm
 
#endif  // CORE_MULTI_SIMULATION_OPTIMIZATION_PARAM_TYPE_LOG_RANGE_PARAM_H_