towards a decent printing system for all matrix like objects

Author: counterclocker

README.md

@@ -13,7 +13,7 @@ This is the multi-types version (last single type version: ce86f8476166f63d0af72
    my own source format...).
 
 2. Cython memoryviews don't seem to be faster in our case. This should be assessed in the optimization phase. For the moment, I'll continue
-   to code with NumPy's C-API. This is **not** the way Cython is intended by my tests indicate that this produces the fastest code. To be continued.
+   to code with NumPy's C-API. This is **not** the way Cython is intended but my tests indicate that this produces the fastest code. To be continued.
 
    Sorry for the fuzz.
 

cysparse/sparse/ll_mat_matrices/ll_mat.cpx

@@ -19,7 +19,7 @@ from cysparse.sparse.csr_mat_matrices.csr_mat_@index@_@type@ cimport MakeCSRSpar
 from cysparse.sparse.csc_mat_matrices.csc_mat_@index@_@type@ cimport MakeCSCSparseMatrix_@index@_@type@
 
 from cysparse.sparse.sparse_utils.generic.generate_indices_@index@ cimport create_c_array_indices_from_python_object_@index@
-from cysparse.sparse.sparse_utils.generic.print_@type@ cimport element_to_string_@type@, empty_to_string_@type@
+from cysparse.sparse.sparse_utils.generic.print_@type@ cimport element_to_string_@type@, conjugated_element_to_string_@type@, empty_to_string_@type@
 
 ########################################################################################################################
 # CySparse include
@@ -2420,6 +2420,33 @@ YOU HAVE TO CAST YOUR NEW TYPE HERE
         # element not found -> return empty cell
         return empty_to_string_@type@(cell_width=cell_width)
 
+    def at_conj_to_string(self, @index@ i, @index@ j, int cell_width=10):
+        """
+        Return a string with a given element if it exists or an "empty" string.
+
+
+        """
+        cdef:
+            @index@ k, next_col
+            @type@ v
+            FLOAT64_t exp
+
+        # EXPLICIT TYPE TESTS
+        k = self.root[i]
+        while k != -1:
+            next_col = self.col[k]
+            if next_col >= j:
+                if next_col == j:
+                    #v = self.val[k]
+                    return conjugated_element_to_string_@type@(self.val[k], cell_width=cell_width)
+                else:  # value not found
+                    break
+
+            k = self.link[k]
+
+        # element not found -> return empty cell
+        return empty_to_string_@type@(cell_width=cell_width)
+
     def print_to(self, OUT, width=9, print_big_matrices=False, transposed=False):
         """
         Print content of matrix to output stream.

cysparse/sparse/ll_mat_matrices/ll_mat_INT32_t_COMPLEX128_t.c

@@ -19,7 +19,7 @@ from cysparse.sparse.csr_mat_matrices.csr_mat_INT32_t_COMPLEX128_t cimport MakeC
 from cysparse.sparse.csc_mat_matrices.csc_mat_INT32_t_COMPLEX128_t cimport MakeCSCSparseMatrix_INT32_t_COMPLEX128_t
 
 from cysparse.sparse.sparse_utils.generic.generate_indices_INT32_t cimport create_c_array_indices_from_python_object_INT32_t
-from cysparse.sparse.sparse_utils.generic.print_COMPLEX128_t cimport element_to_string_COMPLEX128_t, empty_to_string_COMPLEX128_t
+from cysparse.sparse.sparse_utils.generic.print_COMPLEX128_t cimport element_to_string_COMPLEX128_t, conjugated_element_to_string_COMPLEX128_t, empty_to_string_COMPLEX128_t
 
 ########################################################################################################################
 # CySparse include
@@ -2246,6 +2246,33 @@ cdef class LLSparseMatrix_INT32_t_COMPLEX128_t(MutableSparseMatrix_INT32_t_COMPL
         # element not found -> return empty cell
         return empty_to_string_COMPLEX128_t(cell_width=cell_width)
 
+    def at_conj_to_string(self, INT32_t i, INT32_t j, int cell_width=10):
+        """
+        Return a string with a given element if it exists or an "empty" string.
+
+
+        """
+        cdef:
+            INT32_t k, next_col
+            COMPLEX128_t v
+            FLOAT64_t exp
+
+        # EXPLICIT TYPE TESTS
+        k = self.root[i]
+        while k != -1:
+            next_col = self.col[k]
+            if next_col >= j:
+                if next_col == j:
+                    #v = self.val[k]
+                    return conjugated_element_to_string_COMPLEX128_t(self.val[k], cell_width=cell_width)
+                else:  # value not found
+                    break
+
+            k = self.link[k]
+
+        # element not found -> return empty cell
+        return empty_to_string_COMPLEX128_t(cell_width=cell_width)
+
     def print_to(self, OUT, width=9, print_big_matrices=False, transposed=False):
         """
         Print content of matrix to output stream.