version 1.17, 2014/03/12 07:50:37 |
version 1.19, 2014/03/14 05:58:16 |
|
|
/* |
/* |
$OpenXM: OpenXM/src/hgm/mh/src/wmain.c,v 1.16 2014/03/11 05:20:45 takayama Exp $ |
$OpenXM: OpenXM/src/hgm/mh/src/wmain.c,v 1.18 2014/03/14 02:21:40 takayama Exp $ |
License: LGPL |
License: LGPL |
*/ |
*/ |
#include <stdio.h> |
#include <stdio.h> |
Line 84 main(int argc,char *argv[]) { |
|
Line 84 main(int argc,char *argv[]) { |
|
mh_freeWorkArea(); */ |
mh_freeWorkArea(); */ |
mh_main(argc,argv); |
mh_main(argc,argv); |
/* showParam(); */ |
/* showParam(); */ |
|
return(0); |
} |
} |
#endif |
#endif |
struct MH_RESULT *mh_main(int argc,char *argv[]) { |
struct MH_RESULT *mh_main(int argc,char *argv[]) { |
Line 152 struct MH_RESULT *mh_main(int argc,char *argv[]) { |
|
Line 153 struct MH_RESULT *mh_main(int argc,char *argv[]) { |
|
if (MH_Verbose) {for (i=0; i<rank; i++) printf("%lf\n",y0[i]); } |
if (MH_Verbose) {for (i=0; i<rank; i++) printf("%lf\n",y0[i]); } |
rp = (struct MH_RESULT*) mh_malloc(sizeof(struct MH_RESULT)); |
rp = (struct MH_RESULT*) mh_malloc(sizeof(struct MH_RESULT)); |
rp->message = NULL; |
rp->message = NULL; |
|
rp->t_success = 0; |
|
rp->series_error = 1.0e+10; |
|
rp->recommended_abserr = 1.0e-10; |
|
|
*rp=mh_rkmain(x0,y0,xn); |
*rp=mh_rkmain(x0,y0,xn); |
return(rp); |
return(rp); |
} |
} |