  Cryogenic macro-calorimeters instrumented with NTD thermistors have been
developed for several decades. The choice of the optimal bias current is
crucial for a proper operation of these detectors, both in terms of energy
resolution and stability. In this paper we present a set of automatic
measurements and analysis procedures for the characterization and optimization
of the working configuration of the NTD thermistors. The presented procedures
were developed for CUORE, an array of 988 cryogenic macro-calorimeters
instrumented with NTD thermistors that has been taking data since 2017. These
procedures made it possible to characterize a large number of detectors in a
reliable way. They are suitable enough to be used also in other large arrays of
cryogenic detectors, such as CUPID.
