diff options
Diffstat (limited to 'src/algebra/intalg.spad.pamphlet')
| -rw-r--r-- | src/algebra/intalg.spad.pamphlet | 10 |
1 files changed, 5 insertions, 5 deletions
diff --git a/src/algebra/intalg.spad.pamphlet b/src/algebra/intalg.spad.pamphlet index 4bf34804..22c18476 100644 --- a/src/algebra/intalg.spad.pamphlet +++ b/src/algebra/intalg.spad.pamphlet @@ -121,7 +121,7 @@ AlgebraicHermiteIntegration(F,UP,UPUP,R):Exports == Implementation where for k in minColIndex mat .. maxColIndex mat repeat (bc := extendedEuclidean(qelt(mat, j, k), modulus, qelt(vec, i))) case "failed" => return new(0, 0) - qsetelt_!(ans, i, bc.coef1) + qsetelt!(ans, i, bc.coef1) ans sol := particularSolution(map(#1::RF, mat)$MatrixCategoryFunctions2(UP, Vector UP, Vector UP, Matrix UP, RF, @@ -133,7 +133,7 @@ AlgebraicHermiteIntegration(F,UP,UPUP,R):Exports == Implementation where for i in minIndex ans .. maxIndex ans repeat (bc := extendedEuclidean(denom qelt(sol, i), modulus, 1)) case "failed" => return new(0, 0) - qsetelt_!(ans, i, (numer qelt(sol, i) * bc.coef1) rem modulus) + qsetelt!(ans, i, (numer qelt(sol, i) * bc.coef1) rem modulus) ans @ @@ -232,8 +232,8 @@ AlgebraicIntegrate(R0, F, UP, UPUP, R): Exports == Implementation where vf:Vector(QF) := new(n, 0) for i in minIndex v .. maxIndex v repeat r := separate(qelt(v, i) / d)$GP - qsetelt_!(vf, i, r.fracPart) - qsetelt_!(vp, i, r.polyPart) + qsetelt!(vf, i, r.fracPart) + qsetelt!(vp, i, r.polyPart) ff := represents(vf, w := integralBasis()) h := HermiteIntegrate(ff, derivation) p := represents(map(convert(#1)@QF, vp)$VectorFunctions2(GP, QF), w) @@ -262,7 +262,7 @@ AlgebraicIntegrate(R0, F, UP, UPUP, R): Exports == Implementation where -- where the ri's are rational numbers, and fc(z) is arbitrary -- (fc can be linear too) -- la = [b1....,bk] (all rational residues) - la := [- coefficient(q.factor, 0) for q in remove_!(fc::FAC, lf)] + la := [- coefficient(q.factor, 0) for q in remove!(fc::FAC, lf)] -- ld = [D1,...,Dk] where Di is the sum of places where f has residue bi ld := [divisor(rec.num, rec.den, rec.derivden, rec.gd, b::F) for b in la] pp := UPQ2F(fc.factor) |