hs_sparse.cpp File Reference

(* $Id: hs_sparse.cpp 542 2012-07-10 21:04:06Z johnoel $)

#include <fvar.hpp>
#include "hs.h"

Go to the source code of this file.

Defines
#define	XCONST   const
#define	USE_ADJOINT_CODE
Functions
void	report_derivatives (const dvar_vector &x)
static void	xxx (int i)
static void	xxx (ivector &i)
static void	xxxv (ivector &x)
int	varchol (XCONST dvar_hs_smatrix &A, XCONST hs_symbolic &S, dvar_hs_smatrix &L, dcompressed_triplet &sparse_triplet2)
int	tmpxchol (XCONST hs_smatrix &A, XCONST hs_symbolic &S, hs_smatrix &L, ivector &xcount)
hs_smatrix	cs_multiply (const hs_smatrix &A, const hs_smatrix &B)
int	cholnew (XCONST hs_smatrix &A, XCONST hs_symbolic &S, hs_smatrix &L)
double	cs_cumsum (ivector &p, ivector &c, int n)
int	cs_wclear (int mark, int lemax, ivector &w, int n)
void	cs_transpose (XCONST hs_smatrix &_A, int values, hs_smatrix &C)
void	cs_transpose (XCONST dvar_hs_smatrix &_A, int values, dvar_hs_smatrix &C)
dmatrix	make_dmatrix (dcompressed_triplet &M, int n, int m)
dvar_matrix	make_sdvar_matrix (dvar_compressed_triplet &M)
dvar_matrix	make_dvar_matrix (dvar_compressed_triplet &M)
dmatrix	make_dmatrix (dcompressed_triplet &M)
dmatrix	make_sdmatrix (dcompressed_triplet &M)
dvar_matrix	make_dvar_matrix (dvar_compressed_triplet &M, int n, int m)
dvariable	ln_det (XCONST dvar_hs_smatrix &M)
dvariable	ln_det (XCONST dvar_hs_smatrix &_M, int &sgn)
double	ln_det (XCONST hs_smatrix &_M)
double	ln_det (XCONST hs_smatrix &_M, int &sgn)
int	cs_ereach (XCONST hs_smatrix &_A, int k, XCONST ivector &parent, ivector &s, ivector &w)
int	cs_ereach (XCONST dvar_hs_smatrix &_A, int k, XCONST ivector &parent, ivector &s, ivector &w)
void	hs_symperm (XCONST hs_smatrix &_A, XCONST ivector &pinv, hs_smatrix &C)
void	hs_symperm (XCONST dvar_hs_smatrix &_A, XCONST ivector &pinv, dvar_hs_smatrix &C)
void	myacc (int &p, int Lpi1, int ci, const ivector &Li, dvector &x, const dvector &Lx, const double &lki)
int	chol (XCONST hs_smatrix &A, XCONST hs_symbolic &T, hs_smatrix &L)
int	tmpxchol1 (XCONST hs_smatrix &A, XCONST hs_symbolic &T, hs_smatrix &L, ivector &nxcount)
dvector	cs_ipvec (XCONST ivector &p, XCONST dvector &b)
dvector	cs_pvec (XCONST ivector &p, XCONST dvector &b)
dvar_vector	cs_pvec (XCONST ivector &p, XCONST dvar_vector &b)
dvector	cs_lsolve (XCONST hs_smatrix &LL, XCONST dvector &y)
dvar_vector	cs_lsolve (XCONST dvar_hs_smatrix &LL, XCONST dvar_vector &y)
dvar_vector	cs_lsolve (XCONST dvar_hs_smatrix &LL, XCONST dvector &y)
dvector	cs_ltsolve (XCONST hs_smatrix &LL, XCONST dvector &y)
dvar_vector	cs_ltsolve (XCONST dvar_hs_smatrix &LL, XCONST dvar_vector &y)
int	cs_diag (int i, int j, double aij, void *other)
int	cs_diag (int i, int j, const prevariable &aij, void *other)
int	cs_fkeep (hs_smatrix &A, int(*fkeep)(int, int, double, void *), void *other)
int	cs_fkeep (dvar_hs_smatrix &A, int(*fkeep)(int, int, const prevariable &, void *), void *other)
int	cs_scatter (XCONST hs_smatrix &AA, int j, double beta, ivector &w, dvector &x, int mark, hs_smatrix &C, int nz)
int	cs_scatter (XCONST dvar_hs_smatrix &AA, int j, double beta, ivector &w, dvar_vector &x, int mark, dvar_hs_smatrix &C, int nz)
hs_smatrix	cs_add (XCONST hs_smatrix &AA, XCONST hs_smatrix &BB, double alpha, double beta)
dvar_hs_smatrix	cs_add (XCONST dvar_hs_smatrix &AA, XCONST dvar_hs_smatrix &BB, double alpha, double beta)
int	cs_tdfs (int j, int k, ivector &head, XCONST ivector &next, ivector &post, ivector &stack)
ivector	cs_amd (XCONST hs_smatrix &A)
ivector	cs_amd (XCONST dvar_hs_smatrix &A)
ivector	cs_etree (XCONST hs_smatrix &_A)
ivector	cs_etree (XCONST dvar_hs_smatrix &_A)
ivector	cs_post (XCONST ivector &parent, int n)
int	cs_leaf (int i, int j, XCONST ivector &first, ivector &maxfirst, ivector &prevleaf, ivector &ancestor, int *jleaf)
ivector	cs_counts (XCONST hs_smatrix &A, XCONST ivector &parent, XCONST ivector &post)
ivector	cs_counts (XCONST dvar_hs_smatrix &A, XCONST ivector &parent, XCONST ivector &post)
ivector	cs_pinv (XCONST ivector &p, int n)
istream &	operator>> (istream &is, dcompressed_triplet &M)
istream &	operator>> (istream &is, dvar_compressed_triplet &M)
hs_smatrix *	return_choleski_decomp (dcompressed_triplet &st)
dvar_hs_smatrix *	return_choleski_decomp (dvar_compressed_triplet &st)
dvector	return_choleski_decomp_solve (dcompressed_triplet &st, dvector &eps)
dvector	return_choleski_factor_solve (hs_smatrix *PL, dvector &eps)
dvar_vector	return_choleski_factor_solve (dvar_hs_smatrix *PL, dvector &eps)
dvector	solve (dcompressed_triplet &st, dmatrix &Hess, dvector &grad)
dvector	solve (dcompressed_triplet &st, dmatrix &Hess, dvector &grad, hs_symbolic &S)
dvector	solve (const dcompressed_triplet &_st, const dvector &_grad, const hs_symbolic &S)
dvector	solve (const dcompressed_triplet &_st, const dvector &_grad, const hs_symbolic &S, int &ierr)
int	allocated (const dcompressed_triplet &_st)
int	allocated (const dvar_compressed_triplet &_st)
dvariable	ln_det (dvar_compressed_triplet &VM)
dvariable	ln_det (dvar_compressed_triplet &VM, hs_symbolic &S)
dvariable	ln_det (dvar_compressed_triplet &VM, hs_symbolic &S, dcompressed_triplet &s)
double	ln_det (const dcompressed_triplet &VVM, const hs_symbolic &T)
double	ln_det (const dcompressed_triplet &VVM)
static void	dfcholeski_sparse (void)
int	chol (XCONST dvar_hs_smatrix &AA, XCONST hs_symbolic &T, dvar_hs_smatrix &LL)
void	report_dvar_vector_derivatives (void)
void	get_inverse_sparse_hessian (dcompressed_triplet &st, hs_symbolic &S, uostream &ofs1, ofstream &ofs, int usize, int xsize, dvector &u)
hs_smatrix	operator* (const hs_smatrix &A, const hs_smatrix &B)
dcompressed_triplet	make_dcompressed_triplet (const dmatrix &M)
hs_smatrix	make_hs_smatrix (const dmatrix &M)
ostream &	operator<< (const ostream &_ofs, const hs_smatrix &M)
dmatrix	make_dmatrix (const hs_smatrix &M)
Variables
int	check_flag = 0

---

Define Documentation

#define XCONST   const

Definition at line 17 of file hs_sparse.cpp.

#define USE_ADJOINT_CODE

Definition at line 21 of file hs_sparse.cpp.

---

Function Documentation

void report_derivatives

(

const dvar_vector &

x

)

Definition at line 3548 of file hs_sparse.cpp.

static void xxx

(

int

i

)

[static]

Definition at line 24 of file hs_sparse.cpp.

static void xxx

(

ivector &

i

)

[static]

Definition at line 25 of file hs_sparse.cpp.

static void xxxv

(

ivector &

x

)

[static]

Definition at line 26 of file hs_sparse.cpp.

Referenced by tmpxchol1().

int varchol	(	XCONST dvar_hs_smatrix &	A,
		XCONST hs_symbolic &	S,
		dvar_hs_smatrix &	L,
		dcompressed_triplet &	sparse_triplet2
	)

Definition at line 3051 of file hs_sparse.cpp.

Referenced by ln_det().

int tmpxchol	(	XCONST hs_smatrix &	A,
		XCONST hs_symbolic &	S,
		hs_smatrix &	L,
		ivector &	xcount
	)

hs_smatrix cs_multiply	(	const hs_smatrix &	A,
		const hs_smatrix &	B
	)

Definition at line 3588 of file hs_sparse.cpp.

Referenced by operator*().

int cholnew	(	XCONST hs_smatrix &	A,
		XCONST hs_symbolic &	S,
		hs_smatrix &	L
	)

Definition at line 2963 of file hs_sparse.cpp.

double cs_cumsum	(	ivector &	p,
		ivector &	c,
		int	n
	)

Definition at line 58 of file hs_sparse.cpp.

Referenced by cs_transpose(), and hs_symperm().

int cs_wclear	(	int	mark,
		int	lemax,
		ivector &	w,
		int	n
	)

Definition at line 74 of file hs_sparse.cpp.

Referenced by cs_amd().

void cs_transpose	(	XCONST hs_smatrix &	_A,
		int	values,
		hs_smatrix &	C
	)

Definition at line 86 of file hs_sparse.cpp.

Referenced by cs_amd(), and cs_counts().

void cs_transpose	(	XCONST dvar_hs_smatrix &	_A,
		int	values,
		dvar_hs_smatrix &	C
	)

Definition at line 118 of file hs_sparse.cpp.

dmatrix make_dmatrix	(	dcompressed_triplet &	M,
		int	n,
		int	m
	)

Definition at line 152 of file hs_sparse.cpp.

dvar_matrix make_sdvar_matrix

(

dvar_compressed_triplet &

M

)

Definition at line 168 of file hs_sparse.cpp.

dvar_matrix make_dvar_matrix

(

dvar_compressed_triplet &

M

)

Definition at line 186 of file hs_sparse.cpp.

dmatrix make_dmatrix

(

dcompressed_triplet &

M

)

Definition at line 201 of file hs_sparse.cpp.

Referenced by laplace_approximation_calculator::do_newton_raphson_banded().

dmatrix make_sdmatrix

(

dcompressed_triplet &

M

)

Definition at line 215 of file hs_sparse.cpp.

dvar_matrix make_dvar_matrix	(	dvar_compressed_triplet &	M,
		int	n,
		int	m
	)

Definition at line 233 of file hs_sparse.cpp.

dvariable ln_det

(

XCONST dvar_hs_smatrix &

M

)

Definition at line 666 of file hs_sparse.cpp.

dvariable ln_det	(	XCONST dvar_hs_smatrix &	_M,
		int &	sgn
	)

Definition at line 673 of file hs_sparse.cpp.

double ln_det

(

XCONST hs_smatrix &

_M

)

Definition at line 679 of file hs_sparse.cpp.

double ln_det	(	XCONST hs_smatrix &	_M,
		int &	sgn
	)

Definition at line 685 of file hs_sparse.cpp.

int cs_ereach	(	XCONST hs_smatrix &	_A,
		int	k,
		XCONST ivector &	parent,
		ivector &	s,
		ivector &	w
	)

Definition at line 743 of file hs_sparse.cpp.

Referenced by chol(), cholnew(), dfcholeski_sparse(), tmpxchol1(), and varchol().

int cs_ereach	(	XCONST dvar_hs_smatrix &	_A,
		int	k,
		XCONST ivector &	parent,
		ivector &	s,
		ivector &	w
	)

Definition at line 771 of file hs_sparse.cpp.

void hs_symperm	(	XCONST hs_smatrix &	_A,
		XCONST ivector &	pinv,
		hs_smatrix &	C
	)

Definition at line 800 of file hs_sparse.cpp.

Referenced by chol(), cholnew(), dfcholeski_sparse(), tmpxchol1(), and varchol().

void hs_symperm	(	XCONST dvar_hs_smatrix &	_A,
		XCONST ivector &	pinv,
		dvar_hs_smatrix &	C
	)

Definition at line 842 of file hs_sparse.cpp.

void myacc	(	int &	p,
		int	Lpi1,
		int	ci,
		const ivector &	Li,
		dvector &	x,
		const dvector &	Lx,
		const double &	lki
	)

Definition at line 885 of file hs_sparse.cpp.

Referenced by cholnew().

int chol	(	XCONST hs_smatrix &	A,
		XCONST hs_symbolic &	T,
		hs_smatrix &	L
	)

Definition at line 895 of file hs_sparse.cpp.

Referenced by get_inverse_sparse_hessian(), ln_det(), return_choleski_decomp(), return_choleski_decomp_solve(), and solve().

int tmpxchol1	(	XCONST hs_smatrix &	A,
		XCONST hs_symbolic &	T,
		hs_smatrix &	L,
		ivector &	nxcount
	)

Definition at line 976 of file hs_sparse.cpp.

dvector cs_ipvec	(	XCONST ivector &	p,
		XCONST dvector &	b
	)

Definition at line 1105 of file hs_sparse.cpp.

Referenced by get_inverse_sparse_hessian(), return_choleski_decomp_solve(), return_choleski_factor_solve(), and solve().

dvector cs_pvec	(	XCONST ivector &	p,
		XCONST dvector &	b
	)

Definition at line 1120 of file hs_sparse.cpp.

Referenced by get_inverse_sparse_hessian(), return_choleski_decomp_solve(), return_choleski_factor_solve(), and solve().

dvar_vector cs_pvec	(	XCONST ivector &	p,
		XCONST dvar_vector &	b
	)

Definition at line 1134 of file hs_sparse.cpp.

dvector cs_lsolve	(	XCONST hs_smatrix &	LL,
		XCONST dvector &	y
	)

Definition at line 1149 of file hs_sparse.cpp.

Referenced by get_inverse_sparse_hessian(), return_choleski_decomp_solve(), return_choleski_factor_solve(), and solve().

dvar_vector cs_lsolve	(	XCONST dvar_hs_smatrix &	LL,
		XCONST dvar_vector &	y
	)

Definition at line 1174 of file hs_sparse.cpp.

dvar_vector cs_lsolve	(	XCONST dvar_hs_smatrix &	LL,
		XCONST dvector &	y
	)

Definition at line 1199 of file hs_sparse.cpp.

dvector cs_ltsolve	(	XCONST hs_smatrix &	LL,
		XCONST dvector &	y
	)

Definition at line 1225 of file hs_sparse.cpp.

Referenced by get_inverse_sparse_hessian(), and solve().

dvar_vector cs_ltsolve	(	XCONST dvar_hs_smatrix &	LL,
		XCONST dvar_vector &	y
	)

Definition at line 1249 of file hs_sparse.cpp.

int cs_diag	(	int	i,
		int	j,
		double	aij,
		void *	other
	)

Definition at line 1274 of file hs_sparse.cpp.

Referenced by cs_amd().

int cs_diag	(	int	i,
		int	j,
		const prevariable &	aij,
		void *	other
	)

Definition at line 1275 of file hs_sparse.cpp.

int cs_fkeep	(	hs_smatrix &	A,
		int(*)(int, int, double, void *)	fkeep,
		void *	other
	)

Definition at line 1278 of file hs_sparse.cpp.

Referenced by cs_amd().

int cs_fkeep	(	dvar_hs_smatrix &	A,
		int(*)(int, int, const prevariable &, void *)	fkeep,
		void *	other
	)

Definition at line 1304 of file hs_sparse.cpp.

int cs_scatter	(	XCONST hs_smatrix &	AA,
		int	j,
		double	beta,
		ivector &	w,
		dvector &	x,
		int	mark,
		hs_smatrix &	C,
		int	nz
	)

Definition at line 1331 of file hs_sparse.cpp.

Referenced by cs_add(), and cs_multiply().

int cs_scatter	(	XCONST dvar_hs_smatrix &	AA,
		int	j,
		double	beta,
		ivector &	w,
		dvar_vector &	x,
		int	mark,
		dvar_hs_smatrix &	C,
		int	nz
	)

Definition at line 1380 of file hs_sparse.cpp.

hs_smatrix cs_add	(	XCONST hs_smatrix &	AA,
		XCONST hs_smatrix &	BB,
		double	alpha,
		double	beta
	)

Definition at line 1405 of file hs_sparse.cpp.

Referenced by cs_amd().

dvar_hs_smatrix cs_add	(	XCONST dvar_hs_smatrix &	AA,
		XCONST dvar_hs_smatrix &	BB,
		double	alpha,
		double	beta
	)

Definition at line 1450 of file hs_sparse.cpp.

int cs_tdfs	(	int	j,
		int	k,
		ivector &	head,
		XCONST ivector &	next,
		ivector &	post,
		ivector &	stack
	)

Definition at line 1496 of file hs_sparse.cpp.

Referenced by cs_amd(), and cs_post().

ivector cs_amd

(

XCONST hs_smatrix &

A

)

Definition at line 1519 of file hs_sparse.cpp.

ivector cs_amd

(

XCONST dvar_hs_smatrix &

A

)

Definition at line 1852 of file hs_sparse.cpp.

ivector cs_etree

(

XCONST hs_smatrix &

_A

)

Definition at line 2186 of file hs_sparse.cpp.

ivector cs_etree

(

XCONST dvar_hs_smatrix &

_A

)

Definition at line 2218 of file hs_sparse.cpp.

ivector cs_post	(	XCONST ivector &	parent,
		int	n
	)

Definition at line 2251 of file hs_sparse.cpp.

int cs_leaf	(	int	i,
		int	j,
		XCONST ivector &	first,
		ivector &	maxfirst,
		ivector &	prevleaf,
		ivector &	ancestor,
		int *	jleaf
	)

Definition at line 2280 of file hs_sparse.cpp.

Referenced by cs_counts().

ivector cs_counts	(	XCONST hs_smatrix &	A,
		XCONST ivector &	parent,
		XCONST ivector &	post
	)

Definition at line 2301 of file hs_sparse.cpp.

ivector cs_counts	(	XCONST dvar_hs_smatrix &	A,
		XCONST ivector &	parent,
		XCONST ivector &	post
	)

Definition at line 2356 of file hs_sparse.cpp.

ivector cs_pinv	(	XCONST ivector &	p,
		int	n
	)

Definition at line 2413 of file hs_sparse.cpp.

istream& operator>>	(	istream &	is,
		dcompressed_triplet &	M
	)

Definition at line 2603 of file hs_sparse.cpp.

istream& operator>>	(	istream &	is,
		dvar_compressed_triplet &	M
	)

Definition at line 2616 of file hs_sparse.cpp.

hs_smatrix* return_choleski_decomp

(

dcompressed_triplet &

st

)

Definition at line 2629 of file hs_sparse.cpp.

Referenced by calculate_importance_sample_shess().

dvar_hs_smatrix* return_choleski_decomp

(

dvar_compressed_triplet &

st

)

Definition at line 2646 of file hs_sparse.cpp.

dvector return_choleski_decomp_solve	(	dcompressed_triplet &	st,
		dvector &	eps
	)

Definition at line 2663 of file hs_sparse.cpp.

dvector return_choleski_factor_solve	(	hs_smatrix *	PL,
		dvector &	eps
	)

Definition at line 2686 of file hs_sparse.cpp.

Referenced by calculate_importance_sample_shess().

dvar_vector return_choleski_factor_solve	(	dvar_hs_smatrix *	PL,
		dvector &	eps
	)

Definition at line 2703 of file hs_sparse.cpp.

dvector solve	(	dcompressed_triplet &	st,
		dmatrix &	Hess,
		dvector &	grad
	)

Definition at line 2721 of file hs_sparse.cpp.

dvector solve	(	dcompressed_triplet &	st,
		dmatrix &	Hess,
		dvector &	grad,
		hs_symbolic &	S
	)

Definition at line 2751 of file hs_sparse.cpp.

dvector solve	(	const dcompressed_triplet &	_st,
		const dvector &	_grad,
		const hs_symbolic &	S
	)

Definition at line 2786 of file hs_sparse.cpp.

dvector solve	(	const dcompressed_triplet &	_st,
		const dvector &	_grad,
		const hs_symbolic &	S,
		int &	ierr
	)

Definition at line 2812 of file hs_sparse.cpp.

int allocated

(

const dcompressed_triplet &

_st

)

Definition at line 2846 of file hs_sparse.cpp.

int allocated

(

const dvar_compressed_triplet &

_st

)

Definition at line 2851 of file hs_sparse.cpp.

dvariable ln_det

(

dvar_compressed_triplet &

VM

)

Definition at line 2856 of file hs_sparse.cpp.

dvariable ln_det	(	dvar_compressed_triplet &	VM,
		hs_symbolic &	S
	)

Definition at line 2873 of file hs_sparse.cpp.

dvariable ln_det	(	dvar_compressed_triplet &	VM,
		hs_symbolic &	S,
		dcompressed_triplet &	s
	)

Definition at line 2892 of file hs_sparse.cpp.

double ln_det	(	const dcompressed_triplet &	VVM,
		const hs_symbolic &	T
	)

Definition at line 2922 of file hs_sparse.cpp.

double ln_det

(

const dcompressed_triplet &

VVM

)

Definition at line 2944 of file hs_sparse.cpp.

static void dfcholeski_sparse

(

void

)

[static]

Definition at line 3161 of file hs_sparse.cpp.

Referenced by varchol().

int chol	(	XCONST dvar_hs_smatrix &	AA,
		XCONST hs_symbolic &	T,
		dvar_hs_smatrix &	LL
	)

Definition at line 3419 of file hs_sparse.cpp.

void report_dvar_vector_derivatives

(

void

)

Definition at line 3539 of file hs_sparse.cpp.

Referenced by report_derivatives().

void get_inverse_sparse_hessian	(	dcompressed_triplet &	st,
		hs_symbolic &	S,
		uostream &	ofs1,
		ofstream &	ofs,
		int	usize,
		int	xsize,
		dvector &	u
	)

Definition at line 3558 of file hs_sparse.cpp.

Referenced by function_minimizer::hess_routine_noparallel_random_effects().

hs_smatrix operator*	(	const hs_smatrix &	A,
		const hs_smatrix &	B
	)

Definition at line 3637 of file hs_sparse.cpp.

dcompressed_triplet make_dcompressed_triplet

(

const dmatrix &

M

)

Definition at line 3641 of file hs_sparse.cpp.

Referenced by make_hs_smatrix().

hs_smatrix make_hs_smatrix

(

const dmatrix &

M

)

Definition at line 3689 of file hs_sparse.cpp.

ostream& operator<<	(	const ostream &	_ofs,
		const hs_smatrix &	M
	)

Definition at line 3694 of file hs_sparse.cpp.

dmatrix make_dmatrix

(

const hs_smatrix &

M

)

Definition at line 3705 of file hs_sparse.cpp.

---

Variable Documentation

int check_flag = 0

Definition at line 2890 of file hs_sparse.cpp.