Groups | Search | Server Info | Keyboard shortcuts | Login | Register [http] [https] [nntp] [nntps]

Groups > comp.programming > #16274

Re: Another little puzzle

Show all headers | View raw

On 2023-01-08 21:41, Ben Bacarisse wrote:

> But
> the problem is /finding/ a specific average -- the point (or angle) that
> minimises the sum of squares of the distances (or angles) from that
> average point (or angle).

That is not a programming problem. It is a problem space's one. You must 
go to the corresponding part of science or engineering or economics etc 
and formulate it there in terms of that space.

Average is just such formulation of finding a representative member from 
some set having some specific properties. E.g. least squares in physics 
usually has the meaning of least energy etc. So, it goes in this order 
(not in reverse):

1. You need the least energy representative.
2. An average gives you that.
3. You measure the inputs.

After that you ask yourself given these measures how to compute the 
average? And only now it becomes a programming issue.

> The fact that it makes no odds (as everyone knows) whether we consider
> angles (often called central angles in this context) or great circle
> distances is not the issue.  It's finding the average that minimises the
> sum of squares of differences that's the issue.

Minimum of squares (Euclidean norm) is represented by the mathematical mean.

> You say you need a formula, so I'll try...  Let P_n be a collection of n
> unit vectors specifying n points on a unit sphere.  Find the unit vector
> A that minimises
> 
>    Sum_{i=1,n} ( arctan( |A x P_n|  /  A . P_n ) )^2
> 
> (arctan is the "all quadrant" version that is often called atan2 in
> programming languages.)

1. atan2 has two arguments (the adjacent and the opposite legs);
2. What is "A";
3. What operations(?) "x" and "." denote?

>> Programs are written according to the specifications. Numeric programs
>> require a properly stated problem, rather than a bunch of words
>> arbitrarily thrown in a meaningless sentence as above.
> 
> Given the context, I think that's a very biased characterisation of
> what's been said here.
> 
> My first job was as a "numerical analyst", and the very first program I
> was employed to write was for a professor of statistics.  It was to
> calculate a novel kind of fit line.  The specification was just a few
> sentences.  No formulas.  It was perfectly clear, and I could get the
> job done.  I don't think this is unusual.  Words are often enough, and
> they can avoid undue over specification.  For example, the problem in
> question is essentially the same if the points are given by latitude and
> longitude on a non-unit sphere, but the formula would look very
> different.

This just means that you formulated the specifications yourself, being 
able to read your professor's mind. Not everybody's mind is easy read or 
worth reading... (:-))

-- 
Regards,
Dmitry A. Kazakov
http://www.dmitry-kazakov.de

Back to comp.programming | Previous | Next — Previous in thread | Next in thread | Find similar

Thread

Re: Another little puzzle Ben Bacarisse <ben.usenet@bsb.me.uk> - 2022-12-30 01:00 +0000
  Re: Another little puzzle Tim Rentsch <tr.17687@z991.linuxsc.com> - 2022-12-29 23:25 -0800
    Re: Another little puzzle Ben Bacarisse <ben.usenet@bsb.me.uk> - 2022-12-30 14:04 +0000
      Re: Another little puzzle Ben Bacarisse <ben.usenet@bsb.me.uk> - 2022-12-31 00:24 +0000
        Re: Another little puzzle Tim Rentsch <tr.17687@z991.linuxsc.com> - 2022-12-31 06:42 -0800
          Re: Another little puzzle "Dmitry A. Kazakov" <mailbox@dmitry-kazakov.de> - 2022-12-31 17:04 +0100
            Re: Another little puzzle Ben Bacarisse <ben.usenet@bsb.me.uk> - 2023-01-01 01:24 +0000
            Re: Another little puzzle Tim Rentsch <tr.17687@z991.linuxsc.com> - 2023-01-08 07:45 -0800
              Re: Another little puzzle "Dmitry A. Kazakov" <mailbox@dmitry-kazakov.de> - 2023-01-08 17:17 +0100
                Re: Another little puzzle Ben Bacarisse <ben.usenet@bsb.me.uk> - 2023-01-08 20:41 +0000
                Re: Another little puzzle Richard Heathfield <rjh@cpax.org.uk> - 2023-01-08 21:14 +0000
                Re: Another little puzzle "Dmitry A. Kazakov" <mailbox@dmitry-kazakov.de> - 2023-01-08 22:31 +0100
                Re: Another little puzzle Ben Bacarisse <ben.usenet@bsb.me.uk> - 2023-01-09 03:25 +0000
                Re: Another little puzzle "Dmitry A. Kazakov" <mailbox@dmitry-kazakov.de> - 2023-01-09 11:22 +0100
                Re: Another little puzzle Ben Bacarisse <ben.usenet@bsb.me.uk> - 2023-01-09 20:37 +0000
                Re: Another little puzzle "Dmitry A. Kazakov" <mailbox@dmitry-kazakov.de> - 2023-01-10 09:06 +0100
                Re: Another little puzzle Ben Bacarisse <ben.usenet@bsb.me.uk> - 2023-01-11 02:41 +0000
                Re: Another little puzzle "Dmitry A. Kazakov" <mailbox@dmitry-kazakov.de> - 2023-01-11 10:01 +0100
                Re: Another little puzzle Ben Bacarisse <ben.usenet@bsb.me.uk> - 2023-01-12 01:00 +0000
                Re: Another little puzzle Tim Rentsch <tr.17687@z991.linuxsc.com> - 2023-01-10 23:36 -0800
                Re: Another little puzzle Y A <ya00000100000@yahoo.com> - 2023-01-11 02:39 -0800
          Re: Another little puzzle Ben Bacarisse <ben.usenet@bsb.me.uk> - 2023-01-01 01:10 +0000
            Re: Another little puzzle Tim Rentsch <tr.17687@z991.linuxsc.com> - 2023-01-08 07:17 -0800
              Re: Another little puzzle Ben Bacarisse <ben.usenet@bsb.me.uk> - 2023-01-08 19:43 +0000
                Re: Another little puzzle Ben Bacarisse <ben.usenet@bsb.me.uk> - 2023-01-08 19:59 +0000
                Re: Another little puzzle Tim Rentsch <tr.17687@z991.linuxsc.com> - 2023-01-10 23:25 -0800
      Re: Another little puzzle Tim Rentsch <tr.17687@z991.linuxsc.com> - 2022-12-31 06:20 -0800
        Re: Another little puzzle Augǝl <angel0000000001000000000000@mail.ee> - 2022-12-31 10:23 -0800

csiph-web