numpy: эффективное чтение большого массива

У меня есть двоичный файл, содержащий плотную матрицу n * m 32-битных чисел с плавающей запятой. Какой' Является ли наиболее эффективным способом чтения его в упорядоченный по Фортрану массив numpy ?

Размер файла составляет несколько гигабайт. Я могу управлять форматом, но он должен быть компактным (т.е. около 4 * n * m байтов в длину) и должен быть легко создан из кода, отличного от Python.

edit : Крайне важно, чтобы метод создавал матрицу, упорядоченную по Фортрану, напрямую (из-за размера данных я не могу позволить себе создать матрицу с упорядочением по C, а затем преобразовать ее в отдельную копию, упорядоченную по Фортрану).