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