[Octave-patch-tracker] [patch #9997] [octave forge](mapping) areaquad

[Philip Nienhuis]

Follow-up Comment #2, patch #9997 (project octave):

Other things:

1. Unit Sphere
--------------
Your areaquad is based on the wgs84 ellipsoid. Matlab's areaquad.m is based on
a unit sphere (radius = 1).

Since bug #59444 referenceElipsoid.m also has 'Unit Sphere' and it is even the
default, as in Matlab.

In areaquad, entering "Unit Sphere" for the ellipsoid gives NaN as output
because in L.13 'f = 1 ./ (2 .* e)' (where e = E.Eccentricity), variable f
goes to infinity.
In fact, areaquad.m gives NaN for any ellipsoid whose eccentricity 
equals 0 (i.e., is a sphere).

I suppose the easiest (but kludgy) solution is to check for eccentricity < eps
and if that's the case, use another formula.

2. Ellipsoid input
------------------
Matlab's areaquad also allows a vector of the form [semimajor_axis
eccentricity] for the ellipsoid.

Could you please look into these two things?
I'll take care of the other little issues I mentioned in comment #1. In fact I
have already done most of it.

Apologies that I seem to be such a nitpicker with most of your contributions
:-)

    _______________________________________________________

Reply to this item at:

  <https://savannah.gnu.org/patch/?9997>

_______________________________________________
  Message sent via Savannah
  https://savannah.gnu.org/