Shock capturing methods in 1D numerical relativity

A numerical code is presented which uses modern shock capturing methods to evolve spherically symmetric perfect fluid space-times. Harmonic slicing is used to ensure singularity avoidance, which is crucial in strong field situations. Some tests are presented, including an application to the stellar collapse problem.