math_funcs.h

Go to the documentation of this file.

/*****************************************************************************
 ***  File description 
 *****************************************************************************/
 
#ifndef math_funcs_h
#define math_funcs_h
 
#include "array.h"
#include "matpack_concepts.h"
#include "matpack_data.h"
#include "matpack_math.h"
 
Numeric fac(const Index n);
 
Index integer_div(const Index& x, const Index& y);
 
Numeric last(ConstVectorView x);
 
Index last(const ArrayOfIndex& x);
 
void linspace(Vector& x,
              const Numeric start,
              const Numeric stop,
              const Numeric step);
 
void nlinspace(Vector& x,
               const Numeric start,
               const Numeric stop,
               const Index n);
 
void nlinspace(VectorView x,
               const Numeric start,
               const Numeric stop,
               const Index n);
 
void nlogspace(Vector& x,
               const Numeric start,
               const Numeric stop,
               const Index n);
 
Numeric trapz(ConstVectorView x,
              ConstVectorView y);
 
void cumsum(VectorView csum,
            ConstVectorView x);
 
Numeric AngIntegrate_trapezoid(ConstMatrixView Integrand,
                               ConstVectorView za_grid,
                               ConstVectorView aa_grid);
 
Numeric AngIntegrate_trapezoid(ConstVectorView Integrand,
                               ConstVectorView za_grid);
 
Numeric AngIntegrate_trapezoid_opti(ConstMatrixView Integrand,
                                    ConstVectorView za_grid,
                                    ConstVectorView aa_grid,
                                    ConstVectorView grid_stepsize);
 
Numeric sign(const Numeric& x);
 
void mgd(VectorView psd,
         const Vector& x,
         const Numeric& n0,
         const Numeric& mu,
         const Numeric& la,
         const Numeric& ga);
 
void mgd_with_derivatives(VectorView psd,
                          MatrixView jac_data,
                          const Vector& x,
                          const Numeric& n0,
                          const Numeric& mu,
                          const Numeric& la,
                          const Numeric& ga,
                          const bool& do_n0_jac,
                          const bool& do_mu_jac,
                          const bool& do_la_jac,
                          const bool& do_ga_jac);
 
void delanoe_shape_with_derivative(VectorView psd,
                                   MatrixView jac_data,
                                   const Vector& x,
                                   const Numeric& alpha,
                                   const Numeric& beta);
 
Numeric mod_gamma_dist(
    Numeric x, Numeric N0, Numeric Lambda, Numeric mu, Numeric gamma);
 
void unitl(Vector& x);
 
void flat(VectorView x, ConstMatrixView X);
void flat(VectorView x, ConstTensor3View X);
 
void reshape(MatrixView X, ConstVectorView x);
void reshape(Tensor3View X, ConstVectorView x);
 
void calculate_weights_linear(Vector& x, Vector& w, const Index nph);
 
void calculate_int_weights_arbitrary_grid(Vector& w, const Vector& x);
 
 
template <typename MatpackType>
constexpr bool any_negative(const MatpackType& var) noexcept {
  if (var.empty()) return false;
  if (min(var) < 0) return true;
  return false;
}
 
constexpr Index pow_negative_one(Index x) noexcept { return (x % 2) ? -1 : 1; }
 
#endif  // math_funcs_h