| 295 | |
| 296 | namespace log_domain { |
| 297 | double LogSum(const std::vector<double> &xs) { |
| 298 | if (xs.empty()) { |
| 299 | return -1.0 * std::numeric_limits<double>::max(); |
| 300 | } |
| 301 | double sum = xs.front(); |
| 302 | |
| 303 | auto log_add = [](double xa, double xb) { |
| 304 | if (xa > xb) { |
| 305 | std::swap(xa, xb); |
| 306 | } |
| 307 | return xb + std::log1p(std::exp(xa - xb)); |
| 308 | }; |
| 309 | for (int i = 1; i < xs.size(); ++i) { |
| 310 | sum = log_add(sum, xs[i]); |
| 311 | } |
| 312 | return sum; |
| 313 | } |
| 314 | } // namespace log_domain |
| 315 | } // namespace sentencepiece |
no test coverage detected