1. 05 Feb, 2010 1 commit
  2. 27 Dec, 2009 1 commit
  3. 21 Sep, 2009 1 commit
    • Liam Healy's avatar
      QR: More optional arguments, add tests form GSL · 44b9d868
      Liam Healy authored
      Changed required to optional arguments for several of the QR linear
      algebra functions, for returned quantities.  45 new tests added for QR
      that are translated from GSL's tests.  Matrices and solutions that are
      used by more than one of the linear algebra techniques are defined in
      matrix-generation.lisp.
      44b9d868
  4. 20 Sep, 2009 1 commit
    • Liam Healy's avatar
      LU linear algebra tests from GSL · 75f04669
      Liam Healy authored
      Added nine tests derived from the self-tests that GSL has for LU
      decomposition.  All tests except the complex matrix pass within the
      tolerances prescribed by GSL.
      75f04669
  5. 27 Apr, 2009 1 commit
    • Liam Healy's avatar
      Redeclare foreign-allocated variables called by reference · 9848c6ad
      Liam Healy authored
      Variables called using with-foreign-objects that are passed to GSL for
      the purposes of returning computed values (i.e., are not mentioned in
      the CL arglist) need to be declared (:pointer type) because they are
      calls by reference.  Everything now compiles though there are
      execution errors in the tests.
      9848c6ad
  6. 18 Mar, 2009 1 commit
    • Liam Healy's avatar
      LU remultiplication with permutation; unify permutation functions · b7074c43
      Liam Healy authored
      The LU remultiplication tests now permute the answer at the end so
      that it ends up with the original vector with the elements in the
      original order.  The permutation functions #'permute and
      #'permute-inverse now include all arguments, including marrays and
      pointers to C arrays.  Formerely, there were separate functions
      #'permute-vector and #'permute-vector-inverse.
      SBCL64: TOTAL: 1451 assertions passed, 5 failed, 0 execution errors.
      CCL64:  TOTAL: 1452 assertions passed, 4 failed, 0 execution errors.
      b7074c43
  7. 23 Feb, 2009 1 commit
    • Liam Healy's avatar
      Default values for LU functions; added test · c04147d7
      Liam Healy authored
      LU-solve and LU-decomposition now have optional arguments.  For
      LU-solve, the choice of whether to solve in-place or in a new vector
      depends on how this last argument is specified.  Because of the
      different call paths, it is possible that x is a vector or nil.  If it
      is a vector, it will change after the foreign call, which means it
      needs to be added to :outputs.  In order to prevent an error from the
      call to (setf cl-invalid) and (setf c-invalid) when x is nil, these
      two functions have been defined for class T to do nothing but return
      the value.  Also, copy-c-to-cl and copy-cl-to-c check whether the
      argument is a mobject, though this change isn't needed.  
      
      There are two additional tests in 'lu that do a solve and then
      remultiply the results with matrix-product-triangular.  In both SBCL
      and CCL, the answer comes out reversed.  When multiplied by
      matrix-product, it is correct.
      c04147d7
  8. 08 Jan, 2009 1 commit
    • Liam Healy's avatar
      Minor fixes; documentation update · 79a0f72c
      Liam Healy authored
      LU-decomposition now returns signum as a second value.  Remove unused
      *sizet-type*.  Update index.html with clbuild installation procedure
      provided by Glen Henshaw and how to run the tests.
      79a0f72c
  9. 31 Dec, 2008 1 commit
    • Liam Healy's avatar
      Redesign polynomial-solve, invert matrix · 4f581352
      Liam Healy authored
      The function #'polynomial-solve is now directly defined with defmfun
      with a single required argument and two optional arguments: the answer
      vector-complex-double-float and the workspace, both of which default
      to the right thing.  The returned is now the this marray instead of a
      list.  Simplify the definition of #'invert-matrix and return the
      matrix-double-float instead of the CL array.  All tests pass SBCL and
      CCL with trace.
      4f581352
  10. 27 Dec, 2008 1 commit
    • Liam Healy's avatar
      Updated tests; all pass · 64fea6a9
      Liam Healy authored
      Fixes to some definitions 
      lognormal, discrete: use #'mpointer
      lu: new syntax for making array
      quasi: order of arguments fixed in make-quasi-random-number-generator
      sorting: syntax to make combination
      All tests have been regenerated because syntax for arrays etc. has
      changed.  Visually inspected that results match previous check in.
      64fea6a9
  11. 26 Dec, 2008 2 commits
    • Liam Healy's avatar
      Full use of defmobject; elimination of letm · b09dc9d2
      Liam Healy authored
      All GSL objects are full CL objects with garbage collection.  There is
      no need to manually free the associated memory.  All objects may be
      assigned in the normal CL way, e.g. let, so there is not #'letm macro
      anymore.  Everything compiles cleanly in gsll, but nothing has been
      tested.
      b09dc9d2
    • Liam Healy's avatar
      New class foreign-array, new names: marray, make-marray, mobject.lisp · 7c66d57b
      Liam Healy authored
      New class foreign-array represents foreign (C) arrays, without GSL or
      mobject superclass.  New class name 'marray was 'gsl-data,
      representing arrays in GSL.  These objects are made by #'make-marray,
      formerely #'make-array*.  New file names: init/mobject.lisp was
      init/gsl-objects.lisp and data/marray.lisp was data/data.lisp.
      7c66d57b
  12. 08 Dec, 2008 1 commit
    • Liam Healy's avatar
      Array input and output declarations · 1a2ebe96
      Liam Healy authored
      A somewhat comprhensive attempt to make sure :inputs, :outputs, and
      :return arguments in defmfuns match what the function is expecting and
      will return, so that non-native implementations will work correctly
      for copying, and native ones will work correctly for pinning.  The
      :return argument now defaults to :outputs in most circumstances.
      1a2ebe96
  13. 01 Dec, 2008 1 commit
    • Liam Healy's avatar
      GSL array making with make-array* and #m · 83cd8040
      Liam Healy authored
      GSL arrays are garbage-collected objects, and so can be made like any
      other lisp object; they do not need to be make inside a wrapper form.
      #m is a macro which makes them; it expands to a #'make-array* call.
      Examples and usage ported to this new form.  Compiled only; not tested
      yet.  No effort to port to non-native implementation yet; tests as
      expanded in the tests/ directory are not ported yet.
      83cd8040
  14. 29 Nov, 2008 1 commit
  15. 25 Oct, 2008 1 commit
    • Liam Healy's avatar
      Port tests to save-test · 0343f3ce
      Liam Healy authored
      All remaning tests are now defined with save-test, and the
      LISP-UNIT:DEFINE-TEST has been eliminated from the definition files,
      except for arrays.
      0343f3ce
  16. 21 Sep, 2008 1 commit
    • Liam Healy's avatar
      Restoration of tests complete except for arrays · 4e6d5886
      Liam Healy authored
      All the lisp-unit tests defined in the old GSLL have been updated and
      will load by default, with all tests passing on Debian Lenny amd64 and
      GSL 1.11.
      TOTAL: 491 assertions passed, 0 failed, 0 execution errors.
      4e6d5886
  17. 12 Aug, 2008 1 commit
  18. 11 Aug, 2008 1 commit
    • Liam Healy's avatar
      LU decomposition ported to ffa · aa4aa130
      Liam Healy authored
      LU decomposition functions now support real and complex types, and the
      example has been tested and works:
      (invert-matrix
      	     (make-array* '(2 2) 'double-float :initial-contents '(1.0d0 2.0d0 3.0d0 4.0d0)))
      #2A((-1.9999999999999998d0 1.0d0)
          (1.4999999999999998d0 -0.49999999999999994d0))
      aa4aa130
  19. 25 Jul, 2008 4 commits
  20. 18 Feb, 2008 1 commit
  21. 04 Jul, 2006 1 commit
  22. 02 Jul, 2006 1 commit
  23. 29 Apr, 2006 1 commit
  24. 28 Apr, 2006 1 commit
  25. 27 Apr, 2006 1 commit