floating point overflow

[awzq]

int i = 124; int j = 0; "c" и "y" это float с диапазоном от 1 до 2; 
x это long double от 0.01 до 100; 
tempsum - long double; cr - long double;
--------------------------------------------------------------------------
for(; i>0; i--, j++)
{   cr = *c - (stp1.mass[j].y);
            correctmul(&cr, &csp, stp1);
    tempsum += (cr * stp1.mass[j].x * 1000);} //<--- overflow

в указанной строке происходит переполнение
мненужно только 5 знаков после запятой в самом крайнем случае.
с этой ошибкой я довольно часто сталкиваюсь так что мне нужно както решить ее основательно. раньше я просто увеличивал точность(double-> long double).

я разложил цикл получилось что переполнение вызывает
вот это выражение
stp1.mass[j].x * 1000;

[Stilet]

x это long double от 0.01 до 100;

Тут long не нужен.
Покажи весь код, если он небольшой можно его оттрассировать и понять в чем причина.

[awzq]

весь код очень большой это один case
http://www.cyberforum.ru/cpp-beginne...ead646687.html

tempsum = 78.513620529019278;
cr = -0.059998035430908203; - это 0.06 представленная в long double
stp1.mass[j].x = 0.039999999105930328 - это 0.04 представленная в long double
j = 123;
i = 1;

я разложил цикл получилось что переполнение вызывает
вот это выражение
stp1.mass[j].x * 1000;

как вообще бороться с floating point overflow на с++