Остаток от деления double C++

Нужно найти остаток от деления двух чисел типа double с с заданной точностью

long double a, b;
cin >> a >> b;
cout.precision(15);
cout << fixed << "A is: " << a << " B is: " << b << endl;
cout << fixed << fmod(a, b) << endl;

Ввод: 
123123.123123123 3.012345678
Вывод:
A is: 123123.123123123004916 B is: 3.012345678000000 
2.530571907006099

Написал вот такой код, но, видимо точности не хватает, можно ли как то вычислить остаток без fmod более точно?