Direct computation of multivalued phase space solutions for Hamilton-Jacobi equations