Hermitian() wrapper missing a GPU compatible iterator #586

Activity

1. maleadt

   added

   good first issue

   on Feb 10, 2025
2. cpmores commented on Apr 6, 2025

   

   cpmores

   on Apr 6, 2025

   More actions

   Seems that the optimized function 'any' cannot resolve Hermitian{T,S<:AbstractMatrix{<:T}} <: AbstractMatrix{T}.

   'any' in mapreduce.jl looks like:

   Base.any(f::Function, A::AnyGPUArray) = mapreduce(f, |, A)

   But Hermitian or Symmetric is a struct including UP or Down, which cannot go into this optimized function.

   If you use the data instead of the whole struct, any will work properly.

   using CUDA
   using LinearAlgebra
   
   A = CuArray(ones(Float32, 5, 5))
   
   # This works
   S = Symmetric(A)
   println(any(isnan, S.data))
   
   # This works
   H = Hermitian(A)
   println(any(isnan, H.data))
   
   
   B = CuArray(fill(NaN, 5, 5))
   # This works
   SB = Symmetric(B)
   println(any(isnan, SB.data))
   
   # This works
   HB = Hermitian(B)
   println(any(isnan, HB.data))
3. maleadt

   changed the issue type from

   Enhancement

   to

   Bug

   on Oct 19, 2026
4. maleadt

   added

   arraysThe array abstraction: broadcasting, mapreduce, indexing, views, generic linear algebra.The array abstraction: broadcasting, mapreduce, indexing, views, generic linear algebra.

   and removed

   good first issue

   on Oct 19, 2026
5. maleadt commented on Oct 19, 2026

   

   maleadt

   on Oct 19, 2026

   Member

   More actions

   Works on current GPUArrays and CUDA.jl, with correct results for both uplos:

   A = CuArray(ones(Float32, 5, 5))
   any(isnan, Hermitian(A))  # false
6. maleadt

   closed this as completedon Oct 19, 2026