#include "ov.h"
#include "defun.h"
#include "error.h"
#include "dNDArray.h"
#include "fNDArray.h"
#include "lo-specfun.h"

Include dependency graph for psi.cc:

Macros
#define	FLOAT_BRANCH(T, A, M, E)

#define	FLOAT_BRANCH(T, A, M, E)

Functions
OCTAVE_EXPORT octave_value_list	Fpsi (const octave_value_list &args, int) ar

defaults to zero A value of zero computes the digamma a value the trigamma and so on The digamma function is for polygamma	functions (ar{k}higher than 0)

Variables
defaults to zero A value of zero computes the digamma a value the trigamma and so on The digamma function is	ar {z} can have any value real or complex value. However

defaults to zero A value of zero computes the digamma a value the trigamma and so on The digamma function is	defined

defaults to zero A value of zero computes the digamma	function

defaults to zero A value of zero computes the digamma a value	of

Macro Definition Documentation

#define FLOAT_BRANCH	(	T,
		A,
		M,
		E
	)

Value:

if (oct_z.is_ ## T ##_type ())                            \
        {                                                       \
          const A ## NDArray z = oct_z.M ## array_value ();     \
          A ## NDArray psi_z (z.dims ());                       \
                                                                \
          const E* zv = z.data ();                              \
          E* psi_zv = psi_z.fortran_vec ();                     \
          const octave_idx_type n = z.numel ();                 \
          for (octave_idx_type i = 0; i < n; i++)               \
            *psi_zv++ = octave::math::psi (*zv++);              \
                                                                \
          retval = psi_z;                                       \
        }

#define FLOAT_BRANCH	(	T,
		A,
		M,
		E
	)

Value:

if (oct_z.is_ ## T ##_type ())                                    \
        {                                                               \
          const A ## NDArray z = oct_z.M ## array_value ();             \
          A ## NDArray psi_z (z.dims ());                               \
                                                                        \
          const E* zv = z.data ();                                      \
          E* psi_zv = psi_z.fortran_vec ();                             \
          const octave_idx_type n = z.numel ();                         \
          for (octave_idx_type i = 0; i < n; i++)                       \
            {                                                           \
              if (*zv < 0)                                              \
                error ("psi: Z must be non-negative for polygamma (K > 0)"); \
                                                                        \
              *psi_zv++ = octave::math::psi (k, *zv++);                 \
            }                                                           \
          retval = psi_z;                                               \
        }