SCM

SCM Repository

[matrix] Diff of /pkg/Matrix/src/Mutils.h
ViewVC logotype

Diff of /pkg/Matrix/src/Mutils.h

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

revision 2744, Fri Dec 30 18:51:18 2011 UTC revision 2850, Mon Oct 22 13:00:35 2012 UTC
# Line 30  Line 30 
30  # include <alloca.h>  # include <alloca.h>
31  #endif  #endif
32    
33    #ifndef LONG_VECTOR_SUPPORT
34    // notably for  R <= 2.15.x :
35    # define XLENGTH(x) LENGTH(x)
36    # if R_VERSION < R_Version(2,16,0)
37      typedef int R_xlen_t;
38    # endif
39    #endif
40    
41  #define Alloca(n, t)   (t *) alloca( (size_t) ( (n) * sizeof(t) ) )  #define Alloca(n, t)   (t *) alloca( (size_t) ( (n) * sizeof(t) ) )
42    
43  SEXP triangularMatrix_validate(SEXP obj);  SEXP triangularMatrix_validate(SEXP obj);
44  SEXP symmetricMatrix_validate(SEXP obj);  SEXP symmetricMatrix_validate(SEXP obj);
45  SEXP dense_nonpacked_validate(SEXP obj);  SEXP dense_nonpacked_validate(SEXP obj);
46    
47    // La_norm_type() & La_rcond_type()  have been in R_ext/Lapack.h
48    //  but have still not been available to package writers ...
49    char La_norm_type (const char *typstr);
50    char La_rcond_type(const char *typstr);
51    
52  /* enum constants from cblas.h and some short forms */  /* enum constants from cblas.h and some short forms */
53  enum CBLAS_ORDER {CblasRowMajor=101, CblasColMajor=102};  enum CBLAS_ORDER {CblasRowMajor=101, CblasColMajor=102};
54  enum CBLAS_TRANSPOSE {CblasNoTrans=111, CblasTrans=112, CblasConjTrans=113};  enum CBLAS_TRANSPOSE {CblasNoTrans=111, CblasTrans=112, CblasConjTrans=113};
# Line 73  Line 86 
86    
87  void d_packed_getDiag(double *dest, SEXP x, int n);  void d_packed_getDiag(double *dest, SEXP x, int n);
88  void l_packed_getDiag(   int *dest, SEXP x, int n);  void l_packed_getDiag(   int *dest, SEXP x, int n);
89  void tr_d_packed_getDiag(double *dest, SEXP x);  SEXP d_packed_setDiag(double *diag, int l_d, SEXP x, int n);
90  void tr_l_packed_getDiag(   int *dest, SEXP x);  SEXP l_packed_setDiag(   int *diag, int l_d, SEXP x, int n);
91    SEXP d_packed_addDiag(double *diag, int l_d, SEXP x, int n);
92    
93    void tr_d_packed_getDiag(double *dest, SEXP x, int n);
94    void tr_l_packed_getDiag(   int *dest, SEXP x, int n);
95    
96    SEXP tr_d_packed_setDiag(double *diag, int l_d, SEXP x, int n);
97    SEXP tr_l_packed_setDiag(   int *diag, int l_d, SEXP x, int n);
98    SEXP tr_d_packed_addDiag(double *diag, int l_d, SEXP x, int n);
99    
100  SEXP Matrix_getElement(SEXP list, char *nm);  SEXP Matrix_getElement(SEXP list, char *nm);
101    
# Line 122  Line 143 
143    
144  #define uplo_P(_x_) CHAR(STRING_ELT(GET_SLOT(_x_, Matrix_uploSym), 0))  #define uplo_P(_x_) CHAR(STRING_ELT(GET_SLOT(_x_, Matrix_uploSym), 0))
145  #define diag_P(_x_) CHAR(STRING_ELT(GET_SLOT(_x_, Matrix_diagSym), 0))  #define diag_P(_x_) CHAR(STRING_ELT(GET_SLOT(_x_, Matrix_diagSym), 0))
146    #define Diag_P(_x_) (R_has_slot(x, Matrix_diagSym) ?                    \
147                         CHAR(STRING_ELT(GET_SLOT(_x_, Matrix_diagSym), 0)) : " ")
148  #define class_P(_x_) CHAR(asChar(getAttrib(_x_, R_ClassSymbol)))  #define class_P(_x_) CHAR(asChar(getAttrib(_x_, R_ClassSymbol)))
149    
150    
151  // Define this "Cholmod compatible" to some degree  // Define this "Cholmod compatible" to some degree
152  enum x_slot_kind {x_pattern=-1, x_double=0, x_logical=1, x_integer=2, x_complex=3};  enum x_slot_kind {x_pattern=-1, x_double=0, x_logical=1, x_integer=2, x_complex=3};
153  //                  n             d           l            i            z  //                  n             d           l            i            z
# Line 228  Line 252 
252      return FALSE;      return FALSE;
253  }  }
254    
255    SEXP Mmatrix(SEXP args);
256    
257  void make_d_matrix_triangular(double *x, SEXP from);  void make_d_matrix_triangular(double *x, SEXP from);
258  void make_i_matrix_triangular(   int *x, SEXP from);  void make_i_matrix_triangular(   int *x, SEXP from);
# Line 307  Line 332 
332   * These are the ones "everyone" should use -- is() versions, also looking   * These are the ones "everyone" should use -- is() versions, also looking
333   * at super classes:   * at super classes:
334   */   */
 #if R_VERSION < R_Version(2, 15, 0)  
 int Matrix_check_class_etc(SEXP x, const char **valid);  
 #else  
335  # define Matrix_check_class_etc R_check_class_etc  # define Matrix_check_class_etc R_check_class_etc
 #endif  
   
 #if R_VERSION < R_Version(2, 13, 0)  
 int Matrix_check_class_and_super(SEXP x, const char **valid, SEXP rho);  
 #else  
336  # define Matrix_check_class_and_super R_check_class_and_super  # define Matrix_check_class_and_super R_check_class_and_super
 #endif  
   
337    
338  /** Accessing  *sparseVectors :  fast (and recycling)  v[i] for v = ?sparseVector:  /** Accessing  *sparseVectors :  fast (and recycling)  v[i] for v = ?sparseVector:
339   * -> ./sparseVector.c  -> ./t_sparseVector.c :   * -> ./sparseVector.c  -> ./t_sparseVector.c :

Legend:
Removed from v.2744  
changed lines
  Added in v.2850

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