LSST Applications g0f08755f38+9c285cab97,g1635faa6d4+13f3999e92,g1653933729+a8ce1bb630,g1a0ca8cf93+bf6eb00ceb,g28da252d5a+0829b12dee,g29321ee8c0+5700dc9eac,g2bbee38e9b+9634bc57db,g2bc492864f+9634bc57db,g2cdde0e794+c2c89b37c4,g3156d2b45e+41e33cbcdc,g347aa1857d+9634bc57db,g35bb328faa+a8ce1bb630,g3a166c0a6a+9634bc57db,g3e281a1b8c+9f2c4e2fc3,g414038480c+077ccc18e7,g41af890bb2+fde0dd39b6,g5fbc88fb19+17cd334064,g781aacb6e4+a8ce1bb630,g80478fca09+55a9465950,g82479be7b0+d730eedb7d,g858d7b2824+9c285cab97,g9125e01d80+a8ce1bb630,g9726552aa6+10f999ec6a,ga5288a1d22+2a84bb7594,gacf8899fa4+c69c5206e8,gae0086650b+a8ce1bb630,gb58c049af0+d64f4d3760,gc28159a63d+9634bc57db,gcf0d15dbbd+4b7d09cae4,gda3e153d99+9c285cab97,gda6a2b7d83+4b7d09cae4,gdaeeff99f8+1711a396fd,ge2409df99d+5e831397f4,ge79ae78c31+9634bc57db,gf0baf85859+147a0692ba,gf3967379c6+41c94011de,gf3fb38a9a8+8f07a9901b,gfb92a5be7c+9c285cab97,w.2024.46
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: