gh-39273: Fix remaining all imports
    
<!-- ^ Please provide a concise and informative title. -->
<!-- ^ Don't put issue numbers in the title, do this in the PR
description below. -->
<!-- ^ For example, instead of "Fixes #12345" use "Introduce new method
to calculate 1 + 2". -->
<!-- v Describe your changes below in detail. -->
<!-- v Why is this change required? What problem does it solve? -->
<!-- v If this PR resolves an open issue, please link to it here. For
example, "Fixes #12345". -->

This should take care of the remaining imports from `all`. (There are
still a couple of them, but they are okay to leave I think.)


### 📝 Checklist

<!-- Put an `x` in all the boxes that apply. -->

- [ ] The title is concise and informative.
- [ ] The description explains in detail what this PR is about.
- [ ] I have linked a relevant issue or discussion.
- [ ] I have created tests covering the changes.
- [ ] I have updated the documentation and checked the documentation
preview.

### ⌛ Dependencies

<!-- List all open PRs that this PR logically depends on. For example,
-->
<!-- - #12345: short description why this is a dependency -->
<!-- - #34567: ... -->
    
URL: #39273
Reported by: Tobias Diez
Reviewer(s): Dima Pasechnik, Tobias Diez

Author: Release Manager

src/sage/arith/misc.py

@@ -1705,7 +1705,7 @@ def plot(self, xmin=1, xmax=50, k=1, pointsize=30, rgbcolor=(0, 0, 1), join=True
             124.0
         """
         v = [(n, sigma(n, k)) for n in range(xmin, xmax + 1)]
-        from sage.plot.all import list_plot
+        from sage.plot.plot import list_plot
         P = list_plot(v, pointsize=pointsize, rgbcolor=rgbcolor, **kwds)
         if join:
             P += list_plot(v, plotjoined=True, rgbcolor=(0.7, 0.7, 0.7), **kwds)
@@ -3193,7 +3193,7 @@ def plot(self, xmin=1, xmax=50, pointsize=30, rgbcolor=(0, 0, 1),
             46.0
         """
         v = [(n, euler_phi(n)) for n in range(xmin, xmax + 1)]
-        from sage.plot.all import list_plot
+        from sage.plot.plot import list_plot
         P = list_plot(v, pointsize=pointsize, rgbcolor=rgbcolor, **kwds)
         if join:
             P += list_plot(v, plotjoined=True, rgbcolor=(0.7, 0.7, 0.7), **kwds)
@@ -4743,7 +4743,7 @@ def plot(self, xmin=0, xmax=50, pointsize=30, rgbcolor=(0, 0, 1), join=True,
         """
         values = self.range(xmin, xmax + 1)
         v = [(n, values[n - xmin]) for n in range(xmin, xmax + 1)]
-        from sage.plot.all import list_plot
+        from sage.plot.plot import list_plot
         P = list_plot(v, pointsize=pointsize, rgbcolor=rgbcolor, **kwds)
         if join:
             P += list_plot(v, plotjoined=True, rgbcolor=(0.7, 0.7, 0.7), **kwds)
@@ -5271,7 +5271,7 @@ def falling_factorial(x, a):
         (isinstance(a, Expression) and
          a.is_integer())) and a >= 0:
         return prod(((x - i) for i in range(a)), z=x.parent().one())
-    from sage.functions.all import gamma
+    from sage.functions.gamma import gamma
     return gamma(x + 1) / gamma(x - a + 1)
 
 
@@ -5363,7 +5363,7 @@ def rising_factorial(x, a):
         (isinstance(a, Expression) and
          a.is_integer())) and a >= 0:
         return prod(((x + i) for i in range(a)), z=x.parent().one())
-    from sage.functions.all import gamma
+    from sage.functions.gamma import gamma
     return gamma(x + a) / gamma(x)
 
 

src/sage/calculus/desolvers.py

@@ -1365,7 +1365,7 @@ def desolve_rk4_inner(de, dvar):
 
     if not (isinstance(dvar, Expression) and dvar.is_symbol()):
         from sage.symbolic.ring import SR
-        from sage.calculus.all import diff
+        from sage.calculus.functional import diff
         from sage.symbolic.relation import solve
         if isinstance(de, Expression) and de.is_relational():
             de = de.lhs() - de.rhs()

src/sage/calculus/riemann.pyx

@@ -744,7 +744,7 @@ cdef class Riemann_Map:
             sage: m.plot_boundaries(plotjoined=False, rgbcolor=[0,0,1], thickness=6)    # needs sage.plot
             Graphics object consisting of 1 graphics primitive
         """
-        from sage.plot.all import list_plot
+        from sage.plot.plot import list_plot
 
         if rgbcolor is None:
             rgbcolor = [0, 0, 0]
@@ -942,7 +942,8 @@ cdef class Riemann_Map:
             ....:                      thickness=2.0, min_mag=0.1)
         """
         from sage.plot.complex_plot import ComplexPlot
-        from sage.plot.all import list_plot, Graphics
+        from sage.plot.plot import list_plot
+        from sage.plot.graphics import Graphics
 
         cdef int k, i
         if self.exterior:
@@ -1068,7 +1069,7 @@ cdef class Riemann_Map:
             Graphics object consisting of 1 graphics primitive
         """
         from sage.plot.complex_plot import ComplexPlot
-        from sage.plot.all import Graphics
+        from sage.plot.graphics import Graphics
 
         if plot_range is None:
             plot_range = []

src/sage/categories/examples/finite_coxeter_groups.py

@@ -129,7 +129,7 @@ def __contains__(self, x):
 
         (also tested by :meth:`test_an_element` :meth:`test_some_elements`)
         """
-        from sage.structure.all import parent
+        from sage.structure.element import parent
         return parent(x) is self
 
     @cached_method

src/sage/categories/number_fields.py

@@ -185,7 +185,7 @@ def zeta_function(self, prec=53,
                 PARI zeta function associated to Rational Field
             """
             if algorithm == 'gp':
-                from sage.lfunctions.all import Dokchitser
+                from sage.lfunctions.dokchitser import Dokchitser
                 r1, r2 = self.signature()
                 zero = [0]
                 one = [1]

src/sage/categories/sets_cat.py

@@ -36,9 +36,9 @@
 from sage.misc.lazy_attribute import lazy_attribute
 from sage.misc.lazy_import import lazy_import, LazyImport
 from sage.misc.lazy_format import LazyFormat
+# Do not use sage.categories.all here to avoid initialization loop
 from sage.categories.category import Category
 from sage.categories.category_singleton import Category_singleton
-# Do not use sage.categories.all here to avoid initialization loop
 from sage.categories.sets_with_partial_maps import SetsWithPartialMaps
 from sage.categories.subquotients import SubquotientsCategory
 from sage.categories.quotients import QuotientsCategory

src/sage/coding/goppa_code.py

@@ -33,7 +33,7 @@
 from sage.coding.linear_code import AbstractLinearCode
 from sage.coding.encoder import Encoder
 from sage.modules.free_module_element import vector
-from sage.coding.all import codes
+import sage.coding.codes_catalog as codes
 
 
 def _columnize(element):

src/sage/combinat/combinat.py

@@ -122,21 +122,22 @@
 #                  https://www.gnu.org/licenses/
 # ****************************************************************************
 from __future__ import annotations
+
 from collections.abc import Iterator
 
 from sage.arith.misc import bernoulli, factorial
-from sage.rings.integer_ring import ZZ
-from sage.rings.rational_field import QQ
+from sage.combinat.combinat_cython import _stirling_number2
+from sage.misc.cachefunc import cached_function
+from sage.misc.inherit_comparison import InheritComparisonClasscallMetaclass
+from sage.misc.lazy_import import lazy_import
+from sage.misc.misc_c import prod
 from sage.rings.integer import Integer
+from sage.rings.integer_ring import ZZ
 from sage.rings.polynomial.polynomial_element import Polynomial
 from sage.rings.polynomial.polynomial_ring_constructor import PolynomialRing
-from sage.misc.misc_c import prod
-from sage.misc.cachefunc import cached_function
-from sage.structure.sage_object import SageObject
-from sage.misc.lazy_import import lazy_import
-from .combinat_cython import _stirling_number2
-from sage.misc.inherit_comparison import InheritComparisonClasscallMetaclass
+from sage.rings.rational_field import QQ
 from sage.structure.element import Element
+from sage.structure.sage_object import SageObject
 
 lazy_import('sage.interfaces.maxima_lib', 'maxima')
 lazy_import('sage.libs.pari', 'pari')
@@ -339,7 +340,7 @@ def bell_number(n, algorithm='flint', **options) -> Integer:
     if n < 0:
         raise ArithmeticError('Bell numbers not defined for negative indices')
     if algorithm == 'mpmath':
-        from sage.libs.mpmath.all import bell, mp, mag
+        from mpmath import bell, mag, mp
         old_prec = mp.dps
         if 'prec' in options:
             mp.dps = options['prec']
@@ -2052,8 +2053,8 @@ def bell_polynomial(n: Integer, k=None, ordinary=False):
     - [Bel1927]_
     - [Com1974]_
     """
-    from sage.combinat.partition import Partitions
     from sage.arith.misc import multinomial
+    from sage.combinat.partition import Partitions
     if k is None:
         partitions = Partitions(n)
         # We set k = 1 to use the correct ring

src/sage/combinat/designs/block_design.py

@@ -148,7 +148,7 @@ def are_hyperplanes_in_projective_geometry_parameters(v, k, lmbda, return_parame
         ....:         assert are_hyperplanes_in_projective_geometry_parameters(v,k,l+1) is False
         ....:         assert are_hyperplanes_in_projective_geometry_parameters(v,k,l-1) is False
     """
-    import sage.arith.all as arith
+    from sage.arith.misc import gcd
 
     q1 = Integer(v - k)
     q2 = Integer(k - lmbda)
@@ -161,7 +161,7 @@ def are_hyperplanes_in_projective_geometry_parameters(v, k, lmbda, return_parame
     p1,e1 = q1.factor()[0]
     p2,e2 = q2.factor()[0]
 
-    k = arith.gcd(e1,e2)
+    k = gcd(e1,e2)
     d = e1//k
     q = p1**k
     if e2//k != d-1 or lmbda != (q**(d-1)-1)//(q-1):

src/sage/combinat/designs/difference_family.py

@@ -47,11 +47,10 @@
 #                  https://www.gnu.org/licenses/
 # ****************************************************************************
 
-from sage.arith.misc import is_prime_power
+from sage.arith.misc import is_prime_power, is_square, factor
 from sage.misc.cachefunc import cached_function
 
 from sage.categories.sets_cat import EmptySetError
-import sage.arith.all as arith
 from sage.misc.unknown import Unknown
 from sage.rings.finite_rings.integer_mod_ring import Zmod
 from sage.rings.integer import Integer
@@ -572,13 +571,13 @@ def radical_difference_set(K, k, l=1, existence=False, check=True):
         add_zero = True
 
     # q = 4t^2 + 1, t odd
-    elif v % 8 == 5 and k == (v-1)//4 and arith.is_square((v-1)//4):
+    elif v % 8 == 5 and k == (v-1)//4 and is_square((v-1)//4):
         if existence:
             return True
         add_zero = False
 
     # q = 4t^2 + 9, t odd
-    elif v % 8 == 5 and k == (v+3)//4 and arith.is_square((v-9)//4):
+    elif v % 8 == 5 and k == (v+3)//4 and is_square((v-9)//4):
         if existence:
             return True
         add_zero = True
@@ -3823,7 +3822,7 @@ def difference_family(v, k, l=1, existence=False, explain_construction=False, ch
         G = Zmod(v)
         return G, [list(range(1, v))]
 
-    factorization = arith.factor(v)
+    factorization = factor(v)
     if len(factorization) == 1:
         from sage.rings.finite_rings.finite_field_constructor import GF
         K = GF(v,'z')