| 11 | } |
| 12 | |
| 13 | void geod_for(void) { |
| 14 | double s12 = geod_S, lat2, lon2, azi2; |
| 15 | geod_position(&GlobalGeodesicLine, s12, &lat2, &lon2, &azi2); |
| 16 | azi2 += azi2 >= 0 ? -180 : 180; /* Compute back azimuth */ |
| 17 | phi2 = lat2 * DEG_TO_RAD; |
| 18 | lam2 = lon2 * DEG_TO_RAD; |
| 19 | al21 = azi2 * DEG_TO_RAD; |
| 20 | } |
| 21 | |
| 22 | void geod_inv(void) { |
| 23 | double lat1 = phi1 / DEG_TO_RAD, lon1 = lam1 / DEG_TO_RAD, |
no test coverage detected