  We show that orbit-superposition dynamical models (Schwarzschild's method)
provide reliable estimates of nuclear black hole masses and errors when
constructed from adequate orbit libraries and kinematic data. We thus rebut two
recent papers that argue that BH masses obtained from this method are
unreliable. These papers claim to demonstrate that the range of allowable BH
masses derived from a given dataset is artificially too narrow as a result of
an inadequate number of orbits in the library used to construct dynamical
models. This is an elementary error that is easily avoided. We describe a
method to estimate the number and nature of orbits needed for the library. We
provide an example that shows that this prescription is adequate, in the sense
that the range of allowable BH masses is not artificially narrowed by use of
too few orbits. A second point raised by critics is that kinematic data are
generally obtained with insufficient spatial resolution to obtain a reliable
mass. We make the distinction between unreliable determinations and imprecise
ones. We show that there are several different properties of a kinematic
dataset that can lead to imprecise BH determinations, but none of the
attributes we have investigated leads to an unreliable determination. In short,
the degree to which the BH radius of influence is resolved by spectroscopic
observations is already reflected in the BH-mass error envelope, and is not a
hidden source of error. The BH masses published by our group and the Leiden
group are reliable.
