Andrea Pagnani

Simple gpu loop on CUDA does not return

1

On julia 1.7.2 creating a new environment with only the included packages (see below) ``` using Tullio, CUDA, LoopVectorization, CUDAKernels, KernelAbstractions function gpr(N, L) Jseq = rand(Float32, N + 2,...

Possibly unreasonable storage size of Array{Tuple{Int64,Int64,Float64},1}:

12

Dearests, I want to store an array of type `x::Array{Tuple{Int64,Int64,Float64},1}`. Hereafter is the result ``` julia> x 10000-element Array{Tuple{Int64,Int64,Float64},1} ``` now I `@save "x.jld" x` and the size of the...

Should 1^x == 1?

1

Following the [discussion](https://discourse.julialang.org/t/should-1-x-1/77529/5?u=andrea_pagnani) on Discourse on the same issue: do we want to enforce a rule to make 1^x -> 1? There is a concern for x = ±∞ for...

Bad interaction with Enzyme?

6

Hi there! Cross linking here the original [issue](https://github.com/EnzymeAD/Enzyme.jl/issues/1279) on Enzyme following the advice of @wsmoses. Repeating here: I experience a problem with Enzyme interacting with Tullio. The MWE is the...

problem with tullio

5

Hi there, I experience a problem with Enzyme interacting with Tullio. The MWE is the computation of the trace of the product with two matrices $tr(A*B)$. I wrote a method...