#553 closed defect (fixed)
psSpherePrecess - inaccuracies
| Reported by: | Owned by: | Paul Price | |
|---|---|---|---|
| Priority: | high | Milestone: | |
| Component: | PSLib ADD | Version: | unspecified |
| Severity: | normal | Keywords: | |
| Cc: |
Description
In testing psSpherePrecess I've noticed some inaccuracy and I'm not sure where
it's coming from. The following is the output w.r.t. the suggested test cases:
2005:10:15 00:58:59Z|localhost.localdomain|E|testSphereRotPrecess
(tst_psSphereOps.c:341)
Precess r = 1.95818e-05 not equal to expected = 6.23845
2005:10:15 00:58:59Z|localhost.localdomain|E|testSphereRotPrecess
(tst_psSphereOps.c:367)
Precess r = -0.0223402 not equal to expected = 6.26083
2005:10:15 00:58:59Z|localhost.localdomain|E|testSphereRotPrecess
(tst_psSphereOps.c:393)
Precess r = -3.14132 not equal to expected = 3.09662
If you notice, these values are not completely wrong. The first one for
example, has r =~ 0 and expected ~2*pi. I will look into this myself as well,
but I'm hoping that someone has insight into this problem for me. This problem
could be related to, or the same as the previous bug i filed concerning the part
about Ecliptic to ICRS transformation. (bug 552)
Change History (2)
comment:1 by , 21 years ago
| Status: | new → assigned |
|---|
comment:2 by , 21 years ago
| Resolution: | → fixed |
|---|---|
| Status: | assigned → closed |
I updated the functions so that negative values get changed to positive values
by adding 2*pi to the negative number. The values matched after using the new
inversion formula.

Yes, I think it is the same problem as bug 552: the Euler angles aren't being
used correctly. Take a look at the Perl script in bug 461.