  We present results from a numerical study of spherical gravitational collapse
in shift symmetric Einstein dilaton Gauss-Bonnet (EdGB) gravity. This modified
gravity theory has a single coupling parameter that when zero reduces to
general relativity (GR) minimally coupled to a massless scalar field. We first
show results from the weak EdGB coupling limit, where we obtain solutions that
smoothly approach those of the Einstein-Klein-Gordon system of GR. Here, in the
strong field case, though our code does not utilize horizon penetrating
coordinates, we nevertheless find tentative evidence that approaching black
hole formation the EdGB modifications cause the growth of scalar field "hair",
consistent with known static black hole solutions in EdGB gravity. For the
strong EdGB coupling regime, in a companion paper we first showed results that
even in the weak field (i.e. far from black hole formation), the EdGB equations
are of mixed type: evolution of the initially hyperbolic system of partial
differential equations lead to formation of a region where their character
changes to elliptic. Here, we present more details about this regime. In
particular, we show that an effective energy density based on the Misner-Sharp
mass is negative near these elliptic regions, and similarly the null
convergence condition is violated then.
