Opened 13 years ago
Last modified 10 years ago
#1398 reopened defect
Errores en gsl_cdf_beta_Pinv con diferentes versiones de GSL
| Reported by: | Víctor de Buen Remiro | Owned by: | Jorge |
|---|---|---|---|
| Priority: | highest | Milestone: | Mantainance |
| Component: | Math | Version: | head |
| Severity: | blocker | Keywords: | |
| Cc: |
Description
La siguiente línea TOL
Real gsl_cdf_beta_Pinv(0.9999999868417623,24.95643274099755,918.6566331776214);
en windows devuelve 0.06581172725382078 y en linux omitido.
He probado en R y me da lo mismo que en windows por lo que entiendo que está bien.
qbeta(0.9999999868417623,24.95643274099755,918.6566331776214); [1] 0.06581173
En windows estamos usando la versión 1.8 de 2006 que es al que encontré en su día ya precompilada como DLL para llegar y meter en Visual C++. En linux estamos en la 1.13 o la 1.12. Yo he probado a compilar en windows un main.c de un proyecto de ejemplo con la 1.13 y falla como en windows, dando algo más de información que no sé muy bien porqué se está perdiendo pues supuestamente tenemos un hook para recogerla:
gsl: ..\cdf\betainv.c:181: ERROR: inverse failed to converge Default GSL error handler invoked.
Definitivamente hay que retroceder en linux a una GSL que funcione, empezando por la 1.8
Change History (13)
comment:1 Changed 13 years ago by
| Owner: | changed from Víctor de Buen Remiro to Jorge |
|---|---|
| Status: | new → accepted |
comment:2 Changed 13 years ago by
| Resolution: | → fixed |
|---|---|
| Status: | accepted → closed |
Hemos comprobado que con la versión 1.8 esa función retorna los resultados esperados.
comment:4 Changed 13 years ago by
| Summary: | Error en gsl_cdf_beta_Pinv con GSL 1.13 → Errores en gsl_cdf_beta_Pinv con diferentes versiones de GSL |
|---|
Esta otra línea devuelve omitido con la GSL 1.8
Real gsl_cdf_beta_Pinv(0.9446900135371834,0.06873371518266619,6.27892521836533);
En cambio con las GSL 1.13 devuelve 0.058553049991647677 lo mismo que en R.
comment:5 Changed 13 years ago by
| Resolution: | fixed |
|---|---|
| Status: | closed → reopened |
comment:9 Changed 13 years ago by
comment:12 Changed 13 years ago by
| Resolution: | → remind |
|---|---|
| Status: | reopened → closed |
comment:13 Changed 10 years ago by
| Resolution: | remind |
|---|---|
| Status: | closed → reopened |
El test no pasa con TOL-GNU. Véase #1812.
Emite el siguiente mensaje:
Warning: [1] [GSL Error 5] inverse failed to converge At line 181 of file betainv.c
Empiezo por probar la version 1.8 en Linux. He descargado los fuentes de http://gnuwin32.sourceforge.net/downlinks/gsl-src-zip.php ya que el sitio oficial está de "black out"