1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
|
needsPackage "NumericalSchubertCalculus"
debug needsPackage "NumericalSchubertCalculus"
setRandomSeed 21
setDebugOptions("debug"=>1,"verify solutions"=>false)
-- Problem (2,1)^5 in Gr(3,8)
-- a problem with 32 solutions
SchPblm = {
({2,1}, id_(FFF^8)),
({2,1}, random(FFF^8,FFF^8)),
({2,1}, random(FFF^8,FFF^8)),
({2,1}, random(FFF^8,FFF^8)),
({2,1}, random(FFF^8,FFF^8))
}
S = solveSchubertProblem(SchPblm, 3,8);
-- time assert all(S,s->checkIncidenceSolution(s, SchPblm))
print(#S)
end
restart
time load "21e5-G38.m2"
|
