Re: Maxima bug on trigrat function



I would hope that the autoloading capability of the system
would take care of all of the "out-of-core" files if the
autoload properties were set up and the file names and locations
were correct.
  I know of no way to assure that the different files
do not tread on each others' toes, though there would probably
be warning messages for overlaying already-existing program
definitions.  
  Of some interest: trigrat is NOT in the commercial macsyma,
though I suspect the short file could, with a little
effort,  be loaded into it.  I tried it and the reader-dispatch
marker #$  to use the maxima reader  does not work in macsyma.

I'm also cc'ing Bill Gosper in case he wants to know that
there are actually people fussing with the public domain
maxima.
RJF


Camm Maguire wrote:
> 
> Subject: maxima: missing trigrat.lisp
> Package: maxima
> Version: 5.5-beta-4
> Severity: normal
> 
> Greetings, and thank you for your report.  I'm submitting this to the
> Debian Bug tracking system on your behalf to better keep track of this
> issue.  In general, there are many files in the source package which
> are not installed by default in the makefiles, and it appears that the
> installed list has gotten out of date with respect to the
> documentation.  I'm also ccing this message to the maxima list to
> enquire as to whether it would be better to delete the reference to
> trigrat from the documentation, or to include it in the installation.
> 
> Dr. Schelter, this issue has come up before with the mactex.lisp file
> in src/, and I just installed it by hand.  Is there a general wildcard
> file pattern specifying files that we could *safely* install and cover
> any references in the existing documentation?
> 
> Thanks!
> 
> Daniel Martins <dmartins@lcmi.ufsc.br> writes:
> 
> > Dear Camm,
> >
> > As you are the maintainer of the maxima package I am sending in a copy
> > of the problem I found using the testing version of maxima on my
> > debian potato.
> >
> > The function trigrat that is cited in the trigonometric section of the
> > info file with an example silply cannot be accesed
> >
> > In the source code of the maxima 5.4 it appears as trigrat.lisp in the
> > subdirectory maxima5-4/share and is different from the remaining
> > remaing function suffixes. Some of them are .lsp (not .lisp).
> >
> > I dont know how to overcome this problem and I am not a lisp
> > programmer (sorry!)
> >
> > Anyway I am describing below a sample session. Ate the end the
> > comments in the info file regarding trigrat
> >
> >
> > ***************************************************************
> >
> >
> > Thu 4:05pm ~ >maxima
> > GCL (GNU Common Lisp)  Version(2.3) Sun Aug 20 20:05:49 CEST 2000
> > Licensed under GNU Library General Public License
> > Contains Enhancements by W. Schelter
> > Maxima 5.4 Sun Aug 20 20:05:46 CEST 2000 (with enhancements by W. Schelter).
> > Licensed under the GNU Public License (see file COPYING)
> > (C1) trigrat(sin(3*a)/sin(a+%pi/3));
> >
> > Could not find `trigrat' using paths in FILE_SEARCH_LISP,SYSTEM
> >
> >  (combined values: [./###.{o,lsp,lisp},
> >
> > /usr/lib/maxima-5.4/{src,share1,sym}/###.o,
> >
> > /usr/lib/maxima-5.4/{src,share1,sym}/###.o,
> >
> > /usr/lib/maxima-5.4/{src,share1}/###.lisp, /usr/lib/maxima-5.4/{sym}/###.lsp,
> >
> > /usr/lib/maxima-5.4/src/../{src,share,share1,sharem}/foo.{mc,mac}] )
> > #0: TRIGRAT(?_l=[SIN(3*a)/SIN(a+%PI/3)])
> >  -- an error.  Quitting.  To debug this try DEBUGMODE(TRUE);)
> > (C2)
> >
> >
> >
> >
> > ***************************************************************
> >
> >  - Function: TRIGRAT (trigexp)
> >      gives a canonical simplifyed quasilinear form of a trigonometrical
> >      expression; trigexp is a rational fraction of several sin, cos or
> >      tan, the arguments of them are linear forms in some variables (or
> >      kernels) and %pi/n (n integer) with integer coefficients. The
> >      result is a simplifyed fraction with numerator and denominator
> >      linear in sin and cos.  Thus TRIGRAT linearize always when it is
> >      possible.(written by D. Lazard).
> >
> >           (c1) trigrat(sin(3*a)/sin(a+%pi/3));
> >
> >           (d1)                        sqrt(3) sin(2 a) + cos(2 a) - 1
> >
> >      Here is another example (for which the function was intended); see
> >      [Davenport, Siret, Tournier, Calcul Formel, Masson (or in english,
> >      Addison-Wesley), section 1.5.5, Morley theorem). Timings are on
> >      VAX 780.
> >
> >           (c4)   c:%pi/3-a-b;
> >
> >                                                  %pi
> >           (d4)                                 - b - a + ---
> >                                                   3
> >
> >           (c5)   bc:sin(a)*sin(3*c)/sin(a+b);
> >
> >                                    sin(a) sin(3 b + 3 a)
> >           (d5)                             ---------------------
> >                                         sin(b + a)
> >
> >           (c6)   ba:bc,c=a,a=c$
> >
> >           (c7)   ac2:ba^2+bc^2-2*bc*ba*cos(b);
> >
> >               2       2
> >                sin (a) sin (3 b + 3 a)
> >           (d7) -----------------------
> >                     2
> >                  sin (b + a)
> >
> >                                               %pi
> >              2 sin(a) sin(3 a) cos(b) sin(b + a - ---) sin(3 b + 3 a)
> >                                                3
> >            - --------------------------------------------------------
> >                                  %pi
> >                          sin(a - ---) sin(b + a)
> >                                   3
> >
> >                 2             2         %pi
> >              sin (3 a) sin (b + a - ---)
> >                                  3
> >            + ---------------------------
> >                    2     %pi
> >                 sin (a - ---)
> >                           3
> >
> >           (c9)   trigrat(ac2);
> >           Totaltime= 65866 msec.  GCtime= 7716 msec.
> >
> >           (d9)
> >           - (sqrt(3) sin(4 b + 4 a) - cos(4 b + 4 a)
> >
> >           - 2 sqrt(3) sin(4 b + 2 a)
> >
> >           + 2 cos(4 b + 2 a) - 2 sqrt(3) sin(2 b + 4 a) + 2 cos(2 b + 4 a)
> >
> >           + 4 sqrt(3) sin(2 b + 2 a) - 8 cos(2 b + 2 a) - 4 cos(2 b - 2 a)
> >
> >           + sqrt(3) sin(4 b) - cos(4 b) - 2 sqrt(3) sin(2 b) + 10 cos(2 b)
> >
> >           + sqrt(3) sin(4 a) - cos(4 a) - 2 sqrt(3) sin(2 a) + 10 cos(2 a)
> >
> >              - 9)/4
> >
> >
> >
> 
> --
> Camm Maguire                                            camm@enhanced.com
> ==========================================================================
> "The earth is but one country, and mankind its citizens."  --  Baha'u'llah