If your two points on a sphere are within the same hemisphere set one pt as the center of a unit circle in a plane tangential to that point and project the other point onto it then rotate the co-ord ...
Trying to implement cordic32 routine provided in code library, on the following formula. Separate value of sin(), cos etc are ok but when used in formula, the result is weird. i.e.
tmp2 = cos(lat ...