Add helper function ZeroSqrt(DoubleVector) to study intermediate tachyons