- Jan 19, 2010
-
-
Liam Healy authored
Correct default vector size in SV-solve for non-square matrices; patch by Evgeniy Zhemchugov.
-
Liam Healy authored
Define #'cholesky-invert using gsl_cholesky_invert introduced in GSL 1.12. Tests need to be generated. By Evgeniy Zhemchugov.
-
- Jan 17, 2010
-
-
Liam Healy authored
-
Russell Kliese authored
-
- Jan 03, 2010
-
-
ssh://repo.or.cz/srv/git/gsllLiam Healy authored
-
Liam Healy authored
-
- Jan 01, 2010
-
-
Liam Healy authored
-
- Dec 28, 2009
-
-
Liam M. Healy authored
-
- Dec 27, 2009
-
-
Liam Healy authored
Add copyright and license information to each file, include a copy of the GPL under documentation/COPYING.
-
- Dec 26, 2009
-
-
Liam Healy authored
Removed some marrays marked as :input when their values were not used. This failed in CCL because the arrays in those cases were created with no initial values, and no values were set, so that they were marked cl-invalid. Created #'make-marray-or-default to facilitate defaulting of marrays in defmfun, and marking validity appropriately. Test results: SBCL 64: TOTAL: 1627 assertions passed, 6 failed, 0 execution errors. CCL 64: TOTAL: 1620 assertions passed, 13 failed, 0 execution errors.
-
Liam Healy authored
Some symbols that were formerly in the GSLL package and now are in c-array need to imported from that package. This only affects non-native implementation and so was discovered in trying to compile in CCL. There are still too many test failures in CCL.
-
- Dec 25, 2009
-
-
Liam Healy authored
Explain in the documentation the dependence of GSLL on the GSD systems grid and c-array.
-
- Dec 24, 2009
-
-
Liam Healy authored
The grid subdirectory is now the c-array system within the GSD repository, as it has no GSL-dependent definitions and is potentially useful on its own for other C libraries. There is a bare minimum mention of GSD in the documentation, which needs to be enhanced as it is required for GSLL to compile/load.
-
- Dec 23, 2009
-
-
Liam Healy authored
Replace create-matrix and constant-matrix by using the grid functions grid:map-grid and grid:make-grid; eliminate create-vector because vectors can be generated with create-matrix and using nil for the second dimension argument. In the grid:make-grid-data method for make-marray, the first argument, which is the element type, needs to be extracted from the rest spec with a newly-changed grid:spec-scalar-p because compound types like (complex double-float) were not being recognized as element types.
-
- Dec 22, 2009
-
-
Liam Healy authored
The GSLL marrays now function correctly as instances of grid:grid; for example, (grid:transpose *m35*) works. All previous tests pass in SBCL.
-
Liam Healy authored
The grid definitions supply the slot 'grid:data which can be used for CL array storage for marrays. In order for a combination of zero dimension (k) to be made, there is an ignore-errors around the make-affi in initialize-instance :after for foreign-array.
-
- Dec 21, 2009
-
-
Liam Healy authored
Define foreign-array as subclass of grid:grid with access to elements through grid:gref.
-
Liam Healy authored
Moved the definition of the generic function copy-making-destination to copy objects while making a destination object, and the wrapper function #'copy for that and copy-to-destination to grid/ and c-array: package.
-
Liam Healy authored
The generic function definition for copy-to-destination is now in grid/ and the c-array package. This will make it accessible to non-GSLL foreign library interfaces, but the new location is a bit deceptive because a new method needn't be a subclass of foreign-array.
-
Liam Healy authored
Moved the defgeneric and methods for CL arrays and foreign-array for cl-array, maref and (setf maref) to a new file grid/element-reference.lisp. The old file data/maref.lisp remains with the GSL pointer methods.
-
Liam Healy authored
Continuing the separation of definitions useful for any libraries using C arrays, foreign-array.lisp and foreign-friendly.lisp have been moved to the grid subdirectory and placed in the c-array package. All tests pass on SBCL 64, TOTAL: 1627 assertions passed, 6 failed, 0 execution errors.
-
Liam Healy authored
Move number-conversion.lisp from init/ to grid/ and symbol-type portion of init/interface.lisp to grid/symbol-type.lisp. All tests pass TOTAL: 1627 assertions passed, 6 failed, 0 execution errors.
-
- Dec 08, 2009
-
-
Liam Healy authored
Continued separation of GSLL into two packages and eventually systems, the definitions in element-types.lisp are now in the c-array package.
-
Liam Healy authored
Eventually, the definitions used for dealing with types and foreign arrays that aren't specific to GSL will be split off into another system, so that it will be easier to build other foreign library interfaces using these definitions. As a start, the grid/ directory defines a package 'c-array and some of the type definitions formerly in init/. Everything compiles and tests as before, on SBCL TOTAL: 1627 assertions passed, 6 failed, 0 execution errors.
-
- Dec 06, 2009
-
-
Liam Healy authored
-
- Dec 01, 2009
-
- Nov 30, 2009
-
-
Liam M. Healy authored
Function #'number-class finds the CL class from the number type, replacing the previous function #'number-class-from-type in init/defmfun-array.lisp. It is a bit more general.
- Nov 26, 2009
-
-
Sumant Oemrawsingh authored
fft-frequency-vector can be used to make a vector that contains the actual sample frequencies (so you don't have to guess or write your own). fft-shift and fft-inverse-shift shift the zero frequency to the center (for easy plotting), and the inverse (for proper inverse FFT). These functions are comparable to fftfreq, fftshift and ifftshift in numpy (if anyone cares...).
-
- Nov 20, 2009
-
-
Liam Healy authored
Two new functions to write tests in fast-fourier-transforms/example.lisp, #'make-real-noise-fft-test and #'generate-fft-tests. The first one generates, given a size and stride, a pair of forward and inverse real tests, together with the answers. The second one loops over a range of sizes and strides, collecting the sets, and wrapping with a lisp-unit:define-test. Since they rely on #'make-test in test-unit/generate-tests.lisp, they won't normally be compiled. No test file generated yet.
-
- Nov 19, 2009
-
-
Liam Healy authored
The expanded FFT tests with inverse and backward transformations have simplified code. Since the backward and inverse cases must be divided by the length, the complex test will not work in GSL versions prior to 1.12, nor if FSBV is missing. In these cases, the vector is computed by converting to a CL array, dividing by the length, and then converting back.
-
- Nov 18, 2009
-
-
Liam M. Healy authored
The functions elt* and elt+ which used with a scalar use complex, not double, scalars; GSL functions are gsl_vector_complex_scale, gsl_vector_complex_add_constant, gsl_vector_complex_float_scale, gsl_vector_complex_float_add_constant.
-
Sumant Oemrawsingh authored
The backward FFT test for complex numbers doesn't work properly now because of some corruption that happens with elt/. But the forms should be correct.
-
- Nov 17, 2009
-
-
Liam Healy authored
The FFT noise test test-fft-noise calls a new split-out function forward-fft-rc as well as the DFT "answer". In preparation for the generation of the test files, several functions in generate-tests.lisp have been modified to accept a known answer, in lieu of evaluating the form and saving that as the reference.
-
Sumant Oemrawsingh authored
The reason why the complex case didn't work in example.lisp, is that the FFTs (not DFTs) work in-place. for reals it worked because the random (source) vector was copied (or at least, only the real part) before FFT-ing. Now, it's also copied in the complex case.
-
- Nov 16, 2009
-
-
Liam M. Healy authored
Eliminate defparameter *pass-complex-scalar-as-two-reals* and associate comments which predate FSBV.
-
Liam Healy authored
The noise generator make-urand-vector and the noise test, now called test-fft-noise, now work for all four element types, double-float, single-float, (complex double-float), (complex single-float). These are meant to emulate the noise parts of GSL's test_complex_radix2, test_real_radix2, test_complex_float_radix2, test_real_float_radix2, test_complex_float_func, test_complex_func, test_real_float_func, and test_real_func. The real tests pass and the complex tests currently fail.
-
- Nov 15, 2009
-
-
Sumant Oemrawsingh authored
Unpack functions were still not properly handling strides. This is now fixed. The test-real-radix2 function (from GSL's original test_real_radix2) now also works with strides. The function is not yet complete, though, as it does not yet test backward/inverse functions, and also doesn't actually perform comparisons yet.
-
- Nov 14, 2009
-
-
Liam Healy authored
A few overlooked function symbols have been exported and some unused function definitions eliminated. Added an example solve-tridiagonal-example thanks to Mirko Vukovic.
-