求矩阵最大值的回调函数

在编程中，特别是涉及到数值计算和矩阵操作时，常常需要寻找矩阵中的最大值。许多编程语言都提供了库函数或自定义回调函数的方式来实现这个功能。回调函数是一种在特定上下文中执行的函数，通常用于处理某种特定事件或计算过程的结果。

对于矩阵的最大值，回调函数可能看起来像这样（假设使用 Python 和 NumPy 库为例）：

import numpy as np

def find_max(matrix, callback):
    max_val = matrix[0]  # 假设初始化为第一个元素
    max_index = (0, 0)  # 存储最大值的位置

    def _callback(element):  # 内部回调函数
        nonlocal max_val, max_index
        if element > max_val:
            max_val = element
            max_index = np.unravel_index(element.argmax(), matrix.shape)

    np.apply_along_axis(_callback, axis=0, arr=matrix)  # 对每一行应用回调
    np.apply_along_axis(_callback, axis=1, arr=matrix.T)  # 对每一列应用回调（如果矩阵是二维）

    return max_val, max_index

# 使用示例
matrix = np.array([[1, 2, 3], [4, 5, 6], [7, 8, 9]])
max_value, max_index = find_max(matrix, _callback)
print("最大值:", max_value)
print("最大值位置:", max_index)

在这个例子中，`_callback` 是回调函数，它接收矩阵中的每个元素并检查是否大于当前的最大值。如果是，就更新最大值和其索引。`apply_along_axis` 函数会将这个回调函数应用于矩阵的每一维，找到全局最大值。

相关问题：
1. 在其他编程语言中（如 MATLAB 或者 C++），查找矩阵最大值的回调函数实现有何不同？
2. 如果矩阵非常大，如何优化上述代码以提高查找速度？
3. 如何修改回调函数，使其返回最大值的索引而不是仅返回值？