dkahle / algstat Goto Github PK

The following will crash

A=matrix(c(1,1,1),1,3);
b=c(100000);
countFiber(A,b);

with the message

sh: line 1: 25994 Aborted                 (core dumped) '/opt/latte-integrale-1.7.1/bin/count' /tmp/RtmpJWaWTF/2015_06_19_08_45_05/countCode.latte > 'countOut' 2> 'countErr'
Error in file(con, "r") : cannot open the connection
In addition: Warning message:
In file(con, "r") :
  cannot open file 'numOfLatticePoints': No such file or directory

Since

A=matrix(c(1,1,1),1,3);
b=c("100000");
countFiber(A,b)

works just fine, I think the bug is due to the internal handling of 100000 in R which is 1e+05.

rvnorm() along with endgames provides points on the variety (or semi-algebraic set).

Hi,

would it be possible to extend the functionality of count so that it takes matrices A and right-hand sides b as arguments? That is, that something like

A=matrix(c(1,1,1,1),1,4);
b=c(10);
count(A,b);

is possible, where count(A,b) is the number of integer points in the b-fiber of A, i.e., all non-negative integer solutions u of Au=b?

Alternatively, it would also be fine if there would be a method which transforms the set of equations Ax=b, x>=0 into the form demanded by the specifications of count (for the example above: "x1+x2+x3=10","x1>=0","x2>=0","x3>=0").

Thanks!
Tobias

rvnorm() should allow the w argument to be a named list of (min, max) pairs. In particular, this will allow us to force one variable to be non-negative.

Hi, I want to use ggmap without signing up for a google key. Is it possible to do that? If so, which maps are available without a key?
Thanks