JuliaApproximation / JuliaApproximation/FastTransforms.jl
Segmentation fault for spherical harmonic transform
Open
Nobody has claimed this yet.
- Dominant language
- Julia
- Stars
- 281
- Forks
- 27
- PR merge metrics
- No merged PRs in 30d
Description
julia> F = sphrandn(Float64, 1024, 2047); # convenience method
julia> P = plan_sph2fourier(F);
julia> PS = plan_sph_synthesis(F);
julia> PA = plan_sph_analysis(F);
julia> @time G = PS*(P*F);
[9851] signal (11.1): Segmentation fault: 11
in expression starting at REPL[7]:1
Allocations: 8155300 (Pool: 8146430; Big: 8870); GC: 15
zsh: segmentation fault
julia> versioninfo()
Julia Version 1.9.3
Commit bed2cd540a1 (2023-08-24 14:43 UTC)
Build Info:
Official https://julialang.org/ release
Platform Info:
OS: macOS (x86_64-apple-darwin22.4.0)
CPU: 16 × Intel(R) Core(TM) i9-9900K CPU @ 3.60GHz
WORD_SIZE: 64
LIBM: libopenlibm
LLVM: libLLVM-14.0.6 (ORCJIT, skylake)
Threads: 8 on 16 virtual cores
Environment:
JULIA_NUM_THREADS = 8
JULIA_PKG_DEVDIR = /Users/abradley/Dropbox/Julia/Dev
Contributor guide
No contributing guide indexed for this repository
First steps
- Read the whole issue, then the project's contributing guide.
- Comment on the issue to say you are picking it up — it saves two people doing the same work.
- Fork the repository and make your change on a branch.
- Open a pull request that references the issue number.
Research direction
Start by reproducing the provided Julia 1.9.3 example with sphrandn, plan_sph2fourier, plan_sph_synthesis, and plan_sph_analysis. Trace the spherical harmonic transform path involved in PS*(P*F) and verify that the reproduction no longer causes a segmentation fault.
Written by the indexing model from the issue text.
Assessment
- Tech stack
- julia
- Domain
- backend
- Issue type
- Bug
- Difficulty
- 4/5
- Estimated time
- 3-5 days
- Activity status
- Stale
- Clarity
- Needs clarification
- Newbie friendliness
- 30/100