def METHOD_NAME(self, lcs_file, params_file, save_rate=10, **kwargs):
    '''Computes the dataset and iteratively saves the results to a light curves and parameters files.'''
    self.sample_dist()
    self.build_default_bundle()
    lcs = np.nan*np.ones((self.size, self.phases_size))
    params = np.nan*np.ones((self.size, len(self.distributions.keys())))
    for i in range(self.size):
        params[i], lcs[i] = self.compute_model(i, **kwargs)
        if i%save_rate==0:
            np.save(params_file, params)
            np.save(lcs_file, lcs)