public inbox for
 help / color / mirror / Atom feed
From: David Zaslavsky <>
Subject: Re: interp2d discussion
Date: Fri, 10 Jan 2014 06:05:00 -0000	[thread overview]
Message-ID: <> (raw)
In-Reply-To: <>

[-- Attachment #1: Type: text/plain, Size: 2314 bytes --]

Hi Patrick,

On 01/06/2014 11:57 AM, Patrick Alken wrote:
> On 01/05/2014 01:31 PM, David Zaslavsky wrote:
>> On the indexing issue: interp2d treats the x index as the column
>> coordinate and the y index as the row coordinate, so that xsize is the
>> number of columns and ysize is the number of rows. This is done for
>> compatibility with Cartesian coordinates: most people envision x as the
>> horizontal coordinate and y as the vertical one. Given this
>> interpretation, the current implementation of INDEX_2D does use
>> row-major order.
> This makes sense, we can probably leave it the way it is, we'll just
> need to make sure the documentation is very clear on how to construct
> the z array so users aren't confused with the indexing.
Sure, that would be a reasonable option, but I'd like to do a bit of
research (i.e. ask a few more people for opinions) to see what solution
will be the least confusing. After all, if there is a change to be made,
this is the time to do it, not later after the code is incorporated into
GSL. I think backward compatibility will not be too much of a problem
because people will have to make other changes anyway to transition from
the interp2d_* functions to the corresponding gsl_* functions.
>> Also, I should mention that I don't consider interp2d to be
>> feature-complete, nor am I in any particular rush to make it complete.
>> There are some features I'm interested in adding, including additional
>> interpolation schemes, more interpolation methods without boundary
>> checks (for extrapolation), etc. How would this move affect my ability
>> to continue to add features?
> It already has bilinear and bicubic implemented. This is more than
> enough to make it useful to people so I think we can put it in now and
> continue adding features.
> I saw from some comments that you used some code from ALGLIB - can we
> assume you're using the GPL v2 version of ALGLIB? In that case there
> shouldn't be any problem putting it in GSL under GPL v3.
The code was adapted from ALGLIB by Thomas Beutlich, and I'm pretty sure
he told me that he used the GPL v2 version. I'll double-check that, though.

When I have some time I'll fork the GSL repository and work on adding
the interp2d code in a new branch.

:) David

[-- Attachment #2: OpenPGP digital signature --]
[-- Type: application/pgp-signature, Size: 295 bytes --]

  reply	other threads:[~2014-01-10  6:05 UTC|newest]

Thread overview: 11+ messages / expand[flat|nested]  mbox.gz  Atom feed  top
     [not found] <>
2014-01-03 19:00 ` [Help-gsl] Getting an extension listed on GSL page Patrick Alken
2014-01-05 20:31   ` David Zaslavsky
2014-01-06 16:58     ` interp2d discussion Patrick Alken
2014-01-10  6:05       ` David Zaslavsky [this message]
2014-01-10  6:11       ` David Zaslavsky
2014-01-10 16:37         ` Patrick Alken
2014-01-11  0:28           ` jeremy theler
2014-01-11  1:55             ` Patrick Alken
2014-01-14  1:51               ` jeremy theler
2014-01-14  2:23                 ` Patrick Alken
2014-01-11  1:58             ` David Zaslavsky

Reply instructions:

You may reply publicly to this message via plain-text email
using any one of the following methods:

* Save the following mbox file, import it into your mail client,
  and reply-to-all from there: mbox

  Avoid top-posting and favor interleaved quoting:

* Reply using the --to, --cc, and --in-reply-to
  switches of git-send-email(1):

  git send-email \ \ \ \

* If your mail client supports setting the In-Reply-To header
  via mailto: links, try the mailto: link
Be sure your reply has a Subject: header at the top and a blank line before the message body.
This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).