Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Review for DM-5197: new shapelet PSF approximation code #21

Merged
merged 9 commits into from May 13, 2016

Conversation

TallJimbo
Copy link
Member

No description provided.

@@ -238,8 +238,7 @@ struct CModelControl {
LSST_CONTROL_FIELD(
psfName,

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

psfNamePrefix?

Copy link
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Hmm. Maybe. Field names like these are really sequences of namespace-like terms, so I'm not sure whether calling them a "prefix" is warranted (it feels like calling the directory part of a full filename a prefix). Regardless, it's clear that the documentation string below needs to be improved.

@pschella
Copy link

About the range checking. Doesn't Array have a shape?

@pschella
Copy link

pschella commented May 12, 2016

And concerning getCoefficients()[0] =. Even though access to the underlying object is sometimes needed, I still maintain it would be nicer to say setCoefficient(0, value) in this case.

afw::image::Image<Scalar> const & image
) : OptimizerObjective(image.getBBox().getArea(), 4),
_minRadius(ctrl.minRadius),
_maxRadius(ctrl.maxRadiusBoxFraction * std::sqrt(this->dataSize)),

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

You are assuming a square box here?

Copy link
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

No, this is just intended to be a geometric mean of the width and height. But I do think the box will be close enough to square that I'm not concerned about the difference between geometric and arithmetic mean (there really isn't a good reason to pick one or the other).

This will eventually allow us to use analytic derivatives instead
of finite differences.
The current algorithm is not very robust and is sometimes
catastrophically slow.  However, it's able to fit much more general
models than the HSC-side one we'll be adding shortly, so we'd like to
keep it around.
This task is slated for eventual removal, but that's more work than we
want to do on this ticket.  For now, we just need to remove a name
collision.
@TallJimbo
Copy link
Member Author

About the range checking. Doesn't Array have a shape?

Yes, but what matters is that the array shapes agree with the shape of the Image passed to applyPixelFunctor.

@TallJimbo
Copy link
Member Author

And concerning getCoefficients()[0] =. Even though access to the underlying object is sometimes needed, I still maintain it would be nicer to say setCoefficient(0, value) in this case.

My preference would probably be for adding an operator[] overload to ShapeletFunction, actually. But I think either is beyond the scope of this issue.

Many of these passed on my Ubuntu box, but not
on CentOS CI.
@TallJimbo TallJimbo merged commit 244dfb6 into master May 13, 2016
@ktlim ktlim deleted the tickets/DM-5197 branch August 25, 2018 05:49
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
2 participants