From fd79c2cfff5d6d0a462a8b09f25f8093ffd07de7 Mon Sep 17 00:00:00 2001 From: liam <liam@a3d8a0fb-c1db-0310-ace7-a616afeb9e30> Date: Thu, 20 Apr 2006 22:18:20 +0000 Subject: [PATCH] Until vectors with elements other than doubles are implemented, certain functions are not available. git-svn-id: svn+ssh://pop/opt/space/mathematics/gsl/trunk@3042 a3d8a0fb-c1db-0310-ace7-a616afeb9e30 --- README | 38 +++++++++++++++++++++++++++++++------- polynomial.lisp | 5 ++++- 2 files changed, 35 insertions(+), 8 deletions(-) diff --git a/README b/README index bb08f375..268f7ecc 100644 --- a/README +++ b/README @@ -13,11 +13,35 @@ Features: - This library is pure common lisp. No C is used; the only requirement is that the GSL library must be loaded and accessible. -- ??? No actual functions are defined, just the macros for gaining -access to the functions. Defining the functions you need yourself is -very easy. - ====================== -Needs work -bessel-sequence-Jnu in bessel.lisp. -Why do sequence-returns crash on declination? +Needs work: +See Thu Apr 20 2006 + +There are 14 kinds of element types for vectors and matrices: + +gsl_vector_complex_long_double.h +gsl_vector_complex_double.h +gsl_vector_complex_float.h + +gsl_vector_long_double.h +gsl_vector_double.h +gsl_vector_float.h + +gsl_vector_ulong.h +gsl_vector_long.h + +gsl_vector_uint.h +gsl_vector_int.h + +gsl_vector_ushort.h +gsl_vector_short.h + +gsl_vector_uchar.h +gsl_vector_char.h + +I have implemented only one, double. + +This means that the following functions have been put aside: +- #'polynomial-solve (requires vector_complex) +- The index functions in sorting that are not returning permutations, +because the need a vector_int. diff --git a/polynomial.lisp b/polynomial.lisp index 9a12ba68..1db16a15 100644 --- a/polynomial.lisp +++ b/polynomial.lisp @@ -3,7 +3,7 @@ ; description: Polynomials ; date: Tue Mar 21 2006 - 18:33 ; author: Liam M. Healy -; modified: Wed Apr 19 2006 - 09:41 +; modified: Thu Apr 20 2006 - 18:16 ;******************************************************** ;;; $Id: $ @@ -209,12 +209,15 @@ :c-return-value :void) ,workspace)))) +#+future (defun-gsl polynomial-solve-ws ((coefficients (:double n)) (workspace poly-complex-workspace)) "gsl_poly_complex_solve" :return ((gsl-complex (1- n)))) +#+future (export '(polynomial-solve)) +#+future (defmacro polynomial-solve (coefficients &optional workspace) "The roots of the general polynomial @c{$P(x) = a_0 + a_1 x + a_2 x^2 + ... + a_{n-1} x^{n-1}$} -- GitLab