two-to-infinity-loss.R

# rotate Y to align it to X
mypro <- function(X, Y) {
  XY <- crossprod(X, Y)
  s <- svd(XY)
  rotation <- s$v %*% t(s$u)
  Yrot <- Y %*% rotation
  diff <- X - Yrot

  colwise_loss <- Matrix::colSums(diff^2)
  two_infty_loss <- max(sqrt(Matrix::rowSums(diff^2)))

  list(
    rotation = rotation,
    Yrot = Yrot,
    colwise_loss = colwise_loss,
    two_infty_loss = two_infty_loss,
    sum_squared_error = sum(colwise_loss)
  )
}