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 112, Tue Apr 20 17:34:18 2004 UTC revision 113, Tue Apr 20 17:34:49 2004 UTC
# Line 708  Line 708 
708              }              }
709              tmp = Calloc(nr * nci, double); /* scratch storage */              tmp = Calloc(nr * nci, double); /* scratch storage */
710              nr1 = nr + 1;              nr1 = nr + 1;
711                                  /* initialize bVi to zero (cosmetic) */                                  /* initialize bVi to zero */
712              memset(bVi, 0, sizeof(double) * (G2 - G1) * nci);              memset(bVi, 0, sizeof(double) * (G2 - G1) * nci);
713              for (j = G1; j < G2; j += nci) {              for (j = G1; j < G2; j += nci) {
714                  memset(tmp, 0, sizeof(double) * nr * nci);                  memset(tmp, 0, sizeof(double) * nr * nci);
# Line 724  Line 724 
724                          tmp[k * nr + kk] *= DIsqrt[LIi[LIp[j] + kk - 1]];                          tmp[k * nr + kk] *= DIsqrt[LIi[LIp[j] + kk - 1]];
725                      }                      }
726                  }                  }
727                  F77_CALL(dsyrk)("U", "T", &nci, &rr, &one, tmp, &nr,                  F77_CALL(dsyrk)("L", "T", &nci, &rr, &one, tmp, &nr,
728                                  &zero, bVi + (j - G1) * nci, &nci);                                  &zero, bVi + (j - G1) * nci, &nci);
729                  F77_CALL(dpotrf)("U", &nci, bVi + (j - G1) * nci,                  F77_CALL(dpotrf)("L", &nci, bVi + (j - G1) * nci,
730                                   &nci, &kk);                                   &nci, &kk);
731                  if (kk)         /* should never happen */                  if (kk)         /* should never happen */
732                      error(                      error(
# Line 790  Line 790 
790                          tmp[k * nr + kk] *= DIsqrt[ind[kk]];                          tmp[k * nr + kk] *= DIsqrt[ind[kk]];
791                      }                      }
792                  }                  }
793                  F77_CALL(dsyrk)("U", "T", &nci, &nr, &one, tmp, &nr,                  F77_CALL(dsyrk)("L", "T", &nci, &nr, &one, tmp, &nr,
794                                  &zero, mpt + (j - Gp[i])*nci, &nci);                                  &zero, mpt + (j - Gp[i])*nci, &nci);
795                  F77_CALL(dpotrf)("U", &nci, mpt + (j - Gp[i])*nci,                  F77_CALL(dpotrf)("L", &nci, mpt + (j - Gp[i])*nci,
796                                   &nci, &info);                                   &nci, &info);
797                  if (info)       /* should never happen */                  if (info)       /* should never happen */
798                      error(                      error(
# Line 1246  Line 1246 
1246          REML = asLogical(REMLp),          REML = asLogical(REMLp),
1247          cind, i, n = dims[0],          cind, i, n = dims[0],
1248          nf = length(Omega),          nf = length(Omega),
1249          nobs, odind, p, pp1 = dims[1],          nobs, p, pp1 = dims[1],
1250          uncst = asLogical(Uncp);          uncst = asLogical(Uncp);
1251      double      double
1252          *RZX = REAL(RZXsl),          *RZX = REAL(RZXsl),

Legend:
Removed from v.112  
changed lines
  Added in v.113

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