Users want to get rounding mode in `quantileExact` equivalent to "median_low" from Python:
https://docs.python.org/3/library/statistics.html#statistics.median_low
while current implementation works equivalent to "median_high" (?)
DataTypePtr getReturnType() const override
    {
        DataTypePtr res;
        if constexpr (returns_float)
            res = std::make_shared<DataTypeNumber<FloatReturnType>>();
        else
            res = argument_type;
        if constexpr (returns_many)
            return std::make_shared<DataTypeArray>(res);
        else
            return res;
    }
    void add(AggregateDataPtr place, const IColumn ** columns, size_t row_num, Arena *) const override
    {
        auto value = static_cast<const ColVecType &>(*columns[0]).getData()[row_num];
        if constexpr (std::is_same_v<Data, QuantileTiming<Value>>)
        {
            
            if (isNaN(value) || value > std::numeric_limits<Value>::max() || value < std::numeric_limits<Value>::min())
                return;
        }
        if constexpr (has_second_arg)
            this->data(place).add(
                value,
                columns[1]->getUInt(row_num));
        else
            this->data(place).add(value);
    }
    void merge(AggregateDataPtr place, ConstAggregateDataPtr rhs, Arena *) const override