SCM

SCM Repository

[matrix] Diff of /pkg/src/ssclme.c
ViewVC logotype

Diff of /pkg/src/ssclme.c

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 165, Fri May 14 21:00:35 2004 UTC revision 166, Sat May 15 12:59:43 2004 UTC
# Line 1002  Line 1002 
1002    
1003      p = pp1 - 1;      p = pp1 - 1;
1004      b = RZX + p * n;      b = RZX + p * n;
     if (verbose) Rprintf("  EM iterations\n");  
     for (iter = 0; iter <= nEM; iter++) {  
         ssclme_invert(x);  
1005          if (verbose) {          if (verbose) {
1006              SEXP coef = PROTECT(ssclme_coef(x));              SEXP coef = PROTECT(ssclme_coef(x));
1007              int lc = length(coef); double *cc = REAL(coef);              int lc = length(coef); double *cc = REAL(coef);
1008              Rprintf("%3d %.3f", iter, dev[REML ? 1 : 0]);  
1009            Rprintf("  EM iterations\n");
1010            Rprintf("%3d %.3f", 0, dev[REML ? 1 : 0]);
1011              for (i = 0; i < lc; i++) Rprintf(" %#8g", cc[i]);              for (i = 0; i < lc; i++) Rprintf(" %#8g", cc[i]);
1012              Rprintf("\n");              Rprintf("\n");
1013              UNPROTECT(1);              UNPROTECT(1);
1014          }          }
1015        for (iter = 0; iter < nEM; iter++) {
1016            ssclme_invert(x);
1017          for (i = 0; i < nf; i++) {          for (i = 0; i < nf; i++) {
1018              int ki = Gp[i+1] - Gp[i],              int ki = Gp[i+1] - Gp[i],
1019                  nci = nc[i],                  nci = nc[i],
# Line 1046  Line 1047 
1047                        info, i + 1);                        info, i + 1);
1048          }          }
1049          status[0] = status[1] = 0;          status[0] = status[1] = 0;
1050            if (verbose) {
1051                SEXP coef = PROTECT(ssclme_coef(x));
1052                int lc = length(coef); double *cc = REAL(coef);
1053                Rprintf("%3d %.3f", iter + 1, dev[REML ? 1 : 0]);
1054                for (i = 0; i < lc; i++) Rprintf(" %#8g", cc[i]);
1055                Rprintf("\n");
1056                UNPROTECT(1);
1057            }
1058      }      }
1059      ssclme_factor(x);      ssclme_factor(x);
1060      return R_NilValue;      return R_NilValue;

Legend:
Removed from v.165  
changed lines
  Added in v.166

root@r-forge.r-project.org
ViewVC Help
Powered by ViewVC 1.0.0  
Thanks to:
Vienna University of Economics and Business Powered By FusionForge