Finding Unbounded Ray for LP

I am trying to implement the L-Shaped method for two stage stochastic programs with recourse. A sub-problem that arises is to find the unbounded ray if an LP is unbounded. Is there a smart way to detect unboundedness and return the unbounded ray? Thanks!