dyn.load("vecSum.so")
.Call("vecSum", rnorm(10))

vecSum <- function(vec){
  if (!is.vector(vec))
    stop("vec must be a vector")
  if (!is.real(vec)) vec <- as.real(vec)
  .Call("vecSum", vec)
}

vecSum(1:100)
vecSum(matrix(0,2,2))