gbdt.cpp 34.6 KB
Newer Older
1
2
3
4
/*!
 * Copyright (c) 2016 Microsoft Corporation. All rights reserved.
 * Licensed under the MIT License. See LICENSE file in the project root for license information.
 */
Guolin Ke's avatar
Guolin Ke committed
5
6
7
#include "gbdt.h"

#include <LightGBM/metric.h>
Guolin Ke's avatar
Guolin Ke committed
8
#include <LightGBM/network.h>
9
10
11
12
#include <LightGBM/objective_function.h>
#include <LightGBM/prediction_early_stop.h>
#include <LightGBM/utils/common.h>
#include <LightGBM/utils/openmp_wrapper.h>
13
#include <LightGBM/sample_strategy.h>
Guolin Ke's avatar
Guolin Ke committed
14

15
16
#include <chrono>
#include <ctime>
17
#include <queue>
18
19
#include <sstream>

Guolin Ke's avatar
Guolin Ke committed
20
21
namespace LightGBM {

22
23
Common::Timer global_timer;

24
25
26
int LGBM_config_::current_device = lgbm_device_cpu;
int LGBM_config_::current_learner = use_cpu_learner;

27
28
29
GBDT::GBDT()
    : iter_(0),
      train_data_(nullptr),
30
      config_(nullptr),
31
32
      objective_function_(nullptr),
      early_stopping_round_(0),
33
      early_stopping_min_delta_(0.0),
34
35
36
37
38
39
      es_first_metric_only_(false),
      max_feature_idx_(0),
      num_tree_per_iteration_(1),
      num_class_(1),
      num_iteration_for_pred_(0),
      shrinkage_rate_(0.1f),
40
      num_init_iteration_(0) {
Guolin Ke's avatar
Guolin Ke committed
41
  average_output_ = false;
Guolin Ke's avatar
Guolin Ke committed
42
  tree_learner_ = nullptr;
43
  linear_tree_ = false;
44
  data_sample_strategy_.reset(nullptr);
shiyu1994's avatar
shiyu1994 committed
45
46
47
  gradients_pointer_ = nullptr;
  hessians_pointer_ = nullptr;
  boosting_on_gpu_ = false;
Guolin Ke's avatar
Guolin Ke committed
48
49
50
51
52
}

GBDT::~GBDT() {
}

Guolin Ke's avatar
Guolin Ke committed
53
void GBDT::Init(const Config* config, const Dataset* train_data, const ObjectiveFunction* objective_function,
54
                const std::vector<const Metric*>& training_metrics) {
Nikita Titov's avatar
Nikita Titov committed
55
  CHECK_NOTNULL(train_data);
56
  train_data_ = train_data;
57
  if (!config->monotone_constraints.empty()) {
Nikita Titov's avatar
Nikita Titov committed
58
    CHECK_EQ(static_cast<size_t>(train_data_->num_total_features()), config->monotone_constraints.size());
Nikita Titov's avatar
Nikita Titov committed
59
  }
60
  if (!config->feature_contri.empty()) {
Nikita Titov's avatar
Nikita Titov committed
61
    CHECK_EQ(static_cast<size_t>(train_data_->num_total_features()), config->feature_contri.size());
62
  }
63
  iter_ = 0;
wxchan's avatar
wxchan committed
64
  num_iteration_for_pred_ = 0;
65
  max_feature_idx_ = 0;
wxchan's avatar
wxchan committed
66
  num_class_ = config->num_class;
Guolin Ke's avatar
Guolin Ke committed
67
68
  config_ = std::unique_ptr<Config>(new Config(*config));
  early_stopping_round_ = config_->early_stopping_round;
69
  early_stopping_min_delta_ = config->early_stopping_min_delta;
70
  es_first_metric_only_ = config_->first_metric_only;
Guolin Ke's avatar
Guolin Ke committed
71
  shrinkage_rate_ = config_->learning_rate;
72

73
  if (config_->device_type == std::string("cuda")) {
74
    LGBM_config_::current_learner = use_cuda_learner;
75
76
    #ifdef USE_CUDA
    if (config_->device_type == std::string("cuda")) {
77
78
79
      const int gpu_device_id = config_->gpu_device_id >= 0 ? config_->gpu_device_id : 0;
      CUDASUCCESS_OR_FATAL(cudaSetDevice(gpu_device_id));
    }
80
    #endif  // USE_CUDA
81
82
  }

83
  // load forced_splits file
84
85
86
87
88
  if (!config->forcedsplits_filename.empty()) {
    std::ifstream forced_splits_file(config->forcedsplits_filename.c_str());
    std::stringstream buffer;
    buffer << forced_splits_file.rdbuf();
    std::string err;
Guolin Ke's avatar
Guolin Ke committed
89
    forced_splits_json_ = Json::parse(buffer.str(), &err);
90
91
  }

92
93
94
  objective_function_ = objective_function;
  num_tree_per_iteration_ = num_class_;
  if (objective_function_ != nullptr) {
Guolin Ke's avatar
Guolin Ke committed
95
    num_tree_per_iteration_ = objective_function_->NumModelPerIteration();
96
97
98
    if (objective_function_->IsRenewTreeOutput() && !config->monotone_constraints.empty()) {
      Log::Fatal("Cannot use ``monotone_constraints`` in %s objective, please disable it.", objective_function_->GetName());
    }
99
100
  }

101
  data_sample_strategy_.reset(SampleStrategy::CreateSampleStrategy(config_.get(), train_data_, objective_function_, num_tree_per_iteration_));
102
103
  is_constant_hessian_ = GetIsConstHessian(objective_function);

104
105
106
  boosting_on_gpu_ = objective_function_ != nullptr && objective_function_->IsCUDAObjective() &&
                     !data_sample_strategy_->IsHessianChange();  // for sample strategy with Hessian change, fall back to boosting on CPU

107
  tree_learner_ = std::unique_ptr<TreeLearner>(TreeLearner::CreateTreeLearner(config_->tree_learner, config_->device_type,
shiyu1994's avatar
shiyu1994 committed
108
                                                                              config_.get(), boosting_on_gpu_));
109
110
111

  // init tree learner
  tree_learner_->Init(train_data_, is_constant_hessian_);
112
  tree_learner_->SetForcedSplit(&forced_splits_json_);
113
114
115
116
117
118
119
120

  // push training metrics
  training_metrics_.clear();
  for (const auto& metric : training_metrics) {
    training_metrics_.push_back(metric);
  }
  training_metrics_.shrink_to_fit();

121
122
  #ifdef USE_CUDA
  if (config_->device_type == std::string("cuda")) {
shiyu1994's avatar
shiyu1994 committed
123
    train_score_updater_.reset(new CUDAScoreUpdater(train_data_, num_tree_per_iteration_, boosting_on_gpu_));
124
  } else {
125
  #endif  // USE_CUDA
126
    train_score_updater_.reset(new ScoreUpdater(train_data_, num_tree_per_iteration_));
127
  #ifdef USE_CUDA
128
  }
129
  #endif  // USE_CUDA
130
131

  num_data_ = train_data_->num_data();
shiyu1994's avatar
shiyu1994 committed
132

133
134
135
136
137
138
139
  // get max feature index
  max_feature_idx_ = train_data_->num_total_features() - 1;
  // get label index
  label_idx_ = train_data_->label_idx();
  // get feature names
  feature_names_ = train_data_->feature_names();
  feature_infos_ = train_data_->feature_infos();
140
  monotone_constraints_ = config->monotone_constraints;
141
142
  // get parser config file content
  parser_config_str_ = train_data_->parser_config_str();
143

144
145
146
  // check that forced splits does not use feature indices larger than dataset size
  CheckForcedSplitFeatures();

147
  // if need bagging, create buffer
148
149
  data_sample_strategy_->ResetSampleConfig(config_.get(), true);
  ResetGradientBuffers();
150
151
152

  class_need_train_ = std::vector<bool>(num_tree_per_iteration_, true);
  if (objective_function_ != nullptr && objective_function_->SkipEmptyClass()) {
Nikita Titov's avatar
Nikita Titov committed
153
    CHECK_EQ(num_tree_per_iteration_, num_class_);
154
155
    for (int i = 0; i < num_class_; ++i) {
      class_need_train_[i] = objective_function_->ClassNeedTrain(i);
156
157
    }
  }
158
159
160
161

  if (config_->linear_tree) {
    linear_tree_ = true;
  }
wxchan's avatar
wxchan committed
162
163
}

164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
void GBDT::CheckForcedSplitFeatures() {
  std::queue<Json> forced_split_nodes;
  forced_split_nodes.push(forced_splits_json_);
  while (!forced_split_nodes.empty()) {
    Json node = forced_split_nodes.front();
    forced_split_nodes.pop();
    const int feature_index = node["feature"].int_value();
    if (feature_index > max_feature_idx_) {
      Log::Fatal("Forced splits file includes feature index %d, but maximum feature index in dataset is %d",
        feature_index, max_feature_idx_);
    }
    if (node.object_items().count("left") > 0) {
      forced_split_nodes.push(node["left"]);
    }
    if (node.object_items().count("right") > 0) {
      forced_split_nodes.push(node["right"]);
    }
  }
}

wxchan's avatar
wxchan committed
184
void GBDT::AddValidDataset(const Dataset* valid_data,
185
                           const std::vector<const Metric*>& valid_metrics) {
wxchan's avatar
wxchan committed
186
  if (!train_data_->CheckAlign(*valid_data)) {
187
    Log::Fatal("Cannot add validation data, since it has different bin mappers with training data");
188
  }
Guolin Ke's avatar
Guolin Ke committed
189
  // for a validation dataset, we need its score and metric
190
  auto new_score_updater =
191
192
    #ifdef USE_CUDA
    config_->device_type == std::string("cuda") ?
193
194
    std::unique_ptr<CUDAScoreUpdater>(new CUDAScoreUpdater(valid_data, num_tree_per_iteration_,
      objective_function_ != nullptr && objective_function_->IsCUDAObjective())) :
195
    #endif  // USE_CUDA
196
    std::unique_ptr<ScoreUpdater>(new ScoreUpdater(valid_data, num_tree_per_iteration_));
wxchan's avatar
wxchan committed
197
198
  // update score
  for (int i = 0; i < iter_; ++i) {
199
200
201
    for (int cur_tree_id = 0; cur_tree_id < num_tree_per_iteration_; ++cur_tree_id) {
      auto curr_tree = (i + num_init_iteration_) * num_tree_per_iteration_ + cur_tree_id;
      new_score_updater->AddScore(models_[curr_tree].get(), cur_tree_id);
wxchan's avatar
wxchan committed
202
203
    }
  }
Guolin Ke's avatar
Guolin Ke committed
204
  valid_score_updater_.push_back(std::move(new_score_updater));
Guolin Ke's avatar
Guolin Ke committed
205
206
207
208
  valid_metrics_.emplace_back();
  for (const auto& metric : valid_metrics) {
    valid_metrics_.back().push_back(metric);
  }
Guolin Ke's avatar
Guolin Ke committed
209
  valid_metrics_.back().shrink_to_fit();
210

211
212
213
214
215
216
217
  if (early_stopping_round_ > 0) {
    auto num_metrics = valid_metrics.size();
    if (es_first_metric_only_) { num_metrics = 1; }
    best_iter_.emplace_back(num_metrics, 0);
    best_score_.emplace_back(num_metrics, kMinScore);
    best_msg_.emplace_back(num_metrics);
  }
Guolin Ke's avatar
Guolin Ke committed
218
219
}

Guolin Ke's avatar
Guolin Ke committed
220
void GBDT::Boosting() {
221
  Common::FunctionTimer fun_timer("GBDT::Boosting", global_timer);
Guolin Ke's avatar
Guolin Ke committed
222
  if (objective_function_ == nullptr) {
223
    Log::Fatal("No objective function provided");
Guolin Ke's avatar
Guolin Ke committed
224
225
226
  }
  // objective function will calculate gradients and hessians
  int64_t num_score = 0;
227
228
229
230
231
232
233
234
  if (config_->bagging_by_query) {
    data_sample_strategy_->Bagging(iter_, tree_learner_.get(), gradients_.data(), hessians_.data());
    objective_function_->
      GetGradients(GetTrainingScore(&num_score), data_sample_strategy_->num_sampled_queries(), data_sample_strategy_->sampled_query_indices(), gradients_pointer_, hessians_pointer_);
  } else {
    objective_function_->
      GetGradients(GetTrainingScore(&num_score), gradients_pointer_, hessians_pointer_);
  }
Guolin Ke's avatar
Guolin Ke committed
235
236
}

Guolin Ke's avatar
Guolin Ke committed
237
void GBDT::Train(int snapshot_freq, const std::string& model_output_path) {
238
  Common::FunctionTimer fun_timer("GBDT::Train", global_timer);
Guolin Ke's avatar
Guolin Ke committed
239
240
  bool is_finished = false;
  auto start_time = std::chrono::steady_clock::now();
Guolin Ke's avatar
Guolin Ke committed
241
  for (int iter = 0; iter < config_->num_iterations && !is_finished; ++iter) {
Guolin Ke's avatar
Guolin Ke committed
242
243
244
245
    is_finished = TrainOneIter(nullptr, nullptr);
    if (!is_finished) {
      is_finished = EvalAndCheckEarlyStopping();
    }
Guolin Ke's avatar
Guolin Ke committed
246
247
248
249
250
251
252
    auto end_time = std::chrono::steady_clock::now();
    // output used time per iteration
    Log::Info("%f seconds elapsed, finished iteration %d", std::chrono::duration<double,
              std::milli>(end_time - start_time) * 1e-3, iter + 1);
    if (snapshot_freq > 0
        && (iter + 1) % snapshot_freq == 0) {
      std::string snapshot_out = model_output_path + ".snapshot_iter_" + std::to_string(iter + 1);
253
      SaveModelToFile(0, -1, config_->saved_feature_importance_type, snapshot_out.c_str());
Guolin Ke's avatar
Guolin Ke committed
254
255
256
257
    }
  }
}

258
259
260
261
262
void GBDT::RefitTree(const int* tree_leaf_prediction, const size_t nrow, const size_t ncol) {
  CHECK_GT(nrow * ncol, 0);
  CHECK_EQ(static_cast<size_t>(num_data_), nrow);
  CHECK_EQ(models_.size(), ncol);

263
264
  int num_iterations = static_cast<int>(models_.size() / num_tree_per_iteration_);
  std::vector<int> leaf_pred(num_data_);
265
266
  if (linear_tree_) {
    std::vector<int> max_leaves_by_thread = std::vector<int>(OMP_NUM_THREADS(), 0);
267
    #pragma omp parallel for num_threads(OMP_NUM_THREADS()) schedule(static)
268
    for (int i = 0; i < static_cast<int>(nrow); ++i) {
269
      int tid = omp_get_thread_num();
270
271
      for (size_t j = 0; j < ncol; ++j) {
        max_leaves_by_thread[tid] = std::max(max_leaves_by_thread[tid], tree_leaf_prediction[i * ncol + j]);
272
273
274
275
276
277
      }
    }
    int max_leaves = *std::max_element(max_leaves_by_thread.begin(), max_leaves_by_thread.end());
    max_leaves += 1;
    tree_learner_->InitLinear(train_data_, max_leaves);
  }
278

279
280
281
282
  for (int iter = 0; iter < num_iterations; ++iter) {
    Boosting();
    for (int tree_id = 0; tree_id < num_tree_per_iteration_; ++tree_id) {
      int model_index = iter * num_tree_per_iteration_ + tree_id;
283
      #pragma omp parallel for num_threads(OMP_NUM_THREADS()) schedule(static)
284
      for (int i = 0; i < num_data_; ++i) {
285
        leaf_pred[i] = tree_leaf_prediction[i * ncol + model_index];
Nikita Titov's avatar
Nikita Titov committed
286
        CHECK_LT(leaf_pred[i], models_[model_index]->num_leaves());
287
      }
288
      size_t offset = static_cast<size_t>(tree_id) * num_data_;
289
290
      auto grad = gradients_pointer_ + offset;
      auto hess = hessians_pointer_ + offset;
291
292
293
294
295
296
297
      auto new_tree = tree_learner_->FitByExistingTree(models_[model_index].get(), leaf_pred, grad, hess);
      train_score_updater_->AddScore(tree_learner_.get(), new_tree, tree_id);
      models_[model_index].reset(new_tree);
    }
  }
}

Andrew Ziem's avatar
Andrew Ziem committed
298
/* If the custom "average" is implemented it will be used in place of the label average (if enabled)
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
*
* An improvement to this is to have options to explicitly choose
* (i) standard average
* (ii) custom average if available
* (iii) any user defined scalar bias (e.g. using a new option "init_score" that overrides (i) and (ii) )
*
* (i) and (ii) could be selected as say "auto_init_score" = 0 or 1 etc..
*
*/
double ObtainAutomaticInitialScore(const ObjectiveFunction* fobj, int class_id) {
  double init_score = 0.0;
  if (fobj != nullptr) {
    init_score = fobj->BoostFromScore(class_id);
  }
  if (Network::num_machines() > 1) {
    init_score = Network::GlobalSyncUpByMean(init_score);
  }
  return init_score;
}

Guolin Ke's avatar
Guolin Ke committed
319
double GBDT::BoostFromAverage(int class_id, bool update_scorer) {
320
  Common::FunctionTimer fun_timer("GBDT::BoostFromAverage", global_timer);
321
  // boosting from average label; or customized "average" if implemented for the current objective
322
323
324
  if (models_.empty() && !train_score_updater_->has_init_score() && objective_function_ != nullptr) {
    if (config_->boost_from_average || (train_data_ != nullptr && train_data_->num_features() == 0)) {
      double init_score = ObtainAutomaticInitialScore(objective_function_, class_id);
325
      if (std::fabs(init_score) > kEpsilon) {
Guolin Ke's avatar
Guolin Ke committed
326
327
328
329
330
        if (update_scorer) {
          train_score_updater_->AddScore(init_score, class_id);
          for (auto& score_updater : valid_score_updater_) {
            score_updater->AddScore(init_score, class_id);
          }
331
332
333
        }
        Log::Info("Start training from score %lf", init_score);
        return init_score;
Guolin Ke's avatar
Guolin Ke committed
334
      }
335
336
337
    } else if (std::string(objective_function_->GetName()) == std::string("regression_l1")
               || std::string(objective_function_->GetName()) == std::string("quantile")
               || std::string(objective_function_->GetName()) == std::string("mape")) {
338
      Log::Warning("Disabling boost_from_average in %s may cause the slow convergence", objective_function_->GetName());
339
    }
340
  }
Guolin Ke's avatar
Guolin Ke committed
341
342
  return 0.0f;
}
Guolin Ke's avatar
Guolin Ke committed
343

Guolin Ke's avatar
Guolin Ke committed
344
bool GBDT::TrainOneIter(const score_t* gradients, const score_t* hessians) {
345
  Common::FunctionTimer fun_timer("GBDT::TrainOneIter", global_timer);
346
  std::vector<double> init_scores(num_tree_per_iteration_, 0.0);
Guolin Ke's avatar
Guolin Ke committed
347
  // boosting first
Guolin Ke's avatar
Guolin Ke committed
348
  if (gradients == nullptr || hessians == nullptr) {
349
    for (int cur_tree_id = 0; cur_tree_id < num_tree_per_iteration_; ++cur_tree_id) {
Guolin Ke's avatar
Guolin Ke committed
350
      init_scores[cur_tree_id] = BoostFromAverage(cur_tree_id, true);
351
    }
Guolin Ke's avatar
Guolin Ke committed
352
    Boosting();
353
354
355
    gradients = gradients_pointer_;
    hessians = hessians_pointer_;
  } else {
shiyu1994's avatar
shiyu1994 committed
356
    // use customized objective function
357
    // the check below fails unless objective=custom is provided in the parameters on Booster creation
shiyu1994's avatar
shiyu1994 committed
358
    CHECK(objective_function_ == nullptr);
359
    if (data_sample_strategy_->IsHessianChange()) {
shiyu1994's avatar
shiyu1994 committed
360
361
      // need to copy customized gradients when using GOSS
      int64_t total_size = static_cast<int64_t>(num_data_) * num_tree_per_iteration_;
362
      #pragma omp parallel for num_threads(OMP_NUM_THREADS()) schedule(static)
shiyu1994's avatar
shiyu1994 committed
363
364
365
366
367
368
      for (int64_t i = 0; i < total_size; ++i) {
        gradients_[i] = gradients[i];
        hessians_[i] = hessians[i];
      }
      CHECK_EQ(gradients_pointer_, gradients_.data());
      CHECK_EQ(hessians_pointer_, hessians_.data());
369
370
371
372
373
      gradients = gradients_pointer_;
      hessians = hessians_pointer_;
    }
  }

374
  // bagging logic
375
376
377
  if (!config_->bagging_by_query) {
    data_sample_strategy_->Bagging(iter_, tree_learner_.get(), gradients_.data(), hessians_.data());
  }
378
379
380
  const bool is_use_subset = data_sample_strategy_->is_use_subset();
  const data_size_t bag_data_cnt = data_sample_strategy_->bag_data_cnt();
  const std::vector<data_size_t, Common::AlignmentAllocator<data_size_t, kAlignedSize>>& bag_data_indices = data_sample_strategy_->bag_data_indices();
Guolin Ke's avatar
Guolin Ke committed
381

382
383
  if (objective_function_ == nullptr && is_use_subset && bag_data_cnt < num_data_ && !boosting_on_gpu_ && !data_sample_strategy_->IsHessianChange()) {
    ResetGradientBuffers();
shiyu1994's avatar
shiyu1994 committed
384
385
  }

Guolin Ke's avatar
Guolin Ke committed
386
  bool should_continue = false;
387
  for (int cur_tree_id = 0; cur_tree_id < num_tree_per_iteration_; ++cur_tree_id) {
388
    const size_t offset = static_cast<size_t>(cur_tree_id) * num_data_;
389
    std::unique_ptr<Tree> new_tree(new Tree(2, false, false));
390
    if (class_need_train_[cur_tree_id] && train_data_->num_features() > 0) {
391
392
      auto grad = gradients + offset;
      auto hess = hessians + offset;
Guolin Ke's avatar
Guolin Ke committed
393
      // need to copy gradients for bagging subset.
394
395
396
397
      if (is_use_subset && bag_data_cnt < num_data_ && !boosting_on_gpu_) {
        for (int i = 0; i < bag_data_cnt; ++i) {
          gradients_pointer_[offset + i] = grad[bag_data_indices[i]];
          hessians_pointer_[offset + i] = hess[bag_data_indices[i]];
Guolin Ke's avatar
Guolin Ke committed
398
        }
399
400
        grad = gradients_pointer_ + offset;
        hess = hessians_pointer_ + offset;
Guolin Ke's avatar
Guolin Ke committed
401
      }
402
403
      bool is_first_tree = models_.size() < static_cast<size_t>(num_tree_per_iteration_);
      new_tree.reset(tree_learner_->Train(grad, hess, is_first_tree));
404
    }
Guolin Ke's avatar
Guolin Ke committed
405

Guolin Ke's avatar
Guolin Ke committed
406
    if (new_tree->num_leaves() > 1) {
Guolin Ke's avatar
Guolin Ke committed
407
      should_continue = true;
408
      auto score_ptr = train_score_updater_->score() + offset;
409
410
      auto residual_getter = [score_ptr](const label_t* label, int i) {return static_cast<double>(label[i]) - score_ptr[i]; };
      tree_learner_->RenewTreeOutput(new_tree.get(), objective_function_, residual_getter,
411
                                     num_data_, bag_data_indices.data(), bag_data_cnt, train_score_updater_->score());
Guolin Ke's avatar
Guolin Ke committed
412
413
414
      // shrinkage by learning rate
      new_tree->Shrinkage(shrinkage_rate_);
      // update score
415
      UpdateScore(new_tree.get(), cur_tree_id);
416
417
      if (std::fabs(init_scores[cur_tree_id]) > kEpsilon) {
        new_tree->AddBias(init_scores[cur_tree_id]);
Guolin Ke's avatar
Guolin Ke committed
418
      }
419
420
    } else {
      // only add default score one-time
421
      if (models_.size() < static_cast<size_t>(num_tree_per_iteration_)) {
422
423
424
425
426
427
        if (objective_function_ != nullptr && !config_->boost_from_average && !train_score_updater_->has_init_score()) {
          init_scores[cur_tree_id] = ObtainAutomaticInitialScore(objective_function_, cur_tree_id);
          // updates scores
          train_score_updater_->AddScore(init_scores[cur_tree_id], cur_tree_id);
          for (auto& score_updater : valid_score_updater_) {
            score_updater->AddScore(init_scores[cur_tree_id], cur_tree_id);
428
          }
429
        }
430
        new_tree->AsConstantTree(init_scores[cur_tree_id]);
431
432
      }
    }
Guolin Ke's avatar
Guolin Ke committed
433
434
435
    // add model
    models_.push_back(std::move(new_tree));
  }
Guolin Ke's avatar
Guolin Ke committed
436

Guolin Ke's avatar
Guolin Ke committed
437
  if (!should_continue) {
438
    Log::Warning("Stopped training because there are no more leaves that meet the split requirements");
439
440
441
442
    if (models_.size() > static_cast<size_t>(num_tree_per_iteration_)) {
      for (int cur_tree_id = 0; cur_tree_id < num_tree_per_iteration_; ++cur_tree_id) {
        models_.pop_back();
      }
Guolin Ke's avatar
Guolin Ke committed
443
444
445
    }
    return true;
  }
446

Guolin Ke's avatar
Guolin Ke committed
447
448
  ++iter_;
  return false;
Guolin Ke's avatar
Guolin Ke committed
449
}
450

wxchan's avatar
wxchan committed
451
void GBDT::RollbackOneIter() {
452
  if (iter_ <= 0) { return; }
wxchan's avatar
wxchan committed
453
  // reset score
454
  for (int cur_tree_id = 0; cur_tree_id < num_tree_per_iteration_; ++cur_tree_id) {
Guolin Ke's avatar
Guolin Ke committed
455
    auto curr_tree = models_.size() - num_tree_per_iteration_ + cur_tree_id;
wxchan's avatar
wxchan committed
456
    models_[curr_tree]->Shrinkage(-1.0);
457
    train_score_updater_->AddScore(models_[curr_tree].get(), cur_tree_id);
wxchan's avatar
wxchan committed
458
    for (auto& score_updater : valid_score_updater_) {
459
      score_updater->AddScore(models_[curr_tree].get(), cur_tree_id);
wxchan's avatar
wxchan committed
460
461
462
    }
  }
  // remove model
463
  for (int cur_tree_id = 0; cur_tree_id < num_tree_per_iteration_; ++cur_tree_id) {
wxchan's avatar
wxchan committed
464
465
466
467
468
    models_.pop_back();
  }
  --iter_;
}

Guolin Ke's avatar
Guolin Ke committed
469
bool GBDT::EvalAndCheckEarlyStopping() {
470
471
  bool is_met_early_stopping = false;
  // print message for metric
Guolin Ke's avatar
Guolin Ke committed
472
  auto best_msg = OutputMetric(iter_);
Guolin Ke's avatar
Guolin Ke committed
473
474


Guolin Ke's avatar
Guolin Ke committed
475
  is_met_early_stopping = !best_msg.empty();
476
477
  if (is_met_early_stopping) {
    Log::Info("Early stopping at iteration %d, the best iteration round is %d",
478
              iter_, iter_ - early_stopping_round_);
Guolin Ke's avatar
Guolin Ke committed
479
    Log::Info("Output of best iteration round:\n%s", best_msg.c_str());
480
    // pop last early_stopping_round_ models
481
    for (int i = 0; i < early_stopping_round_ * num_tree_per_iteration_; ++i) {
482
483
484
485
      models_.pop_back();
    }
  }
  return is_met_early_stopping;
Guolin Ke's avatar
Guolin Ke committed
486
487
}

488
void GBDT::UpdateScore(const Tree* tree, const int cur_tree_id) {
489
  Common::FunctionTimer fun_timer("GBDT::UpdateScore", global_timer);
Guolin Ke's avatar
Guolin Ke committed
490
  // update training score
491
  if (!data_sample_strategy_->is_use_subset()) {
492
    train_score_updater_->AddScore(tree_learner_.get(), tree, cur_tree_id);
Guolin Ke's avatar
Guolin Ke committed
493

494
    const data_size_t bag_data_cnt = data_sample_strategy_->bag_data_cnt();
Guolin Ke's avatar
Guolin Ke committed
495
    // we need to predict out-of-bag scores of data for boosting
496
    if (num_data_ - bag_data_cnt > 0) {
497
498
      #ifdef USE_CUDA
      if (config_->device_type == std::string("cuda")) {
499
        train_score_updater_->AddScore(tree, data_sample_strategy_->cuda_bag_data_indices().RawData() + bag_data_cnt, num_data_ - bag_data_cnt, cur_tree_id);
500
      } else {
501
      #endif  // USE_CUDA
502
        train_score_updater_->AddScore(tree, data_sample_strategy_->bag_data_indices().data() + bag_data_cnt, num_data_ - bag_data_cnt, cur_tree_id);
503
      #ifdef USE_CUDA
504
      }
505
      #endif  // USE_CUDA
Guolin Ke's avatar
Guolin Ke committed
506
507
    }

Guolin Ke's avatar
Guolin Ke committed
508
  } else {
509
    train_score_updater_->AddScore(tree, cur_tree_id);
Guolin Ke's avatar
Guolin Ke committed
510
  }
Guolin Ke's avatar
Guolin Ke committed
511
512


Guolin Ke's avatar
Guolin Ke committed
513
  // update validation score
Guolin Ke's avatar
Guolin Ke committed
514
  for (auto& score_updater : valid_score_updater_) {
515
    score_updater->AddScore(tree, cur_tree_id);
Guolin Ke's avatar
Guolin Ke committed
516
517
518
  }
}

519
#ifdef USE_CUDA
520
521
522
std::vector<double> GBDT::EvalOneMetric(const Metric* metric, const double* score, const data_size_t num_data) const {
#else
std::vector<double> GBDT::EvalOneMetric(const Metric* metric, const double* score, const data_size_t /*num_data*/) const {
523
524
#endif  // USE_CUDA
  #ifdef USE_CUDA
525
  const bool evaluation_on_cuda = metric->IsCUDAMetric();
shiyu1994's avatar
shiyu1994 committed
526
  if ((boosting_on_gpu_ && evaluation_on_cuda) || (!boosting_on_gpu_ && !evaluation_on_cuda)) {
527
  #endif  // USE_CUDA
528
    return metric->Eval(score, objective_function_);
529
  #ifdef USE_CUDA
shiyu1994's avatar
shiyu1994 committed
530
  } else if (boosting_on_gpu_ && !evaluation_on_cuda) {
531
    const size_t total_size = static_cast<size_t>(num_data) * static_cast<size_t>(num_tree_per_iteration_);
532
533
534
535
536
537
    if (total_size > host_score_.size()) {
      host_score_.resize(total_size, 0.0f);
    }
    CopyFromCUDADeviceToHost<double>(host_score_.data(), score, total_size, __FILE__, __LINE__);
    return metric->Eval(host_score_.data(), objective_function_);
  } else {
538
    const size_t total_size = static_cast<size_t>(num_data) * static_cast<size_t>(num_tree_per_iteration_);
539
540
541
542
543
544
    if (total_size > cuda_score_.Size()) {
      cuda_score_.Resize(total_size);
    }
    CopyFromHostToCUDADevice<double>(cuda_score_.RawData(), score, total_size, __FILE__, __LINE__);
    return metric->Eval(cuda_score_.RawData(), objective_function_);
  }
545
  #endif  // USE_CUDA
Guolin Ke's avatar
Guolin Ke committed
546
547
}

Guolin Ke's avatar
Guolin Ke committed
548
std::string GBDT::OutputMetric(int iter) {
Guolin Ke's avatar
Guolin Ke committed
549
  bool need_output = (iter % config_->metric_freq) == 0;
Guolin Ke's avatar
Guolin Ke committed
550
551
  std::string ret = "";
  std::stringstream msg_buf;
552
  std::vector<std::pair<size_t, size_t>> meet_early_stopping_pairs;
Guolin Ke's avatar
Guolin Ke committed
553
  // print training metric
Guolin Ke's avatar
Guolin Ke committed
554
  if (need_output) {
555
556
    for (auto& sub_metric : training_metrics_) {
      auto name = sub_metric->GetName();
557
      auto scores = EvalOneMetric(sub_metric, train_score_updater_->score(), train_score_updater_->num_data());
Guolin Ke's avatar
Guolin Ke committed
558
      for (size_t k = 0; k < name.size(); ++k) {
Guolin Ke's avatar
Guolin Ke committed
559
560
561
562
563
564
        std::stringstream tmp_buf;
        tmp_buf << "Iteration:" << iter
          << ", training " << name[k]
          << " : " << scores[k];
        Log::Info(tmp_buf.str().c_str());
        if (early_stopping_round_ > 0) {
565
          msg_buf << tmp_buf.str() << '\n';
Guolin Ke's avatar
Guolin Ke committed
566
        }
567
      }
568
    }
Guolin Ke's avatar
Guolin Ke committed
569
570
  }
  // print validation metric
Guolin Ke's avatar
Guolin Ke committed
571
  if (need_output || early_stopping_round_ > 0) {
572
573
    for (size_t i = 0; i < valid_metrics_.size(); ++i) {
      for (size_t j = 0; j < valid_metrics_[i].size(); ++j) {
574
        auto test_scores = EvalOneMetric(valid_metrics_[i][j], valid_score_updater_[i]->score(), valid_score_updater_[i]->num_data());
Guolin Ke's avatar
Guolin Ke committed
575
576
577
578
579
580
581
582
583
584
        auto name = valid_metrics_[i][j]->GetName();
        for (size_t k = 0; k < name.size(); ++k) {
          std::stringstream tmp_buf;
          tmp_buf << "Iteration:" << iter
            << ", valid_" << i + 1 << " " << name[k]
            << " : " << test_scores[k];
          if (need_output) {
            Log::Info(tmp_buf.str().c_str());
          }
          if (early_stopping_round_ > 0) {
585
            msg_buf << tmp_buf.str() << '\n';
586
          }
wxchan's avatar
wxchan committed
587
        }
588
        if (es_first_metric_only_ && j > 0) { continue; }
Guolin Ke's avatar
Guolin Ke committed
589
        if (ret.empty() && early_stopping_round_ > 0) {
590
          auto cur_score = valid_metrics_[i][j]->factor_to_bigger_better() * test_scores.back();
591
          if (cur_score - best_score_[i][j] > early_stopping_min_delta_) {
592
            best_score_[i][j] = cur_score;
593
            best_iter_[i][j] = iter;
Guolin Ke's avatar
Guolin Ke committed
594
            meet_early_stopping_pairs.emplace_back(i, j);
595
          } else {
Guolin Ke's avatar
Guolin Ke committed
596
            if (iter - best_iter_[i][j] >= early_stopping_round_) { ret = best_msg_[i][j]; }
597
          }
wxchan's avatar
wxchan committed
598
599
        }
      }
Guolin Ke's avatar
Guolin Ke committed
600
601
    }
  }
Guolin Ke's avatar
Guolin Ke committed
602
603
604
  for (auto& pair : meet_early_stopping_pairs) {
    best_msg_[pair.first][pair.second] = msg_buf.str();
  }
wxchan's avatar
wxchan committed
605
  return ret;
Guolin Ke's avatar
Guolin Ke committed
606
607
}

608
/*! \brief Get eval result */
609
std::vector<double> GBDT::GetEvalAt(int data_idx) const {
Guolin Ke's avatar
Guolin Ke committed
610
  CHECK(data_idx >= 0 && data_idx <= static_cast<int>(valid_score_updater_.size()));
611
612
  std::vector<double> ret;
  if (data_idx == 0) {
613
    for (auto& sub_metric : training_metrics_) {
614
      auto scores = EvalOneMetric(sub_metric, train_score_updater_->score(), train_score_updater_->num_data());
615
616
617
      for (auto score : scores) {
        ret.push_back(score);
      }
618
    }
619
  } else {
620
621
    auto used_idx = data_idx - 1;
    for (size_t j = 0; j < valid_metrics_[used_idx].size(); ++j) {
622
      auto test_scores = EvalOneMetric(valid_metrics_[used_idx][j], valid_score_updater_[used_idx]->score(), valid_score_updater_[used_idx]->num_data());
623
624
625
      for (auto score : test_scores) {
        ret.push_back(score);
      }
626
627
628
629
630
    }
  }
  return ret;
}

Guolin Ke's avatar
Guolin Ke committed
631
/*! \brief Get training scores result */
632
const double* GBDT::GetTrainingScore(int64_t* out_len) {
633
  *out_len = static_cast<int64_t>(train_score_updater_->num_data()) * num_class_;
Guolin Ke's avatar
Guolin Ke committed
634
  return train_score_updater_->score();
635
636
}

637
void GBDT::PredictContrib(const double* features, double* output) const {
638
  // set zero
Guolin Ke's avatar
Guolin Ke committed
639
640
  const int num_features = max_feature_idx_ + 1;
  std::memset(output, 0, sizeof(double) * num_tree_per_iteration_ * (num_features + 1));
641
642
  const int end_iteration_for_pred = start_iteration_for_pred_ + num_iteration_for_pred_;
  for (int i = start_iteration_for_pred_; i < end_iteration_for_pred; ++i) {
643
644
    // predict all the trees for one iteration
    for (int k = 0; k < num_tree_per_iteration_; ++k) {
Guolin Ke's avatar
Guolin Ke committed
645
      models_[i * num_tree_per_iteration_ + k]->PredictContrib(features, num_features, output + k*(num_features + 1));
646
    }
647
648
649
650
651
652
  }
}

void GBDT::PredictContribByMap(const std::unordered_map<int, double>& features,
                               std::vector<std::unordered_map<int, double>>* output) const {
  const int num_features = max_feature_idx_ + 1;
653
654
  const int end_iteration_for_pred = start_iteration_for_pred_ + num_iteration_for_pred_;
  for (int i = start_iteration_for_pred_; i < end_iteration_for_pred; ++i) {
655
656
657
    // predict all the trees for one iteration
    for (int k = 0; k < num_tree_per_iteration_; ++k) {
      models_[i * num_tree_per_iteration_ + k]->PredictContribByMap(features, num_features, &((*output)[k]));
658
659
660
661
    }
  }
}

Guolin Ke's avatar
Guolin Ke committed
662
663
void GBDT::GetPredictAt(int data_idx, double* out_result, int64_t* out_len) {
  CHECK(data_idx >= 0 && data_idx <= static_cast<int>(valid_score_updater_.size()));
Guolin Ke's avatar
Guolin Ke committed
664

665
  const double* raw_scores = nullptr;
Guolin Ke's avatar
Guolin Ke committed
666
667
  data_size_t num_data = 0;
  if (data_idx == 0) {
wxchan's avatar
wxchan committed
668
    raw_scores = GetTrainingScore(out_len);
Guolin Ke's avatar
Guolin Ke committed
669
670
671
672
673
    num_data = train_score_updater_->num_data();
  } else {
    auto used_idx = data_idx - 1;
    raw_scores = valid_score_updater_[used_idx]->score();
    num_data = valid_score_updater_[used_idx]->num_data();
674
    *out_len = static_cast<int64_t>(num_data) * num_class_;
Guolin Ke's avatar
Guolin Ke committed
675
  }
676
  #ifdef USE_CUDA
677
678
679
680
681
682
  std::vector<double> host_raw_scores;
  if (boosting_on_gpu_) {
    host_raw_scores.resize(static_cast<size_t>(*out_len), 0.0);
    CopyFromCUDADeviceToHost<double>(host_raw_scores.data(), raw_scores, static_cast<size_t>(*out_len), __FILE__, __LINE__);
    raw_scores = host_raw_scores.data();
  }
683
  #endif  // USE_CUDA
Guolin Ke's avatar
Guolin Ke committed
684
  if (objective_function_ != nullptr) {
685
    #pragma omp parallel for num_threads(OMP_NUM_THREADS()) schedule(static)
Guolin Ke's avatar
Guolin Ke committed
686
    for (data_size_t i = 0; i < num_data; ++i) {
Guolin Ke's avatar
Guolin Ke committed
687
      std::vector<double> tree_pred(num_tree_per_iteration_);
688
      for (int j = 0; j < num_tree_per_iteration_; ++j) {
Guolin Ke's avatar
Guolin Ke committed
689
        tree_pred[j] = raw_scores[j * num_data + i];
690
      }
Guolin Ke's avatar
Guolin Ke committed
691
692
      std::vector<double> tmp_result(num_class_);
      objective_function_->ConvertOutput(tree_pred.data(), tmp_result.data());
Guolin Ke's avatar
Guolin Ke committed
693
      for (int j = 0; j < num_class_; ++j) {
694
        out_result[j * num_data + i] = static_cast<double>(tmp_result[j]);
Guolin Ke's avatar
Guolin Ke committed
695
696
      }
    }
697
  } else {
698
    #pragma omp parallel for num_threads(OMP_NUM_THREADS()) schedule(static)
Guolin Ke's avatar
Guolin Ke committed
699
    for (data_size_t i = 0; i < num_data; ++i) {
700
      for (int j = 0; j < num_tree_per_iteration_; ++j) {
Guolin Ke's avatar
Guolin Ke committed
701
        out_result[j * num_data + i] = static_cast<double>(raw_scores[j * num_data + i]);
Guolin Ke's avatar
Guolin Ke committed
702
703
704
705
706
      }
    }
  }
}

707
708
double GBDT::GetUpperBoundValue() const {
  double max_value = 0.0;
Nikita Titov's avatar
Nikita Titov committed
709
  for (const auto &tree : models_) {
710
711
712
713
714
715
716
    max_value += tree->GetUpperBoundValue();
  }
  return max_value;
}

double GBDT::GetLowerBoundValue() const {
  double min_value = 0.0;
Nikita Titov's avatar
Nikita Titov committed
717
  for (const auto &tree : models_) {
718
719
720
721
722
    min_value += tree->GetLowerBoundValue();
  }
  return min_value;
}

Guolin Ke's avatar
Guolin Ke committed
723
724
725
void GBDT::ResetTrainingData(const Dataset* train_data, const ObjectiveFunction* objective_function,
                             const std::vector<const Metric*>& training_metrics) {
  if (train_data != train_data_ && !train_data_->CheckAlign(*train_data)) {
726
    Log::Fatal("Cannot reset training data, since new training data has different bin mappers");
wxchan's avatar
wxchan committed
727
728
  }

Guolin Ke's avatar
Guolin Ke committed
729
  objective_function_ = objective_function;
730
  data_sample_strategy_->UpdateObjectiveFunction(objective_function);
Guolin Ke's avatar
Guolin Ke committed
731
  if (objective_function_ != nullptr) {
Nikita Titov's avatar
Nikita Titov committed
732
    CHECK_EQ(num_tree_per_iteration_, objective_function_->NumModelPerIteration());
733
734
735
    if (objective_function_->IsRenewTreeOutput() && !config_->monotone_constraints.empty()) {
      Log::Fatal("Cannot use ``monotone_constraints`` in %s objective, please disable it.", objective_function_->GetName());
    }
736
  }
737
  is_constant_hessian_ = GetIsConstHessian(objective_function);
738

Guolin Ke's avatar
Guolin Ke committed
739
740
741
742
  // push training metrics
  training_metrics_.clear();
  for (const auto& metric : training_metrics) {
    training_metrics_.push_back(metric);
743
  }
Guolin Ke's avatar
Guolin Ke committed
744
  training_metrics_.shrink_to_fit();
745

746
  #ifdef USE_CUDA
747
748
  boosting_on_gpu_ = objective_function_ != nullptr && objective_function_->IsCUDAObjective() &&
                    !data_sample_strategy_->IsHessianChange();  // for sample strategy with Hessian change, fall back to boosting on CPU
shiyu1994's avatar
shiyu1994 committed
749
  tree_learner_->ResetBoostingOnGPU(boosting_on_gpu_);
750
  #endif  // USE_CUDA
751

Guolin Ke's avatar
Guolin Ke committed
752
753
  if (train_data != train_data_) {
    train_data_ = train_data;
754
    data_sample_strategy_->UpdateTrainingData(train_data);
Guolin Ke's avatar
Guolin Ke committed
755
756
    // not same training data, need reset score and others
    // create score tracker
757
758
    #ifdef USE_CUDA
    if (config_->device_type == std::string("cuda")) {
shiyu1994's avatar
shiyu1994 committed
759
      train_score_updater_.reset(new CUDAScoreUpdater(train_data_, num_tree_per_iteration_, boosting_on_gpu_));
760
    } else {
761
    #endif  // USE_CUDA
762
      train_score_updater_.reset(new ScoreUpdater(train_data_, num_tree_per_iteration_));
763
    #ifdef USE_CUDA
764
    }
765
    #endif  // USE_CUDA
766

Guolin Ke's avatar
Guolin Ke committed
767
768
769
770
771
772
    // update score
    for (int i = 0; i < iter_; ++i) {
      for (int cur_tree_id = 0; cur_tree_id < num_tree_per_iteration_; ++cur_tree_id) {
        auto curr_tree = (i + num_init_iteration_) * num_tree_per_iteration_ + cur_tree_id;
        train_score_updater_->AddScore(models_[curr_tree].get(), cur_tree_id);
      }
773
774
    }

Guolin Ke's avatar
Guolin Ke committed
775
    num_data_ = train_data_->num_data();
776

777
    ResetGradientBuffers();
778

Guolin Ke's avatar
Guolin Ke committed
779
780
781
782
    max_feature_idx_ = train_data_->num_total_features() - 1;
    label_idx_ = train_data_->label_idx();
    feature_names_ = train_data_->feature_names();
    feature_infos_ = train_data_->feature_infos();
783
    parser_config_str_ = train_data_->parser_config_str();
784

785
    tree_learner_->ResetTrainingData(train_data, is_constant_hessian_);
786
    data_sample_strategy_->ResetSampleConfig(config_.get(), true);
787
788
  } else {
    tree_learner_->ResetIsConstantHessian(is_constant_hessian_);
789
  }
790
791
}

Guolin Ke's avatar
Guolin Ke committed
792
793
void GBDT::ResetConfig(const Config* config) {
  auto new_config = std::unique_ptr<Config>(new Config(*config));
794
  if (!config->monotone_constraints.empty()) {
Nikita Titov's avatar
Nikita Titov committed
795
    CHECK_EQ(static_cast<size_t>(train_data_->num_total_features()), config->monotone_constraints.size());
796
797
  }
  if (!config->feature_contri.empty()) {
Nikita Titov's avatar
Nikita Titov committed
798
    CHECK_EQ(static_cast<size_t>(train_data_->num_total_features()), config->feature_contri.size());
799
  }
800
801
802
  if (objective_function_ != nullptr && objective_function_->IsRenewTreeOutput() && !config->monotone_constraints.empty()) {
    Log::Fatal("Cannot use ``monotone_constraints`` in %s objective, please disable it.", objective_function_->GetName());
  }
Guolin Ke's avatar
Guolin Ke committed
803
804
805
  early_stopping_round_ = new_config->early_stopping_round;
  shrinkage_rate_ = new_config->learning_rate;
  if (tree_learner_ != nullptr) {
Guolin Ke's avatar
Guolin Ke committed
806
    tree_learner_->ResetConfig(new_config.get());
807
  }
shiyu1994's avatar
shiyu1994 committed
808

809
810
  boosting_on_gpu_ = objective_function_ != nullptr && objective_function_->IsCUDAObjective() &&
                    !data_sample_strategy_->IsHessianChange();  // for sample strategy with Hessian change, fall back to boosting on CPU
shiyu1994's avatar
shiyu1994 committed
811
812
  tree_learner_->ResetBoostingOnGPU(boosting_on_gpu_);

Guolin Ke's avatar
Guolin Ke committed
813
  if (train_data_ != nullptr) {
814
815
816
817
818
    data_sample_strategy_->ResetSampleConfig(new_config.get(), false);
    if (data_sample_strategy_->NeedResizeGradients()) {
      // resize gradient vectors to copy the customized gradients for goss or bagging with subset
      ResetGradientBuffers();
    }
819
  }
820
  if (config_.get() != nullptr && config_->forcedsplits_filename != new_config->forcedsplits_filename) {
821
822
823
824
825
826
827
    // load forced_splits file
    if (!new_config->forcedsplits_filename.empty()) {
      std::ifstream forced_splits_file(
          new_config->forcedsplits_filename.c_str());
      std::stringstream buffer;
      buffer << forced_splits_file.rdbuf();
      std::string err;
Guolin Ke's avatar
Guolin Ke committed
828
      forced_splits_json_ = Json::parse(buffer.str(), &err);
829
830
831
832
833
834
      tree_learner_->SetForcedSplit(&forced_splits_json_);
    } else {
      forced_splits_json_ = Json();
      tree_learner_->SetForcedSplit(nullptr);
    }
  }
Guolin Ke's avatar
Guolin Ke committed
835
  config_.reset(new_config.release());
Guolin Ke's avatar
Guolin Ke committed
836
837
}

838
839
840
841
void GBDT::ResetGradientBuffers() {
  const size_t total_size = static_cast<size_t>(num_data_) * num_tree_per_iteration_;
  const bool is_use_subset = data_sample_strategy_->is_use_subset();
  const data_size_t bag_data_cnt = data_sample_strategy_->bag_data_cnt();
Guolin Ke's avatar
Guolin Ke committed
842
  if (objective_function_ != nullptr) {
843
844
    #ifdef USE_CUDA
    if (config_->device_type == std::string("cuda") && boosting_on_gpu_) {
845
846
847
      if (cuda_gradients_.Size() < total_size) {
        cuda_gradients_.Resize(total_size);
        cuda_hessians_.Resize(total_size);
Guolin Ke's avatar
Guolin Ke committed
848
      }
849
850
851
      gradients_pointer_ = cuda_gradients_.RawData();
      hessians_pointer_ = cuda_hessians_.RawData();
    } else {
852
    #endif  // USE_CUDA
853
      if (gradients_.size() < total_size) {
shiyu1994's avatar
shiyu1994 committed
854
855
        gradients_.resize(total_size);
        hessians_.resize(total_size);
856
      }
857
858
      gradients_pointer_ = gradients_.data();
      hessians_pointer_ = hessians_.data();
859
    #ifdef USE_CUDA
860
    }
861
    #endif  // USE_CUDA
862
863
864
865
866
867
868
  } else if (data_sample_strategy_->IsHessianChange() || (is_use_subset && bag_data_cnt < num_data_ && !boosting_on_gpu_)) {
    if (gradients_.size() < total_size) {
      gradients_.resize(total_size);
      hessians_.resize(total_size);
    }
    gradients_pointer_ = gradients_.data();
    hessians_pointer_ = hessians_.data();
869
  }
wxchan's avatar
wxchan committed
870
871
}

Guolin Ke's avatar
Guolin Ke committed
872
}  // namespace LightGBM