File: deformed_cyclic.pl

package info (click to toggle)
polymake 4.6-5
  • links: PTS, VCS
  • area: main
  • in suites: bookworm
  • size: 32,288 kB
  • sloc: cpp: 156,233; perl: 42,962; javascript: 30,726; ansic: 2,907; java: 2,654; python: 641; sh: 244; xml: 117; makefile: 61
file content (18 lines) | stat: -rw-r--r-- 509 bytes parent folder | download | duplicates (5) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
 
sub deformed_cyclic{
    my($n,$d) = @_;
    my $m = monomials<Rational,Rational>(1);

    my $vm = new Matrix<PuiseuxFraction<Min>>();
    for (my $k=1; $k<=$n; $k++) {
    	my @l = (1);
    	for (my $i=1; $i<=$d; $i++) {
    	    push @l, $m^($k*$i); #($m+$k)^($i);  #
    	    push @l, $m^(-$k*$i); #1/(($m+$k)^($i));  #
    	}
    	my $v = new Vector<PuiseuxFraction<Min>>(@l);
	$vm /= $v;
    }
    my $um = unit_matrix<PuiseuxFraction<Min>>(2*$d+1);
    $vm /= ($um->minor(~[0],All));
    return $vm;
}