Hi everyone, I am quite confused about the resulting types in below snippet: N = 5 T = Float32 # edit, forgot to mention these variables μ = rand(T, N_states, 1) σ = rand(T, N_states, 1) E₀ = fill(T(0.5), N, 1) # 5×1 Matrix{Float32} U₀_1 = @. ((1 - μ) * σ) * E₀ …

Surprising type result with @

artemsolod March 23, 2024, 9:32am 2

I can’t reproduce on julia 1.10.2 but I guess this is because @. applies to fill as well. Can you try U₀_2 = fill(T(0.5), N, 1) .* (@. (1 - μ) * σ)?

2 Likes

Topic		Replies	Views
Broadcast and .* not the same General Usage	2	967	December 23, 2017
```@.``` throws unexpected error General Usage broadcasting , convert	5	531	May 19, 2020
Broadcast type mismatch General Usage question , broadcast , type	3	549	May 10, 2021
Type Instability with MMatrix General Usage question , package	5	895	November 6, 2017
The Broadcasting API and small Unions General Usage	7	612	May 15, 2018

Surprising type result with @

Related topics