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.

Sign up for GitHub

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

Use similar in creation of DiffResults buffer #95

Merged
merged 5 commits into from
Nov 22, 2022
Merged

Use similar in creation of DiffResults buffer #95

Changes from 2 commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
3d14e53
use similar in creation of diffresults buffer
torfjelde Nov 17, 2022
01d2745
use the input to make the DiffResults buffer
torfjelde Nov 17, 2022
00af713
remove the logdensity argument from _diffresults_buffer
torfjelde Nov 19, 2022
9634308
removed type-annotation
torfjelde Nov 22, 2022
39dd7c1
patch version bump
torfjelde Nov 22, 2022
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
3 changes: 2 additions & 1 deletion src/DiffResults_helpers.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -14,7 +14,7 @@ Allocate a DiffResults buffer for a gradient, taking the element type of `x` int
function _diffresults_buffer(ℓ, x)
T = eltype(x)
S = T <: Real ? float(Real) : Float64 # heuristic
DiffResults.MutableDiffResult(zero(S), (Vector{S}(undef, dimension(ℓ)), ))
DiffResults.MutableDiffResult(zero(S), (similar(x, S), ))
Copy link
Collaborator

Choose a reason for hiding this comment

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

Should this still take into account dimension? I'm not sure but I assume there was a reason why it was used instead of x?

Suggested change
DiffResults.MutableDiffResult(zero(S), (similar(x, S), ))
DiffResults.MutableDiffResult(zero(S), (similar(x, S, dimension(ℓ)), ))

Copy link
Owner

Choose a reason for hiding this comment

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

Currently all it uses from xis its element type, and it takes the dimension from the first argument.

I think it could be done differently, using the dimension from x, and then maybe ℓ would not need to be an argument at all.

Copy link
Contributor Author

Choose a reason for hiding this comment

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

using the dimension from x

So similar basically does this. And AFAIK this is only called from the logdensity_and_gradient method, meaning that x should always have the correct dimensions.

Using dimension from ℓ instead will sometimes result in loss of structure of the original type for x, e.g. ComponentArrays.

Copy link
Owner

Choose a reason for hiding this comment

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

I would recommend something like

function _diffresults_buffer(x)
    T = eltype(x)
    S = T <: Real ? float(Real) : Float64 # heuristic
    DiffResults.MutableDiffResult(zero(S), (similar(x, S), ))
end

with the caller modified accordingly.

Copy link
Contributor Author

Choose a reason for hiding this comment

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

This is what I've impleented, right? Or are my eyes deceiving me?

Copy link
Collaborator

Choose a reason for hiding this comment

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

I think @tpapp's point is that the first argument of the function should be dropped if we do not use it anymore.

end

"""
Expand All @@ -25,5 +25,6 @@ constructed with [`diffresults_buffer`](@ref). Gradient is not copied as caller
vector.
"""
function _diffresults_extract(diffresult::DiffResults.DiffResult)
# NOTE: Is this still needed?
Copy link
Collaborator

Choose a reason for hiding this comment

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

That's not affected by the change above, is it?

Copy link
Owner

Choose a reason for hiding this comment

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

I also don't understand why it would not be needed.

Copy link
Contributor Author

Choose a reason for hiding this comment

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

Nah, it was a note by me so we I could ask why it's there:) As in, can we just remove it?

Copy link
Owner

Choose a reason for hiding this comment

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

What would you use instead?

The purpose was to collect everything related to extraction in a single function. I am open to suggestions, but simply removing it would break the package (the function is used).

Copy link
Contributor Author

Choose a reason for hiding this comment

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

Ah, I thought the ::Real was the reason for it, and was thinking maybe this is now redundant because type-inference is better / DiffResults might have fixed whatever issue required the ::Real in the first place. If that's not the case and it's just for convenience, then I'll remove the comment:)

Copy link
Owner

Choose a reason for hiding this comment

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

I no longer recall the case for ::Real, so I guess you can remove it. We can always put it back if things break and then document why it is there.

Please also remove the comment, and then I am happy to merge.

Copy link
Contributor Author

Choose a reason for hiding this comment

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

Done:)

DiffResults.value(diffresult)::Real, DiffResults.gradient(diffresult)
end