R/softmax.r

In jasonxu90/splitFeas: Multiset Split Feasibility

#' Compute soft-max
#' 
#' \code{softmax} returns the soft maximum of a collection of reals.
#' 
#' @param x input
#' @param a scaling factor
#' @export
#' @examples
#' set.seed(12345)
#' n <- 10
#' x <- rnorm(n)
#' softmax(x)
softmax <- function(x,a=100){
  return(logSumExp(x*a)/a)
}