I think there has been a regression on the issue reported in https://github.com/EnzymeAD/Enzyme.jl/issues/647 on differentiating Bessel functions.
I tried both forward-mode and reverse-mode AD and traced down the error to Base.Math.paynehanek, though I believe the issue also involves Base.Math.DoubleFloat64 since Enzyme could differentiate paynehanek if I commented out the part with DoubleFloat64. I don't know how to read the error messages, but I thought I would open this issue to keep track
Here are the details of my Julia environment
pkg> status Enzyme Bessels
Status `/local/home/lxvm/projects/beams/Project.toml`
[0e736298] Bessels v0.3.0-DEV `~/.julia/dev/Bessels`
[7da242da] Enzyme v0.11.4
julia> versioninfo()
Julia Version 1.9.2
Commit e4ee485e909 (2023-07-05 09:39 UTC)
Platform Info:
OS: Linux (x86_64-linux-gnu)
CPU: 12 × Intel(R) Core(TM) i7-8750H CPU @ 2.20GHz
WORD_SIZE: 64
LIBM: libopenlibm
LLVM: libLLVM-14.0.6 (ORCJIT, skylake)
Threads: 1 on 12 virtual cores
Environment:
LD_PRELOAD = /usr/lib/x86_64-linux-gnu/libstdc++.so.6
I was running Bessels.jl at commit 52b694883c84769c5e51f7fa4d0d5fff823ad098 of the master branch
Hi,
I think there has been a regression on the issue reported in https://github.com/EnzymeAD/Enzyme.jl/issues/647 on differentiating Bessel functions. I tried both forward-mode and reverse-mode AD and traced down the error to
Base.Math.paynehanek
, though I believe the issue also involvesBase.Math.DoubleFloat64
since Enzyme could differentiatepaynehanek
if I commented out the part withDoubleFloat64
. I don't know how to read the error messages, but I thought I would open this issue to keep trackHere are the details of my Julia environment
I was running Bessels.jl at commit 52b694883c84769c5e51f7fa4d0d5fff823ad098 of the master branch
Here is the MWE (taken from https://github.com/EnzymeAD/Enzyme.jl/issues/647) for both forward-mode, which errors, and reverse-mode AD, which segfaults.