CubicSolver.hpp

Go to the documentation of this file.

/*
 * (c) Copyright 2021 CORSIKA Project, corsika-project@lists.kit.edu
 *
 * This software is distributed under the terms of the GNU General Public
 * Licence version 3 (GPL Version 3). See file LICENSE for a full version of
 * the license.
 */

#pragma once

#include <vector>
#include <cmath>
#include <algorithm>
#include <numeric>
#include <complex>

#include <corsika/framework/utility/QuadraticSolver.hpp>
#include <corsika/framework/core/Logging.hpp>

#include <boost/multiprecision/cpp_bin_float.hpp>

namespace corsika {

  namespace andre {
    std::vector<double> solveP3(long double a, long double b, long double c,
                                long double d, double const epsilon = 1e-12);
  } // namespace andre

  inline std::vector<double> solve_cubic_depressed_disciminant_real(
      long double p, long double q, long double const disc, double const epsilon = 1e-12);

  std::vector<double> solve_cubic_depressed_real(long double p, long double q,
                                                 double const epsilon = 1e-12);

  std::vector<double> solve_cubic_real_analytic(long double a, long double b,
                                                long double c, long double d,
                                                double const epsilon = 1e-12);

  template <typename T>
  T cubic_function(T x, T a, T b, T c, T d);

  template <typename T>
  T cubic_function_dfdx(T x, T a, T b, T c);

  template <typename T>
  T cubic_function_d2fd2x(T x, T a, T b);

  std::vector<double> solve_cubic_real(long double a, long double b, long double c,
                                       long double d, double const epsilon = 1e-12);


} // namespace corsika

#include <corsika/detail/framework/utility/CubicSolver.inl>