1. 12 Feb, 2011 1 commit
  2. 30 Jan, 2011 3 commits
    • Liam M. Healy's avatar
      Remove invert-matrix to Antik, optional argument for LU-decomposition · 60c1f52d
      Liam M. Healy authored
      The user-friendly ("higher") functions layered on GSLL are now in
      Antik's math-high system, so the antik directory has been removed, and
      linear-algebra.lisp is enhanced and in Antik. #'LU-invert now takes
      'inverse as an optional argument, with the sensible default.
      60c1f52d
    • Liam M. Healy's avatar
      Define gsll-tests as a system connection · 476d01ee
      Liam M. Healy authored
      Define gsll-tests as a system connection, rather than a separate ASDF
      system.  This means that if lisp-unit and GSLL are loaded (in either
      order), then the tests will automatically be available.
      476d01ee
    • Liam M. Healy's avatar
      In actual-array-class, don't replace real class names · e15aaed2
      Liam M. Healy authored
      In actual-array-class, if 'category is not one of 'vector, 'matrix or
      'both, then assume it's actually a class name and just return it
      without trying to make a class name.  This fixes a problem with
      e.g. elt+ where the second argument specializer is a class name of the
      scalar.
      e15aaed2
  3. 25 Jan, 2011 1 commit
  4. 13 Jan, 2011 1 commit
  5. 12 Jan, 2011 3 commits
  6. 10 Jan, 2011 4 commits
  7. 07 Jan, 2011 1 commit
    • Liam M. Healy's avatar
      Single-float tests have explicit single float numbers (start) · e6380d35
      Liam M. Healy authored
      Change literal floats in tests from an unspecified precision to
      explicit single floats, so that if *read-default-float-format* is
      something other than 'single-float the tests will still pass.  About
      half the failures now fixed.  Function #'cdot is returning a pointer
      instead of a numerical value.
      e6380d35
  8. 01 Jan, 2011 3 commits
    • Liam M. Healy's avatar
      Import symbols from gsll into antik-user · 5d648d16
      Liam M. Healy authored
      Import symbols from the gsll package into the antik-user package,
      managing conflicts.
      5d648d16
    • Liam M. Healy's avatar
      f4bf4c1e
    • Liam M. Healy's avatar
      Methods for arithmetic functions on foreign-arrays · 26271f71
      Liam M. Healy authored
      Depend on Antik instead of GSD.  Add new methods for Antik functions
      on foreign-arrays.  This permits algebra-style computation on arrays,
      e.g.,
      (in-package :antik-user)
      (setf grid:*default-grid-type* 'grid:foreign-array)
      (expt (* (rotate-3d :z #_15_deg) (rotate-3d :z #_45_deg)) -1)
      #m((0.5 0.8660254037844386 0.0) (-0.8660254037844385 0.5000000000000001 -0.0)
         (0.0 0.0 1.0))
      (rotate-3d :z #_-60_deg)
      #m((0.5000000000000001 0.8660254037844386 0.0)
         (-0.8660254037844386 0.5000000000000001 0.0) (0.0 0.0 1.0))
      26271f71
  9. 04 Dec, 2010 1 commit
    • Liam Healy's avatar
      Introduce values-unless-singleton for return values in defmfun expansion · 6010bde3
      Liam Healy authored
      New function values-unless-singleton will wrap the forms in 'values if
      there are more than one, otherwise it just returns the form.  This is
      used in the return value(s) in the defmfun expansion (by way of
      body-expand) so that if the returned form returns multiple values,
      they all come back from the defmfunned function (note that (values
      (values ...)) removes only returns the first value.
      6010bde3
  10. 28 Nov, 2010 1 commit
    • Liam Healy's avatar
      Eliminate calls to grid:complex-to-cl · 75ca6f4e
      Liam Healy authored
      Eliminate calls to grid:complex-to-cl because it has been eliminated;
      rely on fsbv:object to do the conversion.  Use new function
      #'fsbv:converter-defined-p.
      75ca6f4e
  11. 27 Nov, 2010 2 commits
  12. 25 Nov, 2010 1 commit
    • Liam Healy's avatar
      New macros #'maref and #'(setf maref) · 9becff94
      Liam Healy authored
      New macros #'maref and #'(setf maref) get and set elements of the GSL
      array directly from the mpointer.  These are not used by users, but in
      the make-funcallable-form expansion of callbacks when scalarsp = T is
      specified.  This fixes an error in the argument order in the (setf
      get-value) generic function that this replaces.
      9becff94
  13. 24 Nov, 2010 1 commit
    • Liam M. Healy's avatar
      New macro access-value · a2284097
      Liam M. Healy authored
      New macro access-value to call gsl_*_set and gsl_*_get, as yet unused.
      It will eventually replace get-value and (setf get-value) as used in
      reference-foreign-element, but has the arguments in the right order
      and because it's a macro, expands immediately to a
      cffi:foreign-funcall instead of doing a generic function dispatch on
      an eql spcializer of the class name.
      a2284097
  14. 14 Nov, 2010 2 commits
    • Liam Healy's avatar
      Revert the histogram change · a27a80a1
      Liam Healy authored
      Revert the last histogram change because GSD has been changed so that
      the gref compiler macro does not take effect except for foreign-array
      types.
      a27a80a1
    • Liam Healy's avatar
      Definition of grid:gref* for histogram · de942c32
      Liam Healy authored
      Define method grid:gref* for histograms; it was a definition of
      grid:gref, but with the compiler macro that turns that into
      grid:gref*, we need it to be a definition of grid:gref* instead.
      Removed the definition of +foreign-array-pointer+ and
      +foreign-array-type+ to the grid package.
      de942c32
  15. 12 Oct, 2010 1 commit
  16. 03 Oct, 2010 3 commits
  17. 01 Oct, 2010 1 commit
    • Sumant Oemrawsingh's avatar
      Added non-radix-2 tests for radix-2-sized vectors · 53fcec99
      Sumant Oemrawsingh authored
      Effectively, non-radix-2 algorithm tests on power-of-2 vector sizes are added.
      The result is, that some tests fail, and the cause at this point is not yet
      clear. The tests that fail, are all related to non-radix-2 algorithms working
      on power-of-2 vector sizes.
      
      Only tests that are missing, are the so-called bitreverse tests.
      53fcec99
  18. 15 Sep, 2010 3 commits
  19. 04 Sep, 2010 2 commits
  20. 02 Sep, 2010 1 commit
  21. 01 Sep, 2010 1 commit
    • Sumant Oemrawsingh's avatar
      Fix for random errors in tests · 04dfa37f
      Sumant Oemrawsingh authored
      The test functions did not properly initialize newly created vectors to
      contain unique values. This caused comparison of elements that are out of the
      stride to randomly succeed. The new test functions properly initialize the
      vectors just like is done in the GSL test functions. The result is that
      comparison of elements that are out of stride should fail consistently.
      
      This automatically means that the comparison should only check elements which
      are on stride, just like is done in GSL. This is not fixed yet, so the
      fast-fourier-transform unit tests are still left disabled for now.
      04dfa37f
  22. 23 Aug, 2010 1 commit
    • Liam Healy's avatar
      Fix FFT to actually do something, but comment out because of failures · b7859a14
      Liam Healy authored
      The FFT tests were testing nothing of consequence, just that NIL was
      equal to NIL, because it was looking at the result of (reset-urand)
      and then three non-existent multiple values.  This is now fixed, so
      that fft-*-result-check actually binds the results of the
      calculations.  Also, single-float tests are referenced to
      +sgl-epsilon+ and not +dbl-epsilon+.  However, for
      (all-fft-test-forms 9 3 (64 99)), many tests fail:
      FAST-FOURIER-TRANSFORM: 235 assertions passed, 55 failed.
      so the tests are commented out.
      b7859a14
  23. 21 Aug, 2010 1 commit
    • Liam Healy's avatar
      New size-vector-scalar replaces size-vector-real · 43f24411
      Liam Healy authored
      New function size-vector-scalar replaces size-vector-real, and gives
      the size as a double float if the array is real, and as the element
      type if it is complex, in order to be acceptable to the elt* interface
      to GSL's _scale routines.
      43f24411
  24. 19 Aug, 2010 1 commit