Bumping version

spectralDNS · Jan 25, 2019 · 9703302 · 9703302
1 parent f1dbe7a
commit 9703302
Show file tree

Hide file tree

Showing 2 changed files with 18 additions and 1 deletion.
diff --git a/shenfun/__init__.py b/shenfun/__init__.py
@@ -25,7 +25,7 @@
 """
 #pylint: disable=wildcard-import,no-name-in-module
 
-__version__ = '1.2.2'
+__version__ = '1.2.3'
 __author__ = 'Mikael Mortensen'
 
 import numpy as np

diff --git a/shenfun/chebyshev/bases.py b/shenfun/chebyshev/bases.py
@@ -3,6 +3,7 @@
 """
 import functools
 import numpy as np
+import sympy
 from numpy.polynomial import chebyshev as n_cheb
 from scipy.special import eval_chebyt
 from mpi4py_fft import fftw
@@ -88,6 +89,10 @@ def points_and_weights(self, N=None, map_true_domain=False):
     def vandermonde(self, x):
         return n_cheb.chebvander(x, self.N-1)
 
+    def sympy_basis(self, i=0):
+        x = sympy.symbols('x')
+        return sympy.chebyshevt(i, x)
+
     def evaluate_basis(self, x, i=0, output_array=None):
         x = np.atleast_1d(x)
         if output_array is None:
@@ -358,6 +363,10 @@ def _composite_basis(self, V):
         P[:, -1] = (V[:, 0] - V[:, 1])/2
         return P
 
+    def sympy_basis(self, i=0):
+        x = sympy.symbols('x')
+        return sympy.chebyshevt(i, x) - sympy.chebyshevt(i+2, x)
+
     def evaluate_basis(self, x, i=0, output_array=None):
         x = np.atleast_1d(x)
         if output_array is None:
@@ -479,6 +488,10 @@ def _composite_basis(self, V):
         P[:, :-2] = V[:, :-2] - (k[:-2]/(k[:-2]+2))**2*V[:, 2:]
         return P
 
+    def sympy_basis(self, i=0):
+        x = sympy.symbols('x')
+        return sympy.chebyshevt(i, x) - (i/(i+2))**2*sympy.chebyshevt(i+2, x)
+
     def evaluate_basis(self, x, i=0, output_array=None):
         x = np.atleast_1d(x)
         if output_array is None:
@@ -608,6 +621,10 @@ def _composite_basis(self, V):
         P[:, :-4] = V[:, :-4] - (2*(k+2)/(k+3))*V[:, 2:-2] + ((k+1)/(k+3))*V[:, 4:]
         return P
 
+    def sympy_basis(self, i=0):
+        x = sympy.symbols('x')
+        return sympy.chebyshevt(i, x) - (2*(i+2)/(i+3))*sympy.chebyshevt(i+2, x) + (i+1)/(i+3)*sympy.chebyshevt(i+4, x)
+
     def evaluate_basis(self, x, i=0, output_array=None):
         x = np.atleast_1d(x)
         if output_array is None: