In this version, I have rearranged the 4th Do-loop to make it more simple. These simulations actually prove the step from eq. (117) to eq. (118) in the analytical math above.
Oops! I messed that up a little bit. Had to do a trick to get the bivectors to show up in the "cross products vanish" part. It is actually zero times the bivectors but the Clifford package won't do it.
You probably guessed it... Joy's updated local 3-Sphere model using quaternions via the Clifford algebra package. Also with some geometric algebra injected for r0 and Lr0. A cool thing about the Clifford package is that you can extract the pure imaginary part of a quaternion which I did for r1 and r2. But since "i" is used for quaternions it can't be used as an index so have to use "l".