Ring of quasimodular forms as a commutative ring

edit

I'm trying to implement a polynomial ring of two variables, where the coefficient ring is the ring of quasimodular forms (code). I did

QM = QuasiModularForms(1)
PQM.<v1,v2> = QM['v1,v2']

and I got the following error:

TypeError                                 Traceback (most recent call last)
Cell In [9], line 2
      1 QM = QuasiModularForms(Integer(1))
----> 2 PQM = QM['v1,v2']; (v1, v2,) = PQM._first_ngens(2)

File /private/var/tmp/sage-9.8-current/local/var/lib/sage/venv-python3.11.1/lib/python3.11/site-packages/sage/structure/parent.pyx:1274, in sage.structure.parent.Parent.__getitem__ (build/cythonized/sage/structure/parent.c:11540)()
   1272         except AttributeError:
   1273             return self.list()[n]
-> 1274     return meth(n)
   1275 
   1276 #########################################################################

File /private/var/tmp/sage-9.8-current/local/var/lib/sage/venv-python3.11.1/lib/python3.11/site-packages/sage/categories/rings.py:1220, in Rings.ParentMethods.__getitem__(self, arg)
   1217 # 2. Otherwise, try to return a polynomial ring
   1219 from sage.rings.polynomial.polynomial_ring_constructor import PolynomialRing
-> 1220 return PolynomialRing(self, elts)

File /private/var/tmp/sage-9.8-current/local/var/lib/sage/venv-python3.11.1/lib/python3.11/site-packages/sage/rings/polynomial/polynomial_ring_constructor.py:693, in PolynomialRing(base_ring, *args, **kwds)
    690         raise TypeError("variable names specified twice inconsistently: %r and %r" % (names, kwnames))
    692 if multivariate or len(names) != 1:
--> 693     return _multi_variate(base_ring, names, **kwds)
    694 else:
    695     return _single_variate(base_ring, names, **kwds)

File /private/var/tmp/sage-9.8-current/local/var/lib/sage/venv-python3.11.1/lib/python3.11/site-packages/sage/rings/polynomial/polynomial_ring_constructor.py:844, in _multi_variate(base_ring, names, sparse, order, implementation)
    842     else:
    843         constructor = multi_polynomial_ring.MPolynomialRing_polydict
--> 844     R = constructor(base_ring, n, names, order)
    846 if R is None:
    847     raise ValueError("unknown implementation %r for multivariate polynomial rings" % (implementation,))

File /private/var/tmp/sage-9.8-current/local/var/lib/sage/venv-python3.11.1/lib/python3.11/site-packages/sage/rings/polynomial/multi_polynomial_ring.py:130, in MPolynomialRing_polydict.__init__(self, base_ring, n, names, order)
    128 self._gens = tuple(self._gens)
    129 self._zero_tuple = tuple(v)
--> 130 MPolynomialRing_base.__init__(self, base_ring, n, names, order)
    131 self._has_singular = can_convert_to_singular(self)

File /private/var/tmp/sage-9.8-current/local/var/lib/sage/venv-python3.11.1/lib/python3.11/site-packages/sage/rings/polynomial/multi_polynomial_ring_base.pyx:84, in sage.rings.polynomial.multi_polynomial_ring_base.MPolynomialRing_base.__init__ (build/cythonized/sage/rings/polynomial/multi_polynomial_ring_base.c:4101)()
     82 """
     83 if base_ring not in _CommutativeRings:
---> 84     raise TypeError("The base ring %s is not a commutative ring" % base_ring)
     85 
     86 n = int(n)

TypeError: The base ring Ring of Quasimodular Forms for Modular Group SL(2,Z) over Rational Field is not a commutative ring

So first of all, the ring of quasimodular forms is a commutative ring - the default coefficient is over the rationals. I tried to look at the details of the implementation, and here's what I found. The ring QuasiModularForms is essentially a polynomial ring of one variable (E2) with coefficient ModularForms (code), and the ring ModularForms is also a polynomial ring with certain generators. Hence if the base field is just rationals QQ, then everything should be a commutative ring (which is true mathematically). I wonder if there's a way to make it as a commutative ring in Sage. Thanks in advance.