| 342 | } |
| 343 | |
| 344 | void __logistic_trans_compute (MutableClusteredState& state, |
| 345 | const MappedColumnVector& x, const double& y) |
| 346 | { |
| 347 | double sm = 0; |
| 348 | for (int i = 0; i < state.widthOfX; i++) sm += state.coef(i) * x(i); |
| 349 | |
| 350 | double sgn = y > 0 ? -1 : 1; |
| 351 | double t1 = sigma(sgn * sm); |
| 352 | double t2 = sigma(-sgn * sm); |
| 353 | |
| 354 | for (int i = 0; i < state.widthOfX; i++) |
| 355 | state.meat_half(0,i) += t1 * sgn * x(i); |
| 356 | |
| 357 | state.bread += (t1 * t2) * (x * trans(x)); |
| 358 | } |
| 359 | |
| 360 | // ------------------------------------------------------------------------ |
| 361 | |