packages/engine/scram-node/src/expression/numerical.cc

Implementation of various numerical expressions.

Namespaces

Name
scram
scram::mef

Source code

/*
 * Copyright (C) 2014-2018 Olzhas Rakhimov
 * Copyright (C) 2023 OpenPRA ORG Inc.
 *
 * This program is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation; either version 3 of the License, or
 * (at your option) any later version.
 *
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program.  If not, see <http://www.gnu.org/licenses/>.
 */


#include "numerical.h"

#include "src/error.h"

namespace scram::mef {

template <>
void Div::Validate() const {
  auto it = Expression::args().begin();
  for (++it; it != Expression::args().end(); ++it) {
    const auto& expr = *it;
    Interval arg_interval = expr->interval();
    if (expr->value() == 0 || Contains(arg_interval, 0))
      SCRAM_THROW(DomainError("Division by 0."));
  }
}

template <>
void Mod::Validate() const {
  assert(args().size() == 2);
  auto* arg_two = args().back();
  int arg_value = arg_two->value();
  if (arg_value == 0)
    SCRAM_THROW(DomainError("Modulo second operand must not be 0."));
  Interval interval = arg_two->interval();
  int high = interval.upper();
  int low = interval.lower();
  if (high == 0 || low == 0 || (low < 0 && 0 < high)) {
    SCRAM_THROW(
        DomainError("Modulo second operand sample must not contain 0."));
  }
}

template <>
void Pow::Validate() const {
  assert(args().size() == 2);
  auto* arg_one = args().front();
  auto* arg_two = args().back();
  if (arg_one->value() == 0 && arg_two->value() <= 0)
    SCRAM_THROW(DomainError("0 to power 0 or less is undefined."));
  if (Contains(arg_one->interval(), 0) && !IsPositive(arg_two->interval())) {
    SCRAM_THROW(
        DomainError("Power expression 'base' sample range contains 0);"
                    " positive exponent is required."));
  }
}

Mean::Mean(std::vector<Expression*> args) : ExpressionFormula(std::move(args)) {
  if (Expression::args().size() < 2)
    SCRAM_THROW(ValidityError("Expression requires 2 or more arguments."));
}

}  // namespace scram::mef

---

Updated on 2025-11-11 at 16:51:08 +0000