Почему ** 0,5 оказываются более эффективными, чем SQRT()

Question

Я попытался измерить скорость этих двух способов для принятия квадратного корня:

> system.time(expr = replicate(10000, 1:10000 ** (1/2))) 
## user system elapsed 
## 0.027 0.001 0.028 
> system.time(expr = replicate(10000, sqrt(1:10000))) 
## user system elapsed 
## 3.722 0.665 4.494 

Если функция sqrt() не может конкурировать с ** 0.5, почему нам нужна такая функция?

(система OS X Yusemite и R версия 3.1.2)

Answer 1

Вы забыли важные скобки. Вот тайминги после коррекции, что:

system.time(expr = replicate(10000, (1:10000) ** (1/2))) 
#user system elapsed 
#4.76 0.32 5.12 
system.time(expr = replicate(10000, sqrt(1:10000))) 
#user system elapsed 
#2.67 0.57 3.31 

Answer 2

Чтобы добавить @ ответ Роланда, вы попали в Operators precedence «ловушку». ^ предшествует : («** переводится в парсер ^» согласно документации ?"**")

Что произошло на самом деле является

`:`(1, 10000 ** (1/2)) 

Это означает, что вы впервые запустить ** и только потом 1:..

Совет для будущего, попробуйте отладить код перед запуском сложных операций, например, тестирование

1:5 ** (1/2) 
## [1] 1 2 
sqrt(1:5) 
## [1] 1.000000 1.414214 1.732051 2.000000 2.236068 

Выясните проблему.