Average Error: 0.2 → 0.2
Time: 19.1s
Precision: binary64
\[\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{\sin theta \cdot \left(\sin delta \cdot \cos \phi_1\right)}{\frac{\left(\cos delta \cdot \cos delta\right) \cdot \left(\cos delta \cdot \cos delta\right) - \left(\sin \left(\sin^{-1} \left(\cos delta \cdot \sin \phi_1 + \sin delta \cdot \left(\cos \phi_1 \cdot \cos theta\right)\right)\right) \cdot \left(\sin \left(\sin^{-1} \left(\cos delta \cdot \sin \phi_1 + \sin delta \cdot \left(\cos \phi_1 \cdot \cos theta\right)\right)\right) \cdot {\left(\sin \phi_1\right)}^{2}\right)\right) \cdot \left(\sin \left(\sin^{-1} \left(\cos delta \cdot \sin \phi_1 + \sin delta \cdot \left(\cos \phi_1 \cdot \cos theta\right)\right)\right) \cdot \left(\sin \left(\sin^{-1} \left(\cos delta \cdot \sin \phi_1 + \sin delta \cdot \left(\cos \phi_1 \cdot \cos theta\right)\right)\right) \cdot {\left(\sin \phi_1\right)}^{2}\right)\right)}{\left(\cos delta + \sin \phi_1 \cdot \sin \left(\sin^{-1} \left(\cos delta \cdot \sin \phi_1 + \sin delta \cdot \left(\cos \phi_1 \cdot \cos theta\right)\right)\right)\right) \cdot \left(\sin \left(\sin^{-1} \left(\cos delta \cdot \sin \phi_1 + \sin delta \cdot \left(\cos \phi_1 \cdot \cos theta\right)\right)\right) \cdot \left(\sin \left(\sin^{-1} \left(\cos delta \cdot \sin \phi_1 + \sin delta \cdot \left(\cos \phi_1 \cdot \cos theta\right)\right)\right) \cdot {\left(\sin \phi_1\right)}^{2}\right) + {\left(\cos delta\right)}^{2}\right)}}\]
\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{\sin theta \cdot \left(\sin delta \cdot \cos \phi_1\right)}{\frac{\left(\cos delta \cdot \cos delta\right) \cdot \left(\cos delta \cdot \cos delta\right) - \left(\sin \left(\sin^{-1} \left(\cos delta \cdot \sin \phi_1 + \sin delta \cdot \left(\cos \phi_1 \cdot \cos theta\right)\right)\right) \cdot \left(\sin \left(\sin^{-1} \left(\cos delta \cdot \sin \phi_1 + \sin delta \cdot \left(\cos \phi_1 \cdot \cos theta\right)\right)\right) \cdot {\left(\sin \phi_1\right)}^{2}\right)\right) \cdot \left(\sin \left(\sin^{-1} \left(\cos delta \cdot \sin \phi_1 + \sin delta \cdot \left(\cos \phi_1 \cdot \cos theta\right)\right)\right) \cdot \left(\sin \left(\sin^{-1} \left(\cos delta \cdot \sin \phi_1 + \sin delta \cdot \left(\cos \phi_1 \cdot \cos theta\right)\right)\right) \cdot {\left(\sin \phi_1\right)}^{2}\right)\right)}{\left(\cos delta + \sin \phi_1 \cdot \sin \left(\sin^{-1} \left(\cos delta \cdot \sin \phi_1 + \sin delta \cdot \left(\cos \phi_1 \cdot \cos theta\right)\right)\right)\right) \cdot \left(\sin \left(\sin^{-1} \left(\cos delta \cdot \sin \phi_1 + \sin delta \cdot \left(\cos \phi_1 \cdot \cos theta\right)\right)\right) \cdot \left(\sin \left(\sin^{-1} \left(\cos delta \cdot \sin \phi_1 + \sin delta \cdot \left(\cos \phi_1 \cdot \cos theta\right)\right)\right) \cdot {\left(\sin \phi_1\right)}^{2}\right) + {\left(\cos delta\right)}^{2}\right)}}
(FPCore (lambda1 phi1 phi2 delta theta)
 :precision binary64
 (+
  lambda1
  (atan2
   (* (* (sin theta) (sin delta)) (cos phi1))
   (-
    (cos delta)
    (*
     (sin phi1)
     (sin
      (asin
       (+
        (* (sin phi1) (cos delta))
        (* (* (cos phi1) (sin delta)) (cos theta))))))))))
↓
(FPCore (lambda1 phi1 phi2 delta theta)
 :precision binary64
 (+
  lambda1
  (atan2
   (* (sin theta) (* (sin delta) (cos phi1)))
   (/
    (-
     (* (* (cos delta) (cos delta)) (* (cos delta) (cos delta)))
     (*
      (*
       (sin
        (asin
         (+
          (* (cos delta) (sin phi1))
          (* (sin delta) (* (cos phi1) (cos theta))))))
       (*
        (sin
         (asin
          (+
           (* (cos delta) (sin phi1))
           (* (sin delta) (* (cos phi1) (cos theta))))))
        (pow (sin phi1) 2.0)))
      (*
       (sin
        (asin
         (+
          (* (cos delta) (sin phi1))
          (* (sin delta) (* (cos phi1) (cos theta))))))
       (*
        (sin
         (asin
          (+
           (* (cos delta) (sin phi1))
           (* (sin delta) (* (cos phi1) (cos theta))))))
        (pow (sin phi1) 2.0)))))
    (*
     (+
      (cos delta)
      (*
       (sin phi1)
       (sin
        (asin
         (+
          (* (cos delta) (sin phi1))
          (* (sin delta) (* (cos phi1) (cos theta))))))))
     (+
      (*
       (sin
        (asin
         (+
          (* (cos delta) (sin phi1))
          (* (sin delta) (* (cos phi1) (cos theta))))))
       (*
        (sin
         (asin
          (+
           (* (cos delta) (sin phi1))
           (* (sin delta) (* (cos phi1) (cos theta))))))
        (pow (sin phi1) 2.0)))
      (pow (cos delta) 2.0)))))))
double code(double lambda1, double phi1, double phi2, double delta, double theta) {
	return ((double) (lambda1 + ((double) atan2(((double) (((double) (((double) sin(theta)) * ((double) sin(delta)))) * ((double) cos(phi1)))), ((double) (((double) cos(delta)) - ((double) (((double) sin(phi1)) * ((double) sin(((double) asin(((double) (((double) (((double) sin(phi1)) * ((double) cos(delta)))) + ((double) (((double) (((double) cos(phi1)) * ((double) sin(delta)))) * ((double) cos(theta))))))))))))))))));
}
↓
double code(double lambda1, double phi1, double phi2, double delta, double theta) {
	return ((double) (lambda1 + ((double) atan2(((double) (((double) sin(theta)) * ((double) (((double) sin(delta)) * ((double) cos(phi1)))))), (((double) (((double) (((double) (((double) cos(delta)) * ((double) cos(delta)))) * ((double) (((double) cos(delta)) * ((double) cos(delta)))))) - ((double) (((double) (((double) sin(((double) asin(((double) (((double) (((double) cos(delta)) * ((double) sin(phi1)))) + ((double) (((double) sin(delta)) * ((double) (((double) cos(phi1)) * ((double) cos(theta)))))))))))) * ((double) (((double) sin(((double) asin(((double) (((double) (((double) cos(delta)) * ((double) sin(phi1)))) + ((double) (((double) sin(delta)) * ((double) (((double) cos(phi1)) * ((double) cos(theta)))))))))))) * ((double) pow(((double) sin(phi1)), 2.0)))))) * ((double) (((double) sin(((double) asin(((double) (((double) (((double) cos(delta)) * ((double) sin(phi1)))) + ((double) (((double) sin(delta)) * ((double) (((double) cos(phi1)) * ((double) cos(theta)))))))))))) * ((double) (((double) sin(((double) asin(((double) (((double) (((double) cos(delta)) * ((double) sin(phi1)))) + ((double) (((double) sin(delta)) * ((double) (((double) cos(phi1)) * ((double) cos(theta)))))))))))) * ((double) pow(((double) sin(phi1)), 2.0)))))))))) / ((double) (((double) (((double) cos(delta)) + ((double) (((double) sin(phi1)) * ((double) sin(((double) asin(((double) (((double) (((double) cos(delta)) * ((double) sin(phi1)))) + ((double) (((double) sin(delta)) * ((double) (((double) cos(phi1)) * ((double) cos(theta)))))))))))))))) * ((double) (((double) (((double) sin(((double) asin(((double) (((double) (((double) cos(delta)) * ((double) sin(phi1)))) + ((double) (((double) sin(delta)) * ((double) (((double) cos(phi1)) * ((double) cos(theta)))))))))))) * ((double) (((double) sin(((double) asin(((double) (((double) (((double) cos(delta)) * ((double) sin(phi1)))) + ((double) (((double) sin(delta)) * ((double) (((double) cos(phi1)) * ((double) cos(theta)))))))))))) * ((double) pow(((double) sin(phi1)), 2.0)))))) + ((double) pow(((double) cos(delta)), 2.0)))))))))));
}

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 Error: 0.2 bits

    \[\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. SimplifiedError: 0.2 bits

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

  4. Applied flip--Error: 0.2 bits

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

  5. SimplifiedError: 0.2 bits

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

  7. Applied flip--Error: 0.2 bits

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

  8. Applied associate-/l/Error: 0.2 bits

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

  9. SimplifiedError: 0.2 bits

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

  10. Final simplificationError: 0.2 bits

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

Reproduce

herbie shell --seed 2020203 
(FPCore (lambda1 phi1 phi2 delta theta)
  :name "Destination given bearing on a great circle"
  :precision binary64
  (+ lambda1 (atan2 (* (* (sin theta) (sin delta)) (cos phi1)) (- (cos delta) (* (sin phi1) (sin (asin (+ (* (sin phi1) (cos delta)) (* (* (cos phi1) (sin delta)) (cos theta))))))))))