From MAILER-DAEMON Wed Sep 08 12:46:26 2010
Received: from mailman by lists.gnu.org with archive (Exim 4.43)
	id 1OtNn0-0000Js-7p
	for mharc-axiom-developer@gnu.org; Wed, 08 Sep 2010 12:46:26 -0400
Received: from [140.186.70.92] (port=53088 helo=eggs.gnu.org)
	by lists.gnu.org with esmtp (Exim 4.43) id 1OtNmt-0000IC-Ag
	for axiom-developer@nongnu.org; Wed, 08 Sep 2010 12:46:23 -0400
Received: from Debian-exim by eggs.gnu.org with spam-scanned (Exim 4.69)
	(envelope-from <ax87438@martinb.com>) id 1OtNmr-0007gH-VB
	for axiom-developer@nongnu.org; Wed, 08 Sep 2010 12:46:19 -0400
Received: from moutng.kundenserver.de ([212.227.126.171]:57516)
	by eggs.gnu.org with esmtp (Exim 4.69)
	(envelope-from <ax87438@martinb.com>) id 1OtNmr-0007fe-IU
	for axiom-developer@nongnu.org; Wed, 08 Sep 2010 12:46:17 -0400
Received: from suse113.localnet (188-223-141-163.zone14.bethere.co.uk
	[188.223.141.163])
	by mrelayeu.kundenserver.de (node=mrbap0) with ESMTP (Nemesis)
	id 0LuJK3-1OiYF81a5f-011cYx; Wed, 08 Sep 2010 18:46:13 +0200
From: Martin Baker <ax87438@martinb.com>
To: axiom-developer@nongnu.org
Date: Wed, 8 Sep 2010 17:46:29 +0100
User-Agent: KMail/1.13.5 (Linux/2.6.34-12-desktop; KDE/4.4.4; x86_64; ; )
References: <4B5B3779.4060900@axiom-developer.org>
	<4B6436F5.3060201@axiom-developer.org>
	<201001301618.50345.ax87438@martinb.com>
In-Reply-To: <201001301618.50345.ax87438@martinb.com>
MIME-Version: 1.0
Content-Type: Text/Plain;
  charset="iso-8859-1"
Content-Transfer-Encoding: 7bit
Message-Id: <201009081746.29981.ax87438@martinb.com>
X-Provags-ID: V02:K0:9NXPgwm6KUBHsPsvAIV7IruwAP10R5/TiZgaDM2wSiE
	ADs0RRp+mbQSzKeLNNhwS9/juJnKwaJvgVZyKPJfUrNByZhpjz
	5ZR+Pq++1Ly9KnwYWtgy7+9QXJtt5I81TXDWHBn82CtQdyEelI
	KRKWC6ApWoigjH1kNkOHxQ9sJ8zuuskfObqyzJ+IuBEnwNNXyc
	4jL5AKSEox2+uRBikTp3A==
X-detected-operating-system: by eggs.gnu.org: Genre and OS details not
	recognized.
Subject: [Axiom-developer] alternative geometry/graphics framework
X-BeenThere: axiom-developer@nongnu.org
X-Mailman-Version: 2.1.5
Precedence: list
List-Id: Axiom Developers <axiom-developer.nongnu.org>
List-Unsubscribe: <http://lists.nongnu.org/mailman/listinfo/axiom-developer>, 
	<mailto:axiom-developer-request@nongnu.org?subject=unsubscribe>
List-Archive: <http://lists.nongnu.org/archive/html/axiom-developer>
List-Post: <mailto:axiom-developer@nongnu.org>
List-Help: <mailto:axiom-developer-request@nongnu.org?subject=help>
List-Subscribe: <http://lists.nongnu.org/mailman/listinfo/axiom-developer>,
	<mailto:axiom-developer-request@nongnu.org?subject=subscribe>
X-List-Received-Date: Wed, 08 Sep 2010 16:46:24 -0000

I have written an alternative geometry/graphics framework.

I have put links to the code and documentation at the bottom of this
page here:
http://www.euclideanspace.com/maths/standards/program/mycode/graph/

It is based on a scenegraph structure which allows
lines, surfaces and higher dimensional structures to be defined and
modified by transforms, clipping boxes, material (line width, colour and
so on) which control individual nodes or whole branches in the
scenegraph. This means that, for example, transforms can be applied
continuously to a whole branch without altering the nodes in that branch,
this means that we can avoid a cumulative build up of floating point errors
when many transforms are applied.

This framework supports all the 2D and 3D draw and plotting that the
existing framework supports. Use of transforms in the scenegraph means that
many plots can be combined together, either overlayed or next to each other
and combined with various show scales, grids text annotations and so on.

The framework supports:
* Different algebras: vector, matrix, complex numbers and Clifford Algebra.
* Different geometries: Euclidean, projective, conformal (hyperbolic to be
  added later)
* Different coordinate systems: Cartesian coordinates, Argand Plane,
  Spherical, Cylindrical.
* Different number of dimensions: 2,3 (4 and more to be added)

The aim is not just to draw plots but also to allow experimentation with
arbitrary shapes being transformed, for instance to transform a given shape
by reflecting in a circle. To see the sort of thing I am trying to do here
see first chapter of [Dorst,Fontijne & Mann 2007] except this is not just
limited to Clifford algebra and conformal space.

Currently Axiom/FriCAS has domains that represent transforms:
(for example: dhmatrix.spad.pamphlet, moebius.spad.pamphlet and so on) but
they are just standalone domains which don't fit into a wider graphics
framework so it is more difficult to combine them with other things and
to experiment. The aim of this framework is to create a way for geometry
and graphics related entities to interwork.

This scenegraph framework currently exports to file formats like: SVG,X3D,
VRML and Wavefront(obj) but currently there is not the capability to interact
directly in a graphical way, that is, there is no xwindow support. The
existing graphics framework allows user interaction in a very limited way in
that the graphics can be displayed in an xwindow and there is then the
ability to do simple transforms like pan and zoom.
However this is not suitable for the new framework because the existing
framework builds a data structure to represent the graphics
using SPAD and then passes this to 'C' code. The 'C' code then does the pan,
zoom and so on and can then save to a file format.

Apart from being very messy trying to mix different languages, this would not
work for the new framework because:
* We want to use non-linear as well as linear transforms.
* I would like to be able to work with discontinuous (fractal) shapes.
* I want to use different algebras such as vector, matrix, complex numbers
  and Clifford algebras and I would not want to implement these algebras
  using 'C'.

Possible options might be:
* Implement an OpenGL canvas in SPAD which allows SPAD to directly output to
  OpenGL (a thin wrapper for OpenGL calls) and also to read and process mouse
  position.
* Implement Axiom/FriCAS as library code that can be called from other
  languages.
* Implement some form of two-way API which allows other programs to use SPAD
  algebra.

Bill Page mentioned that there is an OpenGL binding for lisp.
http://common-lisp.net/project/cl-opengl/
http://repo.or.cz/w/cl-glfw.git
Would it be possible to create a thin wrapper to allow these functions to be
called from SPAD?
OpenGL consists of a couple of hundred function calls but this would only need
a small subset of those (because transforms would be done by SPAD and not 
OpenGL).
We would also need the ability to read mouse movements from SPAD to allow user
interaction.

Do you think this would be possible?

Martin


From MAILER-DAEMON Wed Sep 08 21:31:44 2010
Received: from mailman by lists.gnu.org with archive (Exim 4.43)
	id 1OtVzL-0004kN-Sf
	for mharc-axiom-developer@gnu.org; Wed, 08 Sep 2010 21:31:43 -0400
Received: from [140.186.70.92] (port=48627 helo=eggs.gnu.org)
	by lists.gnu.org with esmtp (Exim 4.43) id 1OtVzI-0004jj-RZ
	for axiom-developer@nongnu.org; Wed, 08 Sep 2010 21:31:42 -0400
Received: from Debian-exim by eggs.gnu.org with spam-scanned (Exim 4.69)
	(envelope-from <daly@axiom-developer.org>) id 1OtVzG-0007Nz-Hd
	for axiom-developer@nongnu.org; Wed, 08 Sep 2010 21:31:40 -0400
Received: from vs338.rosehosting.com ([209.135.140.38]:45490
	helo=axiom-developer.org) by eggs.gnu.org with esmtp (Exim 4.69)
	(envelope-from <daly@axiom-developer.org>) id 1OtVzG-0007Nn-Cc
	for axiom-developer@nongnu.org; Wed, 08 Sep 2010 21:31:38 -0400
Received: from [192.168.1.100] (dynamic-acs-72-23-235-203.zoominternet.net
	[72.23.235.203])
	by axiom-developer.org (8.12.8/8.12.8) with ESMTP id o891VWYb031624;
	Wed, 8 Sep 2010 20:31:33 -0500
Message-ID: <4C88396D.8040002@axiom-developer.org>
Date: Wed, 08 Sep 2010 21:33:33 -0400
From: Tim Daly <daly@axiom-developer.org>
User-Agent: Thunderbird 2.0.0.21 (Windows/20090302)
MIME-Version: 1.0
To: Martin Baker <ax87438@martinb.com>
Subject: Re: [Axiom-developer] alternative geometry/graphics framework
References: <4B5B3779.4060900@axiom-developer.org>	<4B6436F5.3060201@axiom-developer.org>	<201001301618.50345.ax87438@martinb.com>
	<201009081746.29981.ax87438@martinb.com>
In-Reply-To: <201009081746.29981.ax87438@martinb.com>
Content-Type: text/plain; charset=ISO-8859-1; format=flowed
Content-Transfer-Encoding: 7bit
X-detected-operating-system: by eggs.gnu.org: Linux 2.6? (barebone, rare!)
Cc: axiom-developer@nongnu.org
X-BeenThere: axiom-developer@nongnu.org
X-Mailman-Version: 2.1.5
Precedence: list
List-Id: Axiom Developers <axiom-developer.nongnu.org>
List-Unsubscribe: <http://lists.nongnu.org/mailman/listinfo/axiom-developer>, 
	<mailto:axiom-developer-request@nongnu.org?subject=unsubscribe>
List-Archive: <http://lists.nongnu.org/archive/html/axiom-developer>
List-Post: <mailto:axiom-developer@nongnu.org>
List-Help: <mailto:axiom-developer-request@nongnu.org?subject=help>
List-Subscribe: <http://lists.nongnu.org/mailman/listinfo/axiom-developer>,
	<mailto:axiom-developer-request@nongnu.org?subject=subscribe>
X-List-Received-Date: Thu, 09 Sep 2010 01:31:42 -0000

I'll look into it as soon as I'm able. --Tim

Martin Baker wrote:
> I have written an alternative geometry/graphics framework.
>
> I have put links to the code and documentation at the bottom of this
> page here:
> http://www.euclideanspace.com/maths/standards/program/mycode/graph/
>
> It is based on a scenegraph structure which allows
> lines, surfaces and higher dimensional structures to be defined and
> modified by transforms, clipping boxes, material (line width, colour and
> so on) which control individual nodes or whole branches in the
> scenegraph. This means that, for example, transforms can be applied
> continuously to a whole branch without altering the nodes in that branch,
> this means that we can avoid a cumulative build up of floating point errors
> when many transforms are applied.
>
> This framework supports all the 2D and 3D draw and plotting that the
> existing framework supports. Use of transforms in the scenegraph means that
> many plots can be combined together, either overlayed or next to each other
> and combined with various show scales, grids text annotations and so on.
>
> The framework supports:
> * Different algebras: vector, matrix, complex numbers and Clifford Algebra.
> * Different geometries: Euclidean, projective, conformal (hyperbolic to be
>   added later)
> * Different coordinate systems: Cartesian coordinates, Argand Plane,
>   Spherical, Cylindrical.
> * Different number of dimensions: 2,3 (4 and more to be added)
>
> The aim is not just to draw plots but also to allow experimentation with
> arbitrary shapes being transformed, for instance to transform a given shape
> by reflecting in a circle. To see the sort of thing I am trying to do here
> see first chapter of [Dorst,Fontijne & Mann 2007] except this is not just
> limited to Clifford algebra and conformal space.
>
> Currently Axiom/FriCAS has domains that represent transforms:
> (for example: dhmatrix.spad.pamphlet, moebius.spad.pamphlet and so on) but
> they are just standalone domains which don't fit into a wider graphics
> framework so it is more difficult to combine them with other things and
> to experiment. The aim of this framework is to create a way for geometry
> and graphics related entities to interwork.
>
> This scenegraph framework currently exports to file formats like: SVG,X3D,
> VRML and Wavefront(obj) but currently there is not the capability to interact
> directly in a graphical way, that is, there is no xwindow support. The
> existing graphics framework allows user interaction in a very limited way in
> that the graphics can be displayed in an xwindow and there is then the
> ability to do simple transforms like pan and zoom.
> However this is not suitable for the new framework because the existing
> framework builds a data structure to represent the graphics
> using SPAD and then passes this to 'C' code. The 'C' code then does the pan,
> zoom and so on and can then save to a file format.
>
> Apart from being very messy trying to mix different languages, this would not
> work for the new framework because:
> * We want to use non-linear as well as linear transforms.
> * I would like to be able to work with discontinuous (fractal) shapes.
> * I want to use different algebras such as vector, matrix, complex numbers
>   and Clifford algebras and I would not want to implement these algebras
>   using 'C'.
>
> Possible options might be:
> * Implement an OpenGL canvas in SPAD which allows SPAD to directly output to
>   OpenGL (a thin wrapper for OpenGL calls) and also to read and process mouse
>   position.
> * Implement Axiom/FriCAS as library code that can be called from other
>   languages.
> * Implement some form of two-way API which allows other programs to use SPAD
>   algebra.
>
> Bill Page mentioned that there is an OpenGL binding for lisp.
> http://common-lisp.net/project/cl-opengl/
> http://repo.or.cz/w/cl-glfw.git
> Would it be possible to create a thin wrapper to allow these functions to be
> called from SPAD?
> OpenGL consists of a couple of hundred function calls but this would only need
> a small subset of those (because transforms would be done by SPAD and not 
> OpenGL).
> We would also need the ability to read mouse movements from SPAD to allow user
> interaction.
>
> Do you think this would be possible?
>
> Martin
>
> _______________________________________________
> Axiom-developer mailing list
> Axiom-developer@nongnu.org
> http://lists.nongnu.org/mailman/listinfo/axiom-developer
>
>   


From MAILER-DAEMON Sat Sep 25 02:47:44 2010
Received: from mailman by lists.gnu.org with archive (Exim 4.43)
	id 1OzOXw-0007nr-Ma
	for mharc-axiom-developer@gnu.org; Sat, 25 Sep 2010 02:47:44 -0400
Received: from [140.186.70.92] (port=47381 helo=eggs.gnu.org)
	by lists.gnu.org with esmtp (Exim 4.43) id 1OzOXu-0007nk-3C
	for axiom-developer@nongnu.org; Sat, 25 Sep 2010 02:47:43 -0400
Received: from Debian-exim by eggs.gnu.org with spam-scanned (Exim 4.69)
	(envelope-from <daly@axiom-developer.org>) id 1OzOXs-0004Vj-QT
	for axiom-developer@nongnu.org; Sat, 25 Sep 2010 02:47:41 -0400
Received: from vs338.rosehosting.com ([209.135.140.38]:49193
	helo=axiom-developer.org) by eggs.gnu.org with esmtp (Exim 4.69)
	(envelope-from <daly@axiom-developer.org>) id 1OzOXs-0004VV-Me
	for axiom-developer@nongnu.org; Sat, 25 Sep 2010 02:47:40 -0400
Received: from axiom-developer.org (lincoln.rosehosting.com [127.0.0.1])
	by axiom-developer.org (8.12.8/8.12.8) with ESMTP id o8P6lcYb022781;
	Sat, 25 Sep 2010 01:47:38 -0500
From: daly@axiom-developer.org
Received: (from daly@localhost)
	by axiom-developer.org (8.12.8/8.12.8/Submit) id o8P6lcFO022778;
	Sat, 25 Sep 2010 01:47:38 -0500
Date: Sat, 25 Sep 2010 01:47:38 -0500
Message-Id: <201009250647.o8P6lcFO022778@axiom-developer.org>
To: axiom-developer@nongnu.org, daly@axiom-developer.org
X-detected-operating-system: by eggs.gnu.org: Linux 2.6? (barebone, rare!)
Cc: 
Subject: [Axiom-developer] Axiom September 2010 Release
X-BeenThere: axiom-developer@nongnu.org
X-Mailman-Version: 2.1.5
Precedence: list
List-Id: Axiom Developers <axiom-developer.nongnu.org>
List-Unsubscribe: <http://lists.nongnu.org/mailman/listinfo/axiom-developer>, 
	<mailto:axiom-developer-request@nongnu.org?subject=unsubscribe>
List-Archive: <http://lists.nongnu.org/archive/html/axiom-developer>
List-Post: <mailto:axiom-developer@nongnu.org>
List-Help: <mailto:axiom-developer-request@nongnu.org?subject=help>
List-Subscribe: <http://lists.nongnu.org/mailman/listinfo/axiom-developer>,
	<mailto:axiom-developer-request@nongnu.org?subject=subscribe>
X-List-Received-Date: Sat, 25 Sep 2010 06:47:43 -0000

Axiom has been released for September 2010.
The website has the sources and the ubuntu binary.
Other binaries are being built and will appear when they pass testing.

Tim Daly

=============================================================================

September 2010 Release

This release concentrated on treeshaking the compiler code into book volume 9.
Due to the complexity of this task it will take several releases to complete.

Additional global changes included collecting all of the Axiom-related
latex macros into axiom.sty and removing references to Aldor.

Makefile 
  always run help extractions in parallel
  build src/input quietly

books/bookvol0
   remove references to aldor

books/bookvol1
   remove references to aldor

books/bookvol5
   expose StreamTensor, U32Vector, U32VectorPolynomialOperations
   mark pure common lisp routines
   merge ptrop, varini
   move latex macros to axiom.sty
   remove POLYVEC
   remove compile, duplicated in vol9
   remove memq
   remove references to aldor
   treeshake
   remove $useNewParser

books/bookvol7.1
   move latex macros to axiom.sty
   remove references to aldor
   rewrite \pagehead to \pagetitle

books/bookvol9
   cross-reference functions and variables
   move latex macros to axiom.sty
   remove memq
   treeshake the compiler code
   remove $useNewParser

books/bookvol10
   move latex macros to axiom.sty
   move GOPT0 from bookvol10.3

books/bookvol10.1
   move latex macros to axiom.sty

books/bookvol10.2
   move latex macros to axiom.sty

books/bookvol10.3
   add U32Vector, move GOPT0 from bookvol10.4
   move latex macros to axiom.sty

books/bookvol10.4
   add StreamTensor, U32VectorPolynomialOperations
   fix ScriptTensor regression test
   move latex macros to axiom.sty
   remove POLYVEC
   update Chinese Remainder documentation

books/bookvol10.5
   move latex macros to axiom.sty

books/bookvolbib
   Parnas & Madey [PM95], Parnas & Jin [PJ10], GCL92, AS64, NIST10, RF94, 
   Hamdy [Ham04], Steele [Ste90], Tim Lahey's Sage Integration Test Suite

books/ps/
   v103guessoptionfunctions0, v103u32vector, v104streamtensor,
   v104u32vectorpolynomialoperations, v104u32vectorpolynomialoperations

src/algebra
   Makefile help and test for StreamTensor
   Makefile help and test for U32Vector
   Makefile remove references to aldor
   Makefile test and help for POLYVEC
   Makefile remove POLYVEC
   Makefile add help and test for new algebra
   Makefile handle case-insensitive MAC filesystem
   Makefile remove new algebra scaffolding code

src/doc
   axiom.sty collect all script commands in one place
   axiom.sty consolidate latex macros

src/input
   Makefile add guess.input, manuel.input, risch.input
   guess.input test examples of the GUESS package
   kamke3.input clean up broken tests
   manuel.input add Manuel's integral to test suite
   richlog300-391.input clean up broken tests
   richtrig800-899.input clean up broken tests
   risch.input illustrate the Risch algorithm
   setcmd.input clean up broken tests
   test.input clean up broken tests

src/interp
   Makefile merge varini
   Makefile remove nspadaux, mark, pspad1, pspad2, ptrop, wi1, wi2
   *.lisp remove memq
   treeshake compiler -- br-con, cattable, compiler
   remove nspadaux.lisp, mark.lisp, pspad1.lisp, pspad2.lisp, 
   remove ptrop.lisp, wi1.lisp, wi2.lisp
   add HTMLFormat code i-output.lisp, vmlisp.lisp

src/scripts/tex
   axiom.sty collect all script commands
   axiom.sty consolidate latex macros

src/axiom-website/
   documentation.html add Knuth quote per W. Sit
   download.html add debian, fedora, mandriva, opensuse, slackware, vector
   download.html update ubuntu yum advice





From MAILER-DAEMON Sun Sep 26 09:22:27 2010
Received: from mailman by lists.gnu.org with archive (Exim 4.43)
	id 1OzrBT-0005zz-2O
	for mharc-axiom-developer@gnu.org; Sun, 26 Sep 2010 09:22:27 -0400
Received: from [140.186.70.92] (port=45579 helo=eggs.gnu.org)
	by lists.gnu.org with esmtp (Exim 4.43) id 1OzDHi-0005YO-5O
	for axiom-developer@nongnu.org; Fri, 24 Sep 2010 14:46:15 -0400
Received: from Debian-exim by eggs.gnu.org with spam-scanned (Exim 4.69)
	(envelope-from <d01c@uni-bremen.de>) id 1OzDHg-00056J-UF
	for axiom-developer@nongnu.org; Fri, 24 Sep 2010 14:46:14 -0400
Received: from gabriel-1.zfn.uni-bremen.de ([134.102.20.34]:33904)
	by eggs.gnu.org with esmtp (Exim 4.69)
	(envelope-from <d01c@uni-bremen.de>) id 1OzDHg-000568-Ij
	for axiom-developer@nongnu.org; Fri, 24 Sep 2010 14:46:12 -0400
Received: from [192.168.1.210] (82-198-197-155.briteline.de [82.198.197.155])
	(using TLSv1 with cipher DHE-RSA-AES256-SHA (256/256 bits))
	(No client certificate requested)
	by gabriel-1.zfn.uni-bremen.de (Postfix) with ESMTPSA id 8D951B9FC7
	for <axiom-developer@nongnu.org>; Fri, 24 Sep 2010 20:46:09 +0200 (CEST)
DKIM-Signature: v=1; a=rsa-sha256; c=simple/simple; d=uni-bremen.de; s=dkim;
	t=1285353969; bh=a1iupnHEgvrNn0PECnJhaHHg7CA80vc12JPNizqo2ig=;
	h=Message-ID:Date:From:MIME-Version:To:Content-Type:
	Content-Transfer-Encoding; b=J0JMKtFpOkxFgL/W/qtVT44Dz4pZApsLU/2ss
	joGQYSLyZp9nPDA3OD6+EkW6NIeOupkMlgmdjdIHZ7jpqKCV/kJYK0x3Dtgjb3CDCO3
	0cPhja+/YbenXGOcA5LtUsAB1ppEp5VvwdG+NDXXPN7wgBqPeqsMnBderBCfPPfdmvE
	=
Message-ID: <4C9CF1ED.7020806@uni-bremen.de>
Date: Fri, 24 Sep 2010 20:46:05 +0200
From: "Dr.-Ing. Ingo D. Rullhusen" <d01c@uni-bremen.de>
Organization: =?ISO-8859-15?Q?Universit=E4t_Bremen?=
User-Agent: Mozilla/5.0 (X11; U; Linux x86_64; de;
	rv:1.9.1.13) Gecko/20100916 SeaMonkey/2.0.8
MIME-Version: 1.0
To: axiom-developer@nongnu.org
X-Enigmail-Version: 1.0.1
Content-Type: text/plain; charset=ISO-8859-15
Content-Transfer-Encoding: 7bit
X-detected-operating-system: by eggs.gnu.org: GNU/Linux 2.6 (newer, 2)
X-Mailman-Approved-At: Sun, 26 Sep 2010 09:22:25 -0400
Subject: [Axiom-developer] wrong sign in limit?
X-BeenThere: axiom-developer@nongnu.org
X-Mailman-Version: 2.1.5
Precedence: list
List-Id: Axiom Developers <axiom-developer.nongnu.org>
List-Unsubscribe: <http://lists.nongnu.org/mailman/listinfo/axiom-developer>, 
	<mailto:axiom-developer-request@nongnu.org?subject=unsubscribe>
List-Archive: <http://lists.nongnu.org/archive/html/axiom-developer>
List-Post: <mailto:axiom-developer@nongnu.org>
List-Help: <mailto:axiom-developer-request@nongnu.org?subject=help>
List-Subscribe: <http://lists.nongnu.org/mailman/listinfo/axiom-developer>,
	<mailto:axiom-developer-request@nongnu.org?subject=subscribe>
X-List-Received-Date: Fri, 24 Sep 2010 18:46:15 -0000

Hello,

perhaps i have found a bug in the limit operation. I am using

                        AXIOM Computer Algebra System
                         Version: Axiom (July 2010)
              Timestamp: Wednesday August 18, 2010 at 11:16:49

on a x86_64 with linux. There seems to be an error in the sign, if you
compare the limit result with the de L'Hospital one.

>>>>Example

-- problem with limit

)clear

Z:=sqrt(a^2+h^2)-a
N:=a*sqrt(a^2+h^2) - a^2 - h^2

limit( Z/N, h=0 )

-- this seems to be wrong!
-- because de L'Hospital says

dZ:=differentiate( Z, h, 2 )
dN:=differentiate( N, h, 2 )

eval(dZ/dN,h=0)

-- which seems to be right (see graph): opposite sign!

>>>>EndExample

which generates the output:

                        AXIOM Computer Algebra System
                         Version: Axiom (July 2010)
              Timestamp: Wednesday August 18, 2010 at 11:16:49
-----------------------------------------------------------------------------
   Issue )copyright to view copyright notices.
   Issue )summary for a summary of useful system commands.
   Issue )quit to leave AXIOM and return to shell.
   Visit http://axiom-developer.org for more information
-----------------------------------------------------------------------------
 
   Re-reading compress.daase   Re-reading interp.daase
   Re-reading operation.daase
   Re-reading category.daase
   Re-reading browse.daase
(1) ->
(1) -> )read hos.input

-- problem with limit

)clear
 
   Use )clear all to clear everything in the workspace. Use )clear
      completely to clear everything in the workspace and internal
      tables. Other )clear keyword arguments are
             modes
             operations
             properties
             types
             values
      or abbreviations thereof. Issue )clear? for more information.

Z:=sqrt(a^2+h^2)-a
 

         +-------+
         | 2    2
   (1)  \|h  + a   - a
                                                     Type: Expression
Integer
N:=a*sqrt(a^2+h^2) - a^2 - h^2
 

          +-------+
          | 2    2     2    2
   (2)  a\|h  + a   - h  - a
                                                     Type: Expression
Integer

limit( Z/N, h=0 )
 

        1
   (3)  -
        a
                        Type: Union(OrderedCompletion Expression
Integer,...)

-- this seems to be wrong!
-- because de L'Hospital says

dZ:=differentiate( Z, h, 2 )
 

                  2
                 a
   (4)  -------------------
                  +-------+
          2    2  | 2    2
        (h  + a )\|h  + a
                                                     Type: Expression
Integer
dN:=differentiate( N, h, 2 )
 

                      +-------+
             2     2  | 2    2     3
        (- 2h  - 2a )\|h  + a   + a
   (5)  ----------------------------
                       +-------+
               2    2  | 2    2
             (h  + a )\|h  + a
                                                     Type: Expression
Integer

eval(dZ/dN,h=0)
 

               1
   (6)  - ----------
            +--+
            | 2
          2\|a   - a
                                                     Type: Expression
Integer
(7) ->
(7) ->

(see sign)

Thanks
   Ingo



From MAILER-DAEMON Sun Sep 26 11:47:03 2010
Received: from mailman by lists.gnu.org with archive (Exim 4.43)
	id 1OztRO-00015x-W1
	for mharc-axiom-developer@gnu.org; Sun, 26 Sep 2010 11:47:03 -0400
Received: from [140.186.70.92] (port=35924 helo=eggs.gnu.org)
	by lists.gnu.org with esmtp (Exim 4.43) id 1OztRM-00015i-Li
	for axiom-developer@nongnu.org; Sun, 26 Sep 2010 11:47:01 -0400
Received: from Debian-exim by eggs.gnu.org with spam-scanned (Exim 4.69)
	(envelope-from <daly@axiom-developer.org>) id 1OztRL-0005os-E2
	for axiom-developer@nongnu.org; Sun, 26 Sep 2010 11:47:00 -0400
Received: from vs338.rosehosting.com ([209.135.140.38]:38296
	helo=axiom-developer.org) by eggs.gnu.org with esmtp (Exim 4.69)
	(envelope-from <daly@axiom-developer.org>) id 1OztRL-0005ok-8O
	for axiom-developer@nongnu.org; Sun, 26 Sep 2010 11:46:59 -0400
Received: from [192.168.1.4] (dynamic-acs-72-23-235-203.zoominternet.net
	[72.23.235.203])
	by axiom-developer.org (8.12.8/8.12.8) with ESMTP id o8QFktj5032569;
	Sun, 26 Sep 2010 10:46:56 -0500
Message-ID: <4C9F9507.4010108@axiom-developer.org>
Date: Sun, 26 Sep 2010 14:46:31 -0400
From: Tim Daly <daly@axiom-developer.org>
User-Agent: Mozilla/5.0 (Windows; U; Windows NT 5.1; en-US;
	rv:1.9.2.9) Gecko/20100825 Thunderbird/3.1.3
MIME-Version: 1.0
To: axiom-developer@nongnu.org
Subject: Re: [Axiom-developer] wrong sign in limit?
References: <4C9CF1ED.7020806@uni-bremen.de>
In-Reply-To: <4C9CF1ED.7020806@uni-bremen.de>
Content-Type: text/plain; charset=ISO-8859-15; format=flowed
Content-Transfer-Encoding: 7bit
X-detected-operating-system: by eggs.gnu.org: Linux 2.6? (barebone, rare!)
X-BeenThere: axiom-developer@nongnu.org
X-Mailman-Version: 2.1.5
Precedence: list
List-Id: Axiom Developers <axiom-developer.nongnu.org>
List-Unsubscribe: <http://lists.nongnu.org/mailman/listinfo/axiom-developer>, 
	<mailto:axiom-developer-request@nongnu.org?subject=unsubscribe>
List-Archive: <http://lists.nongnu.org/archive/html/axiom-developer>
List-Post: <mailto:axiom-developer@nongnu.org>
List-Help: <mailto:axiom-developer-request@nongnu.org?subject=help>
List-Subscribe: <http://lists.nongnu.org/mailman/listinfo/axiom-developer>,
	<mailto:axiom-developer-request@nongnu.org?subject=subscribe>
X-List-Received-Date: Sun, 26 Sep 2010 15:47:01 -0000

  Interesting. I'll look at it further. --Tim

On 9/24/2010 2:46 PM, Dr.-Ing. Ingo D. Rullhusen wrote:
> Hello,
>
> perhaps i have found a bug in the limit operation. I am using
>
>                          AXIOM Computer Algebra System
>                           Version: Axiom (July 2010)
>                Timestamp: Wednesday August 18, 2010 at 11:16:49
>
> on a x86_64 with linux. There seems to be an error in the sign, if you
> compare the limit result with the de L'Hospital one.
>
>>>>> Example
> -- problem with limit
>
> )clear
>
> Z:=sqrt(a^2+h^2)-a
> N:=a*sqrt(a^2+h^2) - a^2 - h^2
>
> limit( Z/N, h=0 )
>
> -- this seems to be wrong!
> -- because de L'Hospital says
>
> dZ:=differentiate( Z, h, 2 )
> dN:=differentiate( N, h, 2 )
>
> eval(dZ/dN,h=0)
>
> -- which seems to be right (see graph): opposite sign!
>
>>>>> EndExample
> which generates the output:
>
>                          AXIOM Computer Algebra System
>                           Version: Axiom (July 2010)
>                Timestamp: Wednesday August 18, 2010 at 11:16:49
> -----------------------------------------------------------------------------
>     Issue )copyright to view copyright notices.
>     Issue )summary for a summary of useful system commands.
>     Issue )quit to leave AXIOM and return to shell.
>     Visit http://axiom-developer.org for more information
> -----------------------------------------------------------------------------
>
>     Re-reading compress.daase   Re-reading interp.daase
>     Re-reading operation.daase
>     Re-reading category.daase
>     Re-reading browse.daase
> (1) ->
> (1) ->  )read hos.input
>
> -- problem with limit
>
> )clear
>
>     Use )clear all to clear everything in the workspace. Use )clear
>        completely to clear everything in the workspace and internal
>        tables. Other )clear keyword arguments are
>               modes
>               operations
>               properties
>               types
>               values
>        or abbreviations thereof. Issue )clear? for more information.
>
> Z:=sqrt(a^2+h^2)-a
>
>
>           +-------+
>           | 2    2
>     (1)  \|h  + a   - a
>                                                       Type: Expression
> Integer
> N:=a*sqrt(a^2+h^2) - a^2 - h^2
>
>
>            +-------+
>            | 2    2     2    2
>     (2)  a\|h  + a   - h  - a
>                                                       Type: Expression
> Integer
>
> limit( Z/N, h=0 )
>
>
>          1
>     (3)  -
>          a
>                          Type: Union(OrderedCompletion Expression
> Integer,...)
>
> -- this seems to be wrong!
> -- because de L'Hospital says
>
> dZ:=differentiate( Z, h, 2 )
>
>
>                    2
>                   a
>     (4)  -------------------
>                    +-------+
>            2    2  | 2    2
>          (h  + a )\|h  + a
>                                                       Type: Expression
> Integer
> dN:=differentiate( N, h, 2 )
>
>
>                        +-------+
>               2     2  | 2    2     3
>          (- 2h  - 2a )\|h  + a   + a
>     (5)  ----------------------------
>                         +-------+
>                 2    2  | 2    2
>               (h  + a )\|h  + a
>                                                       Type: Expression
> Integer
>
> eval(dZ/dN,h=0)
>
>
>                 1
>     (6)  - ----------
>              +--+
>              | 2
>            2\|a   - a
>                                                       Type: Expression
> Integer
> (7) ->
> (7) ->
>
> (see sign)
>
> Thanks
>     Ingo
>
>
> _______________________________________________
> Axiom-developer mailing list
> Axiom-developer@nongnu.org
> http://lists.nongnu.org/mailman/listinfo/axiom-developer
>


From MAILER-DAEMON Sun Sep 26 16:52:03 2010
Received: from mailman by lists.gnu.org with archive (Exim 4.43)
	id 1OzyCZ-0001Tj-B4
	for mharc-axiom-developer@gnu.org; Sun, 26 Sep 2010 16:52:03 -0400
Received: from [140.186.70.92] (port=52044 helo=eggs.gnu.org)
	by lists.gnu.org with esmtp (Exim 4.43) id 1OzyCV-0001TF-5E
	for axiom-developer@nongnu.org; Sun, 26 Sep 2010 16:52:00 -0400
Received: from Debian-exim by eggs.gnu.org with spam-scanned (Exim 4.69)
	(envelope-from <ralf@hemmecke.de>) id 1OzyCT-0007Ai-RH
	for axiom-developer@nongnu.org; Sun, 26 Sep 2010 16:51:59 -0400
Received: from mo-p00-ob.rzone.de ([81.169.146.160]:55495)
	by eggs.gnu.org with esmtp (Exim 4.69)
	(envelope-from <ralf@hemmecke.de>) id 1OzyCT-0007AW-Ga
	for axiom-developer@nongnu.org; Sun, 26 Sep 2010 16:51:57 -0400
DKIM-Signature: v=1; a=rsa-sha1; c=relaxed/relaxed; t=1285534315; l=2868;
	s=domk; d=hemmecke.de;
	h=Content-Transfer-Encoding:Content-Type:In-Reply-To:References:
	Subject:To:MIME-Version:From:Date:X-RZG-CLASS-ID:X-RZG-AUTH;
	bh=3pMXYEnfKBvuz/gVPb1Ubp/Wc/Y=;
	b=ghAtZqbjgBC/zROtXEE7iquVZEcg88T7YBDZ9OqoKSb4+O+zLJvUy/PYfNQbYoBlQvd
	uZ3IAN75bjau3t5pgJ9a+UUgTx4QVLukAEQRpuGx69fYWKP4OHzc126tgvpSaVWRXE1Vv
	ha0ATBHZn1aUwBsnyBNJihtma1yEtieES70=
X-RZG-AUTH: :Pm0Ic2CgfvKqpyys4bXn/kzLaHTdYnTFXEqQA90LpSWk1ofrjzFkGjidQA==
X-RZG-CLASS-ID: mo00
Received: from [10.0.0.9] (62-47-55-42.adsl.highway.telekom.at [62.47.55.42])
	by post.strato.de (fruni mo28) (RZmta 23.5)
	with (DHE-RSA-AES128-SHA encrypted) ESMTP id 30022cm8QJZjqa ;
	Sun, 26 Sep 2010 22:51:23 +0200 (MEST)
Message-ID: <4C9FB24A.5000700@hemmecke.de>
Date: Sun, 26 Sep 2010 22:51:22 +0200
From: Ralf Hemmecke <ralf@hemmecke.de>
User-Agent: Mozilla/5.0 (X11; U; Linux i686; en-US;
	rv:1.9.2.9) Gecko/20100915 Thunderbird/3.1.4
MIME-Version: 1.0
To: axiom-developer@nongnu.org
Subject: Re: [Axiom-developer] wrong sign in limit?
References: <4C9CF1ED.7020806@uni-bremen.de>
	<4C9F9507.4010108@axiom-developer.org>
In-Reply-To: <4C9F9507.4010108@axiom-developer.org>
X-Sent-by-Hemmecke: ralf@hemmecke.de
Content-Type: text/plain; charset=ISO-8859-15; format=flowed
Content-Transfer-Encoding: 7bit
X-detected-operating-system: by eggs.gnu.org: Solaris 10 (beta)
X-BeenThere: axiom-developer@nongnu.org
X-Mailman-Version: 2.1.5
Precedence: list
List-Id: Axiom Developers <axiom-developer.nongnu.org>
List-Unsubscribe: <http://lists.nongnu.org/mailman/listinfo/axiom-developer>, 
	<mailto:axiom-developer-request@nongnu.org?subject=unsubscribe>
List-Archive: <http://lists.nongnu.org/archive/html/axiom-developer>
List-Post: <mailto:axiom-developer@nongnu.org>
List-Help: <mailto:axiom-developer-request@nongnu.org?subject=help>
List-Subscribe: <http://lists.nongnu.org/mailman/listinfo/axiom-developer>,
	<mailto:axiom-developer-request@nongnu.org?subject=subscribe>
X-List-Received-Date: Sun, 26 Sep 2010 20:52:00 -0000

I agree that the result is questionable, but what result do you expect?

This one is perhaps better

(0) -> limit(sqrt(a^2+x), x=0)

          +--+
          | 2
    (0)  \|a
                       Type: 
Union(OrderedCompletion(Expression(Integer)),...)

As long as nothing is known about a, that expression cannot be 
simplified further. In fact sqrt stands for *two* solutions.

Look into trigcat.spad.pamphlet.

RadicalCategory(): Category == with
   sqrt   : % -> %
       ++ sqrt(x) returns the square root of x.  The branch cut lies
       ++ along the negative real axis, continuous with quadrant II.

Even if we agree on the common convention that the root symbol for 
positive real arguments denotes a positive value, simplifying the above 
to just a would be wrong, as then plugging in a=-1 would violate the 
convention.

It's a quite subtle result. But it is wrong in the sense that the result 
should be as AXIOM says below.

Maple and Mathematica seem to use similar algorithms as they come up 
with the same questionable result. But look at what FullSimplify gives. 
Mathematica is just not trying hard enough.

I guess, if AXIOM could simplify Z/N, then the result would be better.

Ralf

(1) -> W := sqrt(a^2+h^2)

          +-------+
          | 2    2
    (1)  \|h  + a
                                   Type: Expression(Integer)
(2) -> Z := W -a

          +-------+
          | 2    2
    (2)  \|h  + a   - a
                                   Type: Expression(Integer)
(3) -> N := a*W -W^2

           +-------+
           | 2    2     2    2
    (3)  a\|h  + a   - h  - a
                                   Type: Expression(Integer)
(4) -> W*Z/N

    (4)  - 1
                                   Type: Expression(Integer)
(5) -> Z/N + 1/W

    (5)  0
                                   Type: Expression(Integer)
(6) -> limit(-1/W, h=0)

             1
    (6)  - -----
            +--+
            | 2
           \|a
      Type: Union(OrderedCompletion(Expression(Integer)),...)

==========================================================
BTW, why do you complain? ;-)

Mathematica 7.0 for Linux x86 (32-bit)
Copyright 1988-2008 Wolfram Research, Inc.

In[1]:= w = Sqrt[a^2+h^2]

               2    2
Out[1]= Sqrt[a  + h ]

In[2]:= z = w - a

                    2    2
Out[2]= -a + Sqrt[a  + h ]

In[3]:= n = a*w-w^2

           2    2           2    2
Out[3]= -a  - h  + a Sqrt[a  + h ]

In[4]:= Limit[z/n,h->0]

         1
Out[4]= -
         a

In[5]:= Limit[-1/w,h->0]

                 2
           Sqrt[a ]
Out[5]= -(--------)
               2
              a

In[6]:= Simplify[z/n]

                       2    2
             a - Sqrt[a  + h ]
Out[6]= -------------------------
          2    2           2    2
         a  + h  - a Sqrt[a  + h ]

In[7]:= FullSimplify[z/n]

                 1
Out[7]= -(-------------)
                 2    2
           Sqrt[a  + h ]


==============================

     |\^/|     Maple 11 (IBM INTEL LINUX)
._|\|   |/|_. Copyright (c) Maplesoft, a division of Waterloo Maple Inc. 
2007
  \  MAPLE  /  All rights reserved. Maple is a trademark of
  <____ ____>  Waterloo Maple Inc.
       |       Type ? for help.
 > w := sqrt(a^2+h^2);
                                       2    2 1/2
                                w := (a  + h )

 > z:=w-a;
                                     2    2 1/2
                              z := (a  + h )    - a

 > n:=a*w-w^2;
                                   2    2 1/2    2    2
                          n := a (a  + h )    - a  - h

 > limit(z/n,h=0);
                                       1/a

 > limit(-1/w,h=0);
                                         1
                                    - -------
                                        2 1/2
                                      (a )



From MAILER-DAEMON Tue Sep 28 17:56:13 2010
Received: from mailman by lists.gnu.org with archive (Exim 4.43)
	id 1P0i9l-0000FU-HH
	for mharc-axiom-developer@gnu.org; Tue, 28 Sep 2010 17:56:13 -0400
Received: from [140.186.70.92] (port=44877 helo=eggs.gnu.org)
	by lists.gnu.org with esmtp (Exim 4.43) id 1P0i9h-0000DR-Rh
	for axiom-developer@nongnu.org; Tue, 28 Sep 2010 17:56:11 -0400
Received: from Debian-exim by eggs.gnu.org with spam-scanned (Exim 4.69)
	(envelope-from <camm@maguirefamily.org>) id 1P0i9f-0005i5-Ga
	for axiom-developer@nongnu.org; Tue, 28 Sep 2010 17:56:09 -0400
Received: from tvwna-ip-d-105.princeton.org ([66.180.186.109]:51685
	helo=localhost.m.enhanced.com) by eggs.gnu.org with esmtp (Exim 4.69)
	(envelope-from <camm@maguirefamily.org>) id 1P0i9f-0005hm-Az
	for axiom-developer@nongnu.org; Tue, 28 Sep 2010 17:56:07 -0400
Received: from camm by localhost.m.enhanced.com with local (Exim 4.69)
	(envelope-from <camm@maguirefamily.org>)
	id 1P0i8u-0003QH-Mw; Tue, 28 Sep 2010 17:55:20 -0400
To: daly@axiom-developer.org
Subject: Re: [Axiom-developer] Axiom September 2010 Release
References: <201009250647.o8P6lcFO022778@axiom-developer.org>
From: Camm Maguire <camm@maguirefamily.org>
Date: Tue, 28 Sep 2010 17:55:20 -0400
In-Reply-To: <201009250647.o8P6lcFO022778@axiom-developer.org>
	(daly@axiom-developer.org's message of "Sat\,
	25 Sep 2010 01\:47\:38 -0500")
Message-ID: <87zkv1sfcn.fsf@maguirefamily.org>
User-Agent: Gnus/5.11 (Gnus v5.11) Emacs/22.2 (gnu/linux)
MIME-Version: 1.0
Content-Type: text/plain; charset=us-ascii
X-detected-operating-system: by eggs.gnu.org: GNU/Linux 2.6 (newer, 2)
Cc: axiom-developer@nongnu.org
X-BeenThere: axiom-developer@nongnu.org
X-Mailman-Version: 2.1.5
Precedence: list
List-Id: Axiom Developers <axiom-developer.nongnu.org>
List-Unsubscribe: <http://lists.nongnu.org/mailman/listinfo/axiom-developer>, 
	<mailto:axiom-developer-request@nongnu.org?subject=unsubscribe>
List-Archive: <http://lists.nongnu.org/archive/html/axiom-developer>
List-Post: <mailto:axiom-developer@nongnu.org>
List-Help: <mailto:axiom-developer-request@nongnu.org?subject=help>
List-Subscribe: <http://lists.nongnu.org/mailman/listinfo/axiom-developer>,
	<mailto:axiom-developer-request@nongnu.org?subject=subscribe>
X-List-Received-Date: Tue, 28 Sep 2010 21:56:11 -0000

Greetings!

Please consider this patch, as gcl has a file_exists function and this
prevents image linking when native object relocatin is absent. 

--- ./src/lib/pixmap.c.pamphlet.orig	2010-09-25 05:29:00.000000000 +0000
+++ ./src/lib/pixmap.c.pamphlet	2010-09-28 14:54:03.000000000 +0000
@@ -81,7 +81,7 @@
 /* returns true if the file exists */
 
 int
-file_exists(char *file)
+ax_file_exists(char *file)
 {
     FILE *f;
 
@@ -96,10 +96,10 @@
 {
     char com[512], zfile[512];
 
-    if (file_exists(file))
+    if (ax_file_exists(file))
         return fopen(file, mode);
     sprintf(zfile, "%s.Z", file);
-    if (file_exists(zfile)) {
+    if (ax_file_exists(zfile)) {
         sprintf(com, "gunzip -c %s.Z 2>/dev/null", file);
         return popen(com, mode);
     }



(Also unclear on how -lXpm gets into the link, but it appears to
somehow, so all else is good.  Uploaded into Debian unstable.)

Take care,

daly@axiom-developer.org writes:

> Axiom has been released for September 2010.
> The website has the sources and the ubuntu binary.
> Other binaries are being built and will appear when they pass testing.
>
> Tim Daly
>
> =============================================================================
>
> September 2010 Release
>
> This release concentrated on treeshaking the compiler code into book volume 9.
> Due to the complexity of this task it will take several releases to complete.
>
> Additional global changes included collecting all of the Axiom-related
> latex macros into axiom.sty and removing references to Aldor.
>
> Makefile 
>   always run help extractions in parallel
>   build src/input quietly
>
> books/bookvol0
>    remove references to aldor
>
> books/bookvol1
>    remove references to aldor
>
> books/bookvol5
>    expose StreamTensor, U32Vector, U32VectorPolynomialOperations
>    mark pure common lisp routines
>    merge ptrop, varini
>    move latex macros to axiom.sty
>    remove POLYVEC
>    remove compile, duplicated in vol9
>    remove memq
>    remove references to aldor
>    treeshake
>    remove $useNewParser
>
> books/bookvol7.1
>    move latex macros to axiom.sty
>    remove references to aldor
>    rewrite \pagehead to \pagetitle
>
> books/bookvol9
>    cross-reference functions and variables
>    move latex macros to axiom.sty
>    remove memq
>    treeshake the compiler code
>    remove $useNewParser
>
> books/bookvol10
>    move latex macros to axiom.sty
>    move GOPT0 from bookvol10.3
>
> books/bookvol10.1
>    move latex macros to axiom.sty
>
> books/bookvol10.2
>    move latex macros to axiom.sty
>
> books/bookvol10.3
>    add U32Vector, move GOPT0 from bookvol10.4
>    move latex macros to axiom.sty
>
> books/bookvol10.4
>    add StreamTensor, U32VectorPolynomialOperations
>    fix ScriptTensor regression test
>    move latex macros to axiom.sty
>    remove POLYVEC
>    update Chinese Remainder documentation
>
> books/bookvol10.5
>    move latex macros to axiom.sty
>
> books/bookvolbib
>    Parnas & Madey [PM95], Parnas & Jin [PJ10], GCL92, AS64, NIST10, RF94, 
>    Hamdy [Ham04], Steele [Ste90], Tim Lahey's Sage Integration Test Suite
>
> books/ps/
>    v103guessoptionfunctions0, v103u32vector, v104streamtensor,
>    v104u32vectorpolynomialoperations, v104u32vectorpolynomialoperations
>
> src/algebra
>    Makefile help and test for StreamTensor
>    Makefile help and test for U32Vector
>    Makefile remove references to aldor
>    Makefile test and help for POLYVEC
>    Makefile remove POLYVEC
>    Makefile add help and test for new algebra
>    Makefile handle case-insensitive MAC filesystem
>    Makefile remove new algebra scaffolding code
>
> src/doc
>    axiom.sty collect all script commands in one place
>    axiom.sty consolidate latex macros
>
> src/input
>    Makefile add guess.input, manuel.input, risch.input
>    guess.input test examples of the GUESS package
>    kamke3.input clean up broken tests
>    manuel.input add Manuel's integral to test suite
>    richlog300-391.input clean up broken tests
>    richtrig800-899.input clean up broken tests
>    risch.input illustrate the Risch algorithm
>    setcmd.input clean up broken tests
>    test.input clean up broken tests
>
> src/interp
>    Makefile merge varini
>    Makefile remove nspadaux, mark, pspad1, pspad2, ptrop, wi1, wi2
>    *.lisp remove memq
>    treeshake compiler -- br-con, cattable, compiler
>    remove nspadaux.lisp, mark.lisp, pspad1.lisp, pspad2.lisp, 
>    remove ptrop.lisp, wi1.lisp, wi2.lisp
>    add HTMLFormat code i-output.lisp, vmlisp.lisp
>
> src/scripts/tex
>    axiom.sty collect all script commands
>    axiom.sty consolidate latex macros
>
> src/axiom-website/
>    documentation.html add Knuth quote per W. Sit
>    download.html add debian, fedora, mandriva, opensuse, slackware, vector
>    download.html update ubuntu yum advice
>
>
>
>
> _______________________________________________
> Axiom-developer mailing list
> Axiom-developer@nongnu.org
> http://lists.nongnu.org/mailman/listinfo/axiom-developer
>
>
>
>

-- 
Camm Maguire			     		    camm@maguirefamily.org
==========================================================================
"The earth is but one country, and mankind its citizens."  --  Baha'u'llah


From MAILER-DAEMON Wed Sep 29 11:49:39 2010
Received: from mailman by lists.gnu.org with archive (Exim 4.43)
	id 1P0yuZ-0005j4-KQ
	for mharc-axiom-developer@gnu.org; Wed, 29 Sep 2010 11:49:39 -0400
Received: from [140.186.70.92] (port=49734 helo=eggs.gnu.org)
	by lists.gnu.org with esmtp (Exim 4.43) id 1P0yuV-0005iq-RQ
	for axiom-developer@nongnu.org; Wed, 29 Sep 2010 11:49:37 -0400
Received: from Debian-exim by eggs.gnu.org with spam-scanned (Exim 4.69)
	(envelope-from <camm@maguirefamily.org>) id 1P0yuU-0000Ec-4L
	for axiom-developer@nongnu.org; Wed, 29 Sep 2010 11:49:35 -0400
Received: from tvwna-ip-c-172.princeton.org ([66.180.187.89]:36783
	helo=localhost.m.enhanced.com) by eggs.gnu.org with esmtp (Exim 4.69)
	(envelope-from <camm@maguirefamily.org>) id 1P0yuT-0000EN-6X
	for axiom-developer@nongnu.org; Wed, 29 Sep 2010 11:49:34 -0400
Received: from camm by localhost.m.enhanced.com with local (Exim 4.69)
	(envelope-from <camm@maguirefamily.org>)
	id 1P0yuR-0001jF-JM; Wed, 29 Sep 2010 11:49:31 -0400
To: daly@axiom-developer.org
Subject: Re: [Axiom-developer] Axiom September 2010 Release
References: <201009250647.o8P6lcFO022778@axiom-developer.org>
From: Camm Maguire <camm@maguirefamily.org>
Date: Wed, 29 Sep 2010 11:49:31 -0400
In-Reply-To: <201009250647.o8P6lcFO022778@axiom-developer.org>
	(daly@axiom-developer.org's message of "Sat\,
	25 Sep 2010 01\:47\:38 -0500")
Message-ID: <87iq1opn1w.fsf@maguirefamily.org>
User-Agent: Gnus/5.11 (Gnus v5.11) Emacs/22.2 (gnu/linux)
MIME-Version: 1.0
Content-Type: text/plain; charset=us-ascii
X-detected-operating-system: by eggs.gnu.org: GNU/Linux 2.6 (newer, 2)
Cc: axiom-developer@nongnu.org
X-BeenThere: axiom-developer@nongnu.org
X-Mailman-Version: 2.1.5
Precedence: list
List-Id: Axiom Developers <axiom-developer.nongnu.org>
List-Unsubscribe: <http://lists.nongnu.org/mailman/listinfo/axiom-developer>, 
	<mailto:axiom-developer-request@nongnu.org?subject=unsubscribe>
List-Archive: <http://lists.nongnu.org/archive/html/axiom-developer>
List-Post: <mailto:axiom-developer@nongnu.org>
List-Help: <mailto:axiom-developer-request@nongnu.org?subject=help>
List-Subscribe: <http://lists.nongnu.org/mailman/listinfo/axiom-developer>,
	<mailto:axiom-developer-request@nongnu.org?subject=subscribe>
X-List-Received-Date: Wed, 29 Sep 2010 15:49:37 -0000

Greetings!  How do I enable lisp breaks on error?

Take care,

daly@axiom-developer.org writes:

> Axiom has been released for September 2010.
> The website has the sources and the ubuntu binary.
> Other binaries are being built and will appear when they pass testing.
>
> Tim Daly
>
> =============================================================================
>
> September 2010 Release
>
> This release concentrated on treeshaking the compiler code into book volume 9.
> Due to the complexity of this task it will take several releases to complete.
>
> Additional global changes included collecting all of the Axiom-related
> latex macros into axiom.sty and removing references to Aldor.
>
> Makefile 
>   always run help extractions in parallel
>   build src/input quietly
>
> books/bookvol0
>    remove references to aldor
>
> books/bookvol1
>    remove references to aldor
>
> books/bookvol5
>    expose StreamTensor, U32Vector, U32VectorPolynomialOperations
>    mark pure common lisp routines
>    merge ptrop, varini
>    move latex macros to axiom.sty
>    remove POLYVEC
>    remove compile, duplicated in vol9
>    remove memq
>    remove references to aldor
>    treeshake
>    remove $useNewParser
>
> books/bookvol7.1
>    move latex macros to axiom.sty
>    remove references to aldor
>    rewrite \pagehead to \pagetitle
>
> books/bookvol9
>    cross-reference functions and variables
>    move latex macros to axiom.sty
>    remove memq
>    treeshake the compiler code
>    remove $useNewParser
>
> books/bookvol10
>    move latex macros to axiom.sty
>    move GOPT0 from bookvol10.3
>
> books/bookvol10.1
>    move latex macros to axiom.sty
>
> books/bookvol10.2
>    move latex macros to axiom.sty
>
> books/bookvol10.3
>    add U32Vector, move GOPT0 from bookvol10.4
>    move latex macros to axiom.sty
>
> books/bookvol10.4
>    add StreamTensor, U32VectorPolynomialOperations
>    fix ScriptTensor regression test
>    move latex macros to axiom.sty
>    remove POLYVEC
>    update Chinese Remainder documentation
>
> books/bookvol10.5
>    move latex macros to axiom.sty
>
> books/bookvolbib
>    Parnas & Madey [PM95], Parnas & Jin [PJ10], GCL92, AS64, NIST10, RF94, 
>    Hamdy [Ham04], Steele [Ste90], Tim Lahey's Sage Integration Test Suite
>
> books/ps/
>    v103guessoptionfunctions0, v103u32vector, v104streamtensor,
>    v104u32vectorpolynomialoperations, v104u32vectorpolynomialoperations
>
> src/algebra
>    Makefile help and test for StreamTensor
>    Makefile help and test for U32Vector
>    Makefile remove references to aldor
>    Makefile test and help for POLYVEC
>    Makefile remove POLYVEC
>    Makefile add help and test for new algebra
>    Makefile handle case-insensitive MAC filesystem
>    Makefile remove new algebra scaffolding code
>
> src/doc
>    axiom.sty collect all script commands in one place
>    axiom.sty consolidate latex macros
>
> src/input
>    Makefile add guess.input, manuel.input, risch.input
>    guess.input test examples of the GUESS package
>    kamke3.input clean up broken tests
>    manuel.input add Manuel's integral to test suite
>    richlog300-391.input clean up broken tests
>    richtrig800-899.input clean up broken tests
>    risch.input illustrate the Risch algorithm
>    setcmd.input clean up broken tests
>    test.input clean up broken tests
>
> src/interp
>    Makefile merge varini
>    Makefile remove nspadaux, mark, pspad1, pspad2, ptrop, wi1, wi2
>    *.lisp remove memq
>    treeshake compiler -- br-con, cattable, compiler
>    remove nspadaux.lisp, mark.lisp, pspad1.lisp, pspad2.lisp, 
>    remove ptrop.lisp, wi1.lisp, wi2.lisp
>    add HTMLFormat code i-output.lisp, vmlisp.lisp
>
> src/scripts/tex
>    axiom.sty collect all script commands
>    axiom.sty consolidate latex macros
>
> src/axiom-website/
>    documentation.html add Knuth quote per W. Sit
>    download.html add debian, fedora, mandriva, opensuse, slackware, vector
>    download.html update ubuntu yum advice
>
>
>
>
> _______________________________________________
> Axiom-developer mailing list
> Axiom-developer@nongnu.org
> http://lists.nongnu.org/mailman/listinfo/axiom-developer
>
>
>
>

-- 
Camm Maguire			     		    camm@maguirefamily.org
==========================================================================
"The earth is but one country, and mankind its citizens."  --  Baha'u'llah


From MAILER-DAEMON Wed Sep 29 11:58:20 2010
Received: from mailman by lists.gnu.org with archive (Exim 4.43)
	id 1P0z2y-00006U-IV
	for mharc-axiom-developer@gnu.org; Wed, 29 Sep 2010 11:58:20 -0400
Received: from [140.186.70.92] (port=34763 helo=eggs.gnu.org)
	by lists.gnu.org with esmtp (Exim 4.43) id 1P0z2u-000068-Th
	for axiom-developer@nongnu.org; Wed, 29 Sep 2010 11:58:18 -0400
Received: from Debian-exim by eggs.gnu.org with spam-scanned (Exim 4.69)
	(envelope-from <daly@axiom-developer.org>) id 1P0z2s-0001kY-M0
	for axiom-developer@nongnu.org; Wed, 29 Sep 2010 11:58:16 -0400
Received: from vs338.rosehosting.com ([209.135.140.38]:59552
	helo=axiom-developer.org) by eggs.gnu.org with esmtp (Exim 4.69)
	(envelope-from <daly@axiom-developer.org>) id 1P0z2s-0001k6-IJ
	for axiom-developer@nongnu.org; Wed, 29 Sep 2010 11:58:14 -0400
Received: from [192.168.1.4] (dynamic-acs-72-23-235-203.zoominternet.net
	[72.23.235.203])
	by axiom-developer.org (8.12.8/8.12.8) with ESMTP id o8TFwAj5031623;
	Wed, 29 Sep 2010 10:58:11 -0500
Message-ID: <4CA38C16.9080503@axiom-developer.org>
Date: Wed, 29 Sep 2010 14:57:26 -0400
From: Tim Daly <daly@axiom-developer.org>
User-Agent: Mozilla/5.0 (Windows; U; Windows NT 5.1; en-US;
	rv:1.9.2.9) Gecko/20100825 Thunderbird/3.1.3
MIME-Version: 1.0
To: Camm Maguire <camm@maguirefamily.org>
Subject: Re: [Axiom-developer] Axiom September 2010 Release
References: <201009250647.o8P6lcFO022778@axiom-developer.org>
	<87iq1opn1w.fsf@maguirefamily.org>
In-Reply-To: <87iq1opn1w.fsf@maguirefamily.org>
Content-Type: text/plain; charset=ISO-8859-1; format=flowed
Content-Transfer-Encoding: 7bit
X-detected-operating-system: by eggs.gnu.org: Linux 2.6? (barebone, rare!)
Cc: axiom-developer@nongnu.org
X-BeenThere: axiom-developer@nongnu.org
X-Mailman-Version: 2.1.5
Precedence: list
List-Id: Axiom Developers <axiom-developer.nongnu.org>
List-Unsubscribe: <http://lists.nongnu.org/mailman/listinfo/axiom-developer>, 
	<mailto:axiom-developer-request@nongnu.org?subject=unsubscribe>
List-Archive: <http://lists.nongnu.org/archive/html/axiom-developer>
List-Post: <mailto:axiom-developer@nongnu.org>
List-Help: <mailto:axiom-developer-request@nongnu.org?subject=help>
List-Subscribe: <http://lists.nongnu.org/mailman/listinfo/axiom-developer>,
	<mailto:axiom-developer-request@nongnu.org?subject=subscribe>
X-List-Received-Date: Wed, 29 Sep 2010 15:58:18 -0000

  )set break break


On 9/29/2010 11:49 AM, Camm Maguire wrote:
> Greetings!  How do I enable lisp breaks on error?
>
> Take care,
>
> daly@axiom-developer.org writes:
>
>> Axiom has been released for September 2010.
>> The website has the sources and the ubuntu binary.
>> Other binaries are being built and will appear when they pass testing.
>>
>> Tim Daly
>>
>> =============================================================================
>>
>> September 2010 Release
>>
>> This release concentrated on treeshaking the compiler code into book volume 9.
>> Due to the complexity of this task it will take several releases to complete.
>>
>> Additional global changes included collecting all of the Axiom-related
>> latex macros into axiom.sty and removing references to Aldor.
>>
>> Makefile
>>    always run help extractions in parallel
>>    build src/input quietly
>>
>> books/bookvol0
>>     remove references to aldor
>>
>> books/bookvol1
>>     remove references to aldor
>>
>> books/bookvol5
>>     expose StreamTensor, U32Vector, U32VectorPolynomialOperations
>>     mark pure common lisp routines
>>     merge ptrop, varini
>>     move latex macros to axiom.sty
>>     remove POLYVEC
>>     remove compile, duplicated in vol9
>>     remove memq
>>     remove references to aldor
>>     treeshake
>>     remove $useNewParser
>>
>> books/bookvol7.1
>>     move latex macros to axiom.sty
>>     remove references to aldor
>>     rewrite \pagehead to \pagetitle
>>
>> books/bookvol9
>>     cross-reference functions and variables
>>     move latex macros to axiom.sty
>>     remove memq
>>     treeshake the compiler code
>>     remove $useNewParser
>>
>> books/bookvol10
>>     move latex macros to axiom.sty
>>     move GOPT0 from bookvol10.3
>>
>> books/bookvol10.1
>>     move latex macros to axiom.sty
>>
>> books/bookvol10.2
>>     move latex macros to axiom.sty
>>
>> books/bookvol10.3
>>     add U32Vector, move GOPT0 from bookvol10.4
>>     move latex macros to axiom.sty
>>
>> books/bookvol10.4
>>     add StreamTensor, U32VectorPolynomialOperations
>>     fix ScriptTensor regression test
>>     move latex macros to axiom.sty
>>     remove POLYVEC
>>     update Chinese Remainder documentation
>>
>> books/bookvol10.5
>>     move latex macros to axiom.sty
>>
>> books/bookvolbib
>>     Parnas&  Madey [PM95], Parnas&  Jin [PJ10], GCL92, AS64, NIST10, RF94,
>>     Hamdy [Ham04], Steele [Ste90], Tim Lahey's Sage Integration Test Suite
>>
>> books/ps/
>>     v103guessoptionfunctions0, v103u32vector, v104streamtensor,
>>     v104u32vectorpolynomialoperations, v104u32vectorpolynomialoperations
>>
>> src/algebra
>>     Makefile help and test for StreamTensor
>>     Makefile help and test for U32Vector
>>     Makefile remove references to aldor
>>     Makefile test and help for POLYVEC
>>     Makefile remove POLYVEC
>>     Makefile add help and test for new algebra
>>     Makefile handle case-insensitive MAC filesystem
>>     Makefile remove new algebra scaffolding code
>>
>> src/doc
>>     axiom.sty collect all script commands in one place
>>     axiom.sty consolidate latex macros
>>
>> src/input
>>     Makefile add guess.input, manuel.input, risch.input
>>     guess.input test examples of the GUESS package
>>     kamke3.input clean up broken tests
>>     manuel.input add Manuel's integral to test suite
>>     richlog300-391.input clean up broken tests
>>     richtrig800-899.input clean up broken tests
>>     risch.input illustrate the Risch algorithm
>>     setcmd.input clean up broken tests
>>     test.input clean up broken tests
>>
>> src/interp
>>     Makefile merge varini
>>     Makefile remove nspadaux, mark, pspad1, pspad2, ptrop, wi1, wi2
>>     *.lisp remove memq
>>     treeshake compiler -- br-con, cattable, compiler
>>     remove nspadaux.lisp, mark.lisp, pspad1.lisp, pspad2.lisp,
>>     remove ptrop.lisp, wi1.lisp, wi2.lisp
>>     add HTMLFormat code i-output.lisp, vmlisp.lisp
>>
>> src/scripts/tex
>>     axiom.sty collect all script commands
>>     axiom.sty consolidate latex macros
>>
>> src/axiom-website/
>>     documentation.html add Knuth quote per W. Sit
>>     download.html add debian, fedora, mandriva, opensuse, slackware, vector
>>     download.html update ubuntu yum advice
>>
>>
>>
>>
>> _______________________________________________
>> Axiom-developer mailing list
>> Axiom-developer@nongnu.org
>> http://lists.nongnu.org/mailman/listinfo/axiom-developer
>>
>>
>>
>>


From MAILER-DAEMON Thu Sep 30 09:20:53 2010
Received: from mailman by lists.gnu.org with archive (Exim 4.43)
	id 1P1J49-0007Dp-TK
	for mharc-axiom-developer@gnu.org; Thu, 30 Sep 2010 09:20:53 -0400
Received: from [140.186.70.92] (port=57822 helo=eggs.gnu.org)
	by lists.gnu.org with esmtp (Exim 4.43) id 1P1J44-0007CK-3F
	for axiom-developer@nongnu.org; Thu, 30 Sep 2010 09:20:52 -0400
Received: from Debian-exim by eggs.gnu.org with spam-scanned (Exim 4.69)
	(envelope-from <hebisch@math.uni.wroc.pl>) id 1P1J3x-0005y7-PX
	for axiom-developer@nongnu.org; Thu, 30 Sep 2010 09:20:47 -0400
Received: from hera.math.uni.wroc.pl ([156.17.86.1]:44963)
	by eggs.gnu.org with esmtp (Exim 4.69)
	(envelope-from <hebisch@math.uni.wroc.pl>) id 1P1J3x-0005x6-GE
	for axiom-developer@nongnu.org; Thu, 30 Sep 2010 09:20:41 -0400
Received: from hebisch by hera.math.uni.wroc.pl with local (Exim 4.69)
	(envelope-from <hebisch@math.uni.wroc.pl>) id 1P1J4Z-0000Rp-DS
	for axiom-developer@nongnu.org; Thu, 30 Sep 2010 15:21:19 +0200
Subject: Re: [Axiom-developer] wrong sign in limit?
To: axiom-developer@nongnu.org
Date: Thu, 30 Sep 2010 15:21:19 +0200 (CEST)
In-Reply-To: <4C9CF1ED.7020806@uni-bremen.de>
X-Mailer: ELM [version 2.5 PL8]
MIME-Version: 1.0
Content-Type: text/plain; charset=us-ascii
Content-Transfer-Encoding: 7bit
Message-Id: <E1P1J4Z-0000Rp-DS@hera.math.uni.wroc.pl>
From: Waldek Hebisch <hebisch@math.uni.wroc.pl>
X-detected-operating-system: by eggs.gnu.org: GNU/Linux 2.6 (newer, 3)
X-BeenThere: axiom-developer@nongnu.org
X-Mailman-Version: 2.1.5
Precedence: list
List-Id: Axiom Developers <axiom-developer.nongnu.org>
List-Unsubscribe: <http://lists.nongnu.org/mailman/listinfo/axiom-developer>, 
	<mailto:axiom-developer-request@nongnu.org?subject=unsubscribe>
List-Archive: <http://lists.nongnu.org/archive/html/axiom-developer>
List-Post: <mailto:axiom-developer@nongnu.org>
List-Help: <mailto:axiom-developer-request@nongnu.org?subject=help>
List-Subscribe: <http://lists.nongnu.org/mailman/listinfo/axiom-developer>,
	<mailto:axiom-developer-request@nongnu.org?subject=subscribe>
X-List-Received-Date: Thu, 30 Sep 2010 13:20:52 -0000

Dr.-Ing. Ingo D. Rullhusen wrote:
> 
> Hello,
> 
> perhaps i have found a bug in the limit operation. I am using
> 
>                         AXIOM Computer Algebra System
>                          Version: Axiom (July 2010)
>               Timestamp: Wednesday August 18, 2010 at 11:16:49
> 
> on a x86_64 with linux. There seems to be an error in the sign, if you
> compare the limit result with the de L'Hospital one.
> 
> 
> Z:=sqrt(a^2+h^2)-a
>  
> N:=a*sqrt(a^2+h^2) - a^2 - h^2
> 
> limit( Z/N, h=0 )
>  
> 
>         1
>    (3)  -
>         a
>                         Type: Union(OrderedCompletion Expression
> Integer,...)
> 
> -- this seems to be wrong!
> -- because de L'Hospital says
> 
> dZ:=differentiate( Z, h, 2 )
>  
> dN:=differentiate( N, h, 2 )
>  
> eval(dZ/dN,h=0)
>  
> 
>                1
>    (6)  - ----------
>             +--+
>             | 2
>           2\|a   - a
> 
> (see sign)

Square root has two branches, the result depends on which branch
you choose.  Limit decided that sqrt(a^2) = -a and the result
is correct for this branch (and wrong for the other branch).
de L'Hospital rule can be used only when limit on N is 0, which
hold when sqrt(a^2) = a, but the result is wrong for the other
branch.  If you want result which is correct for both branches
you need something like:

-sqrt(a^2)/a^2

but such expressions cause trouble so limit decided just to pick
one branch...

-- 
                              Waldek Hebisch
hebisch@math.uni.wroc.pl 


