sum = 0;
for (j = 0; j < inIncr; ++j)
{
- if (j==0) dirx= static_cast<double>(*input);
- if (j==1) diry= static_cast<double>(*input);
- if (j==2) dirz= static_cast<double>(*input);
- tmp = static_cast<double>(*input);
+ if (j==0) dirx= static_cast<T>(*input);
+ if (j==1) diry= static_cast<T>(*input);
+ if (j==2) dirz= static_cast<T>(*input);
+ tmp = static_cast<T>(*input);
sum += (tmp * tmp);
++input;
}