Simone Carlo Surace
Simone Carlo Surace
I don't believe that was the case. I was writing some intermediate results from a large calculation to a new file name, and the error occurred intermittently.
I'm trying to come up with one that is more minimal, but I can't right now. The call of `reverse` on `UpperTriangular` is definitely at play, but I don't know...
> Another idea is to use `applicable`. Something like `_tryreverse(x) = applicable(reverse, x) ? reverse(x) : reverse(collect(x))`. It will not give the most natural cotangent type in many cases, but...
> > the issue will probably occur with most matrix wrappers > > > > `reverse(Diagonal([1,2,3]))` is fine, `reverse(Symmetric([1 2; 3 4]))` is an error, again an upstream bug really....
I think that was the suggestion by the participants of the linked discourse thread, but I think the maintainers of Zygote should decide. In case it is moved to Distances.jl,...
Interesting... is the constructor function ever AD'd through? I would have thought that once the object has been constructed, the AD wouldn't be able to see the function that was...
@davidanthoff, this PR would be really handy. Do you think it could be merged before the new year (to set an arbitrary target)?
An attempt at this was started at https://github.com/JuliaGaussianProcesses/AutoGPs.jl
Is this still something worth pursuing?