Create function to compute relative horizontal error ellipse and relative height error between two 3D stations given the variance matrices and corresponding covariance block.
Hi Nic, looks pretty great to me. As this is using a couple of functions that Craig wrote (and he's more familiar with this stuff than I am) I've requested that he review it.
I'd suggest not. I adopt the same strategy as DynAdjust and rotate all vcv components into the local frame based on the position of station 1 (See dnaadjust.cpp line 9958). If the local frame at station 2, or a mid-point, is preferred the user can simply feed in those coordinates as desired.
Create function to compute relative horizontal error ellipse and relative height error between two 3D stations given the variance matrices and corresponding covariance block.