Почему некорректно находит произведение чисел?

Question

[HoShik]

#include <stdio.h>
#include <iostream>

float mult(float count,...)
{
    float comp=1.0;
    for(float *ptr=&count; count>0; count--)
    {
        comp*=*(++ptr);
    }
    return comp;
}
int main()
{
    std::cout<<"Composition of the number 1.0,4.0,94.0 = "<<mult(3.0,1.0,4.0,94.0)<<std::endl;
}

Выводит 0. В чем проблема?

Comments

[vreitech]

В качестве счётчика count использовать переменную типа float? Ну такое... Когда после очередного count-- у вас count окажется равным 0.0000000000000000000002 вместо 0 (как если бы это был int), как думаете, что произойдёт?

[shurshur]

Не надо использовать конструкции типа comp*=*(++ptr);, у них могут быть сайд-эффекты в зависимости от используемого компилятора и его опций. Может умножить как на значение ptr до инкремента, так и на значение после. Я уж не говорю о том, что код от этого становится менее понятным.

До кучи, для разбора аргументов в ... принято использовать cstdarg (stdarg.h). И, кстати, в зависимости от компилятора конечная реализация ... может отличаться от массива, так что подобные хаки с указателями могут и не сработать.

Answer 1

[res2001]

Почему вы в функции с переменным количеством аргументов не используете стандартный механизм для работы с аргументами в подобных функциях (#include <cstdarg>)?

Поставьте отладочный вывод в цикл, где выведите значение count, а заодно и *ptr и comp. То же самое можно быстро посмотреть в отладчике затратив минуты 3 на все. Выводите значение floatов с большим количество знаков после запятой (10-15).
Вы будете удивлены тому, что во float 3.0 это не обязательно равно 3, ну и т.д. То же самое относится к любым типам с плавающей точкой.
Передавайте первым аргументом uint32_t, вторым и следующим float. Тогда все будет нормально. Код, конечно, придется немного переделать.

Comments

[HoShik]

По поводу библиотеки, думал с указателями будет все и так хорошо
float 3.0 я использовал просто для проверки, до этого там стоял инт и в цикле приводил что-то типа *ptr=(float *)&count
Переписал на лад starg - заработало. Спасибо за совет

[Евгений Шатунов]

HoShik, твоя функция умеет домножать на половину числа? Что будет если первым аргументом будет передано число 1.5f?
Ты в курсе что при передаче переменного количества аргументов, сами эти аргументы проходят нормализацию? Внутри mult() у тебя только один первый аргумент имеет тип float, все остальные - double.

[res2001]

HoShik, По идее и с указателями должно работать. Возможно указатель надо уменьшать, а не увеличивать. Стек растет в сторону уменьшения адресов (от старших адресов к младшим). Но на какой-то другой платформе (ОС, ЦПУ) это может быть не так, поэтому стоит использовать промежуточную библиотечную прослойку, которая будет реализована правильно для любой платформы.
Но с int и float - это вам надо вернуть как было. Это правильно. В таком виде как сейчас, возможно, будет работать в отдельно взятом примере, но не как общий случай.
Числа с плавающей запятой это всегда приближение, а в вашем случае это приближение может сказаться на том, что вы возьмете лишний аргумент (или наоборот не доберете).