Skip to contents

rwig 0.2.0

  • Dropped the dependencies on Rcpp and RcppArmadillo. The C++ code talks to R through the native C API (.Call with registered routines) and calls the BLAS/LAPACK libraries shipped with R directly through a small internal matrix layer, which also makes the installed package much smaller. The package now has no compiled-code dependencies at all.
  • Interrupting a long computation (Ctrl-C) now releases worker threads and buffers before returning to R, and surfaces as an R error.
  • Faster log-stabilized sinkhorn() and barycenter(): the soft-min kernels no longer materialize the M x N matrix C - f 1' - 1 g' on every iteration, and worker threads (n_threads) are created once per call instead of once per iteration.
  • Faster wdl() on CPU: the softmax Jacobian is applied in O(N) instead of O(N^2) per column, an unused matrix product per backward step was removed, and the per-batch history buffers use about half the memory.
  • wdl() now draws its random initialization from R’s RNG (rnorm()), so set.seed() reproduces the documented R sequence. Fits with the same seed therefore differ from version 0.1.0, which used Armadillo’s own generator.
  • Fixed: wdl() returned weights, docs_pred and docs_dist in the shuffled training order (the default shuffle = TRUE) next to docs in the input order, so wig() summed document scores into the wrong periods. Per-document outputs are now put back into the input order.
  • Fixed: wig() failed when wig_control was given without group_unit.
  • wdl_specs()/wig_specs() now honor the values passed in barycenter_control (method, max_iter, …) instead of silently overriding them; the wdl_control seed entry is named seed (it was documented as rng_seed but read as seed). verbose defaults to 0 when omitted from a partial sinkhorn_control/barycenter_control.
  • Dropped the dependency on RhpcBLASctl: rwig no longer sets the BLAS thread count to 1 for the whole session when attached. The startup message now explains the recommendation and, if RhpcBLASctl is installed (it is now only suggested), reports the current thread count and the call to change it.
  • Dropped the dependency on lubridate: group_unit is now passed to cut() (?cut.Date), so it accepts “day”, “week”, “month”, “quarter”, “year” and multiples such as “2 months”. Weeks start on Monday.
  • Fixed: threaded (n_threads > 0) log barycenter crashed when the cost matrix had more columns than rows.
  • Fixed: the CUDA build passes an explicit GPU architecture to nvcc (recent CUDA toolkits default to sm_52, which lacks double-precision atomicAdd). Set RWIG_CUDA_ARCH to override the detected flag, or RWIG_NO_CUDA=1 to build without GPU support even when CUDA is installed.

rwig 0.1.0

CRAN release: 2026-04-21