Merge pull request #43 from cneubauerUSGS/m05245

      // Set dalpha to be half the grid spacing for nyquist sampling

      //double dalpha = (PI/180.0)/(2.0*p_demScale);

      double cmin = cos((90.0 - 1.0 / (2.0 * demScale())) * DEG2RAD);

      double dalpha = MAX(cos(g1lat * DEG2RAD), cmin) / (2.0 * demScale() * DEG2RAD);

      double dalpha = MAX(cos(g1lat * DEG2RAD), cmin) / (2.0 * demScale() * RAD2DEG);

      // Previous Sensor version used local version of this method with lat and lon doubles.  Steven said

      // it didn't make a significant difference in speed.

        // put in a test (above) for dd being smaller than the pixel

        // convergence tolerance.  If so the loop exits without an

        // intersection

        dalpha = MAX(cos(g2lat * DEG2RAD), cmin) / (2.0 * demScale() * DEG2RAD);

        dalpha = MAX(cos(g2lat * DEG2RAD), cmin) / (2.0 * demScale() * RAD2DEG);

      } // end while

      SpiceDouble intersectionPoint[3];

  }

  // Do nothing since the DEM intersection was already successful

}

   *                           intersection in intersectSurface() method to prevent early return

   *                           and attempt the iterative method even when the ellipsoid is not

   *                           intersected. Fixes #1438

   *   @history 2018-01-05 Cole Neubauer - Fixed units conversion in intersectSurface so that the

   *                           loop is stepping by radians per pixel, as recommended by Jeff

   *                           Anderson (LROC team). Fixes #5245

   */

  class EquatorialCylindricalShape : public DemShape {

    public:

};

#endif