This mailing list thread discusses reverting the change in 1.21.0 to make unique return a single nan in its output when the input array contains multiple nans. This issue is to track that, needs a decision before 1.22.0rc1 I think.