Floating-point Modulus math.fmod

[Fendercaster via Digitalmars-d-learn]

I'm not quite sure if this is the right forum to ask this 
question:
I've been trying to implement the "floating-point modulus" 
function from the math library. Equivalently that's what I've 
tried in Python too. Problem is - the results are different and 
not even close. Please have a look at both snippets. Maybe 
someone can help me out:

[code = D]
import std.math;
import std.stdio;
immutable real MODULUS = 1e8;

real floatModPow(real base, ulong exp, real mod)
{
    real r = 1.0;

    while (exp) {
        if (exp & 1) r = fmod(r * base, mod);
        exp >>= 1;
        base = fmod(base*base, mod);
    }

    return r;
}

void main()
{
    floatModPow(3.866, 987654321UL, MODULUS).writeln;
}
[/code]

In this case the result is 6.44588e+07.

[code = Python]
import math

MOD = 1e8

def fpmod(b, e, m):
    r = 1.0
    while e:
        if (e & 1):
            r = math.fmod(r * b, m)
        e >>= 1
        b = math.fmod(b * b, m)
    return r

print (fpmod(3.866, 987654321, MOD))
[/code]

In this case the result is 82031250.0.

AFAIK the D.fmod and Python.fmod functions serve the same 
purpose. What am I missing?