4f75512036063f707ac077a22bc028c81edfef27,starry/kepler.py,System,draw,#System#,1123

Before Change

            raise ValueError("Please call `solve()` first.")

        // Number of coefficients
        bodies = [
            body
            for body in [self._primary] + list(self._secondaries)
            if body.map.ydeg > 0
        ]
        N = np.sum([body.map.Ny - 1 for body in bodies])

        // Fast multivariate sampling using the Cholesky factorization
        u = math.cast(np.random.randn(N))

After Change

        // Fast multivariate sampling using the Cholesky factorization
        yhat, cho_ycov = self._solution
        u = math.cast(np.random.randn(N))
        x = yhat + math.dot(cho_ycov, u)

        // Set all the map vectors
        n = 0
        for body in self._solved_bodies:
            inds = slice(n, n + body.map.Ny)
            body.map.amp = x[inds][0]
            body.map[1:, :] = x[inds][1:] / body.map.amp
            n += body.map.Ny

    def lnlike(self, *, design_matrix=None, t=None, woodbury=True):

Instances