JuliaApproximation / JuliaApproximation/FastTransforms.jl
Allocating lmul!
Open
Nobody has claimed this yet.
- Dominant language
- Julia
- Stars
- 281
- Forks
- 27
- PR merge metrics
- No merged PRs in 30d
Description
Hi,
This is not really an issue. I am wondering if the following is expected to allocate:
julia> F = zeros(31,61);
julia> PS = FastSphericalHarmonics.plan_sph_synthesis(copy(F))
FastTransforms plan for FFTW Fourier synthesis on the sphere for 31×61-element array of Float64
julia> @time lmul!(PS,F);
0.000356 seconds (147 allocations: 10.453 KiB)
julia> @time lmul!(PS,F);
0.000368 seconds (147 allocations: 10.453 KiB)
I am trying to optimize FastSphericalHarmonics.sph_evaluate! by caching the plans and stumbled on this allocation.
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 Julia session shown in the issue, including the repeated lmul!(PS,F) calls and allocation counts. Trace the FastSphericalHarmonics.plan_sph_synthesis and lmul! entry points to determine whether the allocations are expected; done means documenting the cause or confirming a specific allocation-related fix.
Written by the indexing model from the issue text.
Assessment
- Tech stack
- julia
- Domain
- performance
- Issue type
- Bug
- Difficulty
- 4/5
- Estimated time
- 3-5 days
- Activity status
- Stale
- Clarity
- Needs clarification
- Newbie friendliness
- 25/100