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 559, Fri Feb 18 14:29:03 2005 UTC revision 565, Tue Feb 22 01:15:45 2005 UTC
# Line 123  Line 123 
123   * where A is a compressed, sparse, blocked matrix and   * where A is a compressed, sparse, blocked matrix and
124   * C is a compressed, sparse, symmetric blocked matrix.   * C is a compressed, sparse, symmetric blocked matrix.
125   *   *
126   * @param uplo 'U' or 'u' for upper triangular storage, else lower.   * @param uplo UPP or LOW for upper or lower
127   * @param trans 'T' or 't' for transpose.   * @param trans TRN or NTR for transpose or no transpose
128   * @param alpha scalar multiplier of outer product   * @param alpha scalar multiplier of outer product
129   * @param A compressed sparse blocked matrix   * @param A compressed sparse blocked matrix
130   * @param beta scalar multiplier of c   * @param beta scalar multiplier of c
# Line 156  Line 156 
156      if (cdims[0] != cdims[1]) error("blocks in C must be square");      if (cdims[0] != cdims[1]) error("blocks in C must be square");
157      if (trans == TRN) {      if (trans == TRN) {
158                                  /* FIXME: Write this part */                                  /* FIXME: Write this part */
159          error("Code for trans == 'T' not yet written");          error("Code for trans == TRN not yet written");
160      } else {      } else {
161          if (adims[0] != cdims[0])          if (adims[0] != cdims[0])
162              error("Inconsistent dimensions: A[%d,%d,%d], C[%d,%d,%d]",              error("Inconsistent dimensions: A[%d,%d,%d], C[%d,%d,%d]",
# Line 335  Line 335 
335   * Perform one of the dgBCMatrix-matrix operations B := alpha*op(A)*B   * Perform one of the dgBCMatrix-matrix operations B := alpha*op(A)*B
336   * or B := alpha*B*op(A)   * or B := alpha*B*op(A)
337   *   *
338   * @param side   * @param side LFT or RGT for left or right
339   * @param uplo   * @param uplo UPP or LOW for upper or lower
340   * @param transa   * @param transa TRN or NTR for transpose or no transpose
341   * @param diag   * @param diag UNT or NUN for unit or non-unit
342     * @param alpha scalar multiplier
343   * @param A pointer to a triangular dgBCMatrix object   * @param A pointer to a triangular dgBCMatrix object
344   * @param B contents of the matrix B   * @param B contents of the matrix B
345   * @param m number of rows in B   * @param m number of rows in B
# Line 374  Line 375 
375   * is a dgBCMatrix triangular matrix and B is a dense matrix.   * is a dgBCMatrix triangular matrix and B is a dense matrix.
376   *   *
377   * @param uplo UPP or LOW   * @param uplo UPP or LOW
378   * @param trans TRN or NTR   * @param transa TRN or NTR
379   * @param diag UNT or NUN   * @param diag UNT or NUN
380     * @param alpha scalar multiplier
381   * @param A pointer to a triangular dgBCMatrix object   * @param A pointer to a triangular dgBCMatrix object
382   * @param m number of rows in B   * @param m number of rows in B
383   * @param n number of columns in B   * @param n number of columns in B
# Line 497  Line 499 
499   * Solve one of the systems op(A)*X = alpha*B or   * Solve one of the systems op(A)*X = alpha*B or
500   * X*op(A) = alpha*B where A dgBCMatrix triangular and B is dgBCMatrix.   * X*op(A) = alpha*B where A dgBCMatrix triangular and B is dgBCMatrix.
501   *   *
502   * @param side 'L' or 'R' for left or right   * @param side LFT or RGT for left or right
503   * @param uplo 'U' or 'L' for upper or lower   * @param uplo UPP or LOW for upper or lower
504   * @param transa 'T' or 'N' for transpose or no transpose   * @param transa TRN or NTR for transpose or no transpose
505   * @param diag 'U' or 'N' for unit diagonal or non-unit   * @param diag UNT or NON for unit or non-unit
506   * @param alpha scalar multiplier   * @param alpha scalar multiplier
507   * @param A pointer to a triangular dgBCMatrix object   * @param A pointer to a triangular dgBCMatrix object
508     * @param Parent parent array for the column blocks of A
509   * @param B pointer to a general dgBCMatrix matrix   * @param B pointer to a general dgBCMatrix matrix
510   */   */
511  void  void
# Line 572  Line 575 
575   *      C := alpha*op(A)*op(B) + beta*C   *      C := alpha*op(A)*op(B) + beta*C
576   * on compressed, sparse, blocked matrices.   * on compressed, sparse, blocked matrices.
577   *   *
578   * @param transa 'T' for transpose of A, else 'N'   * @param transa TRN or NTR for transpose or no transpose of A
579   * @param transb 'T' for transpose of B, else 'N'   * @param transb TRN or NTR for transpose or no transpose of B
580   * @param alpha scalar multiplier   * @param alpha scalar multiplier
581   * @param A pointer to a dgBCMatrix object   * @param A pointer to a dgBCMatrix object
582   * @param B pointer to a dgBCMatrix object   * @param B pointer to a dgBCMatrix object
# Line 641  Line 644 
644      error("Code not yet written");      error("Code not yet written");
645  }  }
646    
647    /**
648     * Coerce a dgBCMatrix to a dgCMatrix
649     *
650     * @param A pointer to a dgBCMatrix object to coerce
651     *
652     * @return pointer to a dgCMatrix
653     */
654  SEXP dgBCMatrix_to_dgCMatrix(SEXP A)  SEXP dgBCMatrix_to_dgCMatrix(SEXP A)
655  {  {
656      SEXP val = PROTECT(NEW_OBJECT(MAKE_CLASS("dgCMatrix"))),      SEXP val = PROTECT(NEW_OBJECT(MAKE_CLASS("dgCMatrix"))),

Legend:
Removed from v.559  
changed lines
  Added in v.565

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