Changes between Version 5 and Version 6 of Ticket #332, comment 2
- Timestamp:
- Jan 7, 2016 10:44:28 AM (9 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
Ticket #332, comment 2
v5 v6 11 11 Perhaps the confusion arose from the different forms of the ellipsoid of revolution. Eq 2.32 gives the ellipsoid of revolution as: 12 12 13 i(h) = ∫,,0,,^pi/2^ Phi^2^(h a √(cos 2 θ + v2 sin2θ) cos θ dθ13 i(h) = ∫,,0,,^pi/2^ Phi^2^(h a √(cos^2^ θ + v^2^ sin^2^ θ) cos θ dθ 14 14 15 15 but the 1-D ellipsoid kernel employs the substitution u = sin theta so it is evaluated as 16 16 17 i(h) = ∫,,0,,^1^ Phi^2^(h a √(1 + u 2(v2- 1)) du17 i(h) = ∫,,0,,^1^ Phi^2^(h a √(1 + u^2^(v^2^ - 1)) du 18 18 19 19 [Note: the u variable above is incorrectly named as cos_alpha in the code.]