Double-differential cross sections of neutron production at angles from 0 to 110 degrees from many reactions induced by light and medium nuclei on targets from 12-C to 208-Pb, at several incident energies from 95 to 600 MeV/A have been measured recently at the Institute of Physical and Chemical Research (RIKEN) Ring Cyclotron in Japan and at the Heavy-Ion Medical Accelerator of the National Institute of Radiological Science in Chiba, Japan using the time-of-flight technique. We have analyzed all these new measurements using the Quantum Molecular Dynamics (QMD) model, the Oak Ridge intranuclear cascade model HIC, the ISABEL intranuclear cascade model from LAHET, and the Los Alamos version of the Quark-Gluon String Model code LAQGSM03. On the whole, all four models used here describe reasonably well most of the measured neutron spectra, although different models agree differently with data from specific reactions and some serious discrepances are observed for some reactions. We present here some illustrative results from our study, discuss possible reasons for some of the observed discrepancies and try to outline ways to further improve the tested codes in order to address these problems.