Repository: sympy/sympy

Issue:
NotImplementedError in limit 
`limit(sin ( sqrt (x))/ sqrt ( sin (x)),x,0,'+')`

Discussion:

Patch:
diff --git a/sympy/core/power.py b/sympy/core/power.py
index a30d0ce0b5..c878fd5ec2 100644
--- a/sympy/core/power.py
+++ b/sympy/core/power.py
@@ -1564,7 +1564,7 @@ def e2int(e):
 
             nuse = n - ei
 
-            if e.is_extended_real and e.is_positive:
+            if e.is_real:
                 lt = b.as_leading_term(x)
 
                 # Try to correct nuse (= m) guess from:
diff --git a/sympy/tensor/array/dense_ndim_array.py b/sympy/tensor/array/dense_ndim_array.py
index 874a029933..c962bb870c 100644
--- a/sympy/tensor/array/dense_ndim_array.py
+++ b/sympy/tensor/array/dense_ndim_array.py
@@ -5,7 +5,7 @@
 from sympy.core.sympify import _sympify
 from sympy.tensor.array.mutable_ndim_array import MutableNDimArray
 from sympy.tensor.array.ndim_array import NDimArray, ImmutableNDimArray
-from sympy.simplify import simplify
+
 
 
 class DenseNDimArray(NDimArray):
@@ -150,8 +150,6 @@ def __setitem__(self, index, value):
     def as_mutable(self):
         return MutableDenseNDimArray(self)
 
-    def _eval_simplify(self, **kwargs):
-        return self.applyfunc(simplify)
 
 class MutableDenseNDimArray(DenseNDimArray, MutableNDimArray):
 


Test patch:
diff --git a/sympy/series/tests/test_limits.py b/sympy/series/tests/test_limits.py
index 73fdefbeee..93eb22b3d8 100644
--- a/sympy/series/tests/test_limits.py
+++ b/sympy/series/tests/test_limits.py
@@ -605,3 +605,6 @@ def test_issue_17431():
 
 def test_issue_17671():
     assert limit(Ei(-log(x)) - log(log(x))/x, x, 1) == EulerGamma
+
+def test_issue_18306():
+    assert limit(sin(sqrt(x))/sqrt(sin(x)), x, 0, '+') == 1
diff --git a/sympy/tensor/array/tests/test_ndim_array.py b/sympy/tensor/array/tests/test_ndim_array.py
index 9a097a3fb8..261f2f87d9 100644
--- a/sympy/tensor/array/tests/test_ndim_array.py
+++ b/sympy/tensor/array/tests/test_ndim_array.py
@@ -1,10 +1,8 @@
 from sympy.testing.pytest import raises
 from sympy import (
     Array, ImmutableDenseNDimArray, ImmutableSparseNDimArray,
-    MutableDenseNDimArray, MutableSparseNDimArray, sin, cos,
-    simplify
+    MutableDenseNDimArray, MutableSparseNDimArray
 )
-from sympy.abc import x, y
 
 
 array_types = [
@@ -33,12 +31,3 @@ def test_array_negative_indices():
         raises(ValueError, lambda: test_array[-3, :])
 
         assert test_array[-1, -1] == 10
-
-
-def test_issue_18361():
-    A = Array([sin(2 * x) - 2 * sin(x) * cos(x)])
-    B = Array([sin(x)**2 + cos(x)**2, 0])
-    C = Array([(x + x**2)/(x*sin(y)**2 + x*cos(y)**2), 2*sin(x)*cos(x)])
-    assert simplify(A) == Array([0])
-    assert simplify(B) == Array([1, 0])
-    assert simplify(C) == Array([x + 1, sin(2*x)])