#include <cmath>
using namespace std;
int main()
{
    cout<<fmod(18.L,0.04L) << endl;     //Ausgabe: 3.659e-19
    cout<<fmod(18.L,(long double)0.04); //Ausgabe: 0.04
}

Die erste Ausgabe ist klar. Aber warum gibt er mir beim 2. etwas anderes aus?
Codepad sagt das gleiche: http://codepad.org/1GE94oYr
Kommt das durch Rundungsfehler beim Umwandeln von "double" in "long double" zustande? Und sollte der Compiler das "double" literal durch den cast nicht zur Compilezeit als "long double" literal auffassen?

Jetzt sollte klar sein, was passiert: 0.04L ist ein long double, der so nah an 0.04 ist wie möglich. static_cast<long double>(0.04) ist ein long double, aber der Wert ist der von dem double-Literal 0.04, welches so nahe wie möglich an 0.04 ist, wie ein double dies erlaubt. Der long double ist zwar potentiell präziser, aber der Cast kann nicht wissen, was hier eigentlich gewollt ist und sieht nur den gerundeten double-Wert. Und dieser wird dann dem long double zugewiesen.

Irgendwie ein bisschen wirr. Ist schwer zu beschreiben. Nochmal mit 1/3:

float foo = 1./3.;  // 0.333333
double bar = 1./3.; // 0.3333333333333
bar = foo;          // 0.333333

Klar geworden?

x % m = y; // wobei 0 <= y < m

und nicht y <= m, zumal bei der ersten rechnung (approximativ) 0 rauskommt (klar weil 0.8*450 = 18, da ist kein rest drinne)

hier wundern mich mehrere dinge....

wie kann denn dabei als ergebnis 0.04 rauskommen?

Die double-Repräsentation von 0.04 in IEEE754 ist 4.00000000000000008326672684689E-2, also ein kleines bisschen größer als 0.04. Das Ergebnis passt schon.

nein, wieso denn?
das modulo hält den werte bereich doch zwischen 0 und eben hier 0.04. da ist es physikalisch doch gar nicht möglich, dass die 0.04 überschritten werden

4.00000000000000008326672684689E-2 > 0.04

Und wo die Physik hier her kommt, versteh ich auch nicht.