LSST Applications g063fba187b+cac8b7c890,g0f08755f38+6aee506743,g1653933729+a8ce1bb630,g168dd56ebc+a8ce1bb630,g1a2382251a+b4475c5878,g1dcb35cd9c+8f9bc1652e,g20f6ffc8e0+6aee506743,g217e2c1bcf+73dee94bd0,g28da252d5a+1f19c529b9,g2bbee38e9b+3f2625acfc,g2bc492864f+3f2625acfc,g3156d2b45e+6e55a43351,g32e5bea42b+1bb94961c2,g347aa1857d+3f2625acfc,g35bb328faa+a8ce1bb630,g3a166c0a6a+3f2625acfc,g3e281a1b8c+c5dd892a6c,g3e8969e208+a8ce1bb630,g414038480c+5927e1bc1e,g41af890bb2+8a9e676b2a,g7af13505b9+809c143d88,g80478fca09+6ef8b1810f,g82479be7b0+f568feb641,g858d7b2824+6aee506743,g89c8672015+f4add4ffd5,g9125e01d80+a8ce1bb630,ga5288a1d22+2903d499ea,gb58c049af0+d64f4d3760,gc28159a63d+3f2625acfc,gcab2d0539d+b12535109e,gcf0d15dbbd+46a3f46ba9,gda6a2b7d83+46a3f46ba9,gdaeeff99f8+1711a396fd,ge79ae78c31+3f2625acfc,gef2f8181fd+0a71e47438,gf0baf85859+c1f95f4921,gfa517265be+6aee506743,gfa999e8aa5+17cd334064,w.2024.51
LSST Data Management Base Package
Loading...
Searching...
No Matches
Public Types | Public Member Functions | Static Public Member Functions | Public Attributes | Static Public Attributes | List of all members
lsst::gauss2d::fit::PriorEvaluation Class Reference

Results from the evaluation of a prior probability function. More...

#include <prior.h>

Inheritance diagram for lsst::gauss2d::fit::PriorEvaluation:
lsst::gauss2d::Object

Public Types

typedef std::map< ParamBaseCRef, std::vector< double > > Jacobians
 

Public Member Functions

double compute_dloglike_dx (const ParamBase &param, bool transformed=true) const
 
std::string repr (bool name_keywords=false, std::string_view namespace_separator=Object::CC_NAMESPACE_SEPARATOR) const override
 Return a full, callable string representation of this.
 
std::string str () const override
 Return a brief, human-readable string representation of this.
 
 PriorEvaluation (double loglike, std::vector< double > residuals={}, Jacobians jacobians={}, bool check_size=true)
 
 ~PriorEvaluation ()
 

Static Public Member Functions

static std::string_view null_str (const std::string_view &namespace_separator)
 

Public Attributes

double loglike
 
std::vector< double > residuals
 
Jacobians jacobians
 

Static Public Attributes

static constexpr std::string_view CC_NAMESPACE_SEPARATOR = "::"
 The C++ namespace separator.
 
static constexpr std::string_view NULL_STR_GENERAL = "None"
 
static constexpr std::string_view PY_NAMESPACE_SEPARATOR = "."
 

Detailed Description

Results from the evaluation of a prior probability function.

Definition at line 16 of file prior.h.

Member Typedef Documentation

◆ Jacobians

Definition at line 18 of file prior.h.

Constructor & Destructor Documentation

◆ PriorEvaluation()

lsst::gauss2d::fit::PriorEvaluation::PriorEvaluation ( double loglike,
std::vector< double > residuals = {},
PriorEvaluation::Jacobians jacobians_ = {},
bool check_size = true )

Definition at line 46 of file prior.cc.

48 : loglike(loglike_), residuals(residuals_), jacobians(jacobians_) {
49 const auto n_resid = residuals.size();
50 if (!((n_resid > 0) || (jacobians.size() == 0))) {
51 throw std::invalid_argument("residuals cannot be empty unless jacobians are too");
52 }
53 if (check_size) {
54 for (const auto& [paramref, values] : jacobians) {
55 if (values.size() != n_resid) {
56 throw std::invalid_argument(paramref.get().str()
57 + " jacobian size=" + std::to_string(values.size())
58 + " != n_resid=" + std::to_string(n_resid));
59 }
60 }
61 } else {
62 }
63}
std::vector< double > residuals
Definition prior.h:21
T size(T... args)
T to_string(T... args)

◆ ~PriorEvaluation()

lsst::gauss2d::fit::PriorEvaluation::~PriorEvaluation ( )

Definition at line 65 of file prior.cc.

65{};

Member Function Documentation

◆ compute_dloglike_dx()

double lsst::gauss2d::fit::PriorEvaluation::compute_dloglike_dx ( const ParamBase & param,
bool transformed = true ) const

Definition at line 67 of file prior.cc.

67 {
68 double dll_dx = 0;
69 auto jacobian = this->jacobians.find(param);
70 if (jacobian != this->jacobians.end()) {
71 // jacobian is dresidual/dx
72 // loglike is -residual^2/2
73 // dloglike/dx = -residual*dresidual/dx
74 const auto& jacobian_values = (*jacobian).second;
75 for (size_t idx = 0; idx < jacobian_values.size(); ++idx) {
76 dll_dx -= jacobian_values[idx] * residuals.at(idx);
77 }
78 }
79 if (!transformed) dll_dx /= param.get_transform_derivative();
80 return dll_dx;
81}
T at(T... args)
T end(T... args)
T find(T... args)
@ jacobian
Compute the model Jacobian.

◆ null_str()

static std::string_view lsst::gauss2d::Object::null_str ( const std::string_view & namespace_separator)
inlinestaticinherited

Definition at line 49 of file object.h.

49 {
50 return namespace_separator == CC_NAMESPACE_SEPARATOR ? "nullptr" : NULL_STR_GENERAL;
51 }
static constexpr std::string_view CC_NAMESPACE_SEPARATOR
The C++ namespace separator.
Definition object.h:45
static constexpr std::string_view NULL_STR_GENERAL
Definition object.h:46

◆ repr()

std::string lsst::gauss2d::fit::PriorEvaluation::repr ( bool name_keywords = false,
std::string_view namespace_separator = Object::CC_NAMESPACE_SEPARATOR ) const
overridevirtual

Return a full, callable string representation of this.

Parameters
name_keywordsWhether to prefix arguments with "{name}=", where name is the arg name in the header (as with keyword arguments in Python).
namespace_separatorThe string to use to delimit namespaces, i.e. :: in C++ and . in Python.
Returns
A callable string representation of this, which should return an an identical object to this.
Note
The representation with name_keywords=false must be callable in C++. The representation with name_keywords=true should be callable in Python, if there are any bindings.

Implements lsst::gauss2d::Object.

Definition at line 83 of file prior.cc.

83 {
84 return type_name_str<PriorEvaluation>(false, namespace_separator) + "("
85 + (name_keywords ? "loglike=" : "") + to_string_float(loglike) + ", "
86 + (name_keywords ? "residuals=" : "") + to_string_float_iter(residuals) + ", "
87 + (name_keywords ? "jacobians=" : "") + str_jacobians(jacobians) + ")";
88}
std::string str_jacobians(const std::map< ParamBaseCRef, std::vector< double > > &jacobians, bool python_style=true)
Definition prior.cc:35
std::string to_string_float_iter(const Container< Value > &container, const int precision=6, const bool scientific=true)
Definition to_string.h:28
std::string to_string_float(const T value, const int precision=6, const bool scientific=true)
Definition to_string.h:15

◆ str()

std::string lsst::gauss2d::fit::PriorEvaluation::str ( ) const
overridevirtual

Return a brief, human-readable string representation of this.

Implements lsst::gauss2d::Object.

Definition at line 90 of file prior.cc.

90 {
91 return type_name_str<PriorEvaluation>(true) + "(loglike=" + to_string_float(loglike) + ", residuals="
92 + to_string_float_iter(residuals) + ", jacobians=" + str_jacobians(jacobians) + ")";
93}

Member Data Documentation

◆ CC_NAMESPACE_SEPARATOR

constexpr std::string_view lsst::gauss2d::Object::CC_NAMESPACE_SEPARATOR = "::"
staticconstexprinherited

The C++ namespace separator.

Definition at line 45 of file object.h.

◆ jacobians

Jacobians lsst::gauss2d::fit::PriorEvaluation::jacobians

Definition at line 22 of file prior.h.

◆ loglike

double lsst::gauss2d::fit::PriorEvaluation::loglike

Definition at line 20 of file prior.h.

◆ NULL_STR_GENERAL

constexpr std::string_view lsst::gauss2d::Object::NULL_STR_GENERAL = "None"
staticconstexprinherited

Definition at line 46 of file object.h.

◆ PY_NAMESPACE_SEPARATOR

constexpr std::string_view lsst::gauss2d::Object::PY_NAMESPACE_SEPARATOR = "."
staticconstexprinherited

Definition at line 47 of file object.h.

◆ residuals

std::vector<double> lsst::gauss2d::fit::PriorEvaluation::residuals

Definition at line 21 of file prior.h.


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