问题描述
|
我想使用GSL进行集成
http://www.gnu.org/software/gsl/manual/html_node/Numerical-Integration.html
但是,我发现没有便捷的方法如何集成功能
(示例http://www.gnu.org/software/gsl/manual/html_node/Numerical-integration-examples.html中的函数f)
可以向集成商报告错误。我想集成一个功能,该功能本身是由可能失败的集成产生的。这是我的示例程序
#include <stdio.h>
#include <math.h>
#include <gsl/gsl_integration.h>
#include <gsl/gsl_errno.h>
double f (double x,void * params) {
GSL_ERROR(\"test error\",GSL_FAILURE);
return 0.0;
}
int main (void)
{
gsl_integration_workspace * w = gsl_integration_workspace_alloc (1000);
double result,error;
gsl_function F;
F.function = &f;
gsl_set_error_handler_off();
int status = gsl_integration_qags (&F,1,1e-7,1000,w,&result,&error);
printf (\"status = %d\\n\",status);
status = GSL_FAILURE;
printf (\"status = %d\\n\",status);
gsl_integration_workspace_free (w);
return 0;
}
产生输出
状态= 0
状态= -1
我认为集成商应该宁愿停止并返回我的错误代码。我该如何实现?
非常感谢您的帮助!!!
2011-04-27:在Brian Gough告诉我之后,我也尝试了这种变体,
#include <stdio.h>
#include <math.h>
#include <gsl/gsl_integration.h>
#include <gsl/gsl_errno.h>
double f (double x,GSL_FAILURE);
return GSL_NAN;
}
int main (void)
{
gsl_integration_workspace * w = gsl_integration_workspace_alloc (1000);
double result,status);
gsl_integration_workspace_free (w);
return 0;
}
它也没有帮助。我现在将填写一个错误报告。解决方法
感谢GSL邮件列表中的Wu Xuebin Wu解决了这个问题:
你好
GSL_ERROR本身是一个宏,看起来像
gsl_error (reason,__FILE__,__LINE__,gsl_errno);
return gsl_errno;
在返回NAN之前,该函数已经返回,因为GSL_ERROR
已被调用。关闭处理程序,仅让第一行执行
没有。默认错误处理程序在打印后中止程序
错误信息。
我不认为这是一个错误。也许您可以编写自己的错误处理程序
解决您的问题。例如,您可以使用\“ goto \”跳出
gsl_integration_qags,或设置一些全局变量以指示
积分结果不正确。
PS:我相信您需要这个宏,
巨集:GSL_ERROR_VAL(原因,gsl_errno,值)
该宏与GSL_ERROR相同,但返回用户定义的值
值而不是错误代码。可用于数学
返回浮点值的函数。
以下示例显示了如何以数学方式返回NaN
使用GSL_ERROR_VAL宏的奇异性
if (x == 0)
{
GSL_ERROR_VAL(\"argument lies on singularity\",GSL_ERANGE,GSL_NAN);
}
所以我根据
#include <stdio.h>
#include <math.h>
#include <gsl/gsl_integration.h>
#include <gsl/gsl_errno.h>
double f (double x,void * params) {
// return GSL_NAN;
GSL_ERROR_VAL (\"argument lies on singularity\",GSL_NAN);
}
int main (void)
{
gsl_integration_workspace * w = gsl_integration_workspace_alloc (1000);
double result,error;
gsl_function F;
F.function = &f;
gsl_set_error_handler_off();
int status = gsl_integration_qags (&F,1,1e-7,1000,w,&result,&error);
printf (\"status = %d\\n\",status);
status = GSL_FAILURE;
printf (\"status = %d\\n\",status);
gsl_integration_workspace_free (w);
return 0;
}
一切都按预期进行...,有点黑,但是我可能会让您的函数存储一些标志。遇到错误时,它将设置标志并为所有后续评估返回零。然后,将其集成后,可以检查该标志以查看结果是否有效。,如何为该函数编写一个包装器,该包装器返回指向结构的指针,其中包含函数结果和错误状态?或者,如果您使用c ++,则可以使用对象...进行这种封装。