I might solve it by myself.
let numerators: [Float] = ...
let denominators: [Float] = ... // actually these are integers in my case
let denominatorsEpsilon = vDSP.add(Float.leastNonzeroMagnitude * 10000000, denominators) // (Float.leastNonzeroMagnitude * 10000000 is minimal number to avoid NaN when dividing.
let divides = vDSP.divide(numerators, denominatorsEpsilon)
let alternativesOfNAN: [Float] = ...
let denominatorsClip = vDSP.clip(denominators, to: 0...1) // denominatorsClip are 0 or 1 since denominators are integers
let result = vDSP.subtract(multiplication: (divides, denominatorsClip), multiplication: (alternativesOfNAN, vDSP.add(-1, denominatorsClip)))
I don't like this code since it is not precise and includes many unnecessary operations, but it is much faster than code checking isNaN.