LM_fit.h File Reference

Levenbaerg-Marquardt non linear data fit routine. More...

#include "tbci/basics.h"
#include "tbci/vector.h"
#include "tbci/matrix.h"
#include "tbci/constants.h"
#include "tbci/solver/gauss_jordan.h"

Include dependency graph for LM_fit.h:

Go to the source code of this file.

Macros
#define	DELTAX   1e-30
#define	DELTA(i, j)
	the kronecker delta
#define	CON2   2.0
#define	CON   SQRT2
#define	BIG   1.0e30
#define	NTAB   10
#define	SAFE   2.0

Functions
	INST (template< typename T > class Vector friend int Basis_Trafo(long int, long int, Vector< T > &);) template< typename T > int Basis_Trafo(long int value
	if (value==0) return 1
	if (n<(int) CSTD__ ceil(MATH__ sqrt(2 *MATH__ log((double) value)/MATH__ log((double) basis)))||basis<=1) return 0
	while (value > 0)
	INST (template< typename T > class Vector friend T Partial_Del(T(*func)(const Vector< T > &, const Vector< T > &), const Vector< T > &, const Vector< T > &, int, T, T &);) template< typename T > T Partial_Del(T(*func)(const Vector< T > &
Matrix< T >	a (10, 10)
Vector< T >	buf_plus (xlen)
Vector< T >	buf_minus (xlen)

Variables
long int	basis
long int Vector< T > &	index
int	i = 0
	return
const Vector< T > const Vector< T > &	x
const Vector< T > const Vector< T > const Vector< T > &	p
const Vector< T > const Vector< T > const Vector< T > int	mu
const Vector< T > const Vector< T > const Vector< T > int T	h
const Vector< T > const Vector< T > const Vector< T > int T T &	err
int	xlen = x.size()

Detailed Description

Levenbaerg-Marquardt non linear data fit routine.

Definition in file LM_fit.h.

Macro Definition Documentation

BIG

#define BIG   1.0e30

Definition at line 88 of file LM_fit.h.

Referenced by buf_minus().

CON

#define CON   SQRT2

Definition at line 87 of file LM_fit.h.

Referenced by buf_minus().

CON2

#define CON2   2.0

Definition at line 86 of file LM_fit.h.

Referenced by buf_minus().

DELTA

#define DELTA	(		i,
			j )

Value:

((i==j)?1:0)

the kronecker delta

Definition at line 53 of file LM_fit.h.

DELTAX

#define DELTAX   1e-30

Definition at line 49 of file LM_fit.h.

NTAB

#define NTAB   10

Definition at line 89 of file LM_fit.h.

Referenced by a(), and buf_minus().

SAFE

#define SAFE   2.0

Definition at line 90 of file LM_fit.h.

Function Documentation

a()

Matrix< T > a	(	10	,
		10	)

References NTAB.

buf_minus()

Vector< T > buf_minus

(

xlen

)

References BIG, buf_minus(), buf_plus(), CON, CON2, err, fac, h, i, MATH__, mu, NTAB, p, REGISTER, STD__, x, and xlen.

Referenced by buf_minus().

buf_plus()

Vector< T > buf_plus

(

xlen

)

References xlen.

Referenced by buf_minus().

if() [1/2]

if

(

)

References basis, CSTD__, and MATH__.

if() [2/2]

if

(

value

= =0

)

INST() [1/2]

INST

(

template< typename T > class Vector friend int Basis_Trafo(long int, long int, Vector< T > &);

)

INST() [2/2]

INST

(

template< typename T > class Vector friend T Partial_Del(T(*func)(const Vector< T > &, const Vector< T > &), const Vector< T > &, const Vector< T > &, int, T, T &);

)

const &

References T.

while()

while	(	value	,
		0	)

Definition at line 76 of file LM_fit.h.

References basis, i, and index.

Variable Documentation

basis

long int basis

Definition at line 66 of file LM_fit.h.

Referenced by basis_trafo(), if(), and while().

err

const Vector<T> const Vector<T> const Vector<T> int T T& err

Initial value:

{
  
  T errt ALIGN(MIN_ALIGN), fac, hh, ans = (T)0

Definition at line 99 of file LM_fit.h.

Referenced by buf_minus(), free_threads(), gaussj(), init_threads(), lina_thread(), lu_decomp(), QMR(), and thread_wait_result().

h

const Vector<T> const Vector<T> const Vector<T> int T h

Definition at line 99 of file LM_fit.h.

Referenced by buf_minus(), Symm_BdMatrix< T >::setval(), ILU0_BdMatrixPreconditioner< T >::solve(), ILU0_Symm_BdMatrixPreconditioner< T >::solve(), sv_decomp(), Update(), and zhseqr_().

i

int i = 0

Definition at line 71 of file LM_fit.h.

Referenced by TBCI::_tbci_fill_fn(), BdMatrix< T >::adddiag(), CRMatrix< T >::allocate(), CSCMatrix< T >::allocate(), Symm_BdMatrix< T >::autoinsert(), basis_trafo(), BdMatrix< T >::BdMatrix(), BdMatrix< T >::BdMatrix(), TVector< T >::BdMatrix< T >, BdMatrixErr::BdMatrixErr(), BiCG(), BiCGSTAB(), bind_threads(), binom(), binomial(), BVector< T >::bubble_sort(), buf_minus(), BVector< T >::BVector(), BVector< T >::BVector(), CTensor< T >::calc_indx(), CTensor< T >::calc_offs(), CTensor< T >::calclayout(), CTensor< T >::calcsize(), cediv(), cediv(), cediv(), cediv(), cediv(), cediv(), cemul(), cemul(), cemul(), cemul(), cemul(), cemul(), CG(), CG2(), CGS(), CHEBY(), BdMatrix< T >::check(), BdMatrix< T >::check_internal(), Tensor< T >::cntrmul(), conj(), conj(), Symm_BdMatrix< T >::construct(), BdMatrix< T >::constructor(), BVector< T >::contains(), FS_Vector< dims, T >::contains(), CRMatrix< T >::copy(), CSCMatrix< T >::copy(), cplx< T >::cplx(), Symm_BdMatrix< T >::create(), CRMatErr::CRMatErr(), CRMatrix< T >::CRMatrix(), CSCMatrix< T >::cscm_vec_mul_exact(), CSCMatErr::CSCMatErr(), CTensor< T >::CTensor(), CTensor< T >::CTensor(), ctrmul(), CRMatrix< T >::destroy(), BdMatrix< T >::div_row(), BdMatrix< T >::div_rows(), do_bdmat_vec_dotmult(), do_bdmat_vec_mult_lnw(), do_bdmat_vec_mult_lnw_opt(), do_bdmat_vec_transmult_lnw(), do_bdmat_vec_transmult_lnw_opt(), BdMatrix< T >::do_copy(), CSCMatrix< T >::do_export(), do_fbdmat_vec_mul(), dot(), ediv(), ediv(), ediv(), ediv(), ediv(), ediv(), eig(), eig(), eig(), emul(), emul(), emul(), emul(), emul(), emul(), exact_sum(), BdMatrix< T >::expand(), expm2(), expm3(), F_BandMatErr::F_BandMatErr(), fac(), CRMatrix< T >::fill(), CSCMatrix< T >::fill(), FS_Vector< dims, T >::fill(), tbci_memalloc_cache< T >::find_by_sz(), fix_condition(), FS_Vector< dims, T >::FS_Vector(), FS_Vector< dims, T >::FS_Vector(), FSVecErr::FSVecErr(), gaussj(), Symm_BdMatrix< T >::get(), TSVector< T >::get(), TVector< T >::get(), BdMatrix< T >::get_col(), F_BandMatrix< T >::get_col(), CTensor< T >::get_lin_idx(), CTensor< T >::get_lin_idx(), BdMatrix< T >::get_row(), CTensor< T >::getcref(), TSVector< T >::getcref(), TVector< T >::getcref(), F_TMatrix< T >::herm(), Index::idx_fill_in1(), Index::idx_fill_in2(), Index::idx_remove1(), Index::idx_remove2(), IdxErr::IdxErr(), if(), imag(), imag(), FS_Vector< dims, T >::incr(), TSVector< T >::incr(), TVector< T >::incr(), Vector< T >::incr(), Index::Index(), CSCMatrix< T >::insert(), inv(), inv(), inv(), IR(), kahan_sum(), CTensor< T >::lin_read(), Index::lin_read(), LU_bkw_subst(), LU_bkw_subst(), lu_decomp(), lu_decomp(), LU_det(), LU_det(), LU_fwd_subst(), LU_fwd_subst(), LU_solve(), LU_solve(), lu_solve(), MatErr::MatErr(), CRMatrix< T >::MatVecMult(), CRMatrix< T >::MatVecMult(), CSCMatrix< T >::MatVecMult(), CSCMatrix< T >::MatVecMult(), metrmul(), CRMatrix< T >::mult(), CSCMatrix< T >::mult(), CSCMatrix< T >::mult(), Tensor< T >::mult(), BdMatrix< T >::mult_row(), BdMatrix< T >::mult_rows(), CSCMatrix< T >::multf(), norm_1(), norm_1(), BdMatrix< T >::operator TMatrix< T >(), operator!=(), BVector< T >::operator()(), BVector< T >::operator()(), CTensor< T >::operator()(), CTensor< T >::operator()(), CTensor< T >::operator()(), CTensor< T >::operator()(), FS_Vector< dims, T >::operator()(), FS_Vector< dims, T >::operator()(), Matrix< T >::operator()(), Matrix< T >::operator()(), Matrix< T >::operator()(), Symm_BdMatrix< T >::operator()(), TMatrix< T >::operator()(), TSVector< T >::operator()(), TVector< T >::operator()(), Vector< T >::operator()(), Vector< T >::operator()(), BdMatrix< T >::operator*(), CRMatrix< T >::operator*(), CRMatrix< T >::operator*(), CRMatrix< T >::operator*(), CSCMatrix< T >::operator*(), CSCMatrix< T >::operator*(), CSCMatrix< T >::operator*(), F_TSMatrix< T >::operator*(), Matrix< T >::operator*(), Symm_BdMatrix< T >::operator*(), Tensor< T >::operator*(), TSMatrix< T >::operator*(), TSMatrix< T >::operator*(), CRMatrix< T >::operator*=(), CSCMatrix< T >::operator*=(), TMatrix< T >::operator*FGD, operator+(), Tensor< T >::operator+(), Tensor< T >::operator+=(), BdMatrix< T >::operator-(), CRMatrix< T >::operator-(), CSCMatrix< T >::operator-(), FS_Vector< dims, T >::operator-(), operator-(), Tensor< T >::operator-(), Tensor< T >::operator-(), Tensor< T >::operator-=(), CRMatrix< T >::operator/(), CSCMatrix< T >::operator/(), FS_Vector< dims, T >::operator/(), Tensor< T >::operator/(), BdMatrix< T >::operator/=(), CRMatrix< T >::operator/=(), CSCMatrix< T >::operator/=(), FS_Vector< dims, T >::operator/=(), operator<<(), operator<<(), operator<<(), operator<<(), operator<<(), operator<<(), operator<<(), operator<<(), operator<<(), operator<<(), Symm_BdMatrix< T >::operator<<FGD, BVector< T >::operator<=(), CRMatrix< T >::operator=(), BdMatrix< T >::operator==(), CRMatrix< T >::operator==(), CSCMatrix< T >::operator==(), operator==(), TSVector< T >::operator==(), TSVector< T >::operator==(), Vector< T >::operator==(), BVector< T >::operator>=(), operator>>(), operator>>(), operator>>(), operator>>(), operator>>(), operator>>(), operator>>(), BVector< T >::operator[](), BVector< T >::operator[](), FS_Vector< dims, T >::operator[](), FS_Vector< dims, T >::operator[](), TSVector< T >::operator[](), TVector< T >::operator[](), Vector< T >::operator[](), Vector< T >::operator[](), own_ev_(), QMR(), real(), real(), BdMatrix< T >::reconfig(), BVector< T >::remove(), BdMatrix< T >::resize(), BdMatrix< T >::resize(), CSCMatrix< T >::resize(), TMatrix< T >::row_expand(), TMatrix< T >::row_expand(), cplx< T >::set(), TVector< T >::set(), F_BandMatrix< T >::set_col(), F_TMatrix< T >::set_ptrs(), TMatrix< T >::set_ptrs(), BdMatrix< T >::set_row(), CRMatrix< T >::setunit(), CSCMatrix< T >::setunit(), F_TMatrix< T >::setunit(), TMatrix< T >::setunit(), CRMatrix< T >::setval(), CRMatrix< T >::setval(), Symm_BdMatrix< T >::setval(), Symm_BdMatrix< T >::setval(), TVector< T >::setval(), TVector< T >::setval(), CRMatrix< T >::size(), DILU_BdMatrixPreconditioner< T >::solve(), ILU0_BdMatrixPreconditioner< T >::solve(), ILU0_Symm_BdMatrixPreconditioner< T >::solve(), sv_decomp(), sv_decomp_backsub(), sv_decomp_backsub(), SymmBdMatrixErr::SymmBdMatrixErr(), TensErr::TensErr(), Tensor< T >::Tensor(), Tensor< T >::Tensor(), BdMatrix< T >::trace(), F_TMatrix< T >::trace(), TMatrix< T >::trace(), F_TMatrix< T >::trans(), CRMatrix< T >::transMult(), CSCMatrix< T >::transMult(), F_BandMatrix< T >::transMult(), CTensor< T >::transpose(), Update(), DiagPreconditioner< T, MatrixType >::update(), DILU_BdMatrixPreconditioner< T >::update(), ILU0_BdMatrixPreconditioner< T >::update(), ILU0_Symm_BdMatrixPreconditioner< T >::update(), VecErr::VecErr(), Vector< T >::Vector(), and while().

index

index

Initial value:

{
   int n = index.size()

Definition at line 66 of file LM_fit.h.

Referenced by basis_trafo(), CTensor< T >::calc_indx(), CTensor< T >::calc_offs(), BdMatrix< T >::check(), BdMatrix< T >::div_row(), Matrix< T >::div_row(), F_BandMatrix< T >::get_col(), TMatrix< T >::get_col(), CRMatrix< T >::get_row(), F_TMatrix< T >::get_row(), Index::idx_remove1(), Index::idx_set1(), BdMatrix< T >::mult_row(), Matrix< T >::mult_row(), BVector< T >::operator()(), BVector< T >::operator()(), CRMatrix< T >::operator()(), F_BandMatrix< T >::operator()(), F_BandMatrix< T >::operator()(), F_Matrix< T >::operator()(), F_TMatrix< T >::operator()(), Matrix< T >::operator()(), Matrix< T >::operator()(), Matrix< T >::operator()(), TMatrix< T >::operator()(), F_BandMatrix< T >::set_col(), F_TMatrix< T >::set_col(), TMatrix< T >::set_col(), TMatrix< T >::set_col_partial(), CRMatrix< T >::set_row(), CRMatrix< T >::set_row(), F_TMatrix< T >::set_row(), TMatrix< T >::set_row(), TMatrix< T >::set_row_partial(), CRMatrix< T >::setval(), CSCMatrix< T >::setval(), F_BandMatrix< T >::setval(), and while().

mu

const Vector<T> const Vector<T> const Vector<T> int mu

Definition at line 99 of file LM_fit.h.

Referenced by buf_minus().

p

const Vector<T> const Vector<T> const Vector<T>& p

Definition at line 98 of file LM_fit.h.

Referenced by BiCG(), BiCGSTAB(), buf_minus(), cback::cback(), CG(), CG2(), CGS(), CHEBY(), expm(), polar(), QMR(), and BdMatrix< T >::setoutopts().

return

return

Definition at line 82 of file LM_fit.h.

Referenced by cgamma_(), chgf_(), and conhyp_().

x

const Vector<T> const Vector<T>& x

Definition at line 98 of file LM_fit.h.

Referenced by arydiv_(), asinh(), bernoulli(), BiCG(), BiCGSTAB(), binom(), binomial(), buf_minus(), CG(), CG2(), cgamma_(), CGS(), CHEBY(), chi_s(), dgbsvx2_(), dgbsvx_(), eig(), eig(), erfc3(), fac(), fact(), IR(), ldgamma(), LU_bkw_subst(), LU_bkw_subst(), LU_fwd_subst(), LU_fwd_subst(), LU_solve(), LU_solve(), lu_solve(), lu_solve(), lu_solve(), lu_solve(), lu_solve(), lu_solve(), lu_solve(), lu_solve_expert(), CSCMatrix< T >::mult(), CSCMatrix< T >::multf(), F_TSMatrix< T >::operator*(), poisson(), pow(), QMR(), TBCI::sign(), store_(), sv_decomp(), sv_decomp_backsub(), sv_decomp_backsub(), svd_solve(), trinomial(), Update(), WEAK(), WEAK(), and WEAK().

xlen

int xlen = x.size()

Definition at line 105 of file LM_fit.h.

Referenced by buf_minus(), and buf_plus().