Переполнение типа с плавающей точкой

#include <stdio.h>
 
int main() {
  float sum1 = 0;
  double sum2 = 0;
 
  for(size_t i=0; i<1000000000; ++i) {
    sum1 += 123456.0;
    sum2 += 123456.0;
  }
 
  printf("sum1 = %15f\nsum2 = %15f\n", sum1, sum2);
}

#include <iostream>
#include <iomanip>
#include <cmath>
 
int main()
{
    const float pi = 3.1415926535897932384626433832795f;
 
    float a = 0;
    std::cout << std::setw( 9 ) << cos( a ) << std::endl; //  1
    a += pi;
    std::cout << std::setw( 9 ) << cos( a ) << std::endl; // -1
 
    a = pi * 1000000;
    std::cout << std::setw( 9 ) << cos( a ) << std::endl; //  0.995356
    a += pi;
    std::cout << std::setw( 9 ) << cos( a ) << std::endl; // -0.979098
 
    a = pi * 10000 * 10000;
    std::cout << std::setw( 9 ) << cos( a ) << std::endl; //  0.210237
    a += pi;
    std::cout << std::setw( 9 ) << cos( a ) << std::endl; //  0.210237
}

#include <cfenv>
#include <iostream>
 
int main() {
    double a = 1000;
    for ( int i = 0; i < 15; ++i ) {
        std::feclearexcept( FE_OVERFLOW );
        a *= a;
        if ( std::fetestexcept( FE_OVERFLOW ) ) {
            std::cout << "Overflow at i = " << i << std::endl;
            break;
        }
    }
}

#include <iostream>
#include <limits>
 
int main()
{
    double max = std::numeric_limits<double>::max();
    std::cout << max << '\n' << max*2;
}