[Bug-gsl] [bug #55687] Bad error handling in gsl_sf_hyperg_1F1_e with NaN arguments

[Simon Byrne]

URL:
  <https://savannah.gnu.org/bugs/?55687>

                 Summary: Bad error handling in gsl_sf_hyperg_1F1_e with NaN
arguments
                 Project: GNU Scientific Library
            Submitted by: simonbyrne
            Submitted on: Sun 10 Feb 2019 04:20:09 AM UTC
                Category: None
                Severity: 3 - Normal
        Operating System: 
                  Status: None
             Assigned to: None
             Open/Closed: Open
                 Release: 2.5
         Discussion Lock: Any

    _______________________________________________________

Details:

The following gives a segfault:

#include <stdio.h>
#include <math.h>
#include <gsl/gsl_errno.h>
#include <gsl/gsl_sf_hyperg.h>

int main (void)
{
  gsl_sf_result result;
  gsl_set_error_handler_off();
  int status = gsl_sf_hyperg_1F1_e(1.0,NAN,-1.0, &result);
  if (status == GSL_SUCCESS) {
    printf ("success: %.18e\n", result.val);
  } else {
    printf ("failure");
  }  
  return 0;
}

>From my understanding of the error handling, the _e functions should return an
invalid error code instead of crashing in such cases.

Downstream issue: https://github.com/JuliaMath/GSL.jl/issues/96




    _______________________________________________________

Reply to this item at:

  <https://savannah.gnu.org/bugs/?55687>

_______________________________________________
  Message sent via Savannah
  https://savannah.gnu.org/