Singular

I would like to see what the function ClosureRees in reesclos.lib actually does by running it on test input,
but it seems not to be exported. What is the best way access this?

Doug

Open the the file reesclos.lib, which is located in your LIB subdirectory,
in a editor and but // in front of the command static.
Compare with the original file:
http://www.singular.uni-kl.de/svn/trunk ... esclos.lib

Save the file and load the the library again.
Now the proc reesclos.lib will be user available.

(If you don't want to modify the code, you can outsmart the behaviour
of the static proc by calling it from an another proc. Wrap it into
Then you may call closurerees instead.)

The appropriate input for ClosureRees is the output from the proc ReesAlgebra. So call, where I is an ideal

More information about the ongoing computation is obtained as follows:

1.) The author of the library has built in some messages which
will be shown when you set the system variable
printlevel to a large positive value. E.g.:
At some stages you have to hit the return key to let the computation
continue. (Reset by printlevel =0;)

2.) If you set TRACE = 3; then the code will be echoed, line by line,
and you have to press each time RETURN. (Reset by TRACE =0;)

Thanks, but....
I copied reesclos.lib to my home directory and put // in front of static as suggeted.
That did indeed let me start computing ClosureRees(ReesAlgebra(i)), which I couldn't do before.
But I quickly got an error message that I don't understand:

LIB "/home/leonada/reesclos.lib";
> ring r=2,(x,y),ds;
> ideal i=x2+x5,x2*x4y,y2+y5;
> list reesi=ReesAlgebra(i);reesi;
[1]:
// characteristic : 2
// number of vars : 5
// block 1 : ordering dp
// : names x y U(1) U(2) U(3)
// block 2 : ordering C
[2]:
// characteristic : 2
// number of vars : 3
// block 1 : ordering dp
// : names x y t
// block 2 : ordering C
> ClosureRees(reesi);
? incompatible type in list assignment
? error occurred in or before /home/leonada/reesclos.lib::ClosureRees line 353: ` map psi=R(1),mapI; // from ReesAlgebra: the map Rees->Kxt`
? expected map-expression. type 'help map;'
? leaving /home/leonada/reesclos.lib::ClosureRees
skipping text from `;` error at token `)`

It is a bit more complicated than I wrote.
Inspecting the code how ClosureRees is called from the proc normalI
one sees that one has to set the new basering first.

So you may proceed in your example:

which gives

Note: In your example you use the ordering ds. Do you really expect to work in the local case? As far as I see, computations are done with global ordering (see the examples and the description). I have chosen above the
ordering dp. But ds will give the same result since internally the
proc switches to the global ordering dp.