#9706: New Version of orthogonal Polynomials
-------------------------------+--------------------------------------------
   Reporter:  maldun           |       Owner:  burcin, maldun                   
       Type:  enhancement      |      Status:  needs_review                     
   Priority:  minor            |   Milestone:  sage-5.0                         
  Component:  symbolics        |    Keywords:  orthogonal polynomials, symbolics
     Author:  Stefan Reiterer  |    Upstream:  N/A                              
   Reviewer:                   |      Merged:                                   
Work_issues:                   |  
-------------------------------+--------------------------------------------
Changes (by maldun):

  * status:  needs_work => needs_review


Old description:

> The current implementation of orthogonal polynomials is just a wrapper
> around maxima. (see http://wiki.sagemath.org/symbolics/)
> There are the following improvements planed:
>
> -using of the pynac class for symbolic functions.
> -faster evaluation in general
> -evaluation of special values
> -mpmath for numeric evaluation

New description:

 The current implementation of orthogonal polynomials is just a wrapper
 around maxima. (see http://wiki.sagemath.org/symbolics/)
 This update holds the following changes:

   *using of the pynac class for symbolic functions.
   *faster evaluation in general
   *evaluation of special values
   *mpmath for numeric evaluation

 '''Remarks:'''
    *The current patch needs scipy-0.8. One has to install it before
 testing (see #9808 for spkg's and installation instructions)
    * Some of the old doctests in the old file don't work any more, due to
 coercion problems with pynac (see #9769)
    * Some doctests in Sage change, due to the fact that new BuiltIn
 functions are added. symbolic.random_test.py had output changes since the
 random expression creation changed of course. The tests in pynac.pyx also
 changed, but this has a strange behavior (see below).

--

Comment:

 @kcrisman thanks for paying attention. I added now an updated patch and
 extended instructions.

 the doctest changes in '''symbolic.random_tests.py''' are easy to explain:
 new functions are involved -> new random expressions. But I had to change
 {{{random_expr(50, nvars=3, coeff_generator=CDF.random_element) }}} to
 {{{random_expr(60, nvars=3, coeff_generator=CDF.random_element) }}} or
 else one gets an expression generated where a division through zero
 occours.

 As mentioned on sage-devel I repaired the doctests in
 '''symbolic.pynac.pyx''', the trick is to enlarge the range of the for
 loop:
 {{{for i in range(get_ginac_serial(), get_ginac_serial()+50):}}}
 changed to
 {{{for i in range(get_ginac_serial(), get_ginac_serial()+100):}}}
 now it works. My explaination: since we have new functions we have longer
 to search, and then we reach our goal. What I can not explain is, that it
 works, when I type it in by hand.

 All doctests pass now, so I think a review would be nice.

 -maldun

-- 
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/9706#comment:28>
Sage <http://www.sagemath.org>
Sage: Creating a Viable Open Source Alternative to Magma, Maple, Mathematica, 
and MATLAB

-- 
You received this message because you are subscribed to the Google Groups 
"sage-trac" group.
To post to this group, send email to [email protected].
To unsubscribe from this group, send email to 
[email protected].
For more options, visit this group at 
http://groups.google.com/group/sage-trac?hl=en.

Reply via email to