Return-Path: X-Original-To: apmail-commons-dev-archive@www.apache.org Delivered-To: apmail-commons-dev-archive@www.apache.org Received: from mail.apache.org (hermes.apache.org [140.211.11.3]) by minotaur.apache.org (Postfix) with SMTP id B3399C287 for ; Mon, 21 May 2012 05:32:39 +0000 (UTC) Received: (qmail 42134 invoked by uid 500); 21 May 2012 05:32:38 -0000 Delivered-To: apmail-commons-dev-archive@commons.apache.org Received: (qmail 42038 invoked by uid 500); 21 May 2012 05:32:37 -0000 Mailing-List: contact dev-help@commons.apache.org; run by ezmlm Precedence: bulk List-Help: List-Unsubscribe: List-Post: List-Id: Reply-To: "Commons Developers List" Delivered-To: mailing list dev@commons.apache.org Received: (qmail 42016 invoked by uid 99); 21 May 2012 05:32:37 -0000 Received: from athena.apache.org (HELO athena.apache.org) (140.211.11.136) by apache.org (qpsmtpd/0.29) with ESMTP; Mon, 21 May 2012 05:32:37 +0000 X-ASF-Spam-Status: No, hits=-2.3 required=5.0 tests=RCVD_IN_DNSWL_MED,SPF_HELO_PASS,SPF_PASS X-Spam-Check-By: apache.org Received-SPF: pass (athena.apache.org: domain of SRS0=Evch=DZ=m4x.org=sebastien.brisard@bounces.m4x.org designates 129.104.30.34 as permitted sender) Received: from [129.104.30.34] (HELO mx1.polytechnique.org) (129.104.30.34) by apache.org (qpsmtpd/0.29) with ESMTP; Mon, 21 May 2012 05:32:30 +0000 Received: from mail-gg0-f171.google.com (mail-gg0-f171.google.com [209.85.161.171]) (using TLSv1 with cipher RC4-MD5 (128/128 bits)) (No client certificate requested) by ssl.polytechnique.org (Postfix) with ESMTPSA id 1CA91140000B5 for ; Mon, 21 May 2012 07:32:08 +0200 (CEST) Received: by ggmi1 with SMTP id i1so5085330ggm.30 for ; Sun, 20 May 2012 22:32:07 -0700 (PDT) MIME-Version: 1.0 Received: by 10.50.207.38 with SMTP id lt6mr5710608igc.34.1337578326764; Sun, 20 May 2012 22:32:06 -0700 (PDT) Received: by 10.64.72.167 with HTTP; Sun, 20 May 2012 22:32:06 -0700 (PDT) In-Reply-To: <4FB96D62.3080308@gmail.com> References: <4FB96D62.3080308@gmail.com> Date: Mon, 21 May 2012 07:32:06 +0200 Message-ID: Subject: Re: [math] numerical stability problems of continued fraction From: =?ISO-8859-1?Q?S=E9bastien_Brisard?= To: Commons Developers List Content-Type: text/plain; charset=ISO-8859-1 Content-Transfer-Encoding: quoted-printable X-AV-Checked: ClamAV using ClamSMTP at svoboda.polytechnique.org (Mon May 21 07:32:08 2012 +0200 (CEST)) X-Org-Mail: sebastien.brisard.1997@polytechnique.org X-Virus-Checked: Checked by ClamAV on apache.org X-Old-Spam-Flag: No, tests=bogofilter, spamicity=0.243106, queueID=5D19814000360 Hi Thomas, > > This seems to solve all the described problems, but the question remains > if this can be included without problems into CM. The algorithm > description in NR is not in source code but on an higher level. > Looking at the relevant chapter in NR, I don't think there is any license issue (no code provided in NR). However, if you really want to play it safe, the paper by Thompson is available here http://www.fresco.org.uk/papers/Thompson-JCP64p490.pdf I guess downloading it is OK, as www.fresco.org.uk seems to be Thompson's own website. What you could do is conform closely to Thompson's notation, instead of NR's notation. S=E9bastien --------------------------------------------------------------------- To unsubscribe, e-mail: dev-unsubscribe@commons.apache.org For additional commands, e-mail: dev-help@commons.apache.org