  We have been constructed a brand-new radiation hydrodynamics solver based
upon Smoothed Particle Hydrodynamics (SPH), which works on parallel computer
system. The code is designed to investigate the formation and evolution of the
first generation objects at $z \gtrsim 10$, where the radiative feedback from
various sources play important roles. The code can compute the fraction of
chemical species e, H$^+$, H,
  H$^-$, H$_2$, and H$_2^+$ by fully implicit time integration. It also can
deal with multiple sources of ionizing radiation, as well as the radiation at
Lyman-Werner band. We compare the results for a few test calculations with the
results of one dimensional simulations, in which we find good agreements with
each other.
  We also evaluate the speedup by parallelization, that is found to be almost
ideal, as far as the number of sources is comparable to the number of
processors.
