SCM

SCM Repository

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

Diff of /pkg/src/dgBCMatrix.c

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

revision 437, Mon Jan 17 12:47:29 2005 UTC revision 438, Mon Jan 17 12:48:32 2005 UTC
# Line 341  Line 341 
341              for (k = Ap[j]; k < k2; k++) {              for (k = Ap[j]; k < k2; k++) {
342                  int ii = Ai[k], K = Ind(ii, ii, Cp, Ci);                  int ii = Ai[k], K = Ind(ii, ii, Cp, Ci);
343    
344                    if (K < 0) error("cscb_syrk: C[%d,%d] not defined", ii, ii);
345                  if (scalar) Cx[K] += alpha * Ax[k] * Ax[k];                  if (scalar) Cx[K] += alpha * Ax[k] * Ax[k];
346                  else F77_CALL(dsyrk)(&uplo, "N", cdims, adims + 1,                  else F77_CALL(dsyrk)(&uplo, "N", cdims, adims + 1,
347                                       &alpha, Ax + k * asz, adims,                                       &alpha, Ax + k * asz, adims,
# Line 348  Line 349 
349    
350                  for (kk = k+1; kk < k2; kk++) {                  for (kk = k+1; kk < k2; kk++) {
351                      int jj = Ai[kk];                      int jj = Ai[kk];
352                      K = (iup) ? Ind(jj, ii, Cp, Ci) : Ind(ii, jj, Cp, Ci);                      K = (iup) ? Ind(ii, jj, Cp, Ci) : Ind(jj, ii, Cp, Ci);
353    
354                        if (K < 0) error("cscb_syrk: C[%d,%d] not defined", ii, jj);
355                      if (scalar) Cx[K] += alpha * Ax[k] * Ax[kk];                      if (scalar) Cx[K] += alpha * Ax[k] * Ax[kk];
356                      else F77_CALL(dgemm)("N", "T", cdims, cdims + 1, adims + 1,                      else F77_CALL(dgemm)("N", "T", cdims, cdims + 1, adims + 1,
357                                           &alpha, Ax + ((iup)?kk:k) * asz, adims,                                           &alpha, Ax + ((iup)?kk:k) * asz, adims,

Legend:
Removed from v.437  
changed lines
  Added in v.438

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