def lu( self ):
        n = self.rows
        L, U = eye(n), Mat(self[:])
        for i in range(n):
            for j in range(i+1,U.rows):
                assert U[i][i] != 0.0, 'LU requires non-zero elements on the diagonal'
                L[j][i] = m = 1.0 * U[j][i] / U[i][i]
                U[j] -= U[i] * m
        assert NPOST or isinstance(L,LowerTri) and isinstance(U,UpperTri) and L*U==self
        return L, U
sets the user's default organization
dataset/ETHPy150Open sunlightlabs/clearspending/completeness/statlib/matfunc.py/Square.lu--crossover with--dataset/ETHPy150Open buildingenergy/buildingenergy-platform/seed/views/accounts.py/set_default_organization