Skip to content
Snippets Groups Projects
Commit 743b14ba authored by Alexander Pace's avatar Alexander Pace
Browse files

Suppressing XLAL Error - <GSL function> (fsolver.c:126): Generic failure

parent a87bcece
No related branches found
No related tags found
No related merge requests found
...@@ -83,6 +83,7 @@ ...@@ -83,6 +83,7 @@
#include <gst/controller/controller.h> #include <gst/controller/controller.h>
#include <math.h> #include <math.h>
#include <string.h> #include <string.h>
#include <gsl/gsl_errno.h>
/* /*
* our own stuff * our own stuff
...@@ -556,16 +557,22 @@ static GstFlowReturn push_nongap(GSTLALItac *element, guint copysamps, guint out ...@@ -556,16 +557,22 @@ static GstFlowReturn push_nongap(GSTLALItac *element, guint copysamps, guint out
{ {
GstBuffer *srcbuf = NULL; GstBuffer *srcbuf = NULL;
GstFlowReturn result = GST_FLOW_OK; GstFlowReturn result = GST_FLOW_OK;
gsl_error_handler_t *old_gsl_error_handler;
union { union {
float complex * as_complex; float complex * as_complex;
double complex * as_double_complex; double complex * as_double_complex;
void * as_void; void * as_void;
} dataptr; } dataptr;
/* make sure the snr threshold is up-to-date */ /* make sure the snr threshold is up-to-date */
element->maxdata->thresh = element->snr_thresh; element->maxdata->thresh = element->snr_thresh;
/* call the peak finding library on a buffer from the adapter if no events are found the result will be a GAP */ /* call the peak finding library on a buffer from the adapter if no events are found the result will be a GAP */
gst_audioadapter_copy_samples(element->adapter, element->data, copysamps, NULL, NULL); gst_audioadapter_copy_samples(element->adapter, element->data, copysamps, NULL, NULL);
/*AEP- 180417 Turning XLAL Errors off*/
old_gsl_error_handler=gsl_set_error_handler_off();
/* put the data pointer one pad length in */ /* put the data pointer one pad length in */
if (element->peak_type == GSTLAL_PEAK_COMPLEX) { if (element->peak_type == GSTLAL_PEAK_COMPLEX) {
...@@ -581,6 +588,10 @@ static GstFlowReturn push_nongap(GSTLALItac *element, guint copysamps, guint out ...@@ -581,6 +588,10 @@ static GstFlowReturn push_nongap(GSTLALItac *element, guint copysamps, guint out
else else
g_assert_not_reached(); g_assert_not_reached();
/*AEP- 180417 Turning XLAL Errors back on*/
gsl_set_error_handler(old_gsl_error_handler);
/* compute \chi^2 values if we can */ /* compute \chi^2 values if we can */
if (element->autocorrelation_matrix) { if (element->autocorrelation_matrix) {
/* compute the chisq norm if it doesn't exist */ /* compute the chisq norm if it doesn't exist */
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment