diff options
Diffstat (limited to 'src/algebra/pleqn.spad.pamphlet')
| -rw-r--r-- | src/algebra/pleqn.spad.pamphlet | 13 |
1 files changed, 3 insertions, 10 deletions
diff --git a/src/algebra/pleqn.spad.pamphlet b/src/algebra/pleqn.spad.pamphlet index e3090a74..0e371880 100644 --- a/src/algebra/pleqn.spad.pamphlet +++ b/src/algebra/pleqn.spad.pamphlet @@ -315,7 +315,6 @@ ParametricLinearEquations(R,Var,Expon,GR): false B1solve (sys:Linsys):Linsoln == - i,j,i1,j1:I rss:L I:=sys.rows nss:L I:=sys.cols k:=sys.rank @@ -349,7 +348,6 @@ ParametricLinearEquations(R,Var,Expon,GR): [p,pbas] regime (y, coef, w, psbf, rk, rkmax, mode) == - i,j:I -- use the y.det nonzero to simplify the groebner basis -- of ideal generated by higher order subdeterminants ydetf:L GR:=factorset y.det @@ -396,8 +394,6 @@ ParametricLinearEquations(R,Var,Expon,GR): if filemode then newfile:=new$FNAME ("",outname,"regime") rksoln:=open$(File Rec3) newfile - y:Rec - k:NNI rrcl:RankConds:= entry? (mode,[1,2,3,7,8,9]$(L I)) => ParCondList (coeff,0) entry? (mode,[4,5,6,10,11,12]$(L I)) => ParCondList (coeff,h) @@ -409,7 +405,7 @@ ParametricLinearEquations(R,Var,Expon,GR): psb:Fgb:= (if rk=rkmax then [] else rrcl.(k+1).fgb) psbf:L L GR:= [factorset x for x in psb] psbf:= minset(psbf) - for y in pc repeat + for y: Rec in pc repeat rec3:Rec3:= regime (y, coeff, w, psbf, rk, rkmax, mode) inconsistent? rec3.wcond => "incompatible system" if filemode then write!(rksoln, rec3) @@ -428,13 +424,11 @@ ParametricLinearEquations(R,Var,Expon,GR): pc:Eqns:=[] npc: Eqns:=[] psbf: Fgb:=[] - rc: Rec done: Boolean := false r:=nrows mat n:=ncols mat maxrk:I:=min(r,n) - k:NNI - for k in min(r,n)..h by -1 until done repeat + for k: NNI in min(r,n)..h by -1 until done repeat pc:= ParCond(mat,k) npc:=[] (empty? pc) and (k >= 1) => maxrk:= k - 1 @@ -443,7 +437,7 @@ ParametricLinearEquations(R,Var,Expon,GR): else zro:L GR:= (if k = maxrk then [] else rcl.1.fgb) covered:Boolean:=false - for rc in pc until covered repeat + for rc: Rec in pc until covered repeat p:GR:= redPol$rp (rc.det, zro) p = 0 => "incompatible or covered subdeterminant" test:=hasoln(zro, [rc.det]) @@ -630,7 +624,6 @@ ParametricLinearEquations(R,Var,Expon,GR): redmat (mat,psb) == - i,j:I r:=nrows(mat) n:=ncols(mat) newmat: M GR:=zero(r,n) |