\start
Date: Sun, 02 Aug 2009 21:41:28 -0400
From: Tim Daly
To: Eduardo Cavazos
Subject: Re: build error

Eduardo,

Where did you get your sources?
One really useful piece of information is available from git. Type:

    git log | head

(that is, git log piped into the 'head' function) at the axiom
top level directory. This will give me the eact commit
information for your sources so I can reproduce the error.

Tim


Eduardo Cavazos wrote:
> Hi!
>
> The build finishes on my system with the following output:
>
> ----------------------------------------------------------------------
> make[4]: Entering directory `/src/axiom/int/input'
> 0 starting algebra regression testing
> make[5]: Entering directory `/src/axiom/int/input'
> make[5]: *** No rule to make target `Any.regress', needed by `all'.  
> Stop.
> make[5]: Leaving directory `/src/axiom/int/input'
> make[4]: *** [all] Error 2
> make[4]: Leaving directory `/src/axiom/int/input'
> make[3]: *** [regress] Error 2
> make[3]: Leaving directory `/src/axiom/src/input'
> make[2]: *** [inputdir] Error 2
> make[2]: Leaving directory `/src/axiom/src'
> make[1]: *** [srcdir] Error 2
> make[1]: Leaving directory `/src/axiom'
> make: *** [all] Error 2
> /src/axiom #
> ----------------------------------------------------------------------
>
> Any guesses?
>
> I'm on a Linux system with gcc 4.2.2. Let me know if you need more 
> particular information.

\start
Date: Mon, 03 Aug 2009 03:44:25 -0500
From: Eduardo Cavazos
To: list
Subject: Re: build error


> Eduardo,
> 
> Where did you get your sources?
> One really useful piece of information is available from git. Type:
> 
>    git log | head
> 
> (that is, git log piped into the 'head' function) at the axiom
> top level directory. This will give me the eact commit
> information for your sources so I can reproduce the error.
> 
> Tim

Hi Tim!

Here ya go:

----------------------------------------------------------------------
/src/axiom # git log | head
commit 10fb9b914bd34291c59928218c3a6d22d19cbcf8
Merge: 0164fcc... f3a6506...
Author: Tim Daly
Date:   Sun Jul 26 02:34:50 2009 -0400

     Merge branch 'master' of /research/silver

commit f3a65067e5ffbb1fb23a09b589b7a4f8b20b138e
Author: Tim Daly
Date:   Sun Jul 26 02:10:20 2009 -0400
/src/axiom #
----------------------------------------------------------------------

\start
Date: Mon, 03 Aug 2009 05:15:18 -0400
From: Tim Daly
To: Eduardo Cavazos
Subject: Re: build error (resend)

Tim Daly wrote:
> Eduardo,
>
> I'm building exactly that version (that's a truly great
> feature of git... I know exactly which sources from the
> hash code).
>
> Try building with the following command line:
>
> make parallel
>
> and let me know if it succeeds. Axiom builds are
> now done in parallel but I may have missed a path.
> I'll try the same thing here.
>
> Tim
>
>
> Eduardo Cavazos wrote:
>> Tim Daly wrote:
>>
>>> Eduardo,
>>>
>>> Where did you get your sources?
>>> One really useful piece of information is available from git. Type:
>>>
>>>    git log | head
>>>
>>> (that is, git log piped into the 'head' function) at the axiom
>>> top level directory. This will give me the eact commit
>>> information for your sources so I can reproduce the error.
>>>
>>> Tim
>>
>> Hi Tim!
>>
>> Here ya go:
>>
>> ----------------------------------------------------------------------
>> /src/axiom # git log | head
>> commit 10fb9b914bd34291c59928218c3a6d22d19cbcf8
>> Merge: 0164fcc... f3a6506...
>> Author: Tim Daly
>> Date:   Sun Jul 26 02:34:50 2009 -0400
>>
>>     Merge branch 'master' of /research/silver
>>
>> commit f3a65067e5ffbb1fb23a09b589b7a4f8b20b138e
>> Author: Tim Daly
>> Date:   Sun Jul 26 02:10:20 2009 -0400
>> /src/axiom #
>> ----------------------------------------------------------------------

\start
Date: Mon, 03 Aug 2009 06:52:47 -0400
From: Tim Daly
To: Eduardo Cavazos
Subject: build error - verified -- bug 7202

Eduardo,

I've confirmed that the problem with Any.regress is due to the
recent parallel build changes. This is bug 7202.

For now, type
   make parallel
and everything should work. I'll fix this shortly.

Tim

Tim Daly wrote:
> Tim Daly wrote:
>> Eduardo,
>>
>> I'm building exactly that version (that's a truly great
>> feature of git... I know exactly which sources from the
>> hash code).
>>
>> Try building with the following command line:
>>
>> make parallel
>>
>> and let me know if it succeeds. Axiom builds are
>> now done in parallel but I may have missed a path.
>> I'll try the same thing here.
>>
>> Tim
>>
>>
>> Eduardo Cavazos wrote:
>>> Tim Daly wrote:
>>>
>>>> Eduardo,
>>>>
>>>> Where did you get your sources?
>>>> One really useful piece of information is available from git. Type:
>>>>
>>>>    git log | head
>>>>
>>>> (that is, git log piped into the 'head' function) at the axiom
>>>> top level directory. This will give me the eact commit
>>>> information for your sources so I can reproduce the error.
>>>>
>>>> Tim
>>>
>>> Hi Tim!
>>>
>>> Here ya go:
>>>
>>> ----------------------------------------------------------------------
>>> /src/axiom # git log | head
>>> commit 10fb9b914bd34291c59928218c3a6d22d19cbcf8
>>> Merge: 0164fcc... f3a6506...
>>> Author: Tim Daly
>>> Date:   Sun Jul 26 02:34:50 2009 -0400
>>>
>>>     Merge branch 'master' of /research/silver
>>>
>>> commit f3a65067e5ffbb1fb23a09b589b7a4f8b20b138e
>>> Author: Tim Daly
>>> Date:   Sun Jul 26 02:10:20 2009 -0400
>>> /src/axiom #
>>> ----------------------------------------------------------------------

\start
Date: Mon, 03 Aug 2009 21:54:32 -0400
From: Tim Daly
To: Eduardo Cavazos
Subject: Parallel make is now the default. fixes bug 7202

Eduardo,

The lastest release is up on the servers. Do
  git pull
to the the changes. You should only need to type
  make
and Axiom should do a parallel make. This fixes bug 7202.

Tim

Tim Daly wrote:
> Eduardo,
>
> I've confirmed that the problem with Any.regress is due to the
> recent parallel build changes. This is bug 7202.
>
> For now, type
>   make parallel
> and everything should work. I'll fix this shortly.
>
> Tim
>
> Tim Daly wrote:
>> Tim Daly wrote:
>>> Eduardo,
>>>
>>> I'm building exactly that version (that's a truly great
>>> feature of git... I know exactly which sources from the
>>> hash code).
>>>
>>> Try building with the following command line:
>>>
>>> make parallel
>>>
>>> and let me know if it succeeds. Axiom builds are
>>> now done in parallel but I may have missed a path.
>>> I'll try the same thing here.
>>>
>>> Tim
>>>
>>>
>>> Eduardo Cavazos wrote:
>>>> Tim Daly wrote:
>>>>
>>>>> Eduardo,
>>>>>
>>>>> Where did you get your sources?
>>>>> One really useful piece of information is available from git. Type:
>>>>>
>>>>>    git log | head
>>>>>
>>>>> (that is, git log piped into the 'head' function) at the axiom
>>>>> top level directory. This will give me the eact commit
>>>>> information for your sources so I can reproduce the error.
>>>>>
>>>>> Tim
>>>>
>>>> Hi Tim!
>>>>
>>>> Here ya go:
>>>>
>>>> ----------------------------------------------------------------------
>>>> /src/axiom # git log | head
>>>> commit 10fb9b914bd34291c59928218c3a6d22d19cbcf8
>>>> Merge: 0164fcc... f3a6506...
>>>> Author: Tim Daly
>>>> Date:   Sun Jul 26 02:34:50 2009 -0400
>>>>
>>>>     Merge branch 'master' of /research/silver
>>>>
>>>> commit f3a65067e5ffbb1fb23a09b589b7a4f8b20b138e
>>>> Author: Tim Daly
>>>> Date:   Sun Jul 26 02:10:20 2009 -0400
>>>> /src/axiom #
>>>> ----------------------------------------------------------------------

\start
Date: Tue, 04 Aug 2009 07:35:26 -0500
From: Gabriel Dos Reis
To: list
Subject: PLMMS 2009: Call for Participation


Apologies for duplicates.
---


      The ACM SIGSAM 2009 International Workshop on Programming
	     Languages for Mechanized Mathematics Systems

		   Munich, Germany; August 21, 2009

   The ACM SIGSAM 2009 International Workshop on Programming Languages
   for Mechanized Mathematics Systems will be co-located with TPHOLs 2009.

   Registration for the PLMMS 2009 workshop is now open.

Important Dates and Deadlines

     * Abstract submission: Now closed
     * Submission deadline: Now closed
     * Author notification: June 22, 2009
     * Camera ready papers: July 15, 2009
     * Registration: Open, http://tphols.in.tum.de/fee.html
     * Workshop: August 21, 2009

General Information

   The scope of this workshop is at the intersection of programming
   languages and mechanized mathematics systems. The latter category
   subsumes present-day computer algebra systems, interactive proof
   assistants, and automated theorem provers, all heading towards
   fully integrated mechanized mathematical assistants.

Tentative Programme

Friday, August 21

   09:00-10:00 Invited talk (Session Chair:)
               Georges Gonthier. Ssreflect: Structured Scripting for
                                    Higher-Order Theorem Proving
=09=09=09=09
   10:00-10:40 Coffee break

   10:40-12:10 Session 1 (Session Chair:)
               Paulo F. Silva, Joost Visser, Jose Oliveira. Galois: A
               Language for Proofs Using Galois Connections and Fork Algebr=
as

               Makarius Wenzel. Parallel Proof Checking in Isabelle/Isar

               Andrea Asperti, Wilmer Riccioti, Claudio Sacerdoti
                         Coen, Enrico Tassi. A New Type for Tactics

   12:10-13:40 Lunch

   13:40-14:40 Invited tutorial (Session Chair:)

               Gabriel Dos Reis. OpenAxiom: A Categorial Platform for
                                    Scientific Computation

   14:40-15:10 Coffee Break

   15:10-16:10 Session 2 (Session Chair:)

               Joe Hurd. OpenTheory: Package Management for Higher
                            Order Logic Theories

               Johannes Holzl. Proving Inequalities Over Reals With
                       Computation in Isabelle/HOL

   16:10-17:00 Business meeting

   Evening     Social Event

Links

     * http://plmms09.cs.tamu.edu/, the PLMMS 2009 workshop web site
     * http://tphols.in.tum.de/, the THOPLs 2009 conference web site
     * http://tphols.in.tum.de/fee.html, registration site

Programme Committee

     * Clemens Ballarin, aicas GmbH
     * Gabriel Dos Reis, Texas A&M University (Co-Chair)
     * Jean-Christophe Filli=E2tre, CNRS Universit=E9 Paris Sud
     * Predrag Janicic, University of Belgrade
     * Jaakko J=E4rvi, Texas A&M University
     * Florina Piroi, Johannes Kepler University
     * Laurent Th=E9ry, INRIA Sophia Antipolis (Co-Chair)
     * Makarius Wenzel, Technische Universit=E4t M=FCnchen

\start
Date: Tue, 4 Aug 2009 20:48:07 +0200
From: Michael Becker
To: list
Subject: unexpected  behaviour of normalize(1-(cos(x))^2)

    Hi,


   Is this (30)  the expected bevaviour of 'normalize' ??


(29) -> normalize ((sin(x))^2+(cos(x))^2)
(29) ->
   (29)  1
                                                     Type: Expression =
Integer



(30) -> normalize (1-(cos(x))^2)
(30) ->
                     x 2
                4tan(-)
                     2
   (30)  ----------------------
             x 4        x 2
         tan(-)  + 2tan(-)  + 1
             2          2
                                                     Type: Expression Integer

\start
Date: Wed, 05 Aug 2009 00:41:12 -0500
From: Eduardo Cavazos
To: list
Subject: Re: Parallel make is now the default. fixes bug 7202

Tim Daly wrote:

> Eduardo,
> 
> The lastest release is up on the servers. Do
>  git pull
> to the the changes. You should only need to type
>  make
> and Axiom should do a parallel make. This fixes bug 7202.
> 
> Tim

Thanks for the prompt support Tim!

Here's where I'm at...

Git says I'm up to date:

~ # cd /src/axiom/
/src/axiom # git pull
Already up-to-date.

I *think* I'm pulling from the right repo:

/src/axiom # git remote show origin
* remote origin
   URL: git://github.com/daly/axiom.git
   Remote branch merged with 'git pull' while on branch master
     master
   Tracked remote branch
     master

Finally, here's the git log again:

/src/axiom # git log | head
commit 10fb9b914bd34291c59928218c3a6d22d19cbcf8
Merge: 0164fcc... f3a6506...
Author: Tim Daly
Date:   Sun Jul 26 02:34:50 2009 -0400

     Merge branch 'master' of /research/silver

commit f3a65067e5ffbb1fb23a09b589b7a4f8b20b138e
Author: Tim Daly
Date:   Sun Jul 26 02:10:20 2009 -0400

However, I ran into what looks like the same error (sorry, I don't have 
it handy). So, what I'll do is try to rebuild from a fresh checkout 
(I'll just clone the local repo) and see how that goes.

Should I have done a "make clean" or anything like that before 
continuing the build?

Side comment/feedback. This build is running on a single-core system so 
'make parallel' reaaaaally causes the system to slow down. :-) Perhaps 
eventually the build process can consider how many cores are available 
and parallelize accordingly. I.e. 4-cores, four "threads". No biggie, 
just mentioning the idea.

\start
Date: Wed, 05 Aug 2009 09:59:59 -0400
From: Tim Daly
To: Michael Becker
Subject: Re: unexpected behaviour of normalize(1-(cos(x))^2)

Michael,

Trig identity substitutions are somewhat problematic in Axiom.
See the src/input/schaum* files for examples.

If the subexpression (1-cos(x)^2) occurs in your expression E you can write:

   sinrule:=rule((1-cos(x)^2) == sin(x)^2)

and then use this rule for your expression E thus

  sinrule(E)

Axiom will not derive several of the trig identities from scratch.

In your expression we have something of the form
    (4a^2) / (a^2 + 1)^2    where a = tan(x/2)
so Axiom needs to show that
   (a^2+1)^2 != 0
   (a^2+1) != 0
   a^2 != -1
   a != i
or, by back-substitution
  tan(x/2) != i
which it does not conclude automatically, even though this
is clearly true in the domain Expression(Integer).



Michael Becker wrote:
>     Hi,
>
>
>    Is this (30)  the expected bevaviour of 'normalize' ??
>
>
> (29) -> normalize ((sin(x))^2+(cos(x))^2)
> (29) ->
>    (29)  1
>                                                      Type: Expression Integer
>
>
>
> (30) -> normalize (1-(cos(x))^2)
> (30) ->
>                      x 2
>                 4tan(-)
>                      2
>    (30)  ----------------------
>              x 4        x 2
>          tan(-)  + 2tan(-)  + 1
>              2          2
>                                                      Type: Expression Integer
\start
Date: Mon, 10 Aug 2009 06:07:37 -0500
From: Tim Daly
To: list
Subject: depsys reduction

The current goal set involves reducing the complexity of the build
process for interpsys. There are two intermediate images to attack,
depsys and bootsys.

The depsys image is a macro-expansion environment for early compiles.
Most of the files in depsys have been merged into vmlisp.lisp.  The
exceptions are two aldor files which will be changed in a different way.

These changes are completed in patches 20090802-20080909.

The net result of this change is to reduce depsys to a single file
with the goal of removing this build completely. This will simplify
the whole build process.

The vmlisp file will be rewritten and restructured later.  It
contains, among other things, macros that I wrote to support the
VMLisp 370 mainframe lisp and the MACLISP lisp styles. When I joined
the Scratchpad group the system was being moved to Common Lisp. The
easy way to make this transition at the time was to provide a macro
layer, which I did.

The next logical step is to remove this macro layer and actually
rewrite the MACLISP/VMLisp primitives into common lisp.

The next actual step, however, is to remove the bootsys and the
boot code.

\start
Date: Mon, 3 Aug 2009 22:13:25 +0000
From: Alex Fabijanic
To: list
Subject: Axiom in Ubuntu 9.04

Dear axiom developers!
I've succesfully installed axiom deb package in Jaunty Jacalope from
Ubuntu repository. But axiom can't solve integrals:
"System error. Couldn't protect."
Elementary operations (such as number calculation) are working good.
Please answer me about probable cause. Alex

\start
Date: Mon, 10 Aug 2009 18:27:04 -0400
From: Tim Daly
To: Alex Fabijanic
Subject: Re: Axiom in Ubuntu 9.04

Alex,

Please do not use the debian version. It is old and out of date.
I cannot get the attention of anyone on the debian list to get it fixed.
See
http://axiom-developer.org/axiom-website/download.html
for newer versions

Tim

alex F wrote:
> Dear axiom developers!
> I've succesfully installed axiom deb package in Jaunty Jacalope from
> Ubuntu repository. But axiom can't solve integrals:
> "System error. Couldn't protect."
> Elementary operations (such as number calculation) are working good.
> Please answer me about probable cause. Alex

\start
Date: Thu, 13 Aug 2009 04:47:09 +0200
From: Michael Becker
To: Tim Daly
Subject: Re: unexpected behaviour of normalize(1-(cos(x))^2)

Am Mittwoch, 5. August 2009 15:59 schrieben Sie:
> Michael,
>
> Trig identity substitutions are somewhat problematic in Axiom.
> See the src/input/schaum* files for examples.
>
> If the subexpression (1-cos(x)^2) occurs in your expression E you can
> write:
>
>    sinrule:=rule((1-cos(x)^2) == sin(x)^2)
>
> and then use this rule for your expression E thus
>
>   sinrule(E)




    Tim,


    this does not always  work (see (6) and (7)) :



(1) -> )set mess auto off
(1) ->  sinrule:=rule((1-cos(x)^2) == sin(x)^2)
(1) ->
                2                   2
   (1)  - cos(x)  + %C + 1 == sin(x)  + %C
                        Type: RewriteRule(Integer,Integer,Expression =
Integer)
(2) -> f:= 1 - cos(x)^2
(2) ->
                2
   (2)  - cos(x)  + 1
                                                     Type: Expression =
Integer
(3) -> sinrule(f)
(3) ->
              2
   (3)  sin(x)
                                                     Type: Expression =
Integer
(4) -> sinrule(f+3)
(4) ->
                2
   (4)  - cos(x)  + 4
                                                     Type: Expression =
Integer
(5) -> sinrule(f+a)
(5) ->
              2
   (5)  sin(x)  + a
                                                     Type: Expression =
Integer
(6) -> sinrule (2*(f+a))
(6) ->
                 2
   (6)  - 2cos(x)  + 2a + 2
                                                     Type: Expression =
Integer
(7) -> sinrule (1/(f+a))
(7) ->
                 1
   (7)  - ---------------
                2
          cos(x)  - a - 1
                                                     Type: Expression =
Integer



   
    - Michael





>
> Axiom will not derive several of the trig identities from scratch.
>
> In your expression we have something of the form
>     (4a^2) / (a^2 + 1)^2    where a = tan(x/2)
> so Axiom needs to show that
>    (a^2+1)^2 != 0
>    (a^2+1) != 0
>    a^2 != -1
>    a != i
> or, by back-substitution
>   tan(x/2) != i
> which it does not conclude automatically, even though this
> is clearly true in the domain Expression(Integer).
>
> Michael Becker wrote:
> >     Hi,
> >
> >
> >    Is this (30)  the expected bevaviour of 'normalize' ??
> >
> >
> > (29) -> normalize ((sin(x))^2+(cos(x))^2)
> > (29) ->
> >    (29)  1
> >                                                      Type: =
Expression
> > Integer
> >
> >
> >
> > (30) -> normalize (1-(cos(x))^2)
> > (30) ->
> >                      x 2
> >                 4tan(-)
> >                      2
> >    (30)  ----------------------
> >              x 4        x 2
> >          tan(-)  + 2tan(-)  + 1
> >              2          2
> >                                                      Type: Expression
> > Integer

\start
Date: Thu, 13 Aug 2009 00:57:22 -0400
From: Bill Page
To: Michael Becker
Subject: Re: unexpected behaviour of normalize(1-(cos(x))^2)

Why not use removeCosSq ?

(3) -> removeCosSq(f)

              2
   (3)  sin(x)
                                                    Type: Expression(Intege=
r)
(4) -> removeCosSq(f+3)

              2
   (4)  sin(x)  + 3
                                                    Type: Expression(Intege=
r)
(5) -> removeCosSq(f+a)

              2
   (5)  sin(x)  + a
                                                    Type: Expression(Intege=
r)

(6) -> removeCosSq(2*(f+a))

               2
   (6)  2sin(x)  + 2a
                                                    Type: Expression(Intege=
r)
(7) -> removeCosSq(1/(f+a))

             1
   (7)  -----------
              2
        sin(x)  + a
                                                    Type: Expression(Intege=
r)

See also:

removeCoshSq
removeSinSq
removeSinhSq


Regards,
Bill Page.

On Wed, Aug 12, 2009 at 10:47 PM, Michael Becker wrote:
> Am Mittwoch, 5. August 2009 15:59 schrieben Sie:
>> Michael,
>>
>> Trig identity substitutions are somewhat problematic in Axiom.
>> See the src/input/schaum* files for examples.
>>
>> If the subexpression (1-cos(x)^2) occurs in your expression E you can
>> write:
>>
>>  sinrule:=rule((1-cos(x)^2) == sin(x)^2)
>>
>> and then use this rule for your expression E thus
>>
>>  sinrule(E)
>
>
>
>
>  Tim,
>
>
>  this does not always work (see (6) and (7)) :
>
>
>
> (1) -> )set mess auto off
> (1) -> sinrule:=rule((1-cos(x)^2) == sin(x)^2)
> (1) ->
>        2          2
>  (1) - cos(x) + %C + 1 == sin(x) + %C
>            Type: RewriteRule(Integer,=
Integer,Expression Integer)
> (2) -> f:= 1 - cos(x)^2
> (2) ->
>        2
>  (2) - cos(x) + 1
>                   =
        Type: Expression Integer
> (3) -> sinrule(f)
> (3) ->
>       2
>  (3) sin(x)
>                   =
        Type: Expression Integer
> (4) -> sinrule(f+3)
> (4) ->
>        2
>  (4) - cos(x) + 4
>                   =
        Type: Expression Integer
> (5) -> sinrule(f+a)
> (5) ->
>       2
>  (5) sin(x) + a
>                   =
        Type: Expression Integer
> (6) -> sinrule (2*(f+a))
> (6) ->
>         2
>  (6) - 2cos(x) + 2a + 2
>                   =
        Type: Expression Integer
> (7) -> sinrule (1/(f+a))
> (7) ->
>         1
>  (7) - ---------------
>        2
>     cos(x) - a - 1
>                   =
        Type: Expression Integer
>
>
>
>
>  - Michael
>
>
>
>
>
>>
>> Axiom will not derive several of the trig identities from scratch.
>>
>> In your expression we have something of the form
>>   (4a^2) / (a^2 + 1)^2  where a = tan(x/2)
>> so Axiom needs to show that
>>  (a^2+1)^2 != 0
>>  (a^2+1) != 0
>>  a^2 != -1
>>  a != i
>> or, by back-substitution
>>  tan(x/2) != i
>> which it does not conclude automatically, even though this
>> is clearly true in the domain Expression(Integer).
>>
>> Michael Becker wrote:
>> >   Hi,
>> >
>> >
>> >  Is this (30) the expected bevaviour of 'normalize' ??
>> >
>> >
>> > (29) -> normalize ((sin(x))^2+(cos(x))^2)
>> > (29) ->
>> >  (29) 1
>> >                  =
         Type: Expression
>> > Integer
>> >
>> >
>> >
>> > (30) -> normalize (1-(cos(x))^2)
>> > (30) ->
>> >           x 2
>> >         4tan(-)
>> >           2
>> >  (30) ----------------------
>> >       x 4    x 2
>> >     tan(-) + 2tan(-) + 1
>> >       2     2
>> >                  =
         Type: Expression
>> > Integer

\start
Date: Sat, 15 Aug 2009 20:42:30 -0400
From: Tim Daly
To: Steven Segletes
Subject: Re: [Axiom-mail] Exponential Integral

Steven,

I found your paper when I was looking for good fit values for E1(x).
I'd be happy to get your newer fit numbers.
I'll rewrite the Axiom code to use them.

Tim Daly
Axiom Lead Developer

Segletes, Steven (Civ,ARL/WMRD) wrote:
>
> Classification: _* UNCLASSIFIED*_
> Caveats: NONE
>  
> I was tickled to find mention of my exponential integral (E1(x)) 
> fitting algorithm at 
> _http://www.axiom-developer.org/axiom-website/patches/20080119.01.tpd.patch_ 
> .  As noted there, the accuracy of my published fit is only 4 decimal 
> places.  Years ago, I  extended that original fit to greater precision 
> (with more parameters of course), if it is of use to anyone, though I 
> have not published the extensions.  I do not subscribe to your 
> maillist, but am happy to provide the fits, should you provide your 
> request to me at Steven Segletes. 
>
> Regards,
> Steven Segletes
> U.S. Army Research Laboratory
> Aberdeen Proving Ground, MD 21005-5066

\start
Date: Sat, 15 Aug 2009 20:47:48 -0400
From: Tim Daly
To: Steven Segletes
Subject: Re: [Axiom-mail] Exponential Integral

If you know of any other fits to interesting functions I'd appreciate
a pointer to the work.  Finding high quality, numerically stable,
smooth piecewise fits is a challenge.

Segletes, Steven (Civ,ARL/WMRD) wrote:
>
> Classification: _* UNCLASSIFIED*_
> Caveats: NONE
>  
> I was tickled to find mention of my exponential integral (E1(x)) 
> fitting algorithm at 
> _http://www.axiom-developer.org/axiom-website/patches/20080119.01.tpd.patch_ 
> .  As noted there, the accuracy of my published fit is only 4 decimal 
> places.  Years ago, I  extended that original fit to greater precision 
> (with more parameters of course), if it is of use to anyone, though I 
> have not published the extensions.  I do not subscribe to your 
> maillist, but am happy to provide the fits, should you provide your 
> request to me at Steven Segletes. 

\start
Date: Tue, 18 Aug 2009 09:49:51 -0400
From: Steven Segletes
To: Tim Daly
Subject: RE: [Axiom-mail] Exponential Integral

Classification:  UNCLASSIFIED
Caveats: NONE

Tim,

I have attached a PDF writeup detailing the more accurate fits to E1(x),
relative to those presented in ARL-TR-1758.  My actual fits are to F1 =
=
[x exp(x) E1(x)] which spans a functional range from 0 to 1.  The best
accuracy I have been yet able to achieve, defined by limiting the value
of [(F1)_fit - F1]/F1 over the domain, is approximately 3.1E-07 with a
12-parameter fit, which unfortunately isn't quite to 32-bit
floating-point accuracy.  Nonetheless, the fit is not a piecewise fit,
but rather a single continuous function over the domain of nonnegative
x, which avoids some of the problems associated with piecewise domain
splicing.

Also, I have not, as yet, had the need to look at other functions that
could stand improved fitting, but you could always suggest a function to
get me thinking about it.

------_=_NextPart_001_01CA200A.C24E8660
	name="E1.pdf"
	filename="E1.pdf"

[encoded data, 93401 characters]
NEM+XQo+PgpzdGFydHhyZWYKNjczNTQKJSVFT0YK

------_=_NextPart_001_01CA200A.C24E8660--

\start
Date: Mon, 24 Aug 2009 04:04:28 -0400
From: Tim Daly
To: sage-devel@googlegroups.com
Subject: Re: [sage-devel] Re: Categories restart: call for reviewers
Cc: Nicolas Thiery

Nicolas

Do you have the inheritance graph of the categories?

Nicolas M. Thiery wrote:
> 	Hi John!
>
> On Sun, Aug 23, 2009 at 12:37:11PM +0100, John Cremona wrote:
>  
>> 2009/8/23 Nicolas M. Thiery Nicolas Thiery:
>>    
>>> So you can focus on looking at the code, doc, and tests in the files
>>> mentioned in:
>>>
>>>        http://trac.sagemath.org/sage_trac/wiki/CategoriesCategoriesReview
>>>
>>> possibly simply by browsing:
>>>
>>>        http://combinat.sagemath.org/hgwebdir.cgi/code/file/tip/sage/categories)
>>>
>>> and make sure they makes sense.
>>>      
>> I don't really know how this type of reviewing differs from normal
>> Sage reviewing of patches attached to trac tickets. 
>>    
>
> Just that you don't need to worry about the technical side (checking
> that the patch applies, pass test, ...). Also setting a positive
> review will be on the wiki page, file by file, rather than on the trac
> ticket. Finally, there is essentially no datastructure/algorithm issue
> in the current code (that will change!)
>
>  
>> But anyway, I looked at one file, finite_fields.py, and the method
>> _call_() in there looks wrong -- it raises an error whatever the
>> input, while the dicstring suggests that it is supposed to try using
>> __call__() and only change the error message if that does not work.
>>    
>
> It's the converse: __call__ tries to do generic stuff (like C(P)
> returns P if P is readilly a parent in the category C), and if that
> does not work calls _call_ (similar to what happens for __mul__ /
> _mul). I actually haven't changed that part; that's how it is in the
> original category framework.
>
> See the doc of _call_ / __call__ in Category (sage.category.category.py=
).
>
> Suggestions to improve the doctests to clarify this are welcome!
>
>  
>> If I have completely misunderstood what is going on, then I will
>> carry on ignoring the category activity until it is finished,
>>    
>
> Please don't!
>
>  
>> at which point I'm sure I will use it all the time!
>>    
>
> Which is why I want as much "early" feedback from you as possible :-)

\start
Date: Fri, 28 Aug 2009 18:53:17 -0400
From: Tim Daly
To: list,
Subject: Interpreter changes

All of the code for the interpreter proper is now pure lisp code.

The code is being rewritten and refactored in detail, documented,
and added to the interpreter book as literate code. The src/interp
directory should be gone in the near future.

There are three side effects of this change.

First, the VMLisp/MACLISP/etc legacy support will be rewritten into
pure common lisp.

Second, the boot translator code is no longer needed and will be
removed in the coming weeks along with the src/boot directory.

Third, the Aldor hooks are being removed since it is clear that
Aldor will never become free and open source. Instead of doing the
translation of this legacy code the effort will be spent to improve
and document the existing compiler.

\start
Date: Fri, 28 Aug 2009 16:53:36 -0700
From: Arthur Ralfs
To: list
Subject: Re: Interpreter changes

Tim Daly wrote:

> Third, the Aldor hooks are being removed since it is clear that
> Aldor will never become free and open source. Instead of doing the
> translation of this legacy code the effort will be spent to improve
> and document the existing compiler.
> 
> Tim
> 

Excellent.

