d85ecc37deb225b960b3c010b892f9ca3f7b889e,pymanopt/manifolds/stiefel.py,Stiefel,proj,#Stiefel#Any#Any#,61

Before Change

    def proj(self, X, U):
        UNew = U - multiprod(
            X, multiprod(multitransp(X), U) + multiprod(multitransp(U), X)) / 2
        return UNew

    egrad2rgrad = proj

    def ehess2rhess(self, X, egrad, ehess, H):

After Change

        return np.tensordot(G, H, axes=G.ndim)

    def proj(self, X, U):
        return U - multiprod(X, multisym(multiprod(multitransp(X), U)))

    egrad2rgrad = proj

Instances