SCM

SCM Repository

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

Diff of /pkg/src/dtpMatrix.c

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

revision 628, Thu Mar 10 17:57:47 2005 UTC revision 629, Thu Mar 10 19:00:57 2005 UTC
# Line 130  Line 130 
130      if (yDim[0] != xDim[1])      if (yDim[0] != xDim[1])
131          error(_("Dimensions of a (%d,%d) and b (%d,%d) do not conform"),          error(_("Dimensions of a (%d,%d) and b (%d,%d) do not conform"),
132                xDim[0], xDim[1], yDim[0], yDim[1]);                xDim[0], xDim[1], yDim[0], yDim[1]);
133      for (j = 0; j < yDim[1]; j++) /* A %*% x  via BLAS 2 DTPMV(.) */      for (j = 0; j < yDim[1]; j++) /* X %*% y[,j]  via BLAS 2 DTPMV(.) */
134          F77_CALL(dtpmv)(uplo, "N", diag, yDim, xx,          F77_CALL(dtpmv)(uplo, "N", diag, yDim, xx,
135                          vx + j * yDim[0], &ione);                          vx + j * yDim[0], &ione);
136      UNPROTECT(1);      UNPROTECT(1);
# Line 154  Line 154 
154                xDim[0], xDim[1], yDim[0], yDim[1]);                xDim[0], xDim[1], yDim[0], yDim[1]);
155      for (i = 0; i < xDim[0]; i++)/* val[i,] := Y' %*% x[i,]  */      for (i = 0; i < xDim[0]; i++)/* val[i,] := Y' %*% x[i,]  */
156          F77_CALL(dtpmv)(uplo, "T", diag, yDim, yx,          F77_CALL(dtpmv)(uplo, "T", diag, yDim, yx,
157                          vx + i * xDim[1], /* incr = */ xDim);                          vx + i, /* incr = */ xDim);
158      UNPROTECT(1);      UNPROTECT(1);
159      return val;      return val;
160  }  }

Legend:
Removed from v.628  
changed lines
  Added in v.629

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