sol() # returns current speed of light
sol(299792458) # use SI units
sol() # speed of light now SI value
eta() # note [t,t] term
u <- as.3vel(c(100,200,300)) # fast terrestrial speed, but not relativistic
boost(u) # boost matrix practically Galilean
is.consistent.boost(boost(u)) # should be TRUE
sol(1) # revert to relativistic units
Run the code above in your browser using DataLab