Header file containing the numerical recipies and GNU Scientific Library routines used in the library.
Definition in file bpm_nr.h.
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <float.h>
#include <string.h>
#include <bpm/bpm_defs.h>
Include dependency graph for bpm_nr.h:
Go to the source code of this file.
Data Structures | |
| struct | lm_fstate |
| struct | gsl_block_struct |
| struct | gsl_matrix |
| struct | _gsl_matrix_view |
| struct | gsl_vector |
| struct | _gsl_vector_view |
| struct | _gsl_vector_const_view |
| struct | complex_t |
Defines | |
| #define | GCF_ITMAX |
| #define | GCF_FPMIN |
| #define | GCF_EPS |
| #define | GSER_EPS |
| #define | GSER_ITMAX |
| #define | RAN1_IA |
| #define | RAN1_IM |
| #define | RAN1_AM |
| #define | RAN1_IQ |
| #define | RAN1_IR |
| #define | RAN1_NTAB |
| #define | RAN1_NDIV |
| #define | RAN1_EPS |
| #define | RAN1_RNMX |
| #define | __LM_BLOCKSZ__ |
| #define | __LM_BLOCKSZ__SQ |
| #define | LINSOLVERS_RETAIN_MEMORY |
| #define | __LM_STATIC__ |
| #define | FABS(x) |
| #define | CNST(x) |
| #define | _LM_POW_ |
| #define | LM_DER_WORKSZ(npar, nmeas) |
| #define | LM_DIF_WORKSZ(npar, nmeas) |
| #define | LM_EPSILON |
| #define | LM_ONE_THIRD |
| #define | LM_OPTS_SZ |
| #define | LM_INFO_SZ |
| #define | LM_INIT_MU |
| #define | LM_STOP_THRESH |
| #define | LM_DIFF_DELTA |
| #define | NR_FFTFORWARD |
| #define | NR_FFTBACKWARD |
| #define | __LM_MEDIAN3(a, b, c) |
| #define | NULL_VECTOR |
| #define | NULL_VECTOR_VIEW |
| #define | NULL_MATRIX |
| #define | NULL_MATRIX_VIEW |
| #define | GSL_DBL_EPSILON |
| #define | OFFSET(N, incX) |
| #define | GSL_MIN(a, b) |
Typedefs | |
| typedef enum CBLAS_TRANSPOSE | CBLAS_TRANSPOSE_t |
| typedef gsl_block_struct | gsl_block |
| typedef _gsl_matrix_view | gsl_matrix_view |
| typedef _gsl_vector_view | gsl_vector_view |
| typedef const _gsl_vector_const_view | gsl_vector_const_view |
Enumerations | |
| enum | CBLAS_TRANSPOSE { CblasNoTrans, CblasTrans, CblasConjTrans } |
| enum | CBLAS_ORDER { CblasRowMajor, CblasColMajor } |
Functions | |
| EXTERN double | nr_gammln (double xx) |
| EXTERN double | nr_gammq (double a, double x) |
| EXTERN int | nr_gcf (double *gammcf, double a, double x, double *gln) |
| EXTERN int | nr_gser (double *gamser, double a, double x, double *gln) |
| EXTERN int | nr_fit (double *x, double y[], int ndata, double sig[], int mwt, double *a, double *b, double *siga, double *sigb, double *chi2, double *q) |
| EXTERN int | nr_is_pow2 (unsigned long n) |
| EXTERN int | nr_four1 (double data[], unsigned long nn, int isign) |
| EXTERN int | nr_realft (double data[], unsigned long n, int isign) |
| EXTERN double | nr_ran1 (long *idum) |
| EXTERN int | nr_seed (long seed) |
| EXTERN double | nr_ranuniform (double lower, double upper) |
| EXTERN double | nr_rangauss (double mean, double std_dev) |
| EXTERN int | nr_lmder (void(*func)(double *p, double *hx, int m, int n, void *adata), void(*jacf)(double *p, double *j, int m, int n, void *adata), double *p, double *x, int m, int n, int itmax, double *opts, double *info, double *work, double *covar, void *adata) |
| EXTERN int | nr_lmdif (void(*func)(double *p, double *hx, int m, int n, void *adata), double *p, double *x, int m, int n, int itmax, double *opts, double *info, double *work, double *covar, void *adata) |
| EXTERN int | nr_lmder_bc (void(*func)(double *p, double *hx, int m, int n, void *adata), void(*jacf)(double *p, double *j, int m, int n, void *adata), double *p, double *x, int m, int n, double *lb, double *ub, int itmax, double *opts, double *info, double *work, double *covar, void *adata) |
| EXTERN int | nr_lmdif_bc (void(*func)(double *p, double *hx, int m, int n, void *adata), double *p, double *x, int m, int n, double *lb, double *ub, int itmax, double *opts, double *info, double *work, double *covar, void *adata) |
| EXTERN void | nr_lmchkjac (void(*func)(double *p, double *hx, int m, int n, void *adata), void(*jacf)(double *p, double *j, int m, int n, void *adata), double *p, int m, int n, void *adata, double *err) |
| EXTERN int | nr_lmcovar (double *JtJ, double *C, double sumsq, int m, int n) |
| EXTERN int | nr_ax_eq_b_LU (double *A, double *B, double *x, int n) |
| EXTERN void | nr_trans_mat_mat_mult (double *a, double *b, int n, int m) |
| EXTERN void | nr_fdif_forw_jac_approx (void(*func)(double *p, double *hx, int m, int n, void *adata), double *p, double *hx, double *hxx, double delta, double *jac, int m, int n, void *adata) |
| EXTERN void | nr_fdif_cent_jac_approx (void(*func)(double *p, double *hx, int m, int n, void *adata), double *p, double *hxm, double *hxp, double delta, double *jac, int m, int n, void *adata) |
| EXTERN double | nr_median (int n, double *arr) |
| EXTERN double | nr_select (int k, int n, double *org_arr) |
| EXTERN gsl_matrix * | gsl_matrix_calloc (const size_t n1, const size_t n2) |
| EXTERN _gsl_vector_view | gsl_matrix_column (gsl_matrix *m, const size_t i) |
| EXTERN _gsl_matrix_view | gsl_matrix_submatrix (gsl_matrix *m, const size_t i, const size_t j, const size_t n1, const size_t n2) |
| EXTERN double | gsl_matrix_get (const gsl_matrix *m, const size_t i, const size_t j) |
| EXTERN void | gsl_matrix_set (gsl_matrix *m, const size_t i, const size_t j, const double x) |
| EXTERN int | gsl_matrix_swap_columns (gsl_matrix *m, const size_t i, const size_t j) |
| EXTERN gsl_matrix * | gsl_matrix_alloc (const size_t n1, const size_t n2) |
| EXTERN _gsl_vector_const_view | gsl_matrix_const_row (const gsl_matrix *m, const size_t i) |
| EXTERN _gsl_vector_view | gsl_matrix_row (gsl_matrix *m, const size_t i) |
| EXTERN _gsl_vector_const_view | gsl_matrix_const_column (const gsl_matrix *m, const size_t j) |
| EXTERN void | gsl_matrix_set_identity (gsl_matrix *m) |
| EXTERN gsl_vector * | gsl_vector_calloc (const size_t n) |
| EXTERN _gsl_vector_view | gsl_vector_subvector (gsl_vector *v, size_t offset, size_t n) |
| EXTERN double | gsl_vector_get (const gsl_vector *v, const size_t i) |
| EXTERN void | gsl_vector_set (gsl_vector *v, const size_t i, double x) |
| EXTERN int | gsl_vector_swap_elements (gsl_vector *v, const size_t i, const size_t j) |
| EXTERN _gsl_vector_const_view | gsl_vector_const_subvector (const gsl_vector *v, size_t i, size_t n) |
| EXTERN void | gsl_vector_free (gsl_vector *v) |
| EXTERN int | gsl_linalg_SV_solve (const gsl_matrix *U, const gsl_matrix *Q, const gsl_vector *S, const gsl_vector *b, gsl_vector *x) |
| EXTERN int | gsl_linalg_bidiag_unpack (const gsl_matrix *A, const gsl_vector *tau_U, gsl_matrix *U, const gsl_vector *tau_V, gsl_matrix *V, gsl_vector *diag, gsl_vector *superdiag) |
| EXTERN int | gsl_linalg_householder_hm (double tau, const gsl_vector *v, gsl_matrix *A) |
| EXTERN int | gsl_linalg_bidiag_unpack2 (gsl_matrix *A, gsl_vector *tau_U, gsl_vector *tau_V, gsl_matrix *V) |
| EXTERN int | gsl_linalg_householder_hm1 (double tau, gsl_matrix *A) |
| EXTERN void | create_givens (const double a, const double b, double *c, double *s) |
| EXTERN double | gsl_linalg_householder_transform (gsl_vector *v) |
| EXTERN int | gsl_linalg_householder_mh (double tau, const gsl_vector *v, gsl_matrix *A) |
| EXTERN void | chop_small_elements (gsl_vector *d, gsl_vector *f) |
| EXTERN void | qrstep (gsl_vector *d, gsl_vector *f, gsl_matrix *U, gsl_matrix *V) |
| EXTERN double | trailing_eigenvalue (const gsl_vector *d, const gsl_vector *f) |
| EXTERN void | create_schur (double d0, double f0, double d1, double *c, double *s) |
| EXTERN void | svd2 (gsl_vector *d, gsl_vector *f, gsl_matrix *U, gsl_matrix *V) |
| EXTERN void | chase_out_intermediate_zero (gsl_vector *d, gsl_vector *f, gsl_matrix *U, size_t k0) |
| EXTERN void | chase_out_trailing_zero (gsl_vector *d, gsl_vector *f, gsl_matrix *V) |
| EXTERN int | gsl_isnan (const double x) |
| EXTERN double | gsl_blas_dnrm2 (const gsl_vector *X) |
| EXTERN double | cblas_dnrm2 (const int N, const double *X, const int incX) |
| EXTERN void | gsl_blas_dscal (double alpha, gsl_vector *X) |
| EXTERN void | cblas_dscal (const int N, const double alpha, double *X, const int incX) |
| EXTERN void | cblas_dgemv (const enum CBLAS_ORDER order, const enum CBLAS_TRANSPOSE TransA, const int M, const int N, const double alpha, const double *A, const int lda, const double *X, const int incX, const double beta, double *Y, const int incY) |
| EXTERN gsl_block * | gsl_block_alloc (const size_t n) |
| EXTERN void | gsl_block_free (gsl_block *b) |
| EXTERN complex_t | complex (double re, double im) |
| EXTERN double | c_real (complex_t z) |
| EXTERN double | c_imag (complex_t z) |
| EXTERN complex_t | c_conj (complex_t z) |
| EXTERN complex_t | c_neg (complex_t z) |
| EXTERN complex_t | c_sum (complex_t z1, complex_t z2) |
| EXTERN complex_t | c_diff (complex_t z1, complex_t z2) |
| EXTERN complex_t | c_mult (complex_t z1, complex_t z2) |
| EXTERN complex_t | c_div (complex_t z1, complex_t z2) |
| EXTERN complex_t | c_scale (double r, complex_t z) |
| EXTERN complex_t | c_sqr (complex_t z) |
| EXTERN complex_t | c_sqrt (complex_t z) |
| EXTERN double | c_norm2 (complex_t z) |
| EXTERN double | c_abs (complex_t z) |
| EXTERN double | c_arg (complex_t z) |
| EXTERN complex_t | c_exp (complex_t z) |
| EXTERN int | c_isequal (complex_t z1, complex_t z2) |
Variables | |
| EXTERN long | bpm_rseed |
1.5.1