File: 2.4.2

package info (click to toggle)
gap-hap 1.66%2Bds-1
  • links: PTS
  • area: main
  • in suites: trixie
  • size: 55,348 kB
  • sloc: xml: 15,368; sh: 216; javascript: 155; makefile: 126; ansic: 57; perl: 36
file content (15 lines) | stat: -rw-r--r-- 377 bytes parent folder | download | duplicates (3) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
 
#This example takes quite some time!
n:=50;; Betti2:=[];; Betti3:=[];;
for p in [10..133]*(1/200) do
 Print("p= ",p,"\n");
 Y:=RandomSimplicialTwoComplex(n,p);;
 K:=CliqueComplex(Y,4);;
 W:=RegularCWComplex(K);;
 W:=ContractedComplex(W);;
 b2:=BettiNumber(W,2);;
 Add(Betti2,[p,b2]);;
 b3:=BettiNumber(W,3);;
 Add(Betti3,[p,b3]);;
od;
ScatterPlot(Betti2);
ScatterPlot(Betti3);