R/sum.r

In float: 32-Bit Floats

#' sum
#' 
#' Sums any combination of float/numeric vector(s)/matri[x|ces].
#' 
#' @details
#' If there are any elements in \code{...}, all elements in the list will first
#' be summed in their native precision, then converted to double precision so
#' they can be combined with \code{base::sum()}.  The final result will be cast
#' to single precision if \code{...} contains only integer and/or float objects.
#' Otherwise, the return will be double precision.
#' 
#' @param x
#' A float matrix.
#' @param ...
#' Additional elements (numeric/float vectors/matrices) to sum.
#' @param na.rm
#' should NA's be removed?
#'  
#' @return 
#' A single value.
#' 
#' @examples
#' library(float)
#' x = flrunif(10, 3)
#' 
#' sum(x)
#' sum(x, 1)
#' 
#' @useDynLib float R_sum_spm
#' @name sum
#' @rdname sum
NULL



sum_float32_ = function(x, na.rm)
{
  ret = .Call(R_sum_spm, DATA(x), na.rm)
  float32(ret)
}

summer = function(x, na.rm)
{
  if (is.float(x))
    dbl(sum_float32_(x, na.rm=na.rm))
  else
    sum(x, na.rm=na.rm)
}

sum_float32 = function(x, ..., na.rm=FALSE)
{
  additional = list(...)
  xsum = sum_float32_(x, na.rm=na.rm)
  
  if (length(additional) > 0)
  {
    additional_sum = base::sum(sapply(additional, summer, na.rm=na.rm))
    if (any(sapply(additional, is.double)))
      dbl(xsum) + additional_sum
    else
      xsum + fl(additional_sum)
  }
  else
    xsum
}



#' @rdname sum
#' @export
setMethod("sum", signature(x="float32"), sum_float32)