Open fgnievinski opened 2 years ago
numerical algorithms will still require a shallow loop over points.
background: https://www.mathworks.com/help/matlab/matlab_prog/vectorization.html
care is needed with elevations below spherical horizon (hint: discard upfront, returning NaN)
push the loop over points deeper into the code, down at the quartic roots solver, replacing roots.m for mroots.m (attached) mroots.zip