u <- r3vel(10)
v <- r3vel(10)
w <- r3vel(10)
x <- as.3vel(c(0.4,0.1,-0.5))
y <- as.3vel(c(0.1,0.2,-0.7))
z <- as.3vel(c(0.2,0.3,-0.1))
gyr(u,v,x)  # gyr[u,v]x
f <- gyrfun(u,v)
g <- gyrfun(v,u)
f(x)
f(r3vel(10))
f(g(x)) - x              # zero, by eqn 9
g(f(x)) - x              # zero, by eqn 9
(x+y) - f(y+x)           # zero by eqn 10
(u+(v+w)) - ((u+v)+f(w)) # zero by eqn 11
# Following taken from Sbitneva 2001:
rbind(x+(y+(x+z))  ,   (x+(y+x))+z)   # left Bol property
rbind((x+y)+(x+y)  ,   x+(y+(y+x)))   # left Bruck property
sol(299792458)   # speed of light in SI
as.3vel(c(1000,3000,1000)) + as.3vel(c(1000,3000,1000))
## should be close to Galilean result
sol(1)   # revert to default c=1
Run the code above in your browser using DataLab