Multi-Band Gutzwiller Wave Functions for Itinerant Ferromagnetism