-
Notifications
You must be signed in to change notification settings - Fork 63
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
Jacobian returns one array #351
Conversation
@YingboMa what's the most useful/common output of the jacobian function in the ecosystem? |
Ok now we are in agreement with ForwardDiff.jl |
@@ -780,6 +787,7 @@ end | |||
adjoint(Duplicated(x, dx), tape) | |||
return dx | |||
end | |||
mapreduce(LinearAlgebra.adjoint, vcat, rows) |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Kind of hacky, but would it be possible to create the zero init'd matrix, then unsafe_array it at the corresponding offset for each index to avoid making a bunch of temporaries
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I will leave that for you as an optimization.
No description provided.