math_funcs.h

Go to the documentation of this file.

/* Copyright (C) 2002-2012
   Patrick Eriksson <Patrick.Eriksson@chalmers.se>
   Stefan Buehler   <sbuehler@ltu.se>
 
   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 2, 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, write to the Free Software
   Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307,
   USA. */
 
/*****************************************************************************
 ***  File description 
 *****************************************************************************/
 
#ifndef math_funcs_h
#define math_funcs_h
 
#include "matpack.h"
#include "matpackI.h"
#include "matpackIII.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);
 
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