如何将numpy linalg求解器用于多维用途？

Question

我正在尝试为超定系统确定解决方案，我正在运行代码，但我认为它看起来很丑陋，因为当我觉得它应该更简单时，我必须做许多for循环，当前看起来像这样：

for s in range(4):
    for m in range(p):
        storing_list = []
        for n in range(tree_length):
            mu.append(tree[n][s][m])
        Q, R = np.linalg.qr(storing_list)
        Qb = np.dot(Q.T, mu)
        x_qr = np.linalg.solve(R, Qb)

由于我是从文件生成数据的，所以我无法真正发布所有代码。如您所见，我必须从树中获取值，将其放入列表中，然后使用该值来求解系统，是否有更干净的方法进行此计算？预先感谢