#include <evaluate.h>

Public Member Functions
	TermsGradient (vecdptr ymc_i=nullptr, double xx_weight_i=0, double xy_weight_i=0, double yy_weight_i=0, double rho_factor_i=0, double sig_x_inv_i=0, double sig_y_inv_i=0, double rho_xy_factor_i=0, double sig_x_src_div_conv_i=0, double sig_y_src_div_conv_i=0, double drho_c_dsig_x_src_i=0, double drho_c_dsig_y_src_i=0, double drho_c_drho_s_i=0, double xmc_t_xy_weight_i=0)

void	set (const Terms &terms, const Ellipse &ell_src, const Ellipse &ell_psf, const Ellipse &ell, vecdptr ymc_)

Public Attributes
vecdptr	ymc = nullptr

double	xx_weight = 0

double	xy_weight = 0

double	yy_weight = 0

double	rho_factor = 0

double	sig_x_inv = 0

double	sig_y_inv = 0

double	rho_xy_factor = 0

double	sig_x_src_div_conv = 0

double	sig_y_src_div_conv = 0

double	drho_c_dsig_x_src = 0

double	drho_c_dsig_y_src = 0

double	drho_c_drho_s = 0

double	xmc_t_xy_weight = 0

Detailed Description

Definition at line 231 of file evaluate.h.

Constructor & Destructor Documentation

◆ TermsGradient()

lsst::gauss2d::detail::TermsGradient::TermsGradient	(	vecdptr	ymc_i = nullptr,
		double	xx_weight_i = 0,
		double	xy_weight_i = 0,
		double	yy_weight_i = 0,
		double	rho_factor_i = 0,
		double	sig_x_inv_i = 0,
		double	sig_y_inv_i = 0,
		double	rho_xy_factor_i = 0,
		double	sig_x_src_div_conv_i = 0,
		double	sig_y_src_div_conv_i = 0,
		double	drho_c_dsig_x_src_i = 0,
		double	drho_c_dsig_y_src_i = 0,
		double	drho_c_drho_s_i = 0,
		double	xmc_t_xy_weight_i = 0 )

inlineexplicit

Definition at line 248 of file evaluate.h.

            : ymc(std::move(ymc_i)),
              xx_weight(xx_weight_i),
              xy_weight(xy_weight_i),
              yy_weight(yy_weight_i),
              rho_factor(rho_factor_i),
              sig_x_inv(sig_x_inv_i),
              sig_y_inv(sig_y_inv_i),
              rho_xy_factor(rho_xy_factor_i),
              sig_x_src_div_conv(sig_x_src_div_conv_i),
              sig_y_src_div_conv(sig_y_src_div_conv_i),
              drho_c_dsig_x_src(drho_c_dsig_x_src_i),
              drho_c_dsig_y_src(drho_c_dsig_y_src_i),
              drho_c_drho_s(drho_c_drho_s_i),
              xmc_t_xy_weight(xmc_t_xy_weight_i) {}

Member Function Documentation

◆ set()

void lsst::gauss2d::detail::TermsGradient::set	(	const Terms &	terms,
		const Ellipse &	ell_src,
		const Ellipse &	ell_psf,
		const Ellipse &	ell,
		vecdptr	ymc_ )

inline

Definition at line 269 of file evaluate.h.

                           {
        double sig_x = ell.get_sigma_x();
        double sig_y = ell.get_sigma_y();
        double rho = ell.get_rho();
 
        this->ymc = std::move(ymc_);
        const double sig_xy = sig_x * sig_y;
        this->xx_weight = terms.xx;
        this->xy_weight = terms.xy;
        this->yy_weight = terms.yy;
        this->rho_factor = rho / (1. - rho * rho);
        this->sig_x_inv = 1 / sig_x;
        this->sig_y_inv = 1 / sig_y;
        this->rho_xy_factor = 1. / (1. - rho * rho) / sig_xy;
        /*
    sigma_conv = sqrt(sigma_src^2 + sigma_psf^2)
    https://www.wolframalpha.com/input/?i=d(sqrt(x%5E2%2By%5E2))%2Fdx
    dsigma_conv/dsigma_src = sigma_src/sigma_conv
    df/dsigma_src = df/d_sigma_conv * dsigma_conv/dsigma_src
 
    rho_conv*sigmaxy_conv = rho_src*sigmaxy_src + rho_psf*sigmaxy_psf
    rho_conv = (rho_src*sigmaxy_src + rho_psf*sigmaxy_psf)/sigmaxy_conv
    drho_conv/drho_src = sigmaxy_src/sigmaxy_conv
 
    drho_conv/dsigmax_src = rho_src*sigmay_src/sigmay_conv*sigmax_psf^2/sigmax_conv^3 (*sigmax_src/sigmax_src)
                          = rho_conv/sigmax_src*(sigmax_psf/sigmax_conv)^2
        + rho_psf*sigmax_psf*sigmay_psf*sigmax_src/sigmay_conv/sigmax_conv^3
 
        */
        double sig_x_src = ell_src.get_sigma_x();
        double sig_y_src = ell_src.get_sigma_y();
        double rho_src = ell_src.get_rho();
        double sig_x_psf = ell_psf.get_sigma_x();
        double sig_y_psf = ell_psf.get_sigma_y();
        double rho_psf = ell_psf.get_rho();
 
        this->sig_x_src_div_conv = sig_x_src / sig_x;
        this->sig_y_src_div_conv = sig_y_src / sig_y;
 
        double covar_psf_dsig_xy = rho_psf * sig_x_psf * sig_y_psf / sig_xy;
        double drho = 0;
        if (sig_x_psf > 0) {
            double sig_p_ratio = sig_x_psf / sig_x;
            drho = rho_src * this->sig_y_src_div_conv * sig_p_ratio * sig_p_ratio / sig_x
                   - covar_psf_dsig_xy * this->sig_x_src_div_conv / sig_x;
        }
        this->drho_c_dsig_x_src = drho;
        if (sig_y_psf > 0) {
            double sig_p_ratio = sig_y_psf / sig_y;
            drho = rho_src * this->sig_x_src_div_conv * sig_p_ratio * sig_p_ratio / sig_y
                   - covar_psf_dsig_xy * this->sig_y_src_div_conv / sig_y;
        }
        this->drho_c_dsig_y_src = drho;
        this->drho_c_drho_s = sig_x_src * sig_y_src / sig_xy;
    }

Member Data Documentation

◆ drho_c_drho_s

double lsst::gauss2d::detail::TermsGradient::drho_c_drho_s = 0

Definition at line 245 of file evaluate.h.

◆ drho_c_dsig_x_src

double lsst::gauss2d::detail::TermsGradient::drho_c_dsig_x_src = 0

Definition at line 243 of file evaluate.h.

◆ drho_c_dsig_y_src

double lsst::gauss2d::detail::TermsGradient::drho_c_dsig_y_src = 0

Definition at line 244 of file evaluate.h.

◆ rho_factor

double lsst::gauss2d::detail::TermsGradient::rho_factor = 0

Definition at line 237 of file evaluate.h.

◆ rho_xy_factor

double lsst::gauss2d::detail::TermsGradient::rho_xy_factor = 0

Definition at line 240 of file evaluate.h.

◆ sig_x_inv

double lsst::gauss2d::detail::TermsGradient::sig_x_inv = 0

Definition at line 238 of file evaluate.h.

◆ sig_x_src_div_conv

double lsst::gauss2d::detail::TermsGradient::sig_x_src_div_conv = 0

Definition at line 241 of file evaluate.h.

◆ sig_y_inv

double lsst::gauss2d::detail::TermsGradient::sig_y_inv = 0

Definition at line 239 of file evaluate.h.

◆ sig_y_src_div_conv

double lsst::gauss2d::detail::TermsGradient::sig_y_src_div_conv = 0

Definition at line 242 of file evaluate.h.

◆ xmc_t_xy_weight

double lsst::gauss2d::detail::TermsGradient::xmc_t_xy_weight = 0

Definition at line 246 of file evaluate.h.

◆ xx_weight

double lsst::gauss2d::detail::TermsGradient::xx_weight = 0

Definition at line 234 of file evaluate.h.

◆ xy_weight

double lsst::gauss2d::detail::TermsGradient::xy_weight = 0

Definition at line 235 of file evaluate.h.

◆ ymc

vecdptr lsst::gauss2d::detail::TermsGradient::ymc = nullptr

Definition at line 233 of file evaluate.h.

◆ yy_weight

double lsst::gauss2d::detail::TermsGradient::yy_weight = 0

Definition at line 236 of file evaluate.h.

The documentation for this class was generated from the following file:

/j/snowflake/release/lsstsw/stack/lsst-scipipe-10.1.0/Linux/gauss2d/gbf99507273+8c5ae1fdc5/include/lsst/gauss2d/evaluate.h

Public Member Functions

Public Attributes