C序与Fortran序区别详解

对于一个文章开发者来说，牢固扎实的基础是十分重要的，golang学习网就来带大家一点点的掌握基础知识点。今天本篇文章带大家了解《NumPy维度顺序解析：C序与Fortran序详解》，主要介绍了，希望对大家的知识积累有所帮助，快点收藏起来吧，否则需要时就找不到了！

理解NumPy多维数组的默认维度顺序 (C-Order)

在使用NumPy创建多维数组时，例如 np.ones((A, B, C))，其默认的维度解释和内存布局遵循C-order（C语言风格的行主序）。这意味着：

1. 逻辑结构: 数组被视为 A 组的 B x C 元素，而每个 B x C 元素又被视为 B 组的 C 元素。直观上，我们可以将其理解为 A 个“层”或“批次”，每个层包含 B 行和 C 列。
2. 内存布局: 在内存中，最后一个维度（即 C 维度）是变化最快的。这意味着，如果一个元素是 x[i, j, k]，那么它在内存中紧邻的下一个元素通常是 x[i, j, k+1]。当 k 达到最大值时，j 会递增，然后 k 从头开始；当 j 达到最大值时，i 会递增，然后 j 和 k 从头开始。

示例解析：np.ones((3, 2, 2))

当您创建 np.ones((3, 2, 2)) 时，NumPy会生成一个三维数组，其形状为 (3, 2, 2)。这表示：

• 3：最外层维度，可以理解为有3个独立的2x2矩阵。
• 2：中间维度，表示每个2x2矩阵有2行。
• 2：最内层维度，表示每个2x2矩阵有2列。

因此，您会得到3个大小为2x2的数组。为了更直观地理解C-order的内存布局，我们可以通过 strides 属性来查看。strides 表示访问数组中每个维度下一个元素所需的字节数。

import numpy as np

# 创建一个形状为 (2, 3, 4) 的数组，使用默认的C-order
arr_c = np.arange(2 * 3 * 4).reshape((2, 3, 4))
print("C-order 数组:\n", arr_c)
print("C-order 数组形状:", arr_c.shape)
print("C-order 数组步长 (bytes/element):\n", arr_c.strides) # 假设元素为4字节整数 (int32)

在上述示例中，如果元素为4字节整数 (int32)，arr_c.strides 的输出将是 (48, 16, 4)。这意味着：

• 要从 arr_c[i, j, k] 移动到 arr_c[i, j, k+1]，需要跳过 4 字节（对应最后一个维度，步长为 1 * sizeof(element)）。
• 要从 arr_c[i, j, k] 移动到 arr_c[i, j+1, k]，需要跳过 16 字节（对应中间维度，步长为 C * sizeof(element)，即 4 * 4）。
• 要从 arr_c[i, j, k] 移动到 arr_c[i+1, j, k]，需要跳过 48 字节（对应第一个维度，步长为 B * C * sizeof(element)，即 3 * 4 * 4）。

这清晰地表明，最后一个维度在内存中是连续的，变化最快。

理解维度与常见映射

许多深度学习框架（如PyTorch）在处理图像数据时，常用的维度顺序可能是 [Channel, Height, Width] 或 [Batch, Channel, Height, Width]。NumPy的C-order默认 (Depth, Height, Width) 或 (Batch, Height, Width, Channel) 在某些情况下可以直接匹配，但在另一些情况下可能需要进行维度转置。

例如，如果您有一个形状为 (3, 256, 256) 的图像数组（3个通道，256x256像素），NumPy的C-order会将其解释为3个256x256的层，每个层内部是行优先存储。这与 [Channel, Row, Columns] 的逻辑是吻合的。然而，如果您的数据源是 (Height, Width, Channel)，您可能需要使用 arr.transpose((2, 0, 1)) 将其转换为 (Channel, Height, Width) 以适应某些模型输入。

Fortran序：另一种内存布局 (F-Order)

NumPy还支持Fortran-order（列主序）的内存布局。您可以通过在创建数组时指定 order='F' 来使用它。

# 创建一个形状为 (2, 3, 4) 的数组，使用Fortran-order
arr_f = np.arange(2 * 3 * 4).reshape((2, 3, 4), order='F')
print("\nFortran-order 数组:\n", arr_f)
print("Fortran-order 数组形状:", arr_f.shape)
print("Fortran-order 数组步长 (bytes/element):\n", arr_f.strides) # 假设元素为4字节整数 (int32)

对于Fortran-order的 (A, B, C) 数组，其内存布局与C-order相反：

• 内存布局: 第一个维度（即 A 维度）是变化最快的。这意味着，如果一个元素是 x[i, j, k]，那么它在内存中紧邻的下一个元素通常是 x[i+1, j, k]。

在上述Fortran-order示例中，如果元素为4字节整数 (int32)，arr_f.strides 的输出将是 (4, 8, 24)。这意味着：

• 要从 arr_f[i, j, k] 移动到 arr_f[i+1, j, k]，需要跳过 4 字节（对应第一个维度，步长为 1 * sizeof(element)）。
• 要从 arr_f[i, j, k] 移动到 arr_f[i, j+1, k]，需要跳过 8 字节（对应中间维度，步长为 `A

今天关于《C序与Fortran序区别详解》的内容就介绍到这里了，是不是学起来一目了然！想要了解更多关于的内容请关注golang学习网公众号！