find package:hmatrix

Find index of elements which satisfy a predicate
>>> find (>0) (ident 3 :: Matrix Double)
[(0,0),(1,1),(2,2)]