Na programação de cálculos de matrizes densas, existe algum motivo para escolher um layout de linha principal do layout maior de colunas?

Eu sei que, dependendo do layout da matriz escolhida, precisamos escrever o código apropriado para usar as memórias de cache efetivamente para fins de velocidade.

O layout da linha principal parece mais natural e mais simples (pelo menos para mim). Mas grandes bibliotecas como LAPACK, que são escritas em Fortran, usam o layout principal da coluna, então deve haver algum motivo para ter feito essa escolha.