From mboxrd@z Thu Jan  1 00:00:00 1970
Return-Path: <gsl-discuss-return-5208-listarch-gsl-discuss=sources.redhat.com@sourceware.org>
Received: (qmail 12855 invoked by alias); 22 Nov 2007 11:33:51 -0000
Received: (qmail 12843 invoked by uid 22791); 22 Nov 2007 11:33:50 -0000
X-Spam-Check-By: sourceware.org
Received: from nf-out-0910.google.com (HELO nf-out-0910.google.com) (64.233.182.191)
    by sourceware.org (qpsmtpd/0.31) with ESMTP; Thu, 22 Nov 2007 11:33:44 +0000
Received: by nf-out-0910.google.com with SMTP id b11so2488651nfh
        for <gsl-discuss@sourceware.org>; Thu, 22 Nov 2007 03:33:41 -0800 (PST)
Received: by 10.86.99.9 with SMTP id w9mr8350081fgb.1195731220877;
        Thu, 22 Nov 2007 03:33:40 -0800 (PST)
Received: from 192.168.0.111 ( [194.44.31.118])
        by mx.google.com with ESMTPS id 12sm940195fgg.2007.11.22.03.33.38
        (version=SSLv3 cipher=RC4-MD5);
        Thu, 22 Nov 2007 03:33:39 -0800 (PST)
Subject: Not clear to understand :)
From: Serhiy Lisovenko <lisovenko.s@gmail.com>
Reply-To: lisovenko.s@gmail.com
To: gsl-discuss@sourceware.org
Content-Type: text/plain
Date: Thu, 22 Nov 2007 11:33:00 -0000
Message-Id: <1195731261.5477.14.camel@localhost.localdomain>
Mime-Version: 1.0
X-Mailer: Evolution 2.0.2 (2.0.2-3) 
Content-Transfer-Encoding: 7bit
Mailing-List: contact gsl-discuss-help@sourceware.org; run by ezmlm
Precedence: bulk
List-Id: <gsl-discuss.sourceware.org>
List-Subscribe: <mailto:gsl-discuss-subscribe@sourceware.org>
List-Archive: <http://sourceware.org/ml/gsl-discuss/>
List-Post: <mailto:gsl-discuss@sourceware.org>
List-Help: <mailto:gsl-discuss-help@sourceware.org>, <http://sourceware.org/ml/#faqs>
Sender: gsl-discuss-owner@sourceware.org
X-SW-Source: 2007-q4/txt/msg00035.txt.bz2

in file fit/linear.c

  for (i = 0; i < n; i++)
    {
      m_x += (x[i * xstride] - m_x) / (i + 1.0);
      m_y += (y[i * ystride] - m_y) / (i + 1.0);
    }

and other same...

code equivalent to

  for (i = 0; i < n; i++)
    {
      m_x +=x[i * xstride];
      m_y +=y[i * ystride];
    }
m_x/=(double)n;
m_y/=(double)n;

But the second code is more clear and some faster (no extra divisions).