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 166, Sat May 15 12:59:43 2004 UTC revision 168, Wed May 19 16:18:48 2004 UTC
# Line 73  Line 73 
73          for (i = 0; i < nf; i++) { /* fill in the rows */          for (i = 0; i < nf; i++) { /* fill in the rows */
74              int j, jj, nci = nc[i], p2 = GpIn[i+1];              int j, jj, nci = nc[i], p2 = GpIn[i+1];
75              for (j = GpIn[i]; j < p2; j++) { /* first col for input col */              for (j = GpIn[i]; j < p2; j++) { /* first col for input col */
76                  int ii = AiIn[j], mj = map[j], ncci = ncc[ii],                  int k, mj = map[j], p3 = ApIn[j+1], pos = ApOut[mj];
77                      pos = ApOut[mj];                  for (k = ApIn[j]; k < p3; k++) {
78                        int ii = AiIn[k], ncci = ncc[ii];
79                  AiOut[pos++] = map[ii];                  AiOut[pos++] = map[ii];
80                  if (ii < j) {   /* above the diagonal */                  if (ii < j) {   /* above the diagonal */
81                      for (jj = 1; jj < ncci; jj++) {                      for (jj = 1; jj < ncci; jj++) {
# Line 90  Line 91 
91                  }                  }
92              }              }
93          }          }
94            }
95          Free(map); Free(ncc);          Free(map); Free(ncc);
96      }      }
97      SET_SLOT(ssc, Matrix_ncSym, allocVector(INTSXP, nf + 2));      SET_SLOT(ssc, Matrix_ncSym, allocVector(INTSXP, nf + 2));
# Line 119  Line 121 
121  }  }
122    
123  SEXP  SEXP
124  ssclme_create(SEXP facs, SEXP ncv, SEXP threshold)  ssclme_create(SEXP facs, SEXP ncv)
125  {  {
126      SEXP ctab, nms, ssc, tmp,      SEXP ctab, nms, ssc, tmp,
127          val = PROTECT(allocVector(VECSXP, 2)),          val = PROTECT(allocVector(VECSXP, 2)),
# Line 1006  Line 1008 
1008          SEXP coef = PROTECT(ssclme_coef(x));          SEXP coef = PROTECT(ssclme_coef(x));
1009          int lc = length(coef); double *cc = REAL(coef);          int lc = length(coef); double *cc = REAL(coef);
1010    
1011            ssclme_factor(x);
1012          Rprintf("  EM iterations\n");          Rprintf("  EM iterations\n");
1013          Rprintf("%3d %.3f", 0, dev[REML ? 1 : 0]);          Rprintf("%3d %.3f", 0, dev[REML ? 1 : 0]);
1014          for (i = 0; i < lc; i++) Rprintf(" %#8g", cc[i]);          for (i = 0; i < lc; i++) Rprintf(" %#8g", cc[i]);
# Line 1050  Line 1053 
1053          if (verbose) {          if (verbose) {
1054              SEXP coef = PROTECT(ssclme_coef(x));              SEXP coef = PROTECT(ssclme_coef(x));
1055              int lc = length(coef); double *cc = REAL(coef);              int lc = length(coef); double *cc = REAL(coef);
1056    
1057                ssclme_factor(x);
1058              Rprintf("%3d %.3f", iter + 1, dev[REML ? 1 : 0]);              Rprintf("%3d %.3f", iter + 1, dev[REML ? 1 : 0]);
1059              for (i = 0; i < lc; i++) Rprintf(" %#8g", cc[i]);              for (i = 0; i < lc; i++) Rprintf(" %#8g", cc[i]);
1060              Rprintf("\n");              Rprintf("\n");

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

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