powfs_utils.h File Reference

Functions
void	genmtch (const parms_t *parms, powfs_t *powfs, const int ipowfs)
void	cog_nea (real *nea, const dmat *ints, const dmat *cogmask, real cogthres, real cogoff, int ntry, rand_t *rstat, real bkgrnd, real bkgrndc, const dmat *bkgrnd2i, const dmat *bkgrnd2ic, real rne)
real	remove_focus_grad (const loc_t *saloc, dmat *grad, real factor)
void	sodium_fit (dmat **sodium, dcell **pgrad, dcell **pi0, dcell **pgx, dcell **pgy, const dcell *i0i, const dccell *sepsf, const dtf_t *dtf, const loc_t *saloc, const dcell *saa, const dcell *srsa, const dcell *srot, const dmat *siglev, const dmat *wvlwts, const dcell *gradoff, real dh, real hs, real htel, real za, real svdthres, int nrep, int save, int use_cache)
void	sodium_fit_wrap (dmat **psodium, dcell **pgrad, dcell **pi0, dcell **pgx, dcell **pgy, const dcell *i0in, const parms_t *parms, powfs_t *powfs, int ipowfs, real r0, real L0, int nrep, int use_cache)

Detailed Description

Contains routines to generate mean short exposure (tip/tilt removed) pixel intensities. Mostly used for LGS pixel intensity for its matched filter. Structure functions from kolmogorov spectrum is used. Not able to take into account outerscale yet.

Todo:

find ways to factor in outerscale effect (use von karman spectrum instead of kolmogorov)

Function Documentation

sodium_fit()

void sodium_fit	(	dmat **	sodium,
		dcell **	pgrad,
		dcell **	pi0,
		dcell **	pgx,
		dcell **	pgy,
		const dcell *	i0i,
		const dccell *	sepsf,
		const dtf_t *	dtf,
		const loc_t *	saloc,
		const dcell *	saa,
		const dcell *	srsa,
		const dcell *	srot,
		const dmat *	siglev,
		const dmat *	wvlwts,
		const dcell *	gradoff,
		real	dh,
		real	hs,
		real	htel,
		real	za,
		real	svdthres,
		int	nrep,
		int	save,
		int	use_cache
	)

Fit i0 to sodium profile using iterative algorithm. The steps are as follows

1. Create sub images for each sodium profile bin
2. Fit such sub-images against i0 to determine the profile
3. Determine subaperture tip/tilt comparing i0 and fitted i0 (using matched filter)
4. Repeat 1-4.

Parameters

sodium	The sodium profile determined by fit
pgrad	The gradients determined by fit.
pi0	The output i0
pgx	The output gx
pgy	The output gy
i0i	The input i0
sepsf	Short exposure PSF
dtf	Detector transfer function
saloc	Saloc
saa	Subaperture area.
srsa	Subaperture to LLT distance
srot	Subaperture to LLT clocking
siglev	Subaperture signal level
wvlwts	Wavelength weights
gradoff	NCPA gradient to be used for pi0,pgx,pgy output.
dh	The sodium profile sampling in meters
hs	LGS focusing height
htel	Telescope hegith
za	Telescope zenith angle
svdthres	SVD threshold
nrep	Number of iterations
save	Save results to file
use_cache	Use cache

sodium_fit_wrap()

void sodium_fit_wrap	(	dmat **	psodium,
		dcell **	pgrad,
		dcell **	pi0,
		dcell **	pgx,
		dcell **	pgy,
		const dcell *	i0in,
		const parms_t *	parms,
		powfs_t *	powfs,
		int	ipowfs,
		real	r0,
		real	L0,
		int	nrep,
		int	use_cache
	)

Fit i0 to sodium profile and replace i0, gx, gy with derived parameters

Parameters

[out]	psodium	sodium profile
[out]	pgrad	estimated actual gradient
[out]	pi0	The output i0
[out]	pgx	The output gx
[out]	pgy	The output gy
[in]	i0in	The input sa intensities. may equal to *pi0
[in]	parms	parms
[in]	powfs	powfs
[in]	ipowfs	ipowfs
[in]	r0	Fried parameter
[in]	L0	outer scale
[in]	nrep	Number of iterations. 1 for mtche, 3 for cog
[in]	use_cache	cache intermediate results.