det, inv, transpose and submatrix

Author: rajithv

ext/symengine/ruby_matrix.c

@@ -229,3 +229,82 @@ VALUE cmatrix_dense_set(VALUE self, VALUE r, VALUE c, VALUE operand)
 
     return self;
 }
+
+
+VALUE cmatrix_dense_det(VALUE self)
+{
+    CDenseMatrix *this;
+    Data_Get_Struct(self, CDenseMatrix, this);
+    
+    basic_struct *cresult;
+    VALUE result;
+
+    cresult = basic_new_heap();
+    
+    dense_matrix_det(cresult, this);
+    result = Data_Wrap_Struct(Klass_of_Basic(cresult), NULL, basic_free_heap,
+                              cresult);
+    return result;
+}
+VALUE cmatrix_dense_inv(VALUE self)
+{
+    CDenseMatrix *this;
+    Data_Get_Struct(self, CDenseMatrix, this);
+    
+    CDenseMatrix *cresult;
+    VALUE result;
+
+    cresult = dense_matrix_new();
+    
+    dense_matrix_inv(cresult, this);
+    result = Data_Wrap_Struct(c_dense_matrix, NULL, dense_matrix_free,
+                              cresult);
+    return result;
+}
+
+VALUE cmatrix_dense_transpose(VALUE self)
+{
+    CDenseMatrix *this;
+    Data_Get_Struct(self, CDenseMatrix, this);
+    
+    CDenseMatrix *cresult;
+    VALUE result;
+
+    cresult = dense_matrix_new();
+    
+    dense_matrix_transpose(cresult, this);
+    result = Data_Wrap_Struct(c_dense_matrix, NULL, dense_matrix_free,
+                              cresult);
+    return result;
+}
+
+VALUE cmatrix_dense_submatrix(VALUE self, VALUE r1, VALUE c1, VALUE r2, VALUE c2, VALUE r_, VALUE c_)
+{
+    CDenseMatrix *this;
+    Data_Get_Struct(self, CDenseMatrix, this);
+    
+    CDenseMatrix *cresult;
+    VALUE result;
+    
+    unsigned long cbasic_r1;
+    cbasic_r1 = NUM2ULONG(r1);
+    unsigned long cbasic_c1;
+    cbasic_c1 = NUM2ULONG(c1);
+    
+    unsigned long cbasic_r2;
+    cbasic_r2 = NUM2ULONG(r2);
+    unsigned long cbasic_c2;
+    cbasic_c2 = NUM2ULONG(c2);
+    
+    unsigned long cbasic_r;
+    cbasic_r = NUM2ULONG(r_);
+    unsigned long cbasic_c;
+    cbasic_c = NUM2ULONG(c_);
+
+    cresult = dense_matrix_new();
+    
+    dense_matrix_submatrix(cresult, this, cbasic_r1, cbasic_c1, cbasic_r2, cbasic_c2, cbasic_r, cbasic_c);
+    result = Data_Wrap_Struct(c_dense_matrix, NULL, dense_matrix_free,
+                              cresult);
+    return result;
+}

ext/symengine/ruby_matrix.h

@@ -14,6 +14,10 @@ VALUE cmatrix_dense_init(VALUE self, VALUE args);
 VALUE cmatrix_dense_to_str(VALUE self);
 VALUE cmatrix_dense_get(VALUE self, VALUE r, VALUE c);
 VALUE cmatrix_dense_set(VALUE self, VALUE r, VALUE c, VALUE operand);
+VALUE cmatrix_dense_det(VALUE self);
+VALUE cmatrix_dense_inv(VALUE self);
+VALUE cmatrix_dense_transpose(VALUE self);
+VALUE cmatrix_dense_submatrix(VALUE self, VALUE r1, VALUE c1, VALUE r2, VALUE c2, VALUE r_, VALUE c_);
 
 void cmatrix_sparse_free(void *ptr);
 VALUE cmatrix_sparse_alloc(VALUE klass);

ext/symengine/symengine.c

@@ -239,6 +239,10 @@ void Init_symengine()
     rb_define_method(c_dense_matrix, "to_s", cmatrix_dense_to_str, 0);
     rb_define_method(c_dense_matrix, "get", cmatrix_dense_get, 2);
     rb_define_method(c_dense_matrix, "set", cmatrix_dense_set, 3);
+    rb_define_method(c_dense_matrix, "det", cmatrix_dense_det, 0);
+    rb_define_method(c_dense_matrix, "inv", cmatrix_dense_inv, 0);
+    rb_define_method(c_dense_matrix, "transpose", cmatrix_dense_transpose, 0);
+    rb_define_method(c_dense_matrix, "submatrix", cmatrix_dense_submatrix, 6);
 
     // SparseMatrix Methods
     rb_define_alloc_func(c_sparse_matrix, cmatrix_sparse_alloc);