e99414f3ea4021ce2816b44b7db75f42bbb34b41,starry/ops/ops.py,OpsReflected,X,#OpsReflected#Any#Any#Any#Any#Any#Any#Any#Any#Any#Any#,460

Before Change

            rTA1Rz = tt.dot(rTA1Rz, A1InvFA1)
        
        // Rotate to the correct phase
        X_rot = tt.set_subtensor(
            tt.zeros((rows, cols))[i_rot], 
            self.dotR(rTA1Rz, inc, obl, theta[i_rot])
        )

        return X_rot

    @autocompile(

After Change

            tt.zeros_like(norm),
            source[:, 0] * norm
        )
        theta_z = tt.arctan2(source[:, 0], source[:, 1])
        rTA1Rz = self.dotRz(rTA1, theta_z)

        // Apply limb darkening?
        if self.filter:
            F = self.F(u, f)
            A1InvFA1 = ts.dot(ts.dot(self.A1Inv, F), self.A1)
            rTA1Rz = tt.dot(rTA1Rz, A1InvFA1)
        
        // Rotate to the correct phase
        X_rot = self.dotR(rTA1Rz, inc, obl, theta)

        // TODO: Implement occultations in reflected light
        // Throw error if there"s an occultation
        X_occ = RaiseValuerErrorIfOp(

Instances