(I tried to post about my problem once before, so my apologies if this
comes through twice.)
I'm receiving a segmentation error when I try to run a procedure in
the Singular console in Sage. I've previously tested this procedure
in Singular (run by itself) and it worked as designed, so I think that
something may be going wrong in the Sage interface.
Here's my procedure:
proc reducePoleOrder (poly pPoly, ideal Jacob, ideal xVariables) {
list divL = division(pPoly, Jacob);
print("division complete");
matrix divM = divL[1];
ideal divI = ideal(divM);
matrix redM = diff(xVariables, divI);
print("reduction complete");
poly redPoly = trace(redM);
print("trace complete");
return(redPoly);
}
If I type:
ring R = (0,t), (a,b,c), dp;
poly Q = 1/3*(a^3+b^3+c^3)-t*a*b*c;
ideal QIdeal = Q;
ideal abcVars = a,b,c;
ideal Jac = diff(abcVars, QIdeal);
ideal JacG = groebner(Jac);
short = 0;
and then try
reducePoleOrder(a^2*b^2*c^2, Jac, abcVars);
the correct result is displayed:
(-2*t^2)/(t^3-1)*a*b*c+(-t)/(t^3-1)*c^3
However, the command
poly testPoly = reducePoleOrder(a^2*b^2*c^2, Jac, abcVars);
testPoly;
returns the error message:
Singular : signal 11 (v: 3032/2007102509):
Segment fault/Bus error occurred at 3f because of 10202 (r:1193943964)
please inform the authors
trying to restart...
I used the Singular monitor("testfile.txt","io"); command to record
the input & output for this series of commands, and it has:
***
poly testPoly = reducePoleOrder(a^2*b^2*c^2, Jac, abcVars);
(null)
print("division complete");
matrix divM = divL[1];
ideal divI = ideal(divM);
matrix redM = diff(xVariables, divI);
print("reduction complete");
poly redPoly = trace(redM);
print("trace complete");
return(redPoly);
;return();
(null)
print("reduction complete");
poly redPoly = trace(redM);
print("trace complete");
return(redPoly);
;return();
(null)
print("trace complete");
return(redPoly);
;return();
***
It looks to me like extra null values and an empty return statement
are being inserted in my procedure, so that a meaningless value is
assigned to testPoly and the program crashes when I try to analyze
testPoly.
UAW
--~--~---------~--~----~------------~-------~--~----~
To post to this group, send email to [email protected]
To unsubscribe from this group, send email to [EMAIL PROTECTED]
For more options, visit this group at
http://groups.google.com/group/sage-support
URLs: http://sage.math.washington.edu/sage/ and http://sage.scipy.org/sage/
-~----------~----~----~----~------~----~------~--~---
