[Octave-patch-tracker] [patch #10022] [octave forge] (mapping) Auxillary Latitudes

[Philip Nienhuis]

Follow-up Comment #10, patch #10022 (project octave):

I pushed rect2geo here:
http://hg.code.sf.net/p/octave/mapping//rev/5d625af9eea8

I edited style (Texinfo & code) a bit.
Note that your rect2geo.m only processed scalars or row vectors properly,
other shapes resulted in scalar (and incorrect) output. Where you had this:

    A = [mu;
        ((3 / 2) * n - (27 / 32) * n3) * sin(2 * mu);
        ((21 / 16) * n2 - (55 / 32) * n4) * sin(4 * mu);
        ((151 / 48) * n3) * sin(6 * mu);
        ((1097 / 512) * n4) * sin(8 * mu)];
    phi = sum (A); 


... I changed it into:

  phi = mu .+ ...
        (((3 / 2) * n - (27 / 32) * n3) * sin (2 * mu)) .+ ...
        (((21 / 16) * n2 - (55 / 32) * n4) * sin (4 * mu)) .+ ...
        (((151 / 48) * n3) * sin (6 * mu)) .+ ...
        (((1097 / 512) * n4) * sin (8 * mu));


... that way rect2geo will accept about any shape of input.

Thanks so far.

    _______________________________________________________

Reply to this item at:

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

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