如何用MathNet.Numerics查找矩阵中符合条件的最小值索引？

你遇到的需求很典型——在Matrix<float>里只针对值≥0的元素，找出所有最小值对应的索引。比如你给出的示例矩阵：

var mat = Matrix<float>.Build.DenseOfArray(new float[,] { {2, 3, 1}, {3, -1, 5}, {1, 4, 3} });

你已经能通过这段代码拿到符合条件的最小值：

var min = mat.Enumerate().Where(x => x >= 0).Min();

但需要定位到两个值为1的元素的索引（[0,2]和[2,0]），而不是被排除的-1的位置。

可行解决方案（附优化建议） #

你找到的方案其实已经很实用了，利用EnumerateIndexed()方法可以直接获取每个元素的行索引、列索引和值，结合筛选就能拿到所有目标索引：

var min = mat.Enumerate().Where(x => x >= 0).Min();
var allMinIndices = mat.EnumerateIndexed()
    .Where(x => x.Item3 >= 0 && x.Item3 == min) // 合并条件让逻辑更闭环
    .Select(x => (Row: x.Item1, Column: x.Item2)) // 转换成更直观的行列结构（可选）
    .ToList();

细节说明：

• EnumerateIndexed()返回的元组包含三个元素：(行索引, 列索引, 元素值)，所以x.Item1对应行号，x.Item2对应列号，x.Item3是元素值。
• 把x >= 0的条件加入筛选更严谨——虽然前面已经通过Where(x => x >=0).Min()确保了min是符合条件的，但合并后逻辑更完整，避免极端情况出现偏差。
• 如果需要更友好的索引格式，可以把结果转换成List<Tuple<int, int>>或者自定义的索引类，方便后续业务使用。

性能优化思路 #

如果是处理超大规模矩阵，你可以考虑合并成一次遍历的写法，减少遍历次数节省开销：

float currentMin = float.MaxValue;
List<(int Row, int Column)> indices = new List<(int Row, int Column)>();

for (int i = 0; i < mat.RowCount; i++)
{
    for (int j = 0; j < mat.ColumnCount; j++)
    {
        float val = mat[i, j];
        if (val >= 0)
        {
            if (val < currentMin)
            {
                currentMin = val;
                indices.Clear();
                indices.Add((i, j));
            }
            else if (val == currentMin)
            {
                indices.Add((i, j));
            }
        }
    }
}

这种写法在超大矩阵场景下内存和时间效率更高，不过对于常规业务场景，之前的LINQ写法已经足够简洁易读，开发效率更高。

内容的提问来源于stack exchange，提问作者MrEighteen