Add pal.polmo

cpal.pxd

@@ -161,6 +161,8 @@ cdef extern from "pal.h":
                  double *ra, double *dec, double *r, int *jstat )
     void palPlantu ( double date, double elong, double phi, double u[13],
                double *ra, double *dec, double *r, int *jstat )
+    void palPolmo( double elongm, double phim, double xp, double yp,
+                   double *elong, double *phi, double *daz )
     void palPm ( double r0, double d0, double pr, double pd,
              double px, double rv, double ep0, double ep1,
              double *r1, double *d1 )

pal.pyx.in

@@ -1489,6 +1489,20 @@ def pa( double ha, double dec, double phi):
      """
      return cpal.palPa( ha, dec, phi )
 
+def polmo( double elongm, double phim, double xp, double yp):
+    """
+    (elong, phi, daz) = polmo( elongm, phim, xp, yp )
+
+    @palPolmo@
+    """
+    cdef double elong
+    cdef double phi
+    cdef double daz
+
+    cpal.palPolmo( elongm, phim, xp, yp, &elong, &phi, &daz )
+
+    return (elong, phi, daz)
+
 def pertel(int jform, double date0, double date1,
             double epoch0, double orbi0, double anode0,
             double perih0, double aorq0, double e0, double am0):

setup.py

@@ -71,7 +71,7 @@
     "palNutc.c", "palOap.c", "palOapqk.c", "palObs.c",
     "palOne2One.c", "palPa.c", "palPertel.c", "palPertue.c",
     "palPlanel.c", "palPlanet.c", "palPlante.c", "palPlantu.c",
-    "palPm.c", "palPrebn.c", "palPrec.c", "palPreces.c",
+    "palPm.c", "palPolmo.c", "palPrebn.c", "palPrec.c", "palPreces.c",
     "palPrenut.c", "palPv2el.c", "palPv2ue.c", "palPvobs.c",
     "palRdplan.c", "palRefco.c", "palRefro.c", "palRefv.c",
     "palRefz.c", "palRverot.c", "palRvgalc.c", "palRvlg.c",

test_pal.py

@@ -895,6 +895,13 @@ def test_pmVector(self):
             self.assertAlmostEqual(r1,r2,12)
             self.assertAlmostEqual(d1,d2,12)
 
+    def test_polmo(self):
+        (elong, phi, daz) = pal.polmo( 0.7, -0.5, 1.0e-6, -2.0e-6)
+
+        self.assertAlmostEqual( elong, 0.7000004837322044, 12 )
+        self.assertAlmostEqual( phi, -0.4999979467222241, 12 )
+        self.assertAlmostEqual( daz, 1.008982781275728e-6, 12 )
+
     def test_prebn(self):
         expected = np.array( [
             [ 9.999257613786738e-1, -1.117444640880939e-2, -4.858341150654265e-3 ],