Path: csiph.com!x330-a1.tempe.blueboxinc.net!usenet.pasdenom.info!gegeweb.org!de-l.enfer-du-nord.net!feeder2.enfer-du-nord.net!feeder.news-service.com!newsfeed.xs4all.nl!newsfeed6.news.xs4all.nl!xs4all!post.news.xs4all.nl!not-for-mail Return-Path: X-Original-To: python-list@python.org Delivered-To: python-list@mail.python.org X-Spam-Status: OK 0.004 X-Spam-Evidence: '*H*': 0.99; '*S*': 0.00; 'subject:Python': 0.04; 'linear': 0.07; 'python)': 0.07; 'terry': 0.07; '>>>': 0.12; 'am,': 0.14; 'wrote:': 0.14; 'already.': 0.16; 'instances.': 0.16; 'reedy': 0.16; 'solver,': 0.16; 'terms).': 0.16; 'tue,': 0.20; 'cc:no real name:2**0': 0.20; 'cc:2**0': 0.20; 'header:In-Reply- To:1': 0.22; 'cc:addr:python-list': 0.22; '(and': 0.22; '(where': 0.23; 'gregory': 0.23; 'replacing': 0.23; '(in': 0.27; 'message- id:@mail.gmail.com': 0.28; 'variables': 0.29; 'cc:addr:python.org': 0.31; 'ewing': 0.31; 'separately': 0.31; 'trick': 0.31; 'actual': 0.31; 'fairly': 0.33; 'constant': 0.35; 'received:209.85.216.46': 0.35; 'received:mail- qw0-f46.google.com': 0.35; 'usually': 0.36; 'received:209.85': 0.37; 'received:google.com': 0.38; 'no.': 0.38; 'used': 0.38; 'explain': 0.39; 'could': 0.39; 'received:209': 0.39; 'solution': 0.40; 'would': 0.40; 'header:Received:5': 0.40; 'simple': 0.60; 'hand': 0.61; '2011': 0.62; 'presented': 0.63; 'formula': 0.84 DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:mime-version:in-reply-to:references:date :message-id:subject:from:to:cc:content-type; bh=jUnNYTDcHaYY3AMXv5f/Is4ZBZzE5tyVVQw5VEDTfjA=; b=uxJvBbsYWvByi2CxQ9JsMKVaCR0joKod6jwAFjcRg8nSsgc5GGdnDVNT7XO17M5rm/ 2gUHYndl5AjWqJ0rMVwjtHONtCyRMvY6467D4iJ8u9vL6uhlbZ2NSCABErHaJLubQ3uV z7JvyIFwRClI9ya43nep04Rb/zqZw0jVkaGdU= DomainKey-Signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=mime-version:in-reply-to:references:date:message-id:subject:from:to :cc:content-type; b=fXPuWuDjWw6R1d/opn8DdAOQaTbOaTJoHr5xwVDKg4APXlXOWu8ARwGtco2ZeeMpcS +e8oa9Avv6FT/ALqY3B7xrESXQZR+WSTL15B4RB0W8qDz0DHTBu1xhio4qkoiJZH/Z6h GZLS76upDs1f9Vmxi2VMS8FkkxM7vAqNVx+OU= MIME-Version: 1.0 In-Reply-To: References: <69c1813d-1a9a-4686-9768-8ec1910a45f8@d19g2000prh.googlegroups.com> <018ur6da31iv7us9gm4dpgl7tfl0i6snb2@4ax.com> <929lmdF59lU1@mid.individual.net> Date: Tue, 3 May 2011 09:51:20 -0700 Subject: Re: Coolest Python recipe of all time From: geremy condra To: Terry Reedy Content-Type: text/plain; charset=ISO-8859-1 Cc: python-list@python.org X-BeenThere: python-list@python.org X-Mailman-Version: 2.1.12 Precedence: list List-Id: General discussion list for the Python programming language List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , Newsgroups: comp.lang.python Message-ID: Lines: 24 NNTP-Posting-Host: 82.94.164.166 X-Trace: 1304441504 news.xs4all.nl 81475 [::ffff:82.94.164.166]:53401 X-Complaints-To: abuse@xs4all.nl Xref: x330-a1.tempe.blueboxinc.net comp.lang.python:4552 On Tue, May 3, 2011 at 8:49 AM, Terry Reedy wrote: > On 5/3/2011 2:29 AM, Gregory Ewing wrote: >> >> Terry Reedy wrote: >>> >>> The trick is that replacing x with j and evaluating therefore causes >>> (in Python) all the coefficients of x (now j) to be added together >>> separately from all the constant terms to reduce the linear equation >>> to a*x+b (= 0 implied). >> >> Hmmm... so if we used quaternions, could we solve systems >> of linear equations in 3 variables? > > Yes and no. The use of 1*j merely collected and added together all the > multipliers of 'x' (and all the constant terms). That is a fairly trivial > matter of constant folding. Systems of linear equations are usually > presented in that form already. The actual solution to the simple equation > is in the formula x = -a/b (where a and b are the sums). The solution > formula for three variables would be far more complex. Or just use a gauss-jordan solver, which has the advantage of being easy to explain and possible to verify by hand on small instances. Geremy Condra