gaussian_process.py 文件源码

def solve_G(self, b):
    """Solve ``Gx=b`` where ``G`` is the Gram matrix of the GP.

    Uses the internally-stored LU decomposition of ``G`` computed in
    ``gp.update()``."""

    assert self.ready
    return linalg.lu_solve((self.LU, self.LU_piv), b, check_finite=True)