![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to wmain.c CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM / src / hgm / mh / src|
Return to wmain.c CVS log | Up to [local] / OpenXM / src / hgm / mh / src |
| version 1.34, 2016/03/02 00:27:02 | version 1.36, 2016/06/06 04:49:11 | ||
|---|---|---|---|
|
|
||
| /* | /* | ||
| $OpenXM: OpenXM/src/hgm/mh/src/wmain.c,v 1.33 2016/03/01 07:13:30 takayama Exp $ | $OpenXM: OpenXM/src/hgm/mh/src/wmain.c,v 1.35 2016/06/04 07:52:14 takayama Exp $ | ||
| License: LGPL | License: LGPL | ||
| */ | */ | ||
| #include <stdio.h> | #include <stdio.h> | ||
|
|
||
| if (MH_strategy) { | if (MH_strategy) { | ||
| if (MH_abserr > SIGDIGIT_DEFAULT*myabs(y0[0])) { | if (MH_abserr > SIGDIGIT_DEFAULT*myabs(y0[0])) { | ||
| MH_success = 0; | MH_success = 0; | ||
| oxprintfe("%%%%Warning, abserr seems not to be small enough, abserr=%lg, y[0]=%lg. Increasing the starting point (q0 or X0g(standalone case)) may or making abserr (err[1] or abserror(standalone case)) smaller will help.\n",MH_abserr,y0[0]); | oxprintfe("%%%%Warning, abserr seems not to be small enough, abserr=%lg, y[0]=%lg. Increasing the starting point (q0 or X0g(standalone case)) may or making abserr (err[1] or abserror(standalone case)) smaller will help, e.g., err=c(%lg,1e-10)\n",MH_abserr,y0[0],y0[0]*(1e-6)); | ||
| }else{ | }else{ | ||
| MH_success = 1; | MH_success = 1; | ||
| } | } |