Permalink
Commits on Sep 11, 2011
  1. REL: set version to 0.10.0b1

    rgommers committed Sep 11, 2011
  2. ENH: special: Merge expit/logit implementation

    Merge branch 'chrisjordansquire/logit'
    pv committed Sep 11, 2011
  3. Added npy_log, npy_exp

    chrisjordansquire committed with pv Sep 1, 2011
  4. ENH: Templating magic

    chrisjordansquire committed with pv Aug 29, 2011
  5. ENH: Py3k compliant

    chrisjordansquire committed with pv Aug 25, 2011
  6. ENH: special: merge ellipk & agm accuracy improvements

    Merge branch 'tecki/elliptic-k-minus1'
    pv committed Sep 11, 2011
  7. TST: special: update ellipk test

    pv committed Sep 11, 2011
  8. Add tecki to THANKS.txt

    pv committed Sep 11, 2011
  9. corrected spelling mistake in docs

    Thanks to Collin Stocks
    tecki committed with pv Aug 17, 2011
  10. add tests for ellipkm1

    The new test values for ellipkm1 are calculated using
    mathematica, using the equality
    
    sqrt(m) K(1 - m) = K(1 - 1 / m)
    
    I used a different naming convention from the boost data,
    in order to show publicly that this data has nothing to do
    with it.
    tecki committed with pv Aug 11, 2011
  11. improved precision of special.agm

    The function special.agm suffered from the lack of precision
    in ellipk. By replacing it with ellipkm1, precision could
    be retained.
    
    Example
    y = agm(1, x)
    
    x        y before  y after
    0.1      0.425     0.425
    1e-10    0.064	   0.064
    1e-16    0.040     0.041
    1e-17    4.4e-309  0.039
    1e-20    4.4e-309  0.033
    tecki committed with pv Aug 11, 2011
  12. Add a more precise version of special.ellipk

    The complete elliptic integral of the first kind K,
    implemented as special.ellipk, diverges around 1.
    This has the unfortunate disadvantage that rounding
    errors in its parameters have tremendous, unwanted
    effects, as 1-1e-17 is rounded to 1 (float64),
    and thus ellipk reports infinte at this point.
    
    This is why the underlying cephes library defines
    a different funktion, K'(m) := K(1 - m). No
    rounding problems here, because the divergence has
    moved to zero, and significantly smaller numbers
    > 0 can be represented with float64.
    
    Unfortunately, scipy re-wrote the cephes library
    to undo just this, killing the precision actually
    supplied by cephes. This patch undoes that, but
    renames the function from ellipk to ellipkm1
    (defined as ellipkm1(m) == ellipk(1 - m)) for
    backward compatibility, and adds a new python
    function ellipk which works as the original function.
    tecki committed with pv Aug 11, 2011
  13. REF: rename io.hb.read/write to hb_read/hb_write, import them in io n…

    …amespace.
    
    Also update the hb_read/hb_write docstrings to note what is returnted and what
    features are currently supported.  Finally, change some imports to fit on one
    line and PEP8-ify the blank lines.
    rgommers committed Sep 11, 2011
  14. Merge branch 'ltidiscrete'

    rgommers committed Sep 11, 2011
  15. ENH: rewrite of the bilinear approximation of the cont2discrete funct…

    …ion.
    
    This implements a more general 'generalized bilinear transformation'.
    Also add:
      - fixes to test case.
      - alpha parameter in cont2discrete now handled for tf and zpk inputs.
      - dlsim now accepts 1-D input vectors.
      - Minor string and doc updates.
    ArmstrongJ committed with rgommers Aug 29, 2011
  16. DOC: various documentation updates supporting discrete LTI functional…

    …ity.
    ArmstrongJ committed with rgommers Aug 18, 2011