rafarob32's profile - activity

BUG: SageMath 9.6 for Ubuntu; inverse() with BUG: SageMath 9.6 for Ubuntu; inverse() A = matrix(IntegerModRing(6), [[2,

Sorry, I don't know how to use mpq_t attribute. Should I import some library for it? How do it? Thanks in advance.

I re-edited with an example. Thanks.

I want to work efficiently with rational numbers under the context %cython but I don`t know how. Can anyone suggest any ideas? Thank you.

An example:

  %cython

    def rational_partitions(n):
        sol = [i/n for i in range(n)]
        for a in sol[0:-1]:
            for b in sol[1:]:
                k=2
                while abs(b-a)/k>1/n:
                    sol.append(abs(b-a)/k)
                    k += 1
        return sol

rational_partitions(10)

There is an inconsistency in the use of the LU decomposition of a matrix. If matrix A is declared in the RDF ring then P, L, U = A.LU () provides PA = LU but if A is declared in the RR ring (or QQ, ZZ,...) then A = PLU. Therefore in a case P is P.transpose () on the other.