Functions isHermitian() and ishpd() return a Boolean.
Function zapim() zaps small imaginary parts of components vector,
returning real if all elements are so zapped.
# NOT RUN {v <- 2^(1:30)
zapim(v+1i*exp(-v))
ishpd(matrix(c(1,0.1i,-0.1i,1),2,2)) # should be TRUEisHermitian(matrix(c(1,3i,-3i,1),2,2)) # should be TRUE# }