Example #1
0
    def init_parameters_from_cfg(self):
        """
        Sets initial parameters in the ParameterSet from the configuration.
        """
        prefix = 'initval_'
        for name in dir(self.cfg):
            if name.startswith(prefix):
                varname = name[len(prefix):]
                value = getattr(self.cfg, name)
                print "Parameter initialization: %015s = %s" % (varname,
                                                                repr(value))
                if isinstance(value, (float, int)):
                    self.ps[varname] = value
                else:
                    self.ps[varname] = post(value)

        prefix = 'initvar_'
        for name in dir(self.cfg):
            if name.startswith(prefix):
                varname = name[len(prefix):]
                variance = getattr(self.cfg, name)
                meanname = 'initmean_' + varname
                mean = getattr(self.cfg, meanname) if meanname in dir(
                    self.cfg) else 0.0
                print "Random parameter initialization: %015s ~ N(%.3f, %.3f)" % (
                    varname, mean, variance)
                self.ps[varname] = post(
                    np.random.normal(mean,
                                     variance,
                                     size=self.ps[varname].shape))
Example #2
0
    def init_parameters(self, var=None, seed=None):
        """Initializes the parameteres of the ParameterSet using a zero-mean normal distribution.
        :param var: Variance of the normal distribution.
                    If None, then the 'initialization_variance' parameter from the configuration is used.
        :param seed: random seed for initialization
        """
        if var is None:
            if self.cfg.has('initialization_variance'):
                var = self.cfg.initialization_variance
            elif self.cfg.has('initvar' ):
                var = self.cfg.initvar
            else:
                var = 0.01
        if seed is None:
            if self.cfg.has('initialization_seed'):
                seed = self.cfg.initialization_seed
            else:
                seed = 1
        print "Random parameter initialization: ALL ~ N(0, %.3f) with seed %d" % (var, seed)
        np.random.seed(seed)
        self.ps.data[:] = post(np.random.normal(0, var, size=self.ps.data.shape))

        if self.cfg.get('positive_weights_init'):
            print "Ensuring positive initialization weights."
            self.ps.data[:] = xp.abs(self.ps.data)

        self.init_parameters_from_cfg()
        self.ps.restore_constants()
Example #3
0
    def init_parameters(self, var=None, seed=None):
        """Initializes the parameteres of the ParameterSet using a zero-mean normal distribution.
        :param var: Variance of the normal distribution.
                    If None, then the 'initialization_variance' parameter from the configuration is used.
        :param seed: random seed for initialization
        """
        if var is None:
            if self.cfg.has('initialization_variance'):
                var = self.cfg.initialization_variance
            elif self.cfg.has('initvar'):
                var = self.cfg.initvar
            else:
                var = 0.01
        if seed is None:
            if self.cfg.has('initialization_seed'):
                seed = self.cfg.initialization_seed
            else:
                seed = 1
        print "Random parameter initialization: ALL ~ N(0, %.3f) with seed %d" % (
            var, seed)
        np.random.seed(seed)
        self.ps.data[:] = post(
            np.random.normal(0, var, size=self.ps.data.shape))

        if self.cfg.get('positive_weights_init'):
            print "Ensuring positive initialization weights."
            self.ps.data[:] = xp.abs(self.ps.data)

        self.init_parameters_from_cfg()
        self.ps.restore_constants()
Example #4
0
    def constants(self, value):
        self._constants = dict(value)

        # build constant selector and zero array for nullifying gradient
        idxs = [np.arange(*self.extents_of_var(var)) for var in self._constants.iterkeys()]
        sel = reduce(lambda a, b: np.concatenate((a, b)), idxs, [])
        sel = np.sort(sel)
        self._constants_selector = post(sel)
        self._constants_zeros = xp.zeros((len(sel),))
Example #5
0
    def constants(self, value):
        self._constants = dict(value)

        # build constant selector and zero array for nullifying gradient
        idxs = [
            np.arange(*self.extents_of_var(var))
            for var in self._constants.iterkeys()
        ]
        sel = reduce(lambda a, b: np.concatenate((a, b)), idxs, [])
        sel = np.sort(sel)
        self._constants_selector = post(sel)
        self._constants_zeros = xp.zeros((len(sel), ))
Example #6
0
 def set_constant_parameters_from_cfg(self):
     """
     Sets constant parameters in the ParameterSet from the configuration.
     """
     prefix = 'const_'
     constants = self.ps.constants
     for name in dir(self.cfg):
         if name.startswith(prefix):
             varname = name[len(prefix):]
             value = getattr(self.cfg, name)
             print "Constant parameter: %015s = %s" % (varname, repr(value))
             constants[varname] = post(value)
     self.ps.constants = constants
Example #7
0
 def set_constant_parameters_from_cfg(self):
     """
     Sets constant parameters in the ParameterSet from the configuration.
     """
     prefix = 'const_'
     constants = self.ps.constants
     for name in dir(self.cfg):
         if name.startswith(prefix):
             varname = name[len(prefix):]
             value = getattr(self.cfg, name)
             print "Constant parameter: %015s = %s" % (varname, repr(value))
             constants[varname] = post(value)
     self.ps.constants = constants
Example #8
0
        def __init__(self, ds, idx, minibatch_size, pad_data, force_cpu):
            self._keys = [k for k in ds.keys() if not k.startswith('meta_')]
            self._minibatch_size = minibatch_size
            self.n_samples = len(idx)
            self.n_bytes = 0

            for key in self._keys:
                data = ds[key][idx, ...]
                if pad_data:
                    data = self._pad(data, minibatch_size)
                self.n_bytes += data.nbytes
                if not force_cpu:
                    data = post(data)
                setattr(self, key, data)
Example #9
0
    def init_parameters_from_cfg(self):
        """
        Sets initial parameters in the ParameterSet from the configuration.
        """
        prefix = 'initval_'
        for name in dir(self.cfg):
            if name.startswith(prefix):
                varname = name[len(prefix):]
                value = getattr(self.cfg, name)
                print "Parameter initialization: %015s = %s" % (varname, repr(value))
                if isinstance(value, (float, int)):
                    self.ps[varname] = value
                else:
                    self.ps[varname] = post(value)

        prefix = 'initvar_'
        for name in dir(self.cfg):
            if name.startswith(prefix):
                varname = name[len(prefix):]
                variance = getattr(self.cfg, name)
                meanname = 'initmean_' + varname
                mean = getattr(self.cfg, meanname) if meanname in dir(self.cfg) else 0.0
                print "Random parameter initialization: %015s ~ N(%.3f, %.3f)" % (varname, mean, variance)
                self.ps[varname] = post(np.random.normal(mean, variance, size=self.ps[varname].shape))
Example #10
0
 def load(cls, state_dir):
     """
     Loads parameters history and best results from the state directory.
     :param state_dir: directory to load from
     :return: a ParameterHistory object loaded from the state directory
     """
     history_filename, results_filename, best_pars_filename = cls._get_result_filenames(state_dir)
     self = cls()
     self.history = np.load(history_filename)['history']
     self.best_pars = post(np.load(best_pars_filename)['best_pars'])
     with open(results_filename, 'rb') as results_file:
         data = json.load(results_file)
         self.best_iter = data['best_iter']
         self.best_val_loss = data['best_val_loss']
         self.best_tst_loss = data['best_tst_loss']
         self.start_time = data['start_time']
         self.end_time = data['end_time']
         self.termination_reason = data['termination_reason']
         if 'cfg' in data:
             self.cfg = data['cfg']
         else:
             self.cfg = {}
     return self
Example #11
0
 def load(cls, state_dir):
     """
     Loads parameters history and best results from the state directory.
     :param state_dir: directory to load from
     :return: a ParameterHistory object loaded from the state directory
     """
     history_filename, results_filename, best_pars_filename = cls._get_result_filenames(
         state_dir)
     self = cls()
     self.history = np.load(history_filename)['history']
     self.best_pars = post(np.load(best_pars_filename)['best_pars'])
     with open(results_filename, 'rb') as results_file:
         data = json.load(results_file)
         self.best_iter = data['best_iter']
         self.best_val_loss = data['best_val_loss']
         self.best_tst_loss = data['best_tst_loss']
         self.start_time = data['start_time']
         self.end_time = data['end_time']
         self.termination_reason = data['termination_reason']
         if 'cfg' in data:
             self.cfg = data['cfg']
         else:
             self.cfg = {}
     return self
Example #12
0
 def load_parameters(self, filename):
     """Loads the parameters of the ParameterSet of the model from the given file."""
     self.ps.data[:] = post(np.load(filename)['ps'])
Example #13
0
 def __setstate__(self, state):
     state['best_pars'] = post(state['best_pars'])
     self.__dict__ = state
Example #14
0
    def generic_training(self,
                         cfg_dir,
                         checkpoint=None,
                         checkpoint_handler=None,
                         loss_record_interval=10,
                         max_missed_val_improvements=200,
                         iteration_gain=1.25,
                         min_improvement=1e-7,
                         reset_termination_criteria=False,
                         desired_loss=None,
                         initialize=True,
                         large_gradient_threshold=0.0,
                         print_gradient_info=False,
                         print_gradient=False,
                         print_parameters=False,
                         log_parameters=[],
                         plot_logged_parameters=True,
                         print_logged_parameters=False,
                         check_gradient_finite=False,
                         log_modules=None):
        """
        Generic training procedure.
        :param cfg_dir: configuration directory
        :param checkpoint: checkpoint
        :param checkpoint_handler: checkpoint handler
        :param loss_record_interval: number of iterations between calculating and recording losses
        :param max_missed_val_improvements: maximum iterations without improvement of loss before training ist stopped
        :param iteration_gain: If not set to 0, then training is performed up to iteration
                               (iteration_gain * iteration_of_last_improvement).
        :param min_improvement: minimum loss change to count as improvement
        :param reset_termination_criteria: resets the termination criteria after loading a checkpoint
        :param desired_loss: if specified, training is terminated with this loss is reached
        :param initialize: if True, the model parameters are initialized using the init_parameters method.
        :param large_gradient_threshold: if specified, a check for large gradient elements that exceed the
                                         given threshold is performed every iteration and they are printed.
        :param print_gradient_info: if True, this function prints diagnostic gradient information
        :param print_gradient: if True, this function prints the full gradient every minibatch
        :param log_parameters: list of parameter names (from self.ps) that should be logged every iteration
                               to file params.out
        :param plot_logged_parameters: if True, logged parameters are plotted to params.png
        :param print_logged_parameters: if True, logged parameters are also printed to standard output
        :param check_gradient_finite: if True, gradient is checked for infs and nans in every iteration
        :param log_modules: list of python modules for which version or latest git commit should be logged
        :return: ParameterHistory object of training
        """

        max_iters = self.cfg.max_iters if self.cfg.has('max_iters') else None

        # build gradient preprocessing chain
        grad_func = self.mb_loss_grad

        # gradient magnitude cap
        if self.cfg.has('gradient_cap'):
            grad_with_cap_orig_grad = grad_func

            def grad_with_cap(pv):
                g = grad_with_cap_orig_grad(pv)
                g_mag = xp.sqrt(xp.sum(g**2))
                if g_mag > self.cfg.gradient_cap:
                    print "gradient magnitude %f is being rescaled" % g_mag
                    g *= self.cfg.gradient_cap / g_mag
                return g

            grad_func = grad_with_cap

        # gradient element cap
        if self.cfg.has('gradient_element_cap'):
            grad_with_element_cap_orig_grad = grad_func

            def grad_with_element_cap(pv):
                g = grad_with_element_cap_orig_grad(pv)
                elems = xp.where(xp.abs(g) > self.cfg.gradient_element_cap)
                g[elems] = xp.sign(g[elems]) * self.cfg.gradient_element_cap
                return g

            grad_func = grad_with_element_cap

        # gradient of constants set to zero
        grad_without_const_orig_grad = grad_func

        def grad_without_const(pv):
            g = grad_without_const_orig_grad(pv)
            return self.ps.nullify_gradient_of_constants(g)

        grad_func = grad_without_const

        # initialize or restore checkpoint, if available
        if not checkpoint:
            itr = 0
            if initialize:
                self.init_parameters()

            his = ParameterHistory(
                cfg=self.cfg,
                state_dir=cfg_dir,
                max_iters=max_iters,
                max_missed_val_improvements=max_missed_val_improvements,
                min_improvement=min_improvement,
                desired_loss=desired_loss,
                iteration_gain=iteration_gain)
            logger = ParameterLogger(out_dir=cfg_dir,
                                     parameters=log_parameters,
                                     plot=plot_logged_parameters,
                                     print_stdout=print_logged_parameters)
            git_log(modules=log_modules, log_dir=cfg_dir)

            # Record initial loss and parameters
            self.record_loss(his, itr)
            logger.log(itr, self.ps)
        else:
            itr = checkpoint['iter']
            self.ps.data[:] = post(checkpoint['data'])
            if 'optimizer_step_rate' in checkpoint:
                self.cfg.optimizer_step_rate = checkpoint[
                    'optimizer_step_rate']

            his = checkpoint['his']
            his.state_dir = cfg_dir
            his.max_missed_val_improvements = max_missed_val_improvements
            his.desired_loss = desired_loss
            his.iteration_gain = iteration_gain

            # start and endtimes in his should have the same length, this is
            # not the case in explicit auto-save checkpoints, therefore set the
            # end to the saved
            if len(his.start_time) != len(his.end_time):
                his.end_time.append(checkpoint['save_time'])
            his.start()

            logger = checkpoint['logger']
            git_log(modules=log_modules, log_dir=cfg_dir, check=True)

        # reset termination criteria if requested
        second_chance_file = join(cfg_dir, "2nd_chance")
        if exists(second_chance_file):
            print "Resetting termination criteria because %s is present" % second_chance_file
            reset_termination_criteria = True
            unlink(second_chance_file)
        if self.cfg.continue_training:
            print "Resetting termination criteria because --continue flag was specified"
            reset_termination_criteria = True
        if reset_termination_criteria:
            his.reset_best()

        if 'step_element_cap' in dir(self.cfg):
            step_element_cap_orig = self.cfg.step_element_cap
        step_element_cap_decrease_iteration = None

        restart = True
        while restart and (max_iters is None or max_iters > 0):
            # create optimizer
            if isinstance(self.cfg.optimizer, dict):

                def wrt_fprime_for_part(partition):
                    wrt_for_part = self.ps.num_partition(partition)

                    def fprime_for_part(pv_part):
                        # we assume that the optimizer updates the ParameterSet inplace and
                        # evaluates the gradient at the current values of the parameters
                        start, stop = self.ps.extents_of_partition(partition)
                        return grad_func(self.ps.num_data)[start:stop]

                    return wrt_for_part, fprime_for_part

                opts_obj = optimizers_from_cfg(self.cfg, wrt_fprime_for_part,
                                               self.mb_loss)
                opts = {
                    part: iter(opt_obj)
                    for part, opt_obj in opts_obj.iteritems()
                }
                partioned_opt = True

                opt_parts = set(opts.keys())
                ps_parts = set(self.ps.partitions)
                if opt_parts != ps_parts:
                    raise ValueError(
                        "optimizer config does not cover all ParameterSet partitions or vice versa: %s"
                        % repr(opt_parts ^ ps_parts))
            else:
                opt = iter(
                    optimizer_from_cfg(self.cfg, self.ps.data, self.mb_loss,
                                       grad_func))
                partioned_opt = False

            # do training
            self.ps.restore_constants()
            last_pars = xp.copy(self.ps.data)
            while not his.should_terminate:
                # call optimizer(s)
                if partioned_opt:
                    for part, opt in opts.iteritems():
                        # print "optimizing %s" % part
                        opt.next()
                else:
                    opt.next()

                # element change cap
                if self.cfg.has('step_element_cap'):
                    d = self.ps.data - last_pars
                    if isinstance(self.cfg.step_element_cap, dict):
                        for par, lim in self.cfg.step_element_cap.iteritems():
                            start, stop = self.ps.extents_of_var(par)
                            dpar = d[start:stop]  # dpar is a subview of d
                            # print "parameter diff for %s is %s (limit is %.4f)" % (par, repr(dpar), lim)
                            elems = xp.where(xp.abs(dpar) > lim)
                            dpar[elems] = xp.sign(dpar[elems]) * lim
                    elif isinstance(self.cfg.step_element_cap, (float, int)):
                        lim = float(self.cfg.step_element_cap)
                        elems = xp.where(xp.abs(d) > lim)
                        d[elems] = xp.sign(d[elems]) * lim
                    else:
                        raise TypeError(
                            "cfg.step_element_cap must either be a dict or a float"
                        )
                    self.ps.data[:] = last_pars + d
                    last_pars = xp.copy(self.ps.data)

                # parameter printout
                if print_parameters:
                    pars = gather(self.ps.data)
                    pars_var = self.ps.split(pars)
                    print "parameters at iteration %d:" % itr
                    for name, value in pars_var.iteritems():
                        print "%10s: %s" % (name, repr(list(value)))

                # obtain gradient if required for debugging operations
                if large_gradient_threshold > 0 or print_gradient_info or print_gradient:
                    gradient = gather(grad_func(self.ps.num_data))
                else:
                    gradient = None

                # check gradient for large elements
                if large_gradient_threshold > 0:
                    lgv = self.ps.find_large_elements(
                        gradient, threshold=large_gradient_threshold)
                    if len(lgv) > 0:
                        print "parameters with large gradient: "
                        for (var, idx), value in lgv.itervalues():
                            print "                                %s[%d] = %.3f" % (
                                var, idx, value)

                # gradient magnitude printout
                if print_gradient_info:
                    gradient_magnitude = np.sqrt(np.sum(gradient**2))
                    print "|gradient| = %.3f" % gradient_magnitude

                # gradient printout
                if print_gradient:
                    gradient_var = self.ps.split(gradient)
                    print "gradient at iteration %d:" % itr
                    for name, value in gradient_var.iteritems():
                        print "%10s: %s" % (name, repr(list(value)))

                # check gradient for NaNs and Infs
                if check_gradient_finite or gradient is not None:
                    if not np.all(np.isfinite(gradient)):
                        his.should_terminate = True
                        his.termination_reason = 'inf_or_nan_gradient'
                        break

                if self.next_minibatch():
                    # iteration finished
                    self.after_iteration(his, itr)

                    itr += 1

                    # log parameters
                    logger.log(itr, self.ps)

                    # calculate losses
                    if itr % loss_record_interval == 0:
                        self.record_loss(his, itr)

                    if step_element_cap_decrease_iteration is not None:
                        if 'step_element_cap_restore_iterations' in dir(
                                self.cfg):
                            restore_itrs = self.cfg.step_element_cap_restore_iterations
                        else:
                            restore_itrs = 100
                        if itr >= step_element_cap_decrease_iteration + restore_itrs:
                            self.cfg.step_element_cap = step_element_cap_orig
                            print "Restored step element cap to %g" % self.cfg.step_element_cap
                            step_element_cap_decrease_iteration = None

                # save checkpoint if necessary
                if checkpoint_handler is not None:
                    if checkpoint_handler.requested:
                        his.stop()
                        checkpoint_handler.save(
                            data=gather(self.ps.data),
                            his=his,
                            iter=itr,
                            logger=logger,
                            optimizer_step_rate=self.cfg.optimizer_step_rate)
                    if his.should_save_checkpoint:
                        checkpoint_handler.save(
                            data=gather(self.ps.data),
                            his=his,
                            iter=itr,
                            logger=logger,
                            optimizer_step_rate=self.cfg.optimizer_step_rate,
                            explicit=True)
                        his.checkpoint_saved()

            # restore best parametes
            self.ps.data[:] = his.best_pars

            # check for retry conditions
            restart = False

            # temporarily reduce step element cap to move over regions with very large gradient
            if (his.should_terminate
                    and his.termination_reason == 'nan_or_inf_loss'
                    and 'step_element_cap' in dir(self.cfg)
                    and 'step_element_cap_min' in dir(self.cfg)
                    and self.cfg.step_element_cap >=
                    self.cfg.step_element_cap_min):
                self.cfg.step_element_cap /= 10.
                step_element_cap_decrease_iteration = itr
                print "Reduced step element cap to %g" % self.cfg.step_element_cap
                his.should_terminate = False
                restart = True

            # advance learning rate schedule
            if (his.should_terminate and his.termination_reason in [
                    'no_improvement', 'nan_or_inf_loss',
                    'user_learning_rate_decrease'
            ] and 'optimizer_step_rate_min' in dir(self.cfg)
                    and self.cfg.optimizer_step_rate / 10. >=
                    self.cfg.optimizer_step_rate_min):
                self.cfg.optimizer_step_rate /= 10.
                print "Decaying optimizer step rate to %g" % self.cfg.optimizer_step_rate
                his.should_terminate = False
                his.last_val_improvement = itr
                restart = True

        # training finished
        self.after_training(his)

        # save results and plot loss
        if checkpoint_handler:
            his.stop()
            checkpoint_handler.save(
                data=gather(self.ps.data),
                his=his,
                iter=itr,
                logger=logger,
                optimizer_step_rate=self.cfg.optimizer_step_rate,
                explicit=True)
        his.finish()
        logger.plot()

        return his
Example #15
0
 def load_parameters(self, filename):
     """Loads the parameters of the ParameterSet of the model from the given file."""
     self.ps.data[:] = post(np.load(filename)['ps'])
Example #16
0
    def generic_training(self, cfg_dir, checkpoint=None, checkpoint_handler=None, loss_record_interval=10,
                         max_missed_val_improvements=200, iteration_gain=1.25, min_improvement=1e-7,
                         reset_termination_criteria=False, desired_loss=None, initialize=True,
                         large_gradient_threshold=0.0, print_gradient_info=False, print_gradient=False,
                         print_parameters=False, log_parameters=[], plot_logged_parameters=True,
                         print_logged_parameters=False, check_gradient_finite=False, log_modules=None):
        """
        Generic training procedure.
        :param cfg_dir: configuration directory
        :param checkpoint: checkpoint
        :param checkpoint_handler: checkpoint handler
        :param loss_record_interval: number of iterations between calculating and recording losses
        :param max_missed_val_improvements: maximum iterations without improvement of loss before training ist stopped
        :param iteration_gain: If not set to 0, then training is performed up to iteration
                               (iteration_gain * iteration_of_last_improvement).
        :param min_improvement: minimum loss change to count as improvement
        :param reset_termination_criteria: resets the termination criteria after loading a checkpoint
        :param desired_loss: if specified, training is terminated with this loss is reached
        :param initialize: if True, the model parameters are initialized using the init_parameters method.
        :param large_gradient_threshold: if specified, a check for large gradient elements that exceed the
                                         given threshold is performed every iteration and they are printed.
        :param print_gradient_info: if True, this function prints diagnostic gradient information
        :param print_gradient: if True, this function prints the full gradient every minibatch
        :param log_parameters: list of parameter names (from self.ps) that should be logged every iteration
                               to file params.out
        :param plot_logged_parameters: if True, logged parameters are plotted to params.png
        :param print_logged_parameters: if True, logged parameters are also printed to standard output
        :param check_gradient_finite: if True, gradient is checked for infs and nans in every iteration
        :param log_modules: list of python modules for which version or latest git commit should be logged
        :return: ParameterHistory object of training
        """

        max_iters = self.cfg.max_iters if self.cfg.has('max_iters') else None

        # build gradient preprocessing chain
        grad_func = self.mb_loss_grad

        # gradient magnitude cap
        if self.cfg.has('gradient_cap'):
            grad_with_cap_orig_grad = grad_func
            def grad_with_cap(pv):
                g = grad_with_cap_orig_grad(pv)
                g_mag = xp.sqrt(xp.sum(g**2))
                if g_mag > self.cfg.gradient_cap:
                    print "gradient magnitude %f is being rescaled" % g_mag
                    g *= self.cfg.gradient_cap / g_mag
                return g
            grad_func = grad_with_cap

        # gradient element cap
        if self.cfg.has('gradient_element_cap'):
            grad_with_element_cap_orig_grad = grad_func
            def grad_with_element_cap(pv):
                g = grad_with_element_cap_orig_grad(pv)
                elems = xp.where(xp.abs(g) > self.cfg.gradient_element_cap)
                g[elems] = xp.sign(g[elems]) * self.cfg.gradient_element_cap
                return g
            grad_func = grad_with_element_cap

        # gradient of constants set to zero
        grad_without_const_orig_grad = grad_func
        def grad_without_const(pv):
            g = grad_without_const_orig_grad(pv)
            return self.ps.nullify_gradient_of_constants(g)
        grad_func = grad_without_const

        # initialize or restore checkpoint, if available
        if not checkpoint:
            itr = 0
            if initialize:
                self.init_parameters()

            his = ParameterHistory(cfg=self.cfg, state_dir=cfg_dir, max_iters=max_iters,
                                   max_missed_val_improvements=max_missed_val_improvements,
                                   min_improvement=min_improvement,
                                   desired_loss=desired_loss, iteration_gain=iteration_gain)
            logger = ParameterLogger(out_dir=cfg_dir, parameters=log_parameters,
                                     plot=plot_logged_parameters, print_stdout=print_logged_parameters)
            git_log(modules=log_modules, log_dir=cfg_dir)

            # Record initial loss and parameters
            self.record_loss(his, itr)
            logger.log(itr, self.ps)
        else:
            itr = checkpoint['iter']
            self.ps.data[:] = post(checkpoint['data'])
            if 'optimizer_step_rate' in checkpoint:
                self.cfg.optimizer_step_rate = checkpoint['optimizer_step_rate']

            his = checkpoint['his']
            his.state_dir = cfg_dir
            his.max_missed_val_improvements = max_missed_val_improvements
            his.desired_loss = desired_loss
            his.iteration_gain = iteration_gain

            # start and endtimes in his should have the same length, this is
            # not the case in explicit auto-save checkpoints, therefore set the
            # end to the saved
            if len(his.start_time) != len(his.end_time):
                his.end_time.append(checkpoint['save_time'])
            his.start()

            logger = checkpoint['logger']
            git_log(modules=log_modules, log_dir=cfg_dir, check=True)

        # reset termination criteria if requested
        second_chance_file = join(cfg_dir, "2nd_chance")
        if exists(second_chance_file):
            print "Resetting termination criteria because %s is present" % second_chance_file
            reset_termination_criteria = True
            unlink(second_chance_file)
        if self.cfg.continue_training:
            print "Resetting termination criteria because --continue flag was specified"
            reset_termination_criteria = True
        if reset_termination_criteria:
            his.reset_best()

        if 'step_element_cap' in dir(self.cfg):
            step_element_cap_orig = self.cfg.step_element_cap
        step_element_cap_decrease_iteration = None

        restart = True
        while restart and (max_iters is None or max_iters > 0):
            # create optimizer
            if isinstance(self.cfg.optimizer, dict):
                def wrt_fprime_for_part(partition):
                    wrt_for_part = self.ps.num_partition(partition)
                    def fprime_for_part(pv_part):
                        # we assume that the optimizer updates the ParameterSet inplace and
                        # evaluates the gradient at the current values of the parameters
                        start, stop = self.ps.extents_of_partition(partition)
                        return grad_func(self.ps.num_data)[start : stop]
                    return wrt_for_part, fprime_for_part
                opts_obj = optimizers_from_cfg(self.cfg, wrt_fprime_for_part, self.mb_loss)
                opts = {part: iter(opt_obj) for part, opt_obj in opts_obj.iteritems()}
                partioned_opt = True

                opt_parts = set(opts.keys())
                ps_parts = set(self.ps.partitions)
                if opt_parts != ps_parts:
                    raise ValueError("optimizer config does not cover all ParameterSet partitions or vice versa: %s" %
                                     repr(opt_parts ^ ps_parts))
            else:
                opt = iter(optimizer_from_cfg(self.cfg, self.ps.data, self.mb_loss, grad_func))
                partioned_opt = False

            # do training
            self.ps.restore_constants()
            last_pars = xp.copy(self.ps.data)
            while not his.should_terminate:
                # call optimizer(s)
                if partioned_opt:
                    for part, opt in opts.iteritems():
                        # print "optimizing %s" % part
                        opt.next()
                else:
                    opt.next()

                # element change cap
                if self.cfg.has('step_element_cap'):
                    d = self.ps.data - last_pars
                    if isinstance(self.cfg.step_element_cap, dict):
                        for par, lim in self.cfg.step_element_cap.iteritems():
                            start, stop = self.ps.extents_of_var(par)
                            dpar = d[start:stop]    # dpar is a subview of d
                            # print "parameter diff for %s is %s (limit is %.4f)" % (par, repr(dpar), lim)
                            elems = xp.where(xp.abs(dpar) > lim)
                            dpar[elems] = xp.sign(dpar[elems]) * lim
                    elif isinstance(self.cfg.step_element_cap, (float, int)):
                        lim = float(self.cfg.step_element_cap)
                        elems = xp.where(xp.abs(d) > lim)
                        d[elems] = xp.sign(d[elems]) * lim
                    else:
                        raise TypeError("cfg.step_element_cap must either be a dict or a float")
                    self.ps.data[:] = last_pars + d
                    last_pars = xp.copy(self.ps.data)

                # parameter printout
                if print_parameters:
                    pars = gather(self.ps.data)
                    pars_var = self.ps.split(pars)
                    print "parameters at iteration %d:" % itr
                    for name, value in pars_var.iteritems():
                        print "%10s: %s" % (name, repr(list(value)))

                # obtain gradient if required for debugging operations
                if large_gradient_threshold > 0 or print_gradient_info or print_gradient:
                    gradient = gather(grad_func(self.ps.num_data))
                else:
                    gradient = None

                # check gradient for large elements
                if large_gradient_threshold > 0:
                    lgv = self.ps.find_large_elements(gradient, threshold=large_gradient_threshold)
                    if len(lgv) > 0:
                        print "parameters with large gradient: "
                        for (var, idx), value in lgv.itervalues():
                            print "                                %s[%d] = %.3f" % (var, idx, value)

                # gradient magnitude printout
                if print_gradient_info:
                    gradient_magnitude = np.sqrt(np.sum(gradient ** 2))
                    print "|gradient| = %.3f" % gradient_magnitude

                # gradient printout
                if print_gradient:
                    gradient_var = self.ps.split(gradient)
                    print "gradient at iteration %d:" % itr
                    for name, value in gradient_var.iteritems():
                        print "%10s: %s" % (name, repr(list(value)))

                # check gradient for NaNs and Infs
                if check_gradient_finite or gradient is not None:
                    if not np.all(np.isfinite(gradient)):
                        his.should_terminate = True
                        his.termination_reason = 'inf_or_nan_gradient'
                        break

                if self.next_minibatch():
                    # iteration finished
                    self.after_iteration(his, itr)

                    itr += 1

                    # log parameters
                    logger.log(itr, self.ps)

                    # calculate losses
                    if itr % loss_record_interval == 0:
                        self.record_loss(his, itr)

                    if step_element_cap_decrease_iteration is not None:
                        if 'step_element_cap_restore_iterations' in dir(self.cfg):
                            restore_itrs = self.cfg.step_element_cap_restore_iterations
                        else:
                            restore_itrs = 100
                        if itr >= step_element_cap_decrease_iteration + restore_itrs:
                            self.cfg.step_element_cap = step_element_cap_orig
                            print "Restored step element cap to %g" % self.cfg.step_element_cap
                            step_element_cap_decrease_iteration = None

                # save checkpoint if necessary
                if checkpoint_handler is not None:
                    if checkpoint_handler.requested:
                        his.stop()
                        checkpoint_handler.save(data=gather(self.ps.data), his=his, iter=itr, logger=logger,
                                                optimizer_step_rate=self.cfg.optimizer_step_rate)
                    if his.should_save_checkpoint:
                        checkpoint_handler.save(data=gather(self.ps.data), his=his, iter=itr, logger=logger,
                                                optimizer_step_rate=self.cfg.optimizer_step_rate,
                                                explicit=True)
                        his.checkpoint_saved()

            # restore best parametes
            self.ps.data[:] = his.best_pars

            # check for retry conditions
            restart = False

            # temporarily reduce step element cap to move over regions with very large gradient
            if (his.should_terminate and his.termination_reason == 'nan_or_inf_loss' and
                    'step_element_cap' in dir(self.cfg) and 'step_element_cap_min' in dir(self.cfg) and
                    self.cfg.step_element_cap >= self.cfg.step_element_cap_min):
                self.cfg.step_element_cap /= 10.
                step_element_cap_decrease_iteration = itr
                print "Reduced step element cap to %g" % self.cfg.step_element_cap
                his.should_terminate = False
                restart = True

            # advance learning rate schedule
            if (his.should_terminate and
                    his.termination_reason in ['no_improvement', 'nan_or_inf_loss', 'user_learning_rate_decrease'] and
                    'optimizer_step_rate_min' in dir(self.cfg) and
                    self.cfg.optimizer_step_rate / 10. >= self.cfg.optimizer_step_rate_min):
                self.cfg.optimizer_step_rate /= 10.
                print "Decaying optimizer step rate to %g" % self.cfg.optimizer_step_rate
                his.should_terminate = False
                his.last_val_improvement = itr
                restart = True

        # training finished
        self.after_training(his)

        # save results and plot loss
        if checkpoint_handler:
            his.stop()
            checkpoint_handler.save(data=gather(self.ps.data), his=his, iter=itr, logger=logger,
                                    optimizer_step_rate=self.cfg.optimizer_step_rate,
                                    explicit=True)
        his.finish()
        logger.plot()

        return his
Example #17
0
 def __setstate__(self, state):
     state['best_pars'] = post(state['best_pars'])
     self.__dict__ = state