function.h

Go to the documentation of this file.

/*
 * -------------------------------------------
 * Copyright (c) 2021 - 2026 Prashant K. Jha
 * -------------------------------------------
 * PeriDEM https://github.com/prashjha/PeriDEM
 *
 * Distributed under the Boost Software License, Version 1.0. (See accompanying
 * file LICENSE)
 */
 
#ifndef UTIL_FUNCTION_H
#define UTIL_FUNCTION_H
 
#include "point.h"              // definition of Point
#include <vector>
 
// tolerance for float comparison
#define COMPARE_EPS 1e-5
 
namespace util {
 
bool isGreater(const double &a, const double &b);
 
bool isLess(const double &a, const double &b);
 
double hatFunction(const double &x, const double &x_min, const double &x_max);
 
double hatFunctionQuick(const double &x, const double &x_min,
                        const double &x_max);
 
double linearStepFunc(const double &x, const double &x1, const double &x2);
 
double gaussian(const double &r, const double &a, const double &beta);
 
double gaussian2d(const util::Point &x, const size_t &dof,
                  const std::vector<double> &params);
 
double doubleGaussian2d(const util::Point &x, const size_t &dof,
                        const std::vector<double> &params);
 
double equivalentMass(const double &m1, const double &m2);
 
double harmonicMean(const double &m1, const double &m2);
 
double normalContactStiffness(const double &K1, const double &K2,
                              const double &horizon, int horizonPower = 5);
 
double selfContactStiffness(const double &K, const double &horizon,
                            int horizonPower = 5);
 
} // namespace util
 
#endif // UTIL_FUNCTION_H