File: deformed_cyclic.pl

package info (click to toggle)
polymake 4.14-2
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid
  • size: 35,888 kB
  • sloc: cpp: 168,933; perl: 43,407; javascript: 31,575; ansic: 3,007; java: 2,654; python: 632; sh: 268; 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;
}