use c++ latex rendering

Author: eendebakpt

symengine/lib/symengine_wrapper.pyx

@@ -3242,32 +3242,10 @@ cdef class DenseMatrixBase(MatrixBase):
         return deref(self.thisptr).__str__().decode("utf-8")
 
     def _repr_latex_(self):
-        MAX_NUMBER_OF_ROWS = 24
-        MAX_NUMBER_OF_COLUMNS = 16
-
-        values = list(self)
-        ncols=self.shape[1]
-        nrows = self.shape[0]
-        if nrows > MAX_NUMBER_OF_ROWS:
-            nrows_display = MAX_NUMBER_OF_ROWS - 2
+        if repr_latex[0]:
+            return latex(deref(self.thisptr)).decode("utf-8")
         else:
-            nrows_display= nrows
-        ncols_display = min(ncols, MAX_NUMBER_OF_COLUMNS)
-        latex = r'$\displaystyle \left[\begin{matrix}'
-
-        newline = r'\\'
-        if ncols_display<ncols:
-            newline = ' & \cdots '  + newline
-        for row in range(nrows_display):
-            vv = values[row*ncols:(row*ncols)+ncols_display]
-            latex += ' & '.join([str(v) for v in vv])
-            if row < self.shape[0]-1:
-                latex += newline
-        if nrows_display < nrows:
-            latex += ' & '.join([r'\vdots' for v in vv])
-
-        latex += r'\end{matrix}\right]$'
-        return latex
+            return None
 
     def __add__(a, b):
         a = _sympify(a, False)
@@ -5440,8 +5418,12 @@ def cse(exprs):
     return (vec_pair_to_list(replacements), vec_basic_to_list(reduced_exprs))
 
 def latex(expr):
-    cdef Basic expr_ = sympify(expr)
-    return symengine.latex(deref(expr_.thisptr)).decode("utf-8")
+    if isinstance(expr, DenseMatrixBase):
+        cdef symengine.DenseMatrix expr_ = sympify(expr)
+        return symengine.latex(deref(expr_.thisptr)).decode("utf-8")
+    else:
+        cdef Basic expr_ = sympify(expr)
+        return symengine.latex(deref(expr_.thisptr)).decode("utf-8")
 
 cdef _flattened_vec(symengine.vec_basic &vec, exprs):
     cdef Basic b