Re: [sage-support] Sage Crash Report

2022-08-22 Thread Robert Parini
nda activate sagetest > sage > > On Sun, Aug 21, 2022 at 6:24 PM Robert Parini wrote: > >> Using conda on macOS 12.4 (with Apple silicon) I get the attached error >> after installing sage with: >> >> conda create -n sagetest sage >> conda activate sagetest

Re: [sage-support] Sage Crash Report

2022-08-22 Thread Robert Parini
: > Can you share the output of the following command? > > conda list -n sagetest > > On Sun, Aug 21, 2022 at 6:24 PM Robert Parini wrote: > >> Using conda on macOS 12.4 (with Apple silicon) I get the attached error >> after installing sage with: >> >&

[sage-support] Sage Crash Report

2022-08-21 Thread Robert Parini
Using conda on macOS 12.4 (with Apple silicon) I get the attached error after installing sage with: conda create -n sagetest sage conda activate sagetest sage -- You received this message because you are subscribed to the Google Groups "sage-support" group. To unsubscribe from this group and s

[sage-support] Re: Bug (inconsistent behaviour) of graph6_string/sparse6_string

2019-11-05 Thread Robert Samal
r 5, 2019 at 1:43:23 AM UTC-8, Robert Samal wrote: > > I noticed the following strange behavior of > graph6_string()/sparse6_string() functions of graphs: > > sage: K2=graphs.CompleteGraph(2) > sage: P=K2.cartesian_product(K2) > > sage: print(P.sparse6_string()) > sage

[sage-support] Bug (weird behaviour) of graph6_string/sparse6_string

2019-11-05 Thread Robert Samal
I noticed the following strange behavior of graph6_string()/sparse6_string() functions of graphs: sage: K2=graphs.CompleteGraph(2) sage: P=K2.cartesian_product(K2) sage: print(P.sparse6_string()) sage: print(Graph(P.graph6_string()).sparse6_string()) :CoKN :Cci To explain: I understand, that

[sage-support] symbolic manipulation -- insufficient simplification

2019-11-01 Thread Robert Samal
I observed the following weird behavior of the symbolic engine. sage: x/x 1 sage: x^2/x x sage: (x^2+x)/x (x^2 + x)/x sage: assume(x>0) sage: assume(x,'real') sage: assumptions() [x > 0, x is real] sage: (x^2+x)/x (x^2 + x)/x To clarify: first, I consider the first two simplifications slightly

Re: [sage-support] Re: solving a linear system of GF(3) -- strange TypeError

2019-10-09 Thread Robert Samal
Indeed it works in Sage 8.4. Thanks! On Wednesday, October 9, 2019 at 8:34:41 AM UTC-7, Dima Pasechnik wrote: > > This got broken in Sage 8.5. > (still works in 8.4) > > > > On Wed, Oct 9, 2019 at 6:09 AM David Joyner > wrote: > >> >> >> On Wed,

[sage-support] Re: solving a linear system of GF(3) -- strange TypeError

2019-10-08 Thread Robert Samal
Sorry, F=GF(3), I made my original example shorter and didn't read it properly. So the full problematic code is B=matrix(GF(3), 2,2,[1,0,1,0], sparse=True) v=vector(GF(3), [1,1]) B.solve_right(v) Thanks, Robert On Tuesday, October 8, 2019 at 5:17:59 PM UTC-7, Robert Samal wrote: >

[sage-support] solving a linear system of GF(3) -- strange TypeError

2019-10-08 Thread Robert Samal
I am trying to solve a rather large linear systems of equations of GF(3). As the matrices are sparse, I thought that adding "sparse=True" to the constructor of the matrix could be of help. However, I ran to a strange error message. B=matrix(GF(3), 2,2,[1,0,1,0], sparse=True) v=vector(F, [1,1])

Re: [sage-support] Re: p.collect(x).collect(y) only collects y

2018-08-19 Thread Robert Webb Jr
Hi, Dima; thanks for taking the time to help out a noob -- again. This is what I did: sage: a,b,y,z=var('a,b,y,z') sage: p=z*x^2+x^2-(x^2+y^2)*(a*x-2*b*y)+z*y^2+y^2 sage: type(p) The book says that I should be able to do this: sage: p.collect(x).collect(y) And get this in return: sage: 2*b*y*

Re: [sage-support] Upgrade to 8.2 hangs at ncurses

2018-05-10 Thread Robert Gross
I must have set some other strange default.--Rob On 5/10/18, Jeroen Demeyer wrote: > On 2018-05-10 19:38, Robert Gross wrote: >> I have aliased cp to "cp -i" and "mv" to "mv -i". > > What do you mean this *exactly*? > > -- > You recei

Re: [sage-support] Upgrade to 8.2 hangs at ncurses

2018-05-09 Thread Robert Gross
On 5/9/18, I wrote: > Hi, > > Mac OS 10.13.4. Upgrade from 8.1 to 8.2. Hangs at > > [ncurses-6.0.p0] config.status: creating include/ncurses_cfg.h > > This happened when I used "sage --upgrade" and it occurs again with "sage > -i ncurses". Same happened building from scratch. The relevant par

[sage-support] Re: Find_root not finding a root

2018-03-15 Thread Robert Dodier
amiliar with, has its own plotting code, dunno about Sage, in any event just reusing the splitting algorithm isn't any big deal. best, Robert Dodier -- You received this message because you are subscribed to the Google Groups "sage-support" group. To unsubscribe from this group and

[sage-support] Re: Error with simple integral, "CEXPT only defined for non-negative integral exponents."

2017-11-09 Thread Robert Dodier
Never use inexact values in integral terms. Well, some experimentation shows that with keepfloat = false (the default), it works OK. But keepfloat = true (as in Sage) it bumps into this error. keepfloat causes trouble in other places too ... every now and then I think we (Maxima project) should du

[sage-support] Re: Integration bug?

2017-10-26 Thread Robert Dodier
ence of 'false' is a bug. If you can make a bug report in the Maxima bug tracker, that would very helpful. https://sourceforge.net/p/maxima/bugs By the way I am working with Maxima 5.40+ (almost 5.41). best, Robert Dodier -- You received this message because you are subscribed to the G

[sage-support] Re: taylor(1/4*x-x^2,x,0.2,3)

2017-07-12 Thread Robert Dodier
prevent that conversion. That causes trouble, given the widespread implicit assumption about exact numbers. It's a bug of course -- perhaps you can submit a bug report to: http://sourceforge.net/p/maxima/bugs but you can work around it by setting keepfloat to false, or writing 1/5 in

[sage-support] Re: Transforming a complex function f:C->C into a function g:R^2->R^2

2017-05-01 Thread Robert Jacobson
uces abs(z) when I ask for the real part of an expression. Consider: sage: y = var('y') sage: assume(x, 'real') sage: assume(y, 'real') sage: log(x+I*y).real_part() log(abs(x + I*y)) Best, Robert -- You received this message because you are subscribed to the Google Gr

[sage-support] Transforming a complex function f:C->C into a function g:R^2->R^2

2017-05-01 Thread Robert Jacobson
eplace abs(z) with sqrt(norm(z))? Best, Robert -- You received this message because you are subscribed to the Google Groups "sage-support" group. To unsubscribe from this group and stop receiving emails from it, send an email to sage-support+unsubscr...@googlegroups.com. To post t

[sage-support] Re: This limit seems to be wrong

2017-01-22 Thread Robert Dodier
/bugs/3280 for a related bug. HTH Robert Dodier -- You received this message because you are subscribed to the Google Groups "sage-support" group. To unsubscribe from this group and stop receiving emails from it, send an email to sage-support+unsubscr...@googlegroups.com. To po

[sage-support] Sage Crash Report

2016-12-27 Thread Robert Samal
-- You received this message because you are subscribed to the Google Groups "sage-support" group. To unsubscribe from this group and stop receiving emails from it, send an email to sage-support+unsubscr...@googlegroups.com. To post to this group, send email to sage-support@googlegroups.com. V

[sage-support] Re: error using simple solve example

2016-08-24 Thread Robert Dodier
ed. Well, Sage is punting 'solve' to Maxima, and it looks like Maxima isn't loaded correctly or something. SET-LOCALE-SUBDIR isn't anything important, it's only setting the subdirectory for the on-line documentation. One could delete that call, but I'm guessing one would then

[sage-support] Re: Easy definite integral

2016-08-19 Thread Robert Dodier
ng to classify the type of problem, and using asksign for that. Maxima can actually get the result without resorting to asksign, but it doesn't try that first. The relevant function, if anybody cares, is METHOD-BY-LIMITS in src/defint.lisp. I dunno if this is a bug. It is needlessly clumsy bu

[sage-support] Re: Possible bug in numerical integration

2016-06-02 Thread Robert Dodier
uming that Sage punts to Maxima for real() here. But integral_numerical is probably not calling Maxima, right? best Robert Dodier -- You received this message because you are subscribed to the Google Groups "sage-support" group. To unsubscribe from this group and stop receiving email

[sage-support] Re: incorrect limit

2016-04-19 Thread Robert Dodier
an experiment) introduces its own problems though. That exposes some bugs in gruntz, and also some results which are different, so it would be necessary to trawl through them and verify that they're correct. best Robert Dodier -- You received this message because you are subscribed to the Go

[sage-support] Re: wrong result for integral(sqrt(cos(x)-cos(x)^3), x, 0, pi/2)

2016-04-14 Thread Robert Dodier
l and domain:complex). Reported as: https://sourceforge.net/p/maxima/bugs/3126/ best Robert Dodier -- You received this message because you are subscribed to the Google Groups "sage-support" group. To unsubscribe from this group and stop receiving emails from it, send an email to sage

[sage-support] Re: External viewer

2016-01-07 Thread Robert Samal
/x86_64-linux-gnu/libgpg-error.so.0 (0x7f627fec5000) On Thursday, January 7, 2016 at 6:38:53 PM UTC+1, Volker Braun wrote: > > Can you tell us more about what eog links to? > > $ sage -sh command -v eog > /usr/bin/eog > $ sage -sh ldd /usr/bin/eog # adjust p

[sage-support] Re: Issue with Solve

2015-02-27 Thread Robert Dodier
e^(-1)] For the record, when I try this directly in Maxima, the equation is solved in both cases. (I assume that Sage calls Maxima to solve equations; is that right?) Heaven knows Maxima has > 1 bug in which results depend on the names of variables, but I guess this isn't one of them. best,

Re: [sage-support] Coercion problem

2015-01-15 Thread Robert Bradshaw
I assume you meant sage: v = P(5) sage: v(oo) A positive finite number This is because the elements of QQ coerce to the parent of oo, which is the "signed infinity ring." This is so we have sage: P. = PolynomialRing(QQ) sage: w = x + 5 sage: v = w - x w(1.0) 6.00 sage: v(1.0) 5.0

[sage-support] Re: sage hardcoded links problem

2015-01-07 Thread Robert McBroom
Fixes to the source files? On Tuesday, January 6, 2015 3:23:10 AM UTC-5, Volker Braun wrote: > > I built a new binary. Harald, can you replace the F21 binary with the new > one on the mirrors? > > buildbot@build:~$ md5sum > binaries/sage-6.4.1-x86_64-Linux-Fedora_21_x86_64.tar.gz > cbdf8a8a1e4ee

[sage-support] sage -upgrade from 6.2 error

2015-01-01 Thread Robert McBroom
sage -upgrade from 6.2 error on Fedora 19. Running in a root window upgrade terminates with a permission error Target: x86_64-redhat-linux Configured with: ../configure --prefix=/usr --mandir=/usr/share/man --infodir=/usr/share/info --with-bugurl=http://bugzilla.redhat.com/bugzilla --enable-b

[sage-support] Re: Manually moving sage data to mathematica?

2014-12-21 Thread Robert Jacobson
I think the Mathematica interface is still broken. I'm looking into what it would take to fix it. What kind of sage object are we talking about? On Sunday, 21 December 2014 17:01:47 UTC-5, Shane Scott wrote: > > Am I correct in thinking the sage-mathematica interface is still broken? > If tha

[sage-support] Status of the Mathematica interface?

2014-12-08 Thread Robert Jacobson
and has been broken for a couple of years. Do I understand correctly that the Mathematica interface is still not working? Best, Robert -- You received this message because you are subscribed to the Google Groups "sage-support" group. To unsubscribe from this group and stop receiv

[sage-support] Re: simple 'solve' error, possibly related to ecl and maxima installation

2014-11-23 Thread Robert Dodier
m working on it, for the record. best Robert Dodier -- You received this message because you are subscribed to the Google Groups "sage-support" group. To unsubscribe from this group and stop receiving emails from it, send an email to sage-support+unsubscr...@googlegroups.com. To post

[sage-support] Re: Why solve(5^( x -1) == (0.04)^(2*x), x) returns empty set?

2014-11-19 Thread Robert Dodier
isn't entirely consistent -- when p isn't decidable to Maxima, you can get a partially-evaluated conditional, but not a partially-evaluated loop (triggers an error), and various programming functions (e.g. length, first, integerp) might act in an unexpected way. This, too, hasn't

[sage-support] Re: Bug in symbolic integral

2014-11-12 Thread Robert Dodier
e(-1,-log(2)) (%i7) %, numer; (%o7) 1.273097216447114 (%i8) quad_qags (foo, x, 2, 3); (%o8) [1.273097216447114,1.413421842285782E-14,21,0] Looks like Maxima handles the definite and indefinite integrals as expected. Or perhaps I have misunderstood the problem? Hope this helps, Robert Dodier -

[sage-support] Re: Having issues with assume

2014-10-28 Thread Robert Dodier
edundatnt as I have already made this clear above. Is this a bug, or am I > missing something? This is a bug. If you have time, can you please report it to the Maxima bug tracker: http://sourceforge.net/p/maxima/bugs best, Robert Dodier -- You received this message because you are subscri

Re: [sage-support] Real Interval Field (RIF) question

2014-10-22 Thread Robert Bradshaw
I do find this behavior quite surprising--diameter should be an alias for either relative or absolute diameter, not depending on the interval. On Wed, Oct 22, 2014 at 2:44 AM, John Cremona wrote: > I am trying to use the Real Interval Field (RIF), which in principle > does exactly what I want. B

[sage-support] Re: matrix log?

2014-10-21 Thread Robert Dodier
mat_function (log, mymatrix); For large, numerical matrices, I'm sure linalg.logm is much faster. But for small or nonnumerical matrices, maybe Maxima is useful. Hope this helps, Robert Dodier -- You received this message because you are subscribed to the Google Groups "sage-support&

[sage-support] Re: solve(-(1/2*sqrt((4*w+1)+1))*t+w==0,w)

2014-10-16 Thread Robert Dodier
an solve a relatively narrow range of problems). But I find that Maxima's 'to_poly_solve' can solve it. Maybe someone here can say how to call it from Sage. > In fact, the solution is: w=t+t^2 Are you sure? Assuming some value for t, plotting the expression doesn't seem to s

[sage-support] Re: Fwd: [sage-cloud] Strange behavior when integrating 1/x^p

2014-10-15 Thread Robert Dodier
would file a bug report. http://sourceforge.net/p/maxima/bugs best, Robert Dodier -- You received this message because you are subscribed to the Google Groups "sage-support" group. To unsubscribe from this group and stop receiving emails from it, send an email to sage-support+unsubscr...@

Re: [sage-support] QEPCAD in sagecell

2014-10-08 Thread Robert Pollak
Am 08.10.2014 um 11:54 schrieb fbotana: > Nevertheless, it seems that QEPCAD is not anymore installed/working. Try > http://sagecell.sagemath.org/?q=xlgebu That's a pity. I was using it [1], but I didn't find the time to learn how to fix (and update) the broken [2] qepcad spkg. R

[sage-support] Unsafe directory at startup

2014-10-06 Thread robert
Installed pre-compiled Sage-6.3.app on Mac OS X 10.9.5 At startup of Terminal Session, there is this ~$ /Applications/Sage-6.3.app/Contents/Resources/sage/sage; exit sys:1: RuntimeWarning: not adding directory '' to sys.path since everybody can write to it. Untrusted users could put files in this

[sage-support] Re: Segfault on integration

2014-09-24 Thread Robert Dodier
On 2014-09-20, Kristoffer Ryhl-Johansen wrote: > f(x)=log(1-x)*log(1+x)/(1+x) >> >> f.integrate(x,0,1) >> > Produces a segfault when I run it on my ubuntu 14.04 computer Fixed by Maxima commit f7921c5265 (bug in Risch code). best Robert Dodier -- You received this

[sage-support] Re: Segfault on integration

2014-09-20 Thread Robert Dodier
ed by one of the other packages loaded by abs_integrate, but that doesn't seem to be the case; so I guess the problem is triggered by abs_integrate itself. I will try to investigate some more. If someone files a bug report, that will help us track it. http://sourceforge.net/p/maxima/bugs b

[sage-support] Re: apparent numerical integration bug in sage

2014-08-29 Thread Robert Dodier
> Harald, thanks. I think now we have six or seven ways to do quadrature in > Sage! Yes, but most of them are QUADPACK, right? best Robert Dodier -- You received this message because you are subscribed to the Google Groups "sage-support" group. To unsubscribe from this group and s

[sage-support] Re: apparent numerical integration bug in sage

2014-08-29 Thread Robert Dodier
functions in Maxima -- I find that quad_qags(f2, x, 1, 10^11) fails (with error=5, "integral is probably divergent or slowly convergent") but quad_qag(f2, x, 1, 10^11, 4) succeeds, likewise quad_qagi(f2, x, 1, inf) succeeds. If Sage is indeed calling QUADPACK, perhaps at least the error num

Re: [sage-support] Why Sage is called Sage ?

2014-08-14 Thread Robert Bradshaw
http://sagemath.blogspot.com/2009/12/mathematical-software-and-me-very.html On Thu, Aug 14, 2014 at 1:14 AM, John Cremona wrote: > When William Stein first started the project it was an acronym SAGE > for (I think) System for Algebra and Geometry Experimentation. But > soon it became a much wide

[sage-support] Re: Simplifying combinations of atomic inequalities

2014-07-18 Thread Robert Pollak
On Wednesday, July 16, 2014 10:25:03 AM UTC+2, I wrote: > > I see the following wrong results: > > sage: x<2 and x<1 > x < 2 > sage: x<2 or x<1 > x < 1 > I found a way to compute these with Sage: sage: qepcad(qepcad_formula.and_(x < 2, x < 1), vars='(x)') x - 1 < 0 sage: qepcad(qepcad_formul

Re: [sage-support] Typo in Documentation of Eigenmatrix

2014-07-17 Thread Robert Bradshaw
https://github.com/sagemath/sage/pull/21 aka http://trac.sagemath.org/ticket/16672 On Thu, Jul 17, 2014 at 9:45 AM, Mahrud Sayrafi wrote: > Hi, > > In this page: > http://www.sagemath.org/doc/constructions/linear_algebra.html#eigenvectors-and-eigenvalues > in the eigenvectors and eigenvalues sect

Re: [sage-support] Re: Simplifying combinations of atomic inequalities

2014-07-17 Thread Robert Pollak
Also, why do I need two steps here?: sage: solve([x==0, x!=1], x) [[x == 0, -1 != 0]] solve([x == 0, -1 != 0], x) [x == 0] -- You received this message because you are subscribed to the Google Groups "sage-support" group. To unsubscribe from this group and stop receiving emails from it, send a

Re: [sage-support] Re: Simplifying combinations of atomic inequalities

2014-07-17 Thread Robert Pollak
Am 17.07.2014 11:32, schrieb Robert Pollak: > In fact, I do not even know how to create the "and" situation. > The following should be "x<=2 and x>=0": > > sage: Polyhedron(ieqs=[(2,-1), (0,0)]).Hrepresentation() > (An inequality (-1) x + 2 >= 0,) Oop

Re: [sage-support] Re: Simplifying combinations of atomic inequalities

2014-07-17 Thread Robert Pollak
Am 16.07.2014 21:06, schrieb slelievre: > Robert Pollak wrote: > I see the following wrong results: > > sage: x<2 and x<1 > x < 2 > sage: x<2 or x<1 > x < 1 > > The best way to manipulate logical combination of inequalities

Re: [sage-support] Re: Simplifying combinations of atomic inequalities

2014-07-17 Thread Robert Pollak
Am 16.07.2014 20:41, schrieb Nils Bruin: > On Wednesday, July 16, 2014 1:25:03 AM UTC-7, robert.pollak wrote: > sage: x<2 and x<1 > x < 2 > sage: x<2 or x<1 > x < 1 > > That's because "and" and "or" are program flow constructs in python, as > they are in C (they have "shortcut eval

[sage-support] Simplifying combinations of atomic inequalities

2014-07-16 Thread Robert Pollak
Hello! I see the following wrong results: sage: x<2 and x<1 x < 2 sage: x<2 or x<1 x < 1 Is this just a syntax problem? How would I enter this correctly? Robert -- You received this message because you are subscribed to the Google Groups "sage-support" group. To

Re: [sage-support] GIT repository and installation

2014-07-15 Thread Robert Bradshaw
The very short answer is to type "make" and wait an hour or three. On Jul 14, 2014 4:54 AM, "Oscar Alberto Castillo Felisola" < o.castillo.felis...@gmail.com> wrote: > Checking it out! Thank you John. > > -- > You received this message because you are subscribed to the Google Groups > "sage-suppor

[sage-support] Re: Newbie qepcad.py syntax problem

2014-07-11 Thread Robert Pollak
> After reading qepcad.py I tried to translate this to a qepcad call: > > qepcad(qepcad_formula.or_(qepcad_formula.and_(x-5 > 0, 3(x-1) <= > x-5), qepcad_formula.and_(x-5 < 0, 3(x-1) >= x-5)), vars='(x)') > > > 3(x-1) tries to call 3 with the argument x-1. That explains the error > y

[sage-support] Newbie qepcad.py syntax problem

2014-07-10 Thread Robert Pollak
Integer(0), Integer(3)(x-Integer(1)) >= x-Integer(5))), vars='(x)') TypeError: 'sage.rings.integer.Integer' object is not callable Can you please help me to compose the correct call? (Because of http://trac.sagemath.org/ticket/16642 I am currently using http://sagecell.sagemath.org/, which has

[sage-support] Re: Won't Integrate

2014-07-05 Thread Robert Dodier
rtainly a drawback. I don't know if e.g. SymPy could solve it; I didn't try. best Robert Dodier -- You received this message because you are subscribed to the Google Groups "sage-support" group. To unsubscribe from this group and stop receiving emails from it, s

[sage-support] Re: Having trouble with syntax and nesting integral into limit

2014-06-05 Thread Robert Dodier
s n increases without bound. So I'm guessing the limit is zero. If you need a proof, maybe you can show the integral is bounded, therefore the limit is zero. Sorry I can't be more helpful, Robert Dodier -- You received this message because you are subscribed to the Google Groups &

[sage-support] Re: ECL says: In function GCD, the value of the second argument is not of the expected type INTEGER

2014-06-03 Thread Robert Dodier
ace 6.02e-23 with 602*10^-25. best Robert Dodier -- You received this message because you are subscribed to the Google Groups "sage-support" group. To unsubscribe from this group and stop receiving emails from it, send an email to sage-support+unsubscr...@googlegroups.com. To post t

[sage-support] Runtime Error: Gap Tables of Marks Library not installed

2014-05-30 Thread Robert Godfroid
host: Windows 8.1 VirtualBox 4.3.12 guest: Ubuntu 14.04 LTS Sage 6.2 Release Date 2014-05-06 Statements that gave rise to the error: A5 = AlternatingGroup(5) A5_sgs = A5.subgroups() len(A5_sgs) => ... RuntimeError: Gap produced error output Error, sorry, the GAP Tables of Marks Library is not ins

[sage-support] Re: Does Sage support content mathML?

2014-05-29 Thread Robert Dodier
MathML and presentation MathML. I tinkered with maximaMathML a couple of months ago and it seemed to work OK (after fixing some bugs). Write me off-list if you're interested. best Robert Dodier -- You received this message because you are subscribed to the Google Groups "sage-support" group

Re: [sage-support] Re: Simplify number

2014-05-29 Thread Robert Bradshaw
What exactly do you mean by "simplify a real number?" On Thu, May 29, 2014 at 8:32 AM, SiL588 . wrote: > Unfortunately I don't know the rules of Phyton language, i just started > using Sage notebook to do linear algebra computation. > I think I did what you said, I assinged m a value that was the

Re: [sage-support] Re: bug in comparison of function field elements

2014-04-30 Thread Robert Bradshaw
On Tue, Apr 29, 2014 at 10:57 PM, Robert Bradshaw wrote: > On Tue, Apr 29, 2014 at 9:07 AM, Volker Braun wrote: >> On Tuesday, April 29, 2014 3:58:14 PM UTC+1, Simon King wrote: >>> >>> Yes there is! The hook is the hash function. >> >> >> CPytho

Re: [sage-support] Re: bug in comparison of function field elements

2014-04-29 Thread Robert Bradshaw
and you'd probably want to cache it. >> ... in some cases only a trivial hash function (such as: hash of the >> parent) should be used. > > > or, better, just 1: set([ZZ(1), QQ(1)]) This is probably a bad idea--it'll lead to very poor and hard-to-diagnose performance

[sage-support] Re: wrong or nonexistent results for various infinite sums

2014-04-23 Thread Robert Dodier
simplify_sum ('sum(x^(3*k)/factorial(2*k),k,0,inf)); => sqrt(%pi)*bessel_i(-1/2,x^(3/2))*x^(3/4)/sqrt(2) Hope this helps, Robert Dodier -- You received this message because you are subscribed to the Google Groups "sage-support" group. To unsubscribe from this group and stop

[sage-support] Re: How to compute sum(log(1-1/n^2),n,2,oo)?

2014-04-16 Thread Robert Dodier
19) log((m+1)/(2*m)) (%i20) limit (%, m, inf); (%o20) -log(2) best, Robert Dodier -- You received this message because you are subscribed to the Google Groups "sage-support" group. To unsubscribe from this group and stop receiving emails from it, send an email to sage-support+unsubscr...@

[sage-support] Re: Integral result discrepancy between SAGE and MAPLE

2014-04-15 Thread Robert Dodier
s and maxima. I > can't recall the precise details right now, but after that I "learned" to > defer substitution to the very end. A different way to handle that is to change the floats to rational numbers, via ratsimp. Then Maxima will be happy with expressions which contain those

Re: [sage-support] Re: Optimizing this function with Cython

2014-03-15 Thread Robert Bradshaw
On Fri, Mar 14, 2014 at 2:03 PM, Georgios Tzanakis wrote: > > On Fri, Mar 14, 2014 at 4:49 PM, Robert Bradshaw > wrote: >> >> Note that >> >> L[i][rows[i]] + j %w == 0: >> >> would probably be just (or nearly) as fast as >> >>

Re: [sage-support] Re: Optimizing this function with Cython

2014-03-14 Thread Robert Bradshaw
Note that L[i][rows[i]] + j %w == 0: would probably be just (or nearly) as fast as (((L[i])[(rows[i])])+j %w)==0 If you're going to be dealing with arrays of ints you might want to look into NumPy and/or memory views for even more speed. On Thu, Mar 13, 2014 at 7:58 PM, Georgios Tzanakis w

[sage-support] Re: Some frustrations with Piecewise()

2014-01-28 Thread Robert Jacobson
to, I see that the patch is 7 months old. I'm curious, is it in active development? I sure would like to see a solid piecewise implementation in an upcoming Sage release. Thanks again. --Robert On Tuesday, 28 January 2014 22:14:39 UTC-5, kcrisman wrote: > > > > On Tuesday, Janu

[sage-support] Some frustrations with Piecewise()

2014-01-28 Thread Robert Jacobson
, 0), 1+x], [(0, 1), 1-x], [(1, 10), 0*x^0]], x) This seems to evaluate correctly within its domain, but I can't plot it, say, with this (or variations thereof): plot(t(x), x, -4, 4) What am I doing wrong? --Robert -- You received this message because you are subscribed to the Google G

Re: [sage-support] Re: README.txt build-from-source buries fact -j8 in "make -j8" ignored -- need export MAKE="make -j8"

2014-01-02 Thread Robert Bradshaw
uld not be enabled in the top-level make, in my opinion. Typically, > "make -jN" makes parallel compiles within the same package (if the package > supports it). The parallel build in sage compiles different packages in > parallel, but each package still compiles as -j1. > >

Re: [sage-support] README.txt build-from-source buries fact -j8 in "make -j8" ignored -- need export MAKE="make -j8"

2014-01-02 Thread Robert Bradshaw
How hard would it be to let "make -jN" actually work from the top-level make? On Tue, Dec 31, 2013 at 4:57 PM, Joseph P. Skudlarek wrote: > This is a request to update the README.txt file used when building from > sources -- the README.txt buries the fact that "-jN" in "make -jN" is > effectively

[sage-support] Re: simplification of symbolic expressions

2013-11-24 Thread Robert Dodier
(%i12) asin(sqrt(3)*sqrt(2)/3) + asin(sqrt(3)/3) - asin(x) + u; (%o12) -asin(x)+u+%pi/2 (%i13) asin(sqrt(3)*sqrt(2)/3) + asin(sqrt(3)/3) + asin(sqrt(1 - u^2)) + asin (u); (%o13) %pi This is just what I got from some half-hearted hacking; I'm sure there are serious limitations. Good

Re: [sage-support] Re: (unknown)

2013-10-23 Thread Robert Bradshaw
sage: Integers(45)['t'] Univariate Polynomial Ring in t over Ring of integers modulo 45 I don't think we have linear algebra over general non-integral-domains, but sage: R = GF(5)['x'] sage: M = random_matrix(R, 4, 4); b = random_vector(R, 4); x = M \ b sage: M*x (4*x^2 + x + 4, x^2 + 2*x + 4, 4*

Re: [sage-support] arrays in Sage's cython

2013-08-21 Thread Robert Bradshaw
Using a Python list is probably the fastest way to iterate over an array of Python objects--it's a PyObject** under the hood and Cython uses the C API calls to get at it. Your "check" might be the bottleneck, especially if it's a Python call. Also, no need to write this as a while loop; just use "

Re: [sage-support] Re: Left and right limits (is default of 'dir' None?)

2013-08-20 Thread Robert Bradshaw
On Mon, Aug 19, 2013 at 6:30 PM, Dima Pasechnik wrote: > On 2013-08-19, Vincent Knight wrote: >> --001a1133aa8653f2ed04e4510b09 >> Content-Type: text/plain; charset=ISO-8859-1 >> >> Thanks for the answer kcrisman but I'm afraid I'm still not sure I >> understand. >> >> If by 'unsigned infinity' y

[sage-support] Re: Integration of formal functions

2013-07-27 Thread Robert Dodier
...)). 'integrate is a formal integral -- it doesn't invoke the code to solve definite integrals -- so it won't bump into that error. best Robert Dodier -- You received this message because you are subscribed to the Google Groups "sage-support" group. To unsubscribe

[sage-support] Re: Error summing inverses of binomial coefficients

2013-07-22 Thread Robert Dodier
On 2013-07-18, Ed Scheinerman wrote: > sage: sum(1/binomial(n,k),k,0,n) > (n + 1)*2^(-n) > > and that answer is wrong. That's a bug in Maxima's simplify_sum -- reported as bug # 2614. https://sourceforge.net/p/maxima/bugs/2614/ best Robert Dodier -- You received this m

[sage-support] Re: Bug in evaluation of Maxima numbers? Or where? Or not a bug at all?

2013-07-19 Thread Robert Dodier
oat called a bigfloat, which is not a CL type. HTH Robert Dodier -- You received this message because you are subscribed to the Google Groups "sage-support" group. To unsubscribe from this group and stop receiving emails from it, send an email to sage-support+unsubscr...@googlegroups

[sage-support] Google OpenID account problems for sagenb.org

2013-05-08 Thread Robert Nendorf
I've recently had issues trying to log in to Sage via sagenb.org I had set up a Sage account using my Google account. Now when I attempt to log in it takes me to a page where I choose which Google account I want to use. So far so good. But when I click on an account, it takes me to the new acc

[sage-support] Re: Substituting a combination of variables in a matrix

2013-04-23 Thread Robert Dodier
nno what might be available in Sage proper. best Robert Dodier -- You received this message because you are subscribed to the Google Groups "sage-support" group. To unsubscribe from this group and stop receiving emails from it, send an email to sage-support+unsubscr...@googlegroups.c

Re: [sage-support] Solve Binary Linear System

2013-04-17 Thread Robert Bradshaw
sage: A = random_matrix(GF(2), 1, 1) sage: A.det() 1 sage: b = random_vector(GF(2), 1) sage: %time x = A \ b CPU times: user 1.61 s, sys: 0.06 s, total: 1.67 s Wall time: 1.67 s sage: A * x == b True On Wed, Apr 17, 2013 at 1:45 PM, Juan Grados wrote: > I have the equation Ax=b where

Re: [sage-support] Inconsistency regarding QQ and ZZ?

2013-03-17 Thread Robert Bradshaw
The syntax "R. = QQ[]" creates a polynomial ring in two variables, with generators A and d (bound to the current session). A^d is not a polynomial in A and d over QQ. sage: R.=QQ[] sage: R Multivariate Polynomial Ring in A, d over Rational Field On Sun, Mar 17, 2013 at 12:51 AM, Rolandb wrote: >

Re: [sage-support] sage vs matlab

2013-03-05 Thread Robert Bradshaw
ists makes more sense, I know a lot of people who have switched from MATLAB to those tools (of which Sage is a superset). It also depends on how much freedom you'll have--if you get a job where everyone uses MATLAB you might not have much choice but if you're doing research than you can use

[sage-support] Re: Bug in solve?

2013-02-24 Thread Robert Dodier
ma source on Linux. (%i1) load (to_poly_solve); Loading maxima-grobner $Revision: 1.6 $ $Date: 2009-06-02 07:49:49 $ (%o1) /home/robert/maxima/maxima-git/maxima-code/share/to_poly_solve/to_poly_s\ olve.mac (%i2) to_poly_solve ([(a*x+b*y)*x*y/c=1,3*log(a + b + c) - log(27*a*b*x*y)],[x,y]); Maxima

[sage-support] zn_poly-09.p09 error on Mac during upgrade to 5.7.

2013-02-23 Thread Robert Gross
Hi, I upgraded to 5.7, and I get an error from zn_poly-0.9.p9 when running the quick self-test: nuss_mul()... FAIL!. This is on a Mac OX 10.6. I did manage to install the rest of sage-5.7, and I can start sage successfully, so I can try to test anything that anyone can suggest to pinpoint the p

Re: [sage-support] using cpickle to store matrix object

2013-02-20 Thread Robert Bradshaw
It's cPickle with a capital P. On Wed, Feb 20, 2013 at 2:30 AM, akhil wrote: > Hello, > > > I want to use cpickle to store a matrix object in a text file. Sample code > is as follows: > > A = matrix(GF(2),2,3) #creating a 2 * 3 matrix having all entries > zero > > import cpickle as pickl

[sage-support] Re: Limit and variable

2013-02-01 Thread Robert Dodier
). So that would work if the exponent is an integer. However, I doubt if Maxima can do anything unless the exponent is a literal integer, so the question seems pointless. I'd have to look at it again before figuring out if the question could be skipped. best Robert Dodier -- You received thi

Re: [sage-support] Link explaining types of decimal numbers

2013-01-23 Thread Robert Bradshaw
1000 bits, not just 53 bits (as it would be if 1.2 was parsed to 53 bits then passed in to RealField(1000)). For all of the above, see RR?, RDF?, etc. for (lots!) more documentation. - Robert On Wed, Jan 23, 2013 at 8:59 AM, LFS wrote: > Hi - I would appreciate if someone could point

[sage-support] parametric_plot and xmin/xmax

2013-01-20 Thread Robert Jacobson
Setting xmin/xmax for parametric_plot doesn't seem to do anything, but ymin/ymax work as expected. What am I doing wrong? t = var('t') parametric_plot( (cos(t), sin(t)), (t, 0, 2*pi), xmin=-2, xmax=2, ymin=-2, ymax=2) -- You received this message because you are subscribed to the Google Groups

[sage-support] Re: solving Diophantine equations in Sage

2012-12-08 Thread Robert Dodier
entations of these approaches in Sage or any upstream project? (e.g. PARI/GP, Singular, I don't know.) best, Robert Dodier -- You received this message because you are subscribed to the Google Groups "sage-support" group. To post to this group, send email to sage-support@googlegrou

[sage-support] solving Diophantine equations in Sage

2012-12-08 Thread Robert Dodier
the basis for an implementation? Pointers to any other resources would be interesting. best Robert Dodier -- You received this message because you are subscribed to the Google Groups "sage-support" group. To post to this group, send email to sage-support@googlegroups.com. To unsubs

[sage-support] Re: Unable to Solve Simple Problem

2012-10-27 Thread Robert Dodier
%if(?%and(-%pi/2 < parg(sqrt(4*c+b^2)+b), parg(sqrt(4*c+b^2)+b) <= %pi/2), [d = (b*sqrt(4*c+b^2)+2*c+b^2)/2],%union())) I didn't check the result; sorry about that. best Robert Dodier -- You received this message because you are su

Re: [sage-support] algebraic number to integer conversion?

2012-10-17 Thread Robert Bradshaw
On Wed, Oct 17, 2012 at 9:37 PM, Eric Kangas wrote: > code: > > b = 11^2 > > a = b^2 > > pri = [int(is_prime(i)) for i in range(a)] > > j = [i for i in range(a)][b+1:a:b] > > k = [i for i in range(a)][(b*2)+1:a:b] > > j.insert(0,0) > > k.insert(0,b) > > m = [matrix(QQ,sqrt(a)/(b/sqrt(b)),pri[j[i]:

Re: [sage-support] algebraic number to integer conversion?

2012-10-17 Thread Robert Bradshaw
What is eval? It looks like a list of eigenvalues, rather than a single one. Could you show how it is computed? On Wed, Oct 17, 2012 at 4:55 PM, Eric Kangas wrote: > Of course complex. > > Here is what happens when I try using the embeddings() function: > > sage: eval[0].embeddings(CC) > >

Re: [sage-support] algebraic number to integer conversion?

2012-10-17 Thread Robert Bradshaw
When you say "plot these values", do you mean as real or complex values? To do so you need to choose an embedding, e.g. sage: K. = QQ[sqrt(5)]; K Number Field in sqrt5 with defining polynomial x^2 - 5 sage: K.embeddings(CC) [ Ring morphism: From: Number Field in sqrt5 with defining polynomial

Re: [sage-support] GPL and code generation using sage

2012-09-27 Thread Robert Bradshaw
On Thu, Sep 27, 2012 at 1:41 PM, Tom Boothby wrote: > On Thu, Sep 27, 2012 at 12:55 AM, Robert Bradshaw > wrote: > >> I would be extremely >> surprised if any Sage developer morally objects to you licensing this >> output as you wish (though opinions may vary w

Re: [sage-support] GPL and code generation using sage

2012-09-27 Thread Robert Bradshaw
On Wed, Sep 26, 2012 at 10:56 PM, Geoffrey Irving wrote: > On Wed, Sep 26, 2012 at 10:42 PM, Robert Bradshaw > wrote: >> On Wed, Sep 26, 2012 at 8:54 PM, Geoffrey Irving wrote: >>> On Wed, Sep 26, 2012 at 6:03 PM, Robert Bradshaw >>> wrote: >>>> On Wed

  1   2   3   4   5   6   7   8   9   10   >