-- -*- M2-comint -*- hash: 4914802418691755791

i1 : kk=ZZ/101;

i2 : setRandomSeed 0
 -- setting random seed to 0

o2 = 0

i3 : g=1;

i4 : rNP=randNicePencil(kk,g);

i5 : qq=rNP.quadraticForm;

i6 : S=rNP.qqRing;

i7 : P=kk[drop(gens S,-2)]

o7 = P

o7 : PolynomialRing

i8 : qs=sub(diff(matrix{{S_(2*g+2), S_(2*g+3)}}, qq), P)

o8 = | x_0y_0-z_1^2 -5x_0^2+12x_0z_1-18y_0z_1-30z_1^2-12x_0z_2+10y_0z_2+6z_1z
     ------------------------------------------------------------------------
     _2-48z_2^2 |

             1      2
o8 : Matrix P  <-- P

i9 : CI=P/ideal qs

o9 = CI

o9 : QuotientRing

i10 : cM=cliffordModule(rNP.matFact1,rNP.matFact2,rNP.baseRing)

o10 = CliffordModule{...6...}

o10 : CliffordModule

i11 : betti (F=cliffordModuleToCIResolution(cM,S,CI))

              0  1  2 3 4 5
o11 = total: 20 16 12 8 5 5
          0: 20 16 12 8 4 1
          1:  .  .  . . 1 4

o11 : BettiTally

i12 : cMu=cliffordModule(rNP.matFactu1,rNP.matFactu2,rNP.baseRing)

o12 = CliffordModule{...6...}

o12 : CliffordModule

i13 : f=cMu.hyperellipticBranchEquation

           3       2 2        3      4
o13 = - 12s t - 50s t  - 16s*t  + 47t

o13 : kk[s, t]

i14 : L=randomLineBundle(0,f);

i15 : betti (FL=cliffordModuleToCIResolution(tensorProduct(cM,L),S,CI))

              0 1 2 3 4  5
o15 = total: 12 8 4 4 8 12
         -1: 12 8 4 . .  .
          0:  . . . 4 8 12

o15 : BettiTally

i16 : betti (FuL=cliffordModuleToCIResolution(tensorProduct(cMu,L),S,CI))

             0 1 2 3 4 5
o16 = total: 7 5 3 2 3 5
         -1: 7 5 3 1 . .
          0: . . . 1 3 5

o16 : BettiTally

i17 : 
