Точность вывода выше, чем двойная точность

Я печатаю некоторые данные из программы C++ для обработки/визуализации с помощью ParaView, но у меня проблема с числами с плавающей запятой. Paraview поддерживает типы данных Float32 и Float64. Float64 эквивалентен double с типичными ограничениями +/-1,7e +/- 308. Но мой код печатает числа вроде 6,5e-318. Это вызывает ошибки в ParaView при чтении данных. Я убедился, что округление этих маленьких чисел до нуля приводит к исчезновению ошибок в ParaView. Я не уверен, почему у меня такой «высокоточный» вывод, возможно, потому, что некоторые числа хранятся с более высокой точностью, чем двойная. Например, следующий код воспроизводит такое же поведение в моей системе:

#include <iostream>
int main(void)
{
  const double var1 = 1.0e-318, var2 = 1.5e-318;
  std::cout << 1.0e-318 << std::endl; 
  std::cout << var1 << std::endl; 
  std::cout << var1 - var2 << std::endl; 
  std::cout.setf(std::ios_base::fixed | std::ios_base::scientific, std::ios_base::floatfield);
  std::cout << 1.0e-318 << std::endl; 
  std::cout << var1 << std::endl; 
  std::cout << var1 - var2 << std::endl; 

  return 0;
}

Мой вывод:

9.99999e-319
9.99999e-319
-4.99999e-319
9.99999e-319
9.99999e-319
-4.99999e-319

Моя система представляет собой Mac OS X Snow Leopard, и я протестировал вышеуказанное с GCC 4.2 и GCC 4.6 с флагами -m32, -m64 и -ffloat-store (не уверен, что это полезно).

На самом деле вывод для меня в порядке, а для ParaView нет. Я просто хочу знать, почему у меня такая разница. Скорее всего, я игнорирую что-то, связанное с числами с плавающей запятой, что может быть важно. Не могли бы вы дать мне некоторое представление об этом выводе/числовом поведении для удвоения?