Make symengine.lib.symengine_wrapper compatible with sympy

symengine/CMakeLists.txt

@@ -2,6 +2,6 @@ add_subdirectory(lib)
 add_subdirectory(tests)
 
 set(PY_PATH ${PYTHON_INSTALL_PATH}/symengine)
-install(FILES __init__.py utilities.py compatibility.py sympy_compat.py
+install(FILES __init__.py utilities.py compatibility.py sympy_compat.py functions.py
     DESTINATION ${PY_PATH}
     )

symengine/__init__.py

@@ -1,12 +1,15 @@
-from .lib.symengine_wrapper import (Symbol, Integer, sympify, SympifyError,
-        Add, Mul, Pow, exp, log, gamma, sqrt, function_symbol, I, E, pi,
-        have_mpfr, have_mpc, have_flint, have_piranha, have_llvm,
-        RealDouble, ComplexDouble, MutableDenseMatrix, ImmutableDenseMatrix,
-        DenseMatrix, Matrix, sin, cos, tan, cot, csc, sec, asin, acos, atan, 
-        acot, acsc, asec, sinh, cosh, tanh, coth, asinh, acosh, atanh, acoth, 
-        Lambdify, LambdifyCSE, DictBasic, series, symarray, diff, zeros, eye, 
-        diag, ones, zeros, add, expand, has_symbol, UndefFunction)
+from .lib.symengine_wrapper import (Symbol, sympify as S, sympify,
+        SympifyError, Add, Mul, Pow, function_symbol, I, E, pi,
+        have_mpfr, have_mpc, have_flint, have_piranha, have_llvm, 
+        Integer, Rational, Float, Number, RealNumber, RealDouble,
+        ComplexDouble, Max, Min, DenseMatrix, Matrix, ImmutableMatrix,
+        ImmutableDenseMatrix, MutableDenseMatrix, MatrixBase, Basic,
+        Lambdify, LambdifyCSE, Lambdify as lambdify, DictBasic, symarray,
+        series, diff, zeros, eye, diag, ones, Derivative, Subs,
+        add, expand, has_symbol, UndefFunction, Function,
+        FunctionSymbol as AppliedUndef)
 from .utilities import var, symbols
+from .functions import *
 
 if have_mpfr:
     from .lib.symengine_wrapper import RealMPFR

symengine/functions.py

@@ -0,0 +1,5 @@
+from .lib.symengine_wrapper import (sin, cos, tan, cot, csc, sec,
+    asin, acos, atan, acot, acsc, asec,
+    sinh, cosh, tanh, coth, sech, csch,
+    asinh, acosh, atanh, acoth, asech, acsch,
+    gamma, log, atan2, sqrt, exp, Abs)

symengine/lib/symengine.pxd

@@ -138,8 +138,7 @@ cdef extern from "<symengine/symengine_rcp.h>" namespace "SymEngine":
     RCP[const Add] rcp_static_cast_Add "SymEngine::rcp_static_cast<const SymEngine::Add>"(RCP[const Basic] &b) nogil
     RCP[const Mul] rcp_static_cast_Mul "SymEngine::rcp_static_cast<const SymEngine::Mul>"(RCP[const Basic] &b) nogil
     RCP[const Pow] rcp_static_cast_Pow "SymEngine::rcp_static_cast<const SymEngine::Pow>"(RCP[const Basic] &b) nogil
-    RCP[const TrigFunction] rcp_static_cast_TrigFunction "SymEngine::rcp_static_cast<const SymEngine::TrigFunction>"(RCP[const Basic] &b) nogil
-    RCP[const HyperbolicFunction] rcp_static_cast_HyperbolicFunction "SymEngine::rcp_static_cast<const SymEngine::HyperbolicFunction>"(RCP[const Basic] &b) nogil
+    RCP[const OneArgFunction] rcp_static_cast_OneArgFunction "SymEngine::rcp_static_cast<const SymEngine::OneArgFunction>"(RCP[const Basic] &b) nogil
     RCP[const FunctionSymbol] rcp_static_cast_FunctionSymbol "SymEngine::rcp_static_cast<const SymEngine::FunctionSymbol>"(RCP[const Basic] &b) nogil
     RCP[const FunctionWrapper] rcp_static_cast_FunctionWrapper "SymEngine::rcp_static_cast<const SymEngine::FunctionWrapper>"(RCP[const Basic] &b) nogil
     RCP[const Abs] rcp_static_cast_Abs "SymEngine::rcp_static_cast<const SymEngine::Abs>"(RCP[const Basic] &b) nogil
@@ -445,9 +444,12 @@ cdef extern from "<symengine/functions.h>" namespace "SymEngine":
     cdef cppclass Function(Basic):
         pass
 
-    cdef cppclass TrigFunction(Function):
+    cdef cppclass OneArgFunction(Function):
         RCP[const Basic] get_arg() nogil
 
+    cdef cppclass TrigFunction(OneArgFunction):
+        pass
+
     cdef cppclass Sin(TrigFunction):
         pass
 
@@ -484,8 +486,8 @@ cdef extern from "<symengine/functions.h>" namespace "SymEngine":
     cdef cppclass ASec(TrigFunction):
         pass
 
-    cdef cppclass HyperbolicFunction(Function):
-        RCP[const Basic] get_arg() nogil
+    cdef cppclass HyperbolicFunction(OneArgFunction):
+        pass
 
     cdef cppclass Sinh(HyperbolicFunction):
         pass
@@ -542,16 +544,16 @@ cdef extern from "<symengine/functions.h>" namespace "SymEngine":
         vec_basic get_variables() nogil
         vec_basic get_point() nogil
 
-    cdef cppclass Abs(Function):
-        RCP[const Basic] get_arg() nogil
+    cdef cppclass Abs(OneArgFunction):
+        pass
 
     cdef cppclass Max(Function):
         pass
 
     cdef cppclass Min(Function):
         pass
 
-    cdef cppclass Gamma(Function):
+    cdef cppclass Gamma(OneArgFunction):
         pass
 
     cdef cppclass ATan2(Function):
@@ -656,17 +658,17 @@ cdef extern from "<symengine/matrix.h>" namespace "SymEngine":
     bool is_a_DenseMatrix "SymEngine::is_a<SymEngine::DenseMatrix>"(const MatrixBase &b) nogil
     DenseMatrix* static_cast_DenseMatrix "static_cast<SymEngine::DenseMatrix*>"(const MatrixBase *a)
     void inverse_FFLU "SymEngine::inverse_fraction_free_LU"(const DenseMatrix &A,
-        DenseMatrix &B) nogil
-    void pivoted_LU (const DenseMatrix &A, DenseMatrix &L, DenseMatrix &U, vector[int] &P) nogil
-    void pivoted_LU_solve (const DenseMatrix &A, const DenseMatrix &b, DenseMatrix &x) nogil
+        DenseMatrix &B) nogil except +
+    void pivoted_LU (const DenseMatrix &A, DenseMatrix &L, DenseMatrix &U, vector[int] &P) nogil except +
+    void pivoted_LU_solve (const DenseMatrix &A, const DenseMatrix &b, DenseMatrix &x) nogil except +
     void inverse_GJ "SymEngine::inverse_gauss_jordan"(const DenseMatrix &A,
-        DenseMatrix &B) nogil
+        DenseMatrix &B) nogil except +
     void FFLU_solve "SymEngine::fraction_free_LU_solve"(const DenseMatrix &A,
-        const DenseMatrix &b, DenseMatrix &x) nogil
+        const DenseMatrix &b, DenseMatrix &x) nogil except +
     void FFGJ_solve "SymEngine::fraction_free_gauss_jordan_solve"(const DenseMatrix &A,
-        const DenseMatrix &b, DenseMatrix &x) nogil
+        const DenseMatrix &b, DenseMatrix &x) nogil except +
     void LDL_solve "SymEngine::LDL_solve"(const DenseMatrix &A, const DenseMatrix &b,
-        DenseMatrix &x) nogil
+        DenseMatrix &x) nogil except +
     void jacobian "SymEngine::sjacobian"(const DenseMatrix &A,
             const DenseMatrix &x, DenseMatrix &result) nogil except +
     void diff "SymEngine::sdiff"(const DenseMatrix &A,

symengine/lib/symengine_wrapper.pxd

@@ -4,159 +4,9 @@ from symengine cimport RCP, map_basic_basic
 cdef class Basic(object):
     cdef RCP[const symengine.Basic] thisptr
 
-cdef class Number(Basic):
-    pass
-
-cdef class Integer(Number):
-    pass
-
-cdef class Rational(Number):
-    pass
-
-cdef class Complex(Number):
-    pass
-
-cdef class RealDouble(Number):
-    pass
-
-cdef class ComplexDouble(Number):
-    pass
-
-cdef class RealMPFR(Number):
-    pass
-
-cdef class ComplexMPC(Number):
-    pass
-
-cdef class PyNumber(Number):
-    pass
-
-cdef class Add(Basic):
-    pass
-
-cdef class Mul(Basic):
-    pass
-
-cdef class Pow(Basic):
-    pass
-
-cdef class Function(Basic):
-    pass
-
-cdef class TrigFunction(Function):
-    pass
-
-cdef class Sin(TrigFunction):
-    pass
-
-cdef class Cos(TrigFunction):
-    pass
-
-cdef class Tan(TrigFunction):
-    pass
-
-cdef class Cot(TrigFunction):
-    pass
-
-cdef class Csc(TrigFunction):
-    pass
-
-cdef class Sec(TrigFunction):
-    pass
-
-cdef class ASin(TrigFunction):
-    pass
-
-cdef class ACos(TrigFunction):
-    pass
-
-cdef class ATan(TrigFunction):
-    pass
-
-cdef class ACot(TrigFunction):
-    pass
-
-cdef class ACsc(TrigFunction):
-    pass
-
-cdef class ASec(TrigFunction):
-    pass
-
-cdef class HyperbolicFunction(Function):
-    pass
-
-cdef class Sinh(HyperbolicFunction):
-    pass
-
-cdef class Cosh(HyperbolicFunction):
-    pass
-
-cdef class Tanh(HyperbolicFunction):
-    pass
-
-cdef class Coth(HyperbolicFunction):
-    pass
-
-cdef class Csch(HyperbolicFunction):
-    pass
-
-cdef class Sech(HyperbolicFunction):
-    pass
-
-cdef class ASinh(HyperbolicFunction):
-    pass
-
-cdef class ACosh(HyperbolicFunction):
-    pass
-
-cdef class ATanh(HyperbolicFunction):
-    pass
-
-cdef class ACoth(HyperbolicFunction):
-    pass
-
-cdef class ACsch(HyperbolicFunction):
-    pass
-
-cdef class ASech(HyperbolicFunction):
-    pass
-
-cdef class FunctionSymbol(Function):
-    pass
-
-cdef class PyFunction(FunctionSymbol):
-    pass
-
-cdef class Abs(Function):
-    pass
-
-cdef class Gamma(Function):
-    pass
-
-cdef class ATan2(Function):
-    pass
-
-cdef class Derivative(Basic):
-    pass
-
-cdef class Subs(Basic):
-    pass
-
 cdef class MatrixBase(object):
     cdef symengine.MatrixBase* thisptr
 
-cdef class DenseMatrixBase(MatrixBase):
-    pass
-
-cdef class MutableDenseMatrix(DenseMatrixBase):
-    pass
-
-cdef class ImmutableDenseMatrix(DenseMatrixBase):
-    pass
-
-cdef class Log(Function):
-    pass
-
 cdef class PyFunctionClass(object):
     cdef RCP[const symengine.PyFunctionClass] thisptr