Funkcja modf - źle działające liczby 1.001 i 2.002

#include <iostream>
#include <iomanip>
#include <cmath>
using namespace std;

int main()

{
    double i = 0, precision, liczba, czescUlamkowa, czescCalkowita;
   
    cout << "precision=";
    cin >> precision;
   
    while( czescUlamkowa != 0 )
    {
        liczba = precision * pow( 10, i );
        czescUlamkowa = modf( liczba, & czescCalkowita );
        i++;
    }
    cout << "Precyzja: " <<( i - 1 ) << " cyfr po przecinku." << endl;
   
    system( "pause" );
    return 0;
}