矩阵求偏迹（约化密度矩阵）python代码

在量子信息和量子计算领域，矩阵的偏迹运算是一项常用操作。不过，在实际编程中，用于量子计算的 Python 库常常与其他 Python 库存在兼容性问题。例如，笔者在使用 PyTorch 和 Cvxpy 中的优化器时，就发现它们与一些量子计算库不太兼容。因此，笔者特意将最常用的偏迹运算代码分享出来。

这段代码的逻辑较为简单，主要是通过对矩阵进行分块，然后不断递归实现的。但目前该函数仅支持 $2^n\times 2^{}$