Average Error: 0.2 → 0.2
Time: 1.5m
Precision: 64
Internal Precision: 576
\[\lambda_1 + \tan^{-1}_* \frac{\left(\sin theta \cdot \sin delta\right) \cdot \cos \phi_1}{\cos delta - \sin \phi_1 \cdot \sin \left(\sin^{-1} \left(\sin \phi_1 \cdot \cos delta + \left(\cos \phi_1 \cdot \sin delta\right) \cdot \cos theta\right)\right)}\]
↓
\[\lambda_1 + \tan^{-1}_* \frac{\left(\cos \phi_1 \cdot \sin delta\right) \cdot \sin theta}{\frac{\cos delta \cdot \cos delta - \left(\sin \left(\sin^{-1} \left(\left(\cos \phi_1 \cdot \sin delta\right) \cdot \cos theta + \sin \phi_1 \cdot \cos delta\right)\right) \cdot \sin \phi_1\right) \cdot \left(\sin \left(\sin^{-1} \left(\left(\cos \phi_1 \cdot \sin delta\right) \cdot \cos theta + \sin \phi_1 \cdot \cos delta\right)\right) \cdot \sin \phi_1\right)}{\cos delta + \sin \left(\sin^{-1} \left(\left(\cos \phi_1 \cdot \sin delta\right) \cdot \cos theta + \sin \phi_1 \cdot \cos delta\right)\right) \cdot \sin \phi_1}}\]

Error

Bits error versus lambda1

Bits error versus phi1

Bits error versus phi2

Bits error versus delta

Bits error versus theta

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Derivation

  1. Initial program 0.2

    \[\lambda_1 + \tan^{-1}_* \frac{\left(\sin theta \cdot \sin delta\right) \cdot \cos \phi_1}{\cos delta - \sin \phi_1 \cdot \sin \left(\sin^{-1} \left(\sin \phi_1 \cdot \cos delta + \left(\cos \phi_1 \cdot \sin delta\right) \cdot \cos theta\right)\right)}\]
  2. Using strategy rm

  3. Applied flip--0.2

    \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\left(\sin theta \cdot \sin delta\right) \cdot \cos \phi_1}{\color{blue}{\frac{\cos delta \cdot \cos delta - \left(\sin \phi_1 \cdot \sin \left(\sin^{-1} \left(\sin \phi_1 \cdot \cos delta + \left(\cos \phi_1 \cdot \sin delta\right) \cdot \cos theta\right)\right)\right) \cdot \left(\sin \phi_1 \cdot \sin \left(\sin^{-1} \left(\sin \phi_1 \cdot \cos delta + \left(\cos \phi_1 \cdot \sin delta\right) \cdot \cos theta\right)\right)\right)}{\cos delta + \sin \phi_1 \cdot \sin \left(\sin^{-1} \left(\sin \phi_1 \cdot \cos delta + \left(\cos \phi_1 \cdot \sin delta\right) \cdot \cos theta\right)\right)}}}\]

  4. Taylor expanded around inf 0.2

    \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\left(\sin theta \cdot \sin delta\right) \cdot \cos \phi_1}{\frac{\cos delta \cdot \cos delta - \left(\sin \phi_1 \cdot \sin \left(\sin^{-1} \left(\sin \phi_1 \cdot \cos delta + \left(\cos \phi_1 \cdot \sin delta\right) \cdot \cos theta\right)\right)\right) \cdot \left(\sin \phi_1 \cdot \sin \left(\sin^{-1} \left(\sin \phi_1 \cdot \cos delta + \left(\cos \phi_1 \cdot \sin delta\right) \cdot \cos theta\right)\right)\right)}{\cos delta + \sin \phi_1 \cdot \sin \left(\sin^{-1} \color{blue}{\left(\cos delta \cdot \sin \phi_1 + \sin delta \cdot \left(\cos theta \cdot \cos \phi_1\right)\right)}\right)}}\]

  5. Applied simplify0.2

    \[\leadsto \color{blue}{\lambda_1 + \tan^{-1}_* \frac{\left(\cos \phi_1 \cdot \sin delta\right) \cdot \sin theta}{\cos delta - \sin \left(\sin^{-1} \left(\left(\cos \phi_1 \cdot \sin delta\right) \cdot \cos theta + \sin \phi_1 \cdot \cos delta\right)\right) \cdot \sin \phi_1}}\]
  6. Using strategy rm

  7. Applied flip--0.2

    \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\left(\cos \phi_1 \cdot \sin delta\right) \cdot \sin theta}{\color{blue}{\frac{\cos delta \cdot \cos delta - \left(\sin \left(\sin^{-1} \left(\left(\cos \phi_1 \cdot \sin delta\right) \cdot \cos theta + \sin \phi_1 \cdot \cos delta\right)\right) \cdot \sin \phi_1\right) \cdot \left(\sin \left(\sin^{-1} \left(\left(\cos \phi_1 \cdot \sin delta\right) \cdot \cos theta + \sin \phi_1 \cdot \cos delta\right)\right) \cdot \sin \phi_1\right)}{\cos delta + \sin \left(\sin^{-1} \left(\left(\cos \phi_1 \cdot \sin delta\right) \cdot \cos theta + \sin \phi_1 \cdot \cos delta\right)\right) \cdot \sin \phi_1}}}\]

Runtime

Time bar (total: 1.5m)Debug logProfile

herbie shell --seed '#(1072743783 989954326 4239155542 3782239461 3602631542 1719177920)' 
(FPCore (lambda1 phi1 phi2 delta theta)
  :name "Destination given bearing on a great circle"
  (+ lambda1 (atan2 (* (* (sin theta) (sin delta)) (cos phi1)) (- (cos delta) (* (sin phi1) (sin (asin (+ (* (sin phi1) (cos delta)) (* (* (cos phi1) (sin delta)) (cos theta))))))))))