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)   ((i==j)?1:0)
	the kronecker delta More...
#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.

CON

#define CON   SQRT2

Definition at line 87 of file LM_fit.h.

CON2

#define CON2   2.0

Definition at line 86 of file LM_fit.h.

DELTA

#define DELTA	(		i,
			j
	)		((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.

SAFE

#define SAFE   2.0

Definition at line 90 of file LM_fit.h.

Function Documentation

a()

Matrix<T> a	(	10	,
		10
	)

buf_minus()

Vector<T> buf_minus

(

xlen

)

buf_plus()

Vector<T> buf_plus

(

xlen

)

if() [1/2]

if

(

value

= =0

)

Referenced by job_vec_sum(), and Vector< unsigned >::Vector().

if() [2/2]

if

(

)

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 &

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(), 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 102 of file LM_fit.h.

Referenced by 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 97 of file LM_fit.h.

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

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(), BiCG(), BiCGSTAB(), bind_threads(), binom(), binomial(), BVector< unsigned >::bubble_sort(), BVector< unsigned >::BVector(), CTensor< T >::calc_indx(), CTensor< T >::calc_offs(), CTensor< T >::calclayout(), CTensor< T >::calcsize(), cediv(), cemul(), CG(), CG2(), CGS(), CHEBY(), BdMatrix< T >::check(), BdMatrix< T >::check_internal(), Tensor< T >::cntrmul(), conj(), Symm_BdMatrix< T >::construct(), BdMatrix< T >::constructor(), FS_Vector< dims, T >::contains(), BVector< unsigned >::contains(), CRMatrix< T >::copy(), CSCMatrix< T >::copy(), Symm_BdMatrix< T >::create(), CRMatrix< T >::CRMatrix(), CSCMatrix< T >::cscm_vec_mul_exact(), 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(), eig(), emul(), exact_sum(), BdMatrix< T >::expand(), expm2(), expm3(), fac(), CRMatrix< T >::fill(), FS_Vector< dims, T >::fill(), CSCMatrix< T >::fill(), tbci_memalloc_cache< T >::find_by_sz(), fix_condition(), FS_Vector< dims, T >::FS_Vector(), gaussj(), Symm_BdMatrix< T >::get(), TVector< T >::get(), TSVector< T >::get(), F_BandMatrix< T >::get_col(), BdMatrix< T >::get_col(), CTensor< T >::get_lin_idx(), BdMatrix< T >::get_row(), CTensor< T >::getcref(), TVector< T >::getcref(), TSVector< T >::getcref(), F_TMatrix< T >::herm(), Index::idx_fill_in1(), Index::idx_fill_in2(), Index::idx_remove1(), Index::idx_remove2(), if(), imag(), FS_Vector< dims, T >::incr(), TVector< T >::incr(), TSVector< T >::incr(), Vector< unsigned >::incr(), Index::Index(), CSCMatrix< T >::insert(), inv(), IR(), kahan_sum(), Index::lin_read(), CTensor< T >::lin_read(), LU_bkw_subst(), lu_decomp(), LU_det(), LU_fwd_subst(), lu_solve(), LU_solve(), CRMatrix< T >::MatVecMult(), CSCMatrix< T >::MatVecMult(), metrmul(), CRMatrix< T >::mult(), CSCMatrix< T >::mult(), Tensor< T >::mult(), BdMatrix< T >::mult_row(), BdMatrix< T >::mult_rows(), CSCMatrix< T >::multf(), norm_1(), BdMatrix< T >::operator TMatrix< T >(), operator!=(), FS_Vector< dims, T >::operator()(), Symm_BdMatrix< T >::operator()(), CTensor< T >::operator()(), BVector< unsigned >::operator()(), TVector< T >::operator()(), TSVector< T >::operator()(), Vector< unsigned >::operator()(), Matrix< T >::operator()(), Symm_BdMatrix< T >::operator*(), CRMatrix< T >::operator*(), CSCMatrix< T >::operator*(), BdMatrix< T >::operator*(), Tensor< T >::operator*(), F_TSMatrix< T >::operator*(), TSMatrix< T >::operator*(), Matrix< T >::operator*(), CRMatrix< T >::operator*=(), CSCMatrix< T >::operator*=(), Tensor< T >::operator+(), operator+(), Tensor< T >::operator+=(), CRMatrix< T >::operator-(), FS_Vector< dims, T >::operator-(), CSCMatrix< T >::operator-(), BdMatrix< T >::operator-(), Tensor< T >::operator-(), operator-(), Tensor< T >::operator-=(), FS_Vector< dims, T >::operator/(), CRMatrix< T >::operator/(), CSCMatrix< T >::operator/(), Tensor< T >::operator/(), FS_Vector< dims, T >::operator/=(), CRMatrix< T >::operator/=(), CSCMatrix< T >::operator/=(), BdMatrix< T >::operator/=(), operator<<(), BVector< unsigned >::operator<=(), CRMatrix< T >::operator=(), CRMatrix< T >::operator==(), CSCMatrix< T >::operator==(), operator==(), BdMatrix< T >::operator==(), TSVector< T >::operator==(), Vector< unsigned >::operator==(), BVector< unsigned >::operator>=(), operator>>(), FS_Vector< dims, T >::operator[](), own_ev_(), QMR(), real(), BdMatrix< T >::reconfig(), BVector< unsigned >::remove(), CSCMatrix< T >::resize(), BdMatrix< T >::resize(), 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(), Symm_BdMatrix< T >::setval(), TVector< T >::setval(), CRMatrix< T >::size(), ILU0_Symm_BdMatrixPreconditioner< T >::solve(), DILU_BdMatrixPreconditioner< T >::solve(), ILU0_BdMatrixPreconditioner< T >::solve(), sv_decomp(), sv_decomp_backsub(), Tensor< T >::Tensor(), F_TMatrix< T >::trace(), BdMatrix< T >::trace(), TMatrix< T >::trace(), F_TMatrix< T >::trans(), F_BandMatrix< T >::transMult(), CRMatrix< T >::transMult(), CSCMatrix< T >::transMult(), CTensor< T >::transpose(), Update(), DiagPreconditioner< T, MatrixType >::update(), ILU0_Symm_BdMatrixPreconditioner< T >::update(), DILU_BdMatrixPreconditioner< T >::update(), ILU0_BdMatrixPreconditioner< T >::update(), Vector< unsigned >::Vector(), and while().

index

index

Initial value:

{
   int n = index.size()

Definition at line 69 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(), F_BandMatrix< T >::operator()(), CRMatrix< T >::operator()(), BVector< unsigned >::operator()(), F_TMatrix< T >::operator()(), TMatrix< T >::operator()(), F_Matrix< T >::operator()(), Matrix< 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(), F_TMatrix< T >::set_row(), TMatrix< T >::set_row(), TMatrix< T >::set_row_partial(), F_BandMatrix< T >::setval(), CRMatrix< T >::setval(), CSCMatrix< T >::setval(), and while().

mu

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

Definition at line 97 of file LM_fit.h.

p

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

Definition at line 97 of file LM_fit.h.

Referenced by BiCG(), BiCGSTAB(), 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 97 of file LM_fit.h.

Referenced by arydiv_(), asinh(), bernoulli(), BiCG(), BiCGSTAB(), binom(), binomial(), CG(), CG2(), cgamma_(), CGS(), CHEBY(), chi_s(), FS_Vector< dims, T >::contains(), eig(), erfc3(), fac(), fact(), FS_Vector< dims, T >::fill(), FS_Vector< dims, T >::FS_Vector(), IR(), ldgamma(), LU_bkw_subst(), LU_fwd_subst(), lu_solve(), LU_solve(), lu_solve_expert(), FS_Vector< dims, T >::max(), FS_Vector< dims, T >::min(), CSCMatrix< T >::mult(), CSCMatrix< T >::multf(), FS_Vector< dims, T >::operator BVector< T >(), FS_Vector< dims, T >::operator()(), F_TSMatrix< T >::operator*(), FS_Vector< dims, T >::operator-(), FS_Vector< dims, T >::operator/(), FS_Vector< dims, T >::operator/=(), FS_Vector< dims, T >::operator=(), FS_Vector< dims, T >::operator==(), FS_Vector< dims, T >::operator[](), poisson(), pow(), QMR(), TBCI::sign(), store_(), FS_Vector< dims, T >::sum(), sv_decomp(), sv_decomp_backsub(), svd_solve(), trinomial(), Update(), and WEAK().

xlen

int xlen = x.size()

Definition at line 105 of file LM_fit.h.