sage.plot: Update # needs

Author: Matthias Koeppe

src/sage/plot/graphics.py

@@ -64,7 +64,7 @@ def is_Graphics(x):
         sage: from sage.plot.graphics import is_Graphics
         sage: is_Graphics(1)
         False
-        sage: is_Graphics(disk((0.0, 0.0), 1, (0, pi/2)))
+        sage: is_Graphics(disk((0.0, 0.0), 1, (0, pi/2)))                               # needs sage.symbolic
         True
     """
     return isinstance(x, Graphics)
@@ -1359,6 +1359,7 @@ def _set_scale(self, subplot, scale=None, base=None):
 
         EXAMPLES::
 
+            sage: # needs sage.symbolic
             sage: p = plot(x, 1, 10)
             sage: fig = p.matplotlib()
             sage: ax = fig.get_axes()[0]
@@ -1373,6 +1374,7 @@ def _set_scale(self, subplot, scale=None, base=None):
 
         TESTS::
 
+            sage: # needs sage.symbolic
             sage: p._set_scale(ax, 'log')
             Traceback (most recent call last):
             ...

src/sage/plot/histogram.py

@@ -274,8 +274,8 @@ def histogram(datalist, **options):
 
         sage: nv = normalvariate
         sage: H = histogram([nv(0, 1) for _ in range(1000)], bins=20, density=True, range=[-5, 5])
-        sage: P = plot(1/sqrt(2*pi)*e^(-x^2/2), (x, -5, 5), color='red', linestyle='--')
-        sage: H+P
+        sage: P = plot(1/sqrt(2*pi)*e^(-x^2/2), (x, -5, 5), color='red', linestyle='--')    # needs sage.symbolic
+        sage: H + P                                                                         # needs sage.symbolic
         Graphics object consisting of 2 graphics primitives
 
     .. PLOT::

src/sage/plot/line.py

@@ -590,6 +590,7 @@ def line2d(points, **options):
 
     A red, blue, and green "cool cat"::
 
+        sage: # needs sage.symbolic
         sage: G = plot(-cos(x), -2, 2, thickness=5, rgbcolor=(0.5,1,0.5))
         sage: P = polygon([[1,2], [5,6], [5,0]], rgbcolor=(1,0,0))
         sage: Q = polygon([(-x,y) for x,y in P[0]], rgbcolor=(0,0,1))

src/sage/plot/plot3d/base.pyx

@@ -1998,9 +1998,10 @@ end_scene""".format(
 
         This works when faces have more then 3 sides::
 
-            sage: P = polytopes.dodecahedron()                                          # needs sage.geometry.polyhedron
-            sage: Q = P.plot().all[-1]                                                  # needs sage.geometry.polyhedron
-            sage: print(Q.stl_binary()[:40].decode('ascii'))                            # needs sage.geometry.polyhedron
+            sage: # needs sage.geometry.polyhedron sage.groups
+            sage: P = polytopes.dodecahedron()
+            sage: Q = P.plot().all[-1]
+            sage: print(Q.stl_binary()[:40].decode('ascii'))
             STL binary file / made by SageMath / ###
         """
         import struct
@@ -2060,9 +2061,10 @@ end_scene""".format(
 
         Now works when faces have more then 3 sides::
 
-            sage: P = polytopes.dodecahedron()                                          # needs sage.geometry.polyhedron
-            sage: Q = P.plot().all[-1]                                                  # needs sage.geometry.polyhedron
-            sage: print(Q.stl_ascii_string().splitlines()[:7])                          # needs sage.geometry.polyhedron
+            sage: # needs sage.geometry.polyhedron sage.groups
+            sage: P = polytopes.dodecahedron()
+            sage: Q = P.plot().all[-1]
+            sage: print(Q.stl_ascii_string().splitlines()[:7])
             ['solid surface',
              'facet normal 0.0 0.5257311121191338 0.8506508083520399',
              '    outer loop',

src/sage/plot/plot3d/revolution_plot3d.py

@@ -1,4 +1,4 @@
-# sage.doctest: needs sage.plot
+# sage.doctest: needs sage.plot sage.symbolic
 """
 Surfaces of revolution
 

src/sage/plot/plot3d/transform.pyx

@@ -174,11 +174,11 @@ def rotate_arbitrary(v, double theta):
 
         sage: rotate_arbitrary((1,2,3), -1).det()
         1.0000000000000002
-        sage: rotate_arbitrary((1,1,1), 2*pi/3) * vector(RDF, (1,2,3))  # rel tol 2e-15
+        sage: rotate_arbitrary((1,1,1), 2*pi/3) * vector(RDF, (1,2,3))  # rel tol 2e-15     # needs sage.symbolic
         (1.9999999999999996, 2.9999999999999996, 0.9999999999999999)
         sage: rotate_arbitrary((1,2,3), 5) * vector(RDF, (1,2,3))  # rel tol 2e-15
         (1.0000000000000002, 2.0, 3.000000000000001)
-        sage: rotate_arbitrary((1,1,1), pi/7)^7  # rel tol 2e-15
+        sage: rotate_arbitrary((1,1,1), pi/7)^7  # rel tol 2e-15                            # needs sage.symbolic
         [-0.33333333333333337   0.6666666666666671   0.6666666666666665]
         [  0.6666666666666665 -0.33333333333333337   0.6666666666666671]
         [  0.6666666666666671   0.6666666666666667 -0.33333333333333326]
@@ -194,7 +194,7 @@ def rotate_arbitrary(v, double theta):
 
         Setup some variables::
 
-            sage: vx,vy,vz,theta = var('x y z theta')
+            sage: vx,vy,vz,theta = var('x y z theta')                                       # needs sage.symbolic
 
         Symbolic rotation matrices about X and Y axis::
 
@@ -205,37 +205,37 @@ def rotate_arbitrary(v, double theta):
         way to tell Maxima that `x^2+y^2+z^2=1` which would make for
         a much cleaner calculation::
 
-            sage: vy = sqrt(1-vx^2-vz^2)
+            sage: vy = sqrt(1-vx^2-vz^2)                                                    # needs sage.symbolic
 
         Now we rotate about the `x`-axis so `v` is in the `xy`-plane::
 
-            sage: t = arctan(vy/vz)+pi/2
-            sage: m = rotX(t)
-            sage: new_y = vy*cos(t) - vz*sin(t)
+            sage: t = arctan(vy/vz)+pi/2                                                    # needs sage.symbolic
+            sage: m = rotX(t)                                                               # needs sage.symbolic
+            sage: new_y = vy*cos(t) - vz*sin(t)                                             # needs sage.symbolic
 
         And rotate about the `z` axis so `v` lies on the `x` axis::
 
-            sage: s = arctan(vx/new_y) + pi/2
-            sage: m = rotZ(s) * m
+            sage: s = arctan(vx/new_y) + pi/2                                               # needs sage.symbolic
+            sage: m = rotZ(s) * m                                                           # needs sage.symbolic
 
         Rotating about `v` in our old system is the same as rotating
         about the `x`-axis in the new::
 
-            sage: m = rotX(theta) * m
+            sage: m = rotX(theta) * m                                                       # needs sage.symbolic
 
         Do some simplifying here to avoid blow-up::
 
-            sage: m = m.simplify_rational()
+            sage: m = m.simplify_rational()                                                 # needs sage.symbolic
 
         Now go back to the original coordinate system::
 
-            sage: m = rotZ(-s) * m
-            sage: m = rotX(-t) * m
+            sage: m = rotZ(-s) * m                                                          # needs sage.symbolic
+            sage: m = rotX(-t) * m                                                          # needs sage.symbolic
 
         And simplify every single entry (which is more effective that simplify
         the whole matrix like above)::
 
-            sage: m.parent()([x.simplify_full() for x in m._list()])  # long time; random
+            sage: m.parent()([x.simplify_full() for x in m._list()])  # random  # long time, needs sage.symbolic
             [                                       -(cos(theta) - 1)*x^2 + cos(theta)              -(cos(theta) - 1)*sqrt(-x^2 - z^2 + 1)*x + sin(theta)*abs(z)      -((cos(theta) - 1)*x*z^2 + sqrt(-x^2 - z^2 + 1)*sin(theta)*abs(z))/z]
             [             -(cos(theta) - 1)*sqrt(-x^2 - z^2 + 1)*x - sin(theta)*abs(z)                           (cos(theta) - 1)*x^2 + (cos(theta) - 1)*z^2 + 1 -((cos(theta) - 1)*sqrt(-x^2 - z^2 + 1)*z*abs(z) - x*z*sin(theta))/abs(z)]
             [     -((cos(theta) - 1)*x*z^2 - sqrt(-x^2 - z^2 + 1)*sin(theta)*abs(z))/z -((cos(theta) - 1)*sqrt(-x^2 - z^2 + 1)*z*abs(z) + x*z*sin(theta))/abs(z)                                        -(cos(theta) - 1)*z^2 + cos(theta)]