William,
I think I can reproduce one of your Itanium GAP bug; the workspace
filename gets mangled by SaveWorkspace.
(so the workspace gets saved, but with a horribly wrong name...)
Should be next to trivial to fix...
Dima
PS. Let's move this discussion over to sage-devel, as already
suggested here
I could be wrong but that problem might relate to the fact that
plotting is often done in floats, which can't handle quantities like
15^1024. Other types in Sage can handle such things, so you might
have to work around that limitation by plotting the log of the
function or something similar.
-M.
Thanks, indeed this solved the problem in the example.
Unfortunately, there is still a problem, if the degree of both
polynomials U and V increases to, say d = 1024. Note that the degree
of the rational function P = U/V is still 0 and both poles (0 and 1)
are far enough outside of the range wher
On Sat, Jan 9, 2010 at 11:36 AM, Dr. David Kirkby
wrote:
> Dima Pasechnik wrote:
>
>> William,
>> I emailed you few weeks back asking what exactly you mean by Itanium
>> environment, as I was unable to reproduce your problems on an Itanium
>> cluster I have access to--- at least not with gcc.
>> A
On Sat, Jan 9, 2010 at 1:11 PM, Shing Hing Man wrote:
> Hi,
> First, I define the ring of polynomial over the rationals, S.
> Then define a polynomial g in S.
> Is there a simple way to convert g to type
> 'sage.symbolic.expression.Expression' ?
>
> sage: S. = PolynomialRing(QQ);S
> Univ
Hi,
First, I define the ring of polynomial over the rationals, S.
Then define a polynomial g in S.
Is there a simple way to convert g totype
'sage.symbolic.expression.Expression' ?
sage: S. = PolynomialRing(QQ);S
Univariate Polynomial Ring in x over Rational Field
sage: g = x^3 - 11*
Somehow the nested form of the polynomials is causing problems. Using
expand, rather than full_simplify, seems to solve the problem:
d = 128
R = PolynomialRing(ZZ, x)
U = x^d + R.random_element(d-1)
V = x^d - x^(d-1)
U = expand(U)
V = expand(V)
P = U/V
G = P.plot(2, 15)
G.show()
Hope that helps,
Dima Pasechnik wrote:
William,
I emailed you few weeks back asking what exactly you mean by Itanium
environment, as I was unable to reproduce your problems on an Itanium
cluster I have access to--- at least not with gcc.
At least not in a stand-alone build of GAP.
Intel compilers showed to be tr
Hi,
I have two monic polynomials U, V of equal degree d with integer
coefficients. Furthermore the second one is of particularly simple
form, namely with roots only at 0 and 1. I want to plot the value of
their quotient in the range from 2 to 15.
Evaluating the quotient for any specific x in th
