[BACK]Return to mh.c CVS log [TXT][DIR] Up to [local] / OpenXM / src / hgm / mh / src

Diff for /OpenXM/src/hgm/mh/src/mh.c between version 1.6 and 1.7

version 1.6, 2013/03/05 06:00:52 version 1.7, 2013/03/05 06:35:54
Line 1 
Line 1 
 /* $OpenXM: OpenXM/src/hgm/mh/src/mh.c,v 1.5 2013/02/25 12:11:23 takayama Exp $ */  /* $OpenXM: OpenXM/src/hgm/mh/src/mh.c,v 1.6 2013/03/05 06:00:52 takayama Exp $ */
 #include <stdio.h>  #include <stdio.h>
 #include "sfile.h"  #include "sfile.h"
 #include "mh.h"  #include "mh.h"
Line 144  struct cWishart *mh_cwishart_hgm(int m,int n,double be
Line 144  struct cWishart *mh_cwishart_hgm(int m,int n,double be
 main() {  main() {
   double beta[5]={1.0,2.0,3.0,4.0,5.0};    double beta[5]={1.0,2.0,3.0,4.0,5.0};
   struct cWishart *cw;    struct cWishart *cw;
     struct SFILE *sfp;
     char *s;
     char str[1024];
     double x;
   cw=mh_cwishart_hgm(3,5,beta,0.3,7,  0.01,1,10);    cw=mh_cwishart_hgm(3,5,beta,0.3,7,  0.01,1,10);
   if (cw != NULL) {    if (cw != NULL) {
     printf("x=%lf, y=%lf\n",cw->x,(cw->f)[0]);      printf("x=%lf, y=%lf\n",cw->x,(cw->f)[0]);
Line 152  main() {
Line 156  main() {
   cw=mh_cwishart_hgm(4,5,beta,0.3,7,  0.01,1,10);    cw=mh_cwishart_hgm(4,5,beta,0.3,7,  0.01,1,10);
   if (cw != NULL) {    if (cw != NULL) {
     printf("x=%lf, y=%lf\n",cw->x,(cw->f)[0]);      printf("x=%lf, y=%lf\n",cw->x,(cw->f)[0]);
       s = (char *)cw->aux;
     /* printf("%s",(char *)cw->aux); */      /* printf("%s",(char *)cw->aux); */
       sfp = mh_fopen(s,"r",0);
       while (mh_fgets(str,1024,sfp)) {
         sscanf(str,"%lg",&x); printf("%lg\n",x);
       }
       mh_fclose(sfp);
   }    }
 }  }
 main1() {  main1() {

Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7

FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>