[bug #65728] Add sparse functionalities

2024-05-12 Thread Patrick Alken
mments: --- Date: Sun 12 May 2024 01:09:04 PM UTC By: Patrick Alken from Brijesh Upadhaya I don't know if this is necessary, but the scaling of a complex (sparse) matrix with a real vector, complex sparse matrix-vector multiplication, and unpacking(in place) a complex vector

[bug #63679] C99 compatibility fix for the configure script

2024-05-10 Thread Patrick Alken
Update of bug #63679 (group gsl): Status:None => Fixed Open/Closed:Open => Closed ___ Follow-up Comment #2: Patch applied in commit

Re: Typo in Mersenne Twister PRNG's leading comment

2023-06-26 Thread Patrick Alken
Thank you David, I fixed it on the git On 6/26/23 10:47, David Grimm wrote: Hello, as of commit b6c4b6e077679a4ce9b4a6e0889d0e59672acc3a of the git repository git.savannah.gnu.org/gsl.git, in the file rng/mt.c, lines 43-44, the title of the paper describing the Mersenne Twister algorithm is

Re: -lcblas needed with -latlas

2022-04-25 Thread Patrick Alken
If you're using the ATLAS library, you need both -latlas and -lcblas together. If you are using the GSL cblas, you only need -lgslcblas On 4/25/22 15:15, Crawford, Christopher wrote: https://www.gnu.org/software/gsl/doc/html/usage.html#an-example-program I was looking through this

Re: Segmentation fault in gsl_integration_fixed (GSL Integration)

2022-04-19 Thread Patrick Alken
I don't think this is a bug. For Hermite integration, b must be > 0 (see documentation). For your input parameters, the weight function is 1, so you could just use Legendre integration On 4/19/22 18:57, Zhoulai Fu wrote: A segmentation fault is raised in gsl_integration_fixed function in GSL

Re: Potential security bug: Buffer overflow in gsl_stats_quantile_from_sorted_data (of library Statistics)

2022-04-16 Thread Patrick Alken
Hi Zhoulai,   I appreciate you following up on this, and for your patience. I have just committed a fix to the git repository which should address the issue. I added a check on the input parameter 'f' which returns an error if it is out of range. Best regards, Patrick On 4/16/22 08:46,

[bug #59624] Buffer overflow in gsl_stats_quantile_from_sorted_data

2022-04-16 Thread Patrick Alken
Update of bug #59624 (project gsl): Status:None => Fixed Open/Closed:Open => Closed ___ Follow-up Comment #2: fixed in commit

Re: Fwd: bug #59781

2022-03-27 Thread Patrick Alken
. The website says that anyone interested in contributing should start with working on the issues mentioned in the bug tracker. Being acquainted with the QR algorithm I decided to work on bug #59781. Thank you very much to Patrick Alken for including well written commentary in the code. Without

[bug #59781] Eigenvalue failure

2022-03-27 Thread Patrick Alken
Update of bug #59781 (project gsl): Status:None => Fixed Open/Closed:Open => Closed ___ Follow-up Comment #2: fixed in commit

Re: Fwd: bug #59781

2022-03-24 Thread Patrick Alken
to Patrick Alken for including well written commentary in the code. Without the guidance to Matrix Computations and DLAHQR, working with the code would have been substantially more difficult. Part of algorithm 7.5.2 from Matrix Computations involves looking for subdiagonal entries that are sufficiently

[bug #59781] Eigenvalue failure

2022-03-24 Thread Patrick Alken
Follow-up Comment #1, bug #59781 (project gsl): emails from M. Dunlap mdunlap1 =at= email =dot= arizona =dot= edu Part of algorithm 7.5.2 from Matrix Computations involves looking for subdiagonal entries that are sufficiently close to zero. If any are found we can break the problem up into two

[bug #39473] more efficient algorithm for 3j, 6j, 9j calculations (gsl_sf_coupling_{3j, 6j, 9j}_e

2021-11-07 Thread Patrick Alken
Additional Item Attachment, bug #39473 (project gsl): File name: report.pdf Size:33 KB File name: coupling3j.patch Size:15 KB

[bug #61094] gsl_linalg_LU_decomp fails to compute LU decomposition for some singular matrix

2021-08-31 Thread Patrick Alken
Update of bug #61094 (project gsl): Status:None => Fixed Open/Closed:Open => Closed ___ Follow-up Comment #1: Thank you for the

Re: missing declarations in Basis Splines

2021-07-03 Thread Patrick Alken
Hello,   I mistakenly updated the doc website with some experimental docs for some new B-spline features. These features are not in v2.7, or on the git yet. I have now updated the website with the correct 2.7 documentation. Sorry for the trouble. Thanks, Patrick On 7/3/21 3:01 PM, florian

Re: Undefined reference to `gsl_spmatrix_uchar_norm1'

2021-06-28 Thread Patrick Alken
gsl_spmatrix_uchar.h#n175> there should be a function gsl_spmatrix_uchar_norm1(), so swig gets confused and creates a wrapper for a non-existing function. Maybe that line 175 could be commented out? Håkon On Fri, Jun 25, 2021 at 6:46 PM Patrick Alken <mailto:al...@colorado.edu>> wrote:

Re: Undefined reference to `gsl_spmatrix_uchar_norm1'

2021-06-25 Thread Patrick Alken
Hello,   Unsigned versions of this function were excluded (uchar, uint, ushort), because there is the possibility of adding two large uchar values together during the 1-norm calculation and ending up with a smaller value due to it wrapping back to 0. This didn't make much sense to me when I

Re: Inaccurate Results for Lambert W function

2021-06-06 Thread Patrick Alken
Thank you, I have logged the issue into the bug tracker On 6/6/21 10:56 AM, Daming Zou wrote: Hi, The gsl_sf_lambert_W0(x) function gives inaccurate results for inputs near 0. I attached a code snippet to demonstrate it at the end of this mail. The details are as following: [Error

[bug #60741] Inaccurate Results for Lambert W function

2021-06-06 Thread Patrick Alken
URL: Summary: Inaccurate Results for Lambert W function Project: GNU Scientific Library Submitted by: psa Submitted on: Mon 07 Jun 2021 03:49:16 AM UTC Category: Runtime error

Re: Describe the error

2021-05-19 Thread Patrick Alken
Thanks, I fixed this bug on the git On 5/19/21 9:13 AM, Mark Galassi wrote: >> The descriptions of function gsl_histogram_max_bin and >> function gsl_histogram_min_bin have something in common. >> Is that a problem? > I think you mean to say "the description of gsl_histogram_min_bin has an >

Re: Problem with Schur decomposition

2021-05-18 Thread Patrick Alken
Note that the documentation says: "If T is desired, it is stored in the upper portion of A on output". This means that the gsl_eigen_nonsymm / gsl_eigen_nonsymm_Z functions do not zero out the lower portion of the output matrix. The lower portion of A is used as temporary workspace by the

[bug #60635] physical constants may need updating

2021-05-18 Thread Patrick Alken
URL: Summary: physical constants may need updating Project: GNU Scientific Library Submitted by: psa Submitted on: Wed 19 May 2021 02:23:28 AM UTC Category: Accuracy problem

Re: Bug in the complex no header file of the gsl_complex.h .

2021-05-13 Thread Patrick Alken
Try this instead: gsl_complex z; GSL_SET_COMPLEX(,5.7,8.0); printf("output %f",gsl_complex_arg(z)); On 5/13/21 12:50 AM, Chetan wrote: >    error: expected identifier or ‘(’ before ‘do’ >        95 | #define GSL_SET_COMPLEX(zp,x,y) do {(zp)->dat[0]=(x); >    (zp)->dat[1]=(y);} while(0) > >   

Re: [bug #60457] Feature request: complex tridiagonal solvers

2021-04-26 Thread Patrick Alken
Hello, certainly I can add that function. It would be best if you send me a patch with 'git diff' Patrick On 4/26/21 12:23 PM, Ed wrote: > URL: > > > Summary: Feature request: complex tridiagonal solvers > Project: GNU

[bug #55412] Bessel function returns incorrect value

2021-04-22 Thread Patrick Alken
Update of bug #55412 (project gsl): Status:None => Fixed Open/Closed:Open => Closed ___ Follow-up Comment #1: Note that the

[bug #60370] Interpolation inconsistencies

2021-04-11 Thread Patrick Alken
Update of bug #60370 (project gsl): Status:None => Duplicate Open/Closed:Open => Closed ___ Follow-up Comment #1: This is a duplicate,

[bug #60371] Interpolation domain error handling

2021-04-10 Thread Patrick Alken
URL: Summary: Interpolation domain error handling Project: GNU Scientific Library Submitted by: psa Submitted on: Sun 11 Apr 2021 03:59:07 AM UTC Category: Runtime error

[bug #60370] Interpolation inconsistencies

2021-04-10 Thread Patrick Alken
URL: Summary: Interpolation inconsistencies Project: GNU Scientific Library Submitted by: psa Submitted on: Sun 11 Apr 2021 03:54:57 AM UTC Category: Documentation

[bug #60335] Test failure in spmatrix

2021-04-04 Thread Patrick Alken
Update of bug #60335 (project gsl): Status:None => Fixed Open/Closed:Open => Closed ___ Follow-up Comment #1: spmatrix test failure

[bug #60335] Test failure in spmatrix

2021-04-03 Thread Patrick Alken
URL: Summary: Test failure in spmatrix Project: GNU Scientific Library Submitted by: psa Submitted on: Sun 04 Apr 2021 05:13:07 AM UTC Category: Runtime error

[bug #36577] vector/matrix _add, _scale routines take wrong scalar type

2021-02-18 Thread Patrick Alken
Update of bug #36577 (project gsl): Status: Confirmed => Fixed Open/Closed:Open => Closed ___ Follow-up Comment #1: fixed in commit

[bug #60026] Incorporate MIXMAX random number extension into GSL

2021-02-08 Thread Patrick Alken
URL: Summary: Incorporate MIXMAX random number extension into GSL Project: GNU Scientific Library Submitted by: psa Submitted on: Tue 09 Feb 2021 04:09:05 AM UTC Category: None

Re: gsl_ran_gaussian_ziggurat doesn't work with gsl_rng_slatec and gsl_rng_uni

2021-02-05 Thread Patrick Alken
Hello,   I would welcome your help with fixing rng issues. We don't currently have a GSL developer with expertise in rngs (I myself have never worked on them). The best approach would be to clone the git repository (instructions on the GSL webpage), make your changes, and then send me a git diff

[bug #45521] Erroneous use of GSL_ERROR_NULL instead of GSL_ERROR macro

2021-01-21 Thread Patrick Alken
Update of bug #45521 (project gsl): Status:None => Fixed Open/Closed:Open => Closed ___ Follow-up Comment #1: fixed in commit

[bug #59910] Update information for Python wrapper

2021-01-21 Thread Patrick Alken
Update of bug #59910 (project gsl): Status:None => Fixed Open/Closed:Open => Closed ___ Follow-up Comment #1: added link for

[bug #59494] Bug in gsl_rstat

2021-01-20 Thread Patrick Alken
Update of bug #59494 (project gsl): Status:None => Fixed Open/Closed:Open => Closed ___ Follow-up Comment #1: The issue with n=5 has

[bug #59758] Doc bug

2021-01-19 Thread Patrick Alken
Update of bug #59758 (project gsl): Status:None => Fixed Open/Closed:Open => Closed ___ Follow-up Comment #1: fixed in commit

[bug #59710] Undefined behavior error in Nonlinear Least-Squares Fitting

2021-01-19 Thread Patrick Alken
Update of bug #59710 (project gsl): Status:None => Fixed ___ Reply to this item at: ___

[bug #58202] Wrong Result Median for gsl_rstat_quantile_get with fix

2021-01-19 Thread Patrick Alken
Update of bug #58202 (project gsl): Status:None => Fixed Open/Closed:Open => Closed ___ Follow-up Comment #1: fixed in commit

[bug #59914] Native build of GSL-2.5 on windows 10

2021-01-19 Thread Patrick Alken
URL: Summary: Native build of GSL-2.5 on windows 10 Project: GNU Scientific Library Submitted by: psa Submitted on: Wed 20 Jan 2021 02:43:53 AM UTC Category: None

[bug #59913] gsl 2.3.0 problem in gsl_integration_cquad

2021-01-19 Thread Patrick Alken
URL: Summary: gsl 2.3.0 problem in gsl_integration_cquad Project: GNU Scientific Library Submitted by: psa Submitted on: Wed 20 Jan 2021 02:41:50 AM UTC Category: None

[bug #59912] gsl_permutation header files

2021-01-19 Thread Patrick Alken
URL: Summary: gsl_permutation header files Project: GNU Scientific Library Submitted by: psa Submitted on: Wed 20 Jan 2021 02:38:43 AM UTC Category: Documentation

[bug #59911] Problem with qagui 1D integrator

2021-01-19 Thread Patrick Alken
URL: Summary: Problem with qagui 1D integrator Project: GNU Scientific Library Submitted by: psa Submitted on: Wed 20 Jan 2021 02:36:50 AM UTC Category: None

[bug #59910] Update information for Python wrapper

2021-01-19 Thread Patrick Alken
URL: Summary: Update information for Python wrapper Project: GNU Scientific Library Submitted by: psa Submitted on: Wed 20 Jan 2021 02:34:41 AM UTC Category: None

[bug #59900] Add truncated normal distribution

2021-01-17 Thread Patrick Alken
URL: Summary: Add truncated normal distribution Project: GNU Scientific Library Submitted by: psa Submitted on: Sun 17 Jan 2021 10:26:25 PM UTC Category: None

[bug #59624] Buffer overflow in gsl_stats_quantile_from_sorted_data

2021-01-17 Thread Patrick Alken
Follow-up Comment #1, bug #59624 (project gsl): from msunet =at= shellblade =dot= net This one is a "fix" for bug 59624, simply checking that |f| is in the expected range. Though I don't think there is ultimately a way to check in C whether the given array has the right size, so a buggy

[bug #59834] Misaligned memory access error in deque.c

2021-01-06 Thread Patrick Alken
URL: Summary: Misaligned memory access error in deque.c Project: GNU Scientific Library Submitted by: psa Submitted on: Wed 06 Jan 2021 03:50:58 PM UTC Category: Runtime error

[bug #59781] Eigenvalue failure

2020-12-29 Thread Patrick Alken
URL: Summary: Eigenvalue failure Project: GNU Scientific Library Submitted by: psa Submitted on: Tue 29 Dec 2020 11:38:25 PM UTC Category: Runtime error Severity: 3 -

Re: gsl: francis.c:209: ERROR: maximum iterations reached without finding all eigenvalues

2020-12-29 Thread Patrick Alken
0, 0, 0, 0, 0, 0, 0, 0, 0, 0 }; > > Wolfram Mathematica: > In[4]:= Eigenvalues[data] > Out[4]= {-9, -9, -9, -9, 9, 9, 9, 9, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0} > > > With best regards, > Dmitry Cheshkov > > 29.12.2020 20:11, Patrick Alken пишет: >> Dear Dmitry, >> >>

Re: gsl: francis.c:209: ERROR: maximum iterations reached without finding all eigenvalues

2020-12-29 Thread Patrick Alken
Dear Dmitry,   I successfully found the eigenvalues of your matrix using GSL 2.6 (see attached code). The output of the program is: --- $ ./eig eval[1] = 0.e+00 + i 0.e+00 eval[2] = 0.e+00 + i 0.e+00 eval[3] = 0.e+00 + i

[bug #59759] spmatrix test fails on x86_64

2020-12-23 Thread Patrick Alken
Follow-up Comment #1, bug #59759 (project gsl): from noloader =at= gmail =dot= com ___ Reply to this item at: ___ Message sent via Savannah

[bug #59759] spmatrix test fails on x86_64

2020-12-23 Thread Patrick Alken
:47 PM Patrick Alken wrote: > > Hello, I just committed a fix to the git repository. Would you be able > to test the git version to see if its resolved on your platform? Hi Patrick. Here's the result for the tip of master: make check-TESTS make[2]: Entering directory '/home/jwalton/Bu

[bug #59758] Doc bug

2020-12-23 Thread Patrick Alken
URL: Summary: Doc bug Project: GNU Scientific Library Submitted by: psa Submitted on: Wed 23 Dec 2020 10:01:14 PM UTC Category: Documentation Severity: 3 - Normal

[bug #59710] Undefined behavior error in Nonlinear Least-Squares Fitting

2020-12-17 Thread Patrick Alken
Update of bug #59710 (project gsl): Open/Closed:Open => Closed ___ Follow-up Comment #1: Thanks for your report. The problem is that on the first iteration, the lambda value jumps to

[bug #59710] Undefined behavior error in Nonlinear Least-Squares Fitting

2020-12-17 Thread Patrick Alken
URL: Summary: Undefined behavior error in Nonlinear Least-Squares Fitting Project: GNU Scientific Library Submitted by: psa Submitted on: Thu 17 Dec 2020 04:17:55 PM UTC Category:

Re: Undefined behavior error in Nonlinear Least-Squares Fitting

2020-12-13 Thread Patrick Alken
Hello,   Thanks for your report. The problem is that on the first iteration, the lambda value jumps to -12.8, and the value of exp(12.8 * 60) is too large to represent in double precision, so infs and nans appear in the residual vector and Jacobian matrix. For this particular problem, the value

[bug #59624] Buffer overflow in gsl_stats_quantile_from_sorted_data

2020-12-04 Thread Patrick Alken
URL: Summary: Buffer overflow in gsl_stats_quantile_from_sorted_data Project: GNU Scientific Library Submitted by: psa Submitted on: Fri 04 Dec 2020 10:06:37 PM UTC Category: Runtime

[bug #59494] Bug in gsl_rstat

2020-11-20 Thread Patrick Alken
URL: Summary: Bug in gsl_rstat Project: GNU Scientific Library Submitted by: psa Submitted on: Fri 20 Nov 2020 04:21:29 PM UTC Category: Runtime error Severity: 3 -

Re: bug in gsl_rstat_workspace

2020-11-20 Thread Patrick Alken
Thanks, I've logged the issue into the bug tracker, and will investigate when I find time On 11/20/20 3:40 AM, Wilhelm Braunschober wrote: > Yours Wilhelm Braunschober

Re: spmatrix tests fail on x86_64

2020-08-28 Thread Patrick Alken
Hello, I just committed a fix to the git repository. Would you be able to test the git version to see if its resolved on your platform? Patrick On 8/28/20 6:26 PM, Jeffrey Walton wrote: > Hi Everyone, > > I'm building the GSL 2.6 release tarball on Ubuntu 18.04 LTS x86_64 > fully patched. 'make

[bug #58977] Any plans to make new release?

2020-08-18 Thread Patrick Alken
Update of bug #58977 (project gsl): Open/Closed:Open => Closed ___ Reply to this item at: ___

Re: Observed large error from spherical Bessel functions at high orders

2020-07-14 Thread Patrick Alken
Hello,   Unfortunately some of the special functions do not produce the expected accuracy for all inputs, and it seems you've found another example. It would be great if you could look into algorithms for computing these functions and possibly find a way to improve the accuracy for these inputs.

Re: Extrapolation returns NaN

2020-06-19 Thread Patrick Alken
On 6/19/20 6:13 AM, Yury Pakhomov wrote: > Hello All! > > gsl_spline_eval returns NaN values when x is out of range. > > Since GSL-1.15 this situation remains unresolved > http://git.savannah.gnu.org/cgit/gsl.git/tree/NEWS#n470 > > BW, Yuri > This is the documented behavior ((see gsl_interp_eval

[bug #52540] Complex QR decomposition

2020-06-07 Thread Patrick Alken
Update of bug #52540 (project gsl): Status:None => Fixed ___ Reply to this item at: ___

[bug #52540] Complex QR decomposition

2020-06-07 Thread Patrick Alken
Update of bug #52540 (project gsl): Open/Closed:Open => Closed ___ Follow-up Comment #2: Complex QR has been added to GSL ___

Re: 'syrk' very slow compared to 'gemm'

2020-02-19 Thread Patrick Alken
Hello,   The GSL CBLAS library is solely intended so that users can compile applications "out of the box" without needing to install an optimized BLAS library. The GSL CBLAS is not intended for serious work, and should never be used for performance critical applications. We always recommend

Re: Inverse Gamma Issue

2020-01-28 Thread Patrick Alken
Thanks Robert,   There have been many bug reports sent in recently - I need to add them all the tracker. Hopefully I can find time to look at your work, though it may take a while. Thanks, Patrick On 1/28/20 4:09 AM, Janes, Robert (Columbus) via Bug reports for the GNU Scientific Library

Re: [bug #57173] Feature request: zeta function for complex arguments

2019-11-05 Thread Patrick Alken
Jungman wrote: > On 11/5/19 11:26 AM, Mark Galassi wrote: >>>>>>> "Patrick" == Patrick Alken writes: >> >> >> Patrick> I don't think any progress was made on this >> >> One update to the discussion from back then is that C++1

Re: [bug #57173] Feature request: zeta function for complex arguments

2019-11-05 Thread Patrick Alken
I don't think any progress was made on this On 11/5/19 9:49 AM, Mark Shoulson wrote: > URL: > > > Summary: Feature request: zeta function for complex > arguments > Project: GNU Scientific Library >

Re: [Bug-gsl] Possible error handling/return bug in vector/copy_source.c

2019-09-22 Thread Patrick Alken
Yes you're right, thanks for catching this. I put the vector length checks first now on the git. This code was just recently changed for v2.6 to add the BLAS calls, so I overlooked that. Patrick On 9/22/19 12:55 PM, Diego via Bug-gsl wrote: > I've been playing with a wrapper for GSL for Chicken

Re: [Bug-gsl] About the web page of GSL - GNU Scientific Library

2019-09-13 Thread Patrick Alken
Thanks, I have updated the link On 8/20/19 9:20 PM, mrkmh...@tmu.ac.jp wrote: > In the web page of >GSL - GNU Scientific Library > URL > > You write > >> A Japanese translation is also available online (may not be the most > recent version). >>GSL

Re: [Bug-gsl] the multifit_nlinear subpackage is mucked up

2019-08-08 Thread Patrick Alken
On 8/8/19 3:30 AM, JohnT wrote: > It fails to test ok. > > > > Can you give details, or post the output from your test*.log files?

Re: [Bug-gsl] Hang caused by Inf inputs in hermite special function

2019-07-16 Thread Patrick Alken
Hello, some work was done on the hermite functions in the current branch of GSL on the git repository. Could you test this issue by checking out the latest git version? Thanks, Patrick On 7/15/19 7:20 PM, Jackson Vanover wrote: > Using version 2.5 of GSL from

Re: [Bug-gsl] Building DLLs on Windows

2019-06-19 Thread Patrick Alken
Philip, I forwarded your email to Brian (cc'd) who has done a lot of work on windows builds of GSL. His response is pasted below. To build GSL as it is now requires that WIN32 is defined somewhere in the GSL code. The _WIN32 define is always defined by the MSVC compiler so it this were

Re: [Bug-gsl] Minor inaccuracy in comment

2019-02-20 Thread Patrick Alken
Thank you, I have corrected this on the git On 2/19/19 1:46 AM, Chris Hiszpanski wrote: > Hello, > > Noticed a comment in rng/taus113.c is inaccurate -- z_1 is computed from z_0, > not the other way around. Patch is provided below. > > Chris > > diff --git a/rng/taus113.c b/rng/taus113.c > index

Re: [Bug-gsl] gsl 2.3.0 bug?

2018-11-20 Thread Patrick Alken
Hello,   I converted your main.cpp to C (I couldn't make it compile with my g++). Also I linked it against GSL v2.5, but the CQUAD code hasn't changed since 2.3 as far as I know. The program works as expected (C code attached). $ ./main result = 4.390157623791e-164 abserr =

[Bug-gsl] [bug #54925] Reduce cache misses for source_gemm_r

2018-10-31 Thread Patrick Alken
URL: Summary: Reduce cache misses for source_gemm_r Project: GNU Scientific Library Submitted by: psa Submitted on: Wed 31 Oct 2018 01:57:19 PM UTC Category: Performance

[Bug-gsl] [bug #54140] Updating email

2018-10-31 Thread Patrick Alken
Update of bug #54140 (project gsl): Status:None => Fixed Open/Closed:Open => Closed ___ Follow-up Comment #1: fixed in commit

[Bug-gsl] [bug #54921] delete keyword clash in gsl_movstat.h

2018-10-31 Thread Patrick Alken
Update of bug #54921 (project gsl): Status:None => Fixed Open/Closed:Open => Closed ___ Follow-up Comment #1: Fixed in commit

Re: [Bug-gsl] bug with gsl_integration

2018-10-19 Thread Patrick Alken
Thank you for finding a solution Johannes. I've never worked with mex/GSL myself. If the default handler is used, then the program will abort/segfault anyway, so I don't think it should matter if you are calling from mex or from a command line. On 10/18/2018 11:26 PM, Johannes Hendriks wrote: >

Re: [Bug-gsl] [bug #54140] Updating email

2018-06-20 Thread Patrick Alken
Hello, it is done On 06/18/2018 03:22 AM, Nicolas Darnis wrote: > URL: > > > Summary: Updating email > Project: GNU Scientific Library > Submitted by: ndarnis > Submitted on: Mon 18 Jun 2018 11:22:21

Re: [Bug-gsl] Small documentation bug

2018-06-14 Thread Patrick Alken
Hello, thanks for spotting this. I've fixed it on the git On 06/14/2018 08:40 AM, Damien Desfontaines wrote: The diagram for gsl_ran_binomial_pdf looks wrong: https://www.gnu.org/software/gsl/doc/html/randist.html#c.gsl_ran_binomial_pdf It looks like this is a binomial distribution with

[Bug-gsl] [bug #53919] handle large values correctly in (log_)erf(c) functions

2018-05-15 Thread Patrick Alken
Follow-up Comment #1, bug #53919 (project gsl): The patch uses HUGE_VAL (C99) and isnan() checks - there is probably a cleaner way to handle the problem, but will require further investigation ___ Reply to this item at:

[Bug-gsl] [bug #53919] handle large values correctly in (log_)erf(c) functions

2018-05-15 Thread Patrick Alken
URL: Summary: handle large values correctly in (log_)erf(c) functions Project: GNU Scientific Library Submitted by: psa Submitted on: Wed 16 May 2018 02:55:06 AM UTC Category: Runtime

Re: [Bug-gsl] Bug in Hypergeometric function

2018-05-14 Thread Patrick Alken
Hello,   I have logged your bug report into our bug tracker. It may take some time to develop a fix for this. Thanks, Patrick On 05/01/2018 08:51 AM, Salehi Vaziri Kamran wrote: Dear Sir or Madam, Let me first thank you for your amazing package that is very well-written with a nice

[Bug-gsl] [bug #53904] Bug gsl_matrix_complex_set

2018-05-14 Thread Patrick Alken
URL: Summary: Bug gsl_matrix_complex_set Project: GNU Scientific Library Submitted by: psa Submitted on: Mon 14 May 2018 09:39:47 PM UTC Category: None Severity: 3 -

[Bug-gsl] [bug #53903] Test failure with gsl_sf_synchrotron_1_e on x86

2018-05-14 Thread Patrick Alken
URL: Summary: Test failure with gsl_sf_synchrotron_1_e on x86 Project: GNU Scientific Library Submitted by: psa Submitted on: Mon 14 May 2018 09:38:11 PM UTC Category: Runtime error

Re: [Bug-gsl] git source: doc/examples wants gaussfilt.c, but it does not exist

2018-05-10 Thread Patrick Alken
Hello,   It should work now Patrick On 05/10/2018 10:46 AM, Hans-Bernhard Bröker wrote: > Hello all, > > as of yesterday's commit 4cbb37d11189bc5bab3c8243bcfa45ff2ef75d22 > the git version of doc/examples/Makefile.am requests a non-existing > source-file which was apparently supposed to be

Re: [Bug-gsl] typo error

2018-04-21 Thread Patrick Alken
Thanks, I have fixed it on the git On 04/21/2018 05:54 PM, Alex Henrie wrote: > On Sat, Apr 21, 2018 at 4:42 PM Patrick Alken <al...@colorado.edu> wrote: >> Err, sorry I was looking at the wrong PDF. On page 534, I see a >> discussion of Autoconf macros, still no

Re: [Bug-gsl] typo error

2018-04-21 Thread Patrick Alken
Err, sorry I was looking at the wrong PDF. On page 534, I see a discussion of Autoconf macros, still no equation for 'e' On 04/21/2018 04:39 PM, Patrick Alken wrote: > Hello, can you give the chapter/section number? On page 534, I see an > example program for sparse linear algebra (sectio

Re: [Bug-gsl] typo error

2018-04-21 Thread Patrick Alken
Hello, can you give the chapter/section number? On page 534, I see an example program for sparse linear algebra (section 43.3) and I don't see any equation for 'e' on this page. On 04/21/2018 02:05 PM, Србислав Нешић wrote: > In PDF document > >

Re: [Bug-gsl] incorrect convergence check result

2018-03-22 Thread Patrick Alken
positive (converged) result. Regards, Edwin On 3/21/18, 12:26 PM, "Patrick Alken" <al...@colorado.edu> wrote: Hello, Can you provide a minimal working example program which illustrates the issue? Thanks, Patrick On 03/

Re: [Bug-gsl] incorrect convergence check result

2018-03-21 Thread Patrick Alken
Hello,   Can you provide a minimal working example program which illustrates the issue? Thanks, Patrick On 03/21/2018 08:02 AM, Sarkissian, Edwin (398E) wrote: Dear Madam/Sir, Sometimes, the function gsl_multifit_nlinear_test sets the argument info to 2 (converged based on gradient check)

Re: [Bug-gsl] please check lib gsl_multifit_nlinear.h for bugs

2018-03-20 Thread Patrick Alken
Hello,   The website documentation and examples are for GSL v2.4, so they will likely not work with v2.1 since the nonlinear least squares code has changed significantly in recent versions. Can you please install v2.4 and try the example program again? Also please send along any example code

Re: [Bug-gsl] error in multifit/multiwlinear.c in gsl-2.3 and gsl-2.4

2018-02-11 Thread Patrick Alken
Hello,   I responded to this yesterday. The issue was the gsl_multifit_wlinear_svd function is deprecated (replaced by gsl_multifit_wlinear_tsvd), and I hadn't tested the old function. The issue should be fixed now in the git repository. Could you checkout the latest git code and test your

Re: [Bug-gsl] multiwlinear.c:73: ERROR: size of workspace does not match size of observation matrix

2018-02-10 Thread Patrick Alken
Vlad,   I have uploaded a fix for this issue to the git repository. Now when I run your original program I get this output: $ ./test solution: x=0.5166  y=0.5166   (rank=1, chi2=0.004) Thanks, Patrick On 02/10/2018 10:52 AM, Patrick Alken wrote: > So as it turns

Re: [Bug-gsl] multiwlinear.c:73: ERROR: size of workspace does not match size of observation matrix

2018-02-10 Thread Patrick Alken
: gsl_multifit_wlinear_tsvd(A, w, y, tol, x, cov, , , wsp); This should make your program run correctly. In the meantime I will figure out what to do with gsl_multifit_wlinear_svd - maybe just make it a wrapper function for the new tsvd call. Sorry for the trouble, Patrick On 02/10/2018 10:28 AM, Patrick Alken

[Bug-gsl] [bug #52927] make check fails on Bessel j2 test

2018-01-18 Thread Patrick Alken
s [407] On Thu, Jan 18, 2018 at 3:01 PM, Patrick Alken <al...@colorado.edu> wrote: > Hello, > > It looks like its calculating the value correctly, but the test > tolerance needs to be relaxed slightly. Can you locate this line in > specfun/test_bessel.c: > > 186

Re: [Bug-gsl] FFT not generated

2018-01-07 Thread Patrick Alken
Hi,   Did the library pass 'make check'. What happens if you do: cd fft ; make test ; ./test On 01/05/2018 05:25 PM, Peter Schacht wrote: > Hello, > > > > I installed the GSL in win10 64bit -system with minGw and msys using the > cofigure script. The installaion seemed to work fine and many

Re: [Bug-gsl] Minor improvement to gsl_spmatrix.h suggested

2018-01-04 Thread Patrick Alken
Manuel,   Sounds like a good suggestion to me - I've made the change on the git and hopefully didn't break anything. If you're currently using gsl_spmatrix in your work and can test the update I'd appreciate it. Thanks, Patrick On 12/22/2017 09:44 AM, Schmitz Manuel (LBC) wrote: Dear Sirs,

[Bug-gsl] [bug #52570] Inaccuracy of the Airy function due to invocation of GSL's cosine function with large input parameters

2017-12-01 Thread Patrick Alken
URL: Summary: Inaccuracy of the Airy function due to invocation of GSL's cosine function with large input parameters Project: GNU Scientific Library Submitted by: psa Submitted on: Fri 01 Dec 2017

  1   2   3   4   >