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 75, Wed Apr 14 20:43:06 2004 UTC revision 76, Wed Apr 14 20:43:40 2004 UTC
# Line 1191  Line 1191 
1191                              &alpha, REAL(VECTOR_ELT(bVar, i)), &nci,                              &alpha, REAL(VECTOR_ELT(bVar, i)), &nci,
1192                              &one, vali, &nci);                              &one, vali, &nci);
1193              if (REML) {              if (REML) {
1194                  int mp = mi * p;                  int j;
1195                  F77_CALL(dsyrk)("U", "N", &nci, &mp,                  for (j = 0; j < p; j++) {
1196                                  &alpha, RZX + Gp[i], &nci,                      F77_CALL(dsyrk)("U", "N", &nci, &mi,
1197                                    &alpha, RZX + Gp[i] + j*n, &nci,
1198                                  &one, vali, &nci);                                  &one, vali, &nci);
1199              }              }
1200                }
1201              F77_CALL(dpotrf)("U", &nci, vali, &nci, &info);              F77_CALL(dpotrf)("U", &nci, vali, &nci, &info);
1202              if (info) error("DPOTRF returned error code %d", info);              if (info) error("DPOTRF returned error code %d", info);
1203              F77_CALL(dpotri)("U", &nci, vali, &nci, &info);              F77_CALL(dpotri)("U", &nci, vali, &nci, &info);

Legend:
Removed from v.75  
changed lines
  Added in v.76

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