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

[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,

[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,

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:

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

2018-05-10 Thread Robert Gross
me other strange default.--Rob On 5/10/18, Jeroen Demeyer <j.deme...@ugent.be> 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 receiv

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

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

2018-03-15 Thread Robert Dodier
iliar 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 s

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

2017-11-09 Thread Robert Dodier
k we (Maxima project) should dump it. 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 pos

[sage-support] Re: Integration bug?

2017-10-26 Thread Robert Dodier
n(3*x)+3*sin(x))/24 Obviously the presence 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 yo

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

2017-07-12 Thread Robert Dodier
e, to 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 instead

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

2017-05-01 Thread Robert Jacobson
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 Groups "sage-support" group. To unsubscribe from th

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

2017-05-01 Thread Robert Jacobson
))? 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 to this group, send email to sa

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

2017-01-22 Thread Robert Dodier
ceforge.net/p/maxima/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. T

[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.

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

2016-08-24 Thread Robert Dodier
run into a similar error with the next call to a Maxima function. Sorry I can't be more helpful, 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 ema

[sage-support] Re: Easy definite integral

2016-08-19 Thread Robert Dodier
msy but not incorrect. HTH in some way. 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.

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

2016-06-02 Thread Robert Dodier
I guess I'm also assuming 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 an

[sage-support] Re: incorrect limit

2016-04-19 Thread Robert Dodier
I tried as 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 subscrib

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

2016-04-14 Thread Robert Dodier
e same for domain:real 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

[sage-support] Re: External viewer

2016-01-07 Thread Robert Samal
f627fec5000) 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 path as necessary > > Also, are you overridin

[sage-support] Re: Issue with Solve

2015-02-27 Thread Robert Dodier
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, Robert Dodier -- You received this message because

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.x = PolynomialRing(QQ) sage: w = x + 5 sage: v = w - x w(1.0) 6.00 sage: v(1.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

[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

[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

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

2014-11-23 Thread Robert Dodier
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 to this group, send

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

2014-11-19 Thread Robert Dodier
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 caused trouble, from what I remember. FWIW Robert Dodier -- You

[sage-support] Re: Bug in symbolic integral

2014-11-12 Thread Robert Dodier
(-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

[sage-support] Re: Having issues with assume

2014-10-28 Thread Robert Dodier
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 subscribed to the Google Groups sage-support

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 john.crem...@gmail.com wrote: I am trying to use the Real Interval Field (RIF), which in principle does

[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 group. To unsubscribe

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

2014-10-16 Thread Robert Dodier
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 show a solution at w = t + t^2. best Robert Dodier

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

2014-10-15 Thread Robert Dodier
://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...@googlegroups.com. To post

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. Robert [1] https

[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

[sage-support] Re: Segfault on integration

2014-09-24 Thread Robert Dodier
On 2014-09-20, Kristoffer Ryhl-Johansen kristofferr...@gmail.com 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 message

[sage-support] Re: Segfault on integration

2014-09-20 Thread Robert Dodier
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 best Robert Dodier -- You

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

2014-08-29 Thread Robert Dodier
(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 number can be reported? For what it's worth, Robert

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

2014-08-29 Thread Robert Dodier
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 stop receiving emails from it, send an email to sage-support+unsubscr

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 john.crem...@gmail.com wrote: When William Stein first started the project it was an acronym SAGE for (I think) System for Algebra and Geometry Experimentation. But soon

[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: x2 and x1 x 2 sage: x2 or x1 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_formula.or_(x 2, x

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: x2 and x1 x 2 sage: x2 or x1 x 1 That's because and and or are program flow constructs in python, as they are in C (they have shortcut evaluation behaviour).

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: x2 and x1 x 2 sage: x2 or x1 x 1 The best way to manipulate logical combination of inequalities might be to use polyhedra. Looking at the documentation [1] I do

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,) Oops, mistake! This should be sage: Polyhedron(ieqs=[(2,-1

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 an

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 sayraf...@gmail.com wrote: Hi, In this page: http://www.sagemath.org/doc/constructions/linear_algebra.html#eigenvectors-and-eigenvalues in the eigenvectors and

[sage-support] Simplifying combinations of atomic inequalities

2014-07-16 Thread Robert Pollak
Hello! I see the following wrong results: sage: x2 and x1 x 2 sage: x2 or x1 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 unsubscribe from this group and stop

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-support group.

[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 you're

[sage-support] Newbie qepcad.py syntax problem

2014-07-10 Thread Robert Pollak
: '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 'Version B 1.50, 22 May 2008' of qepcad.) Best regards, Robert Pollak -- You received

[sage-support] Re: Won't Integrate

2014-07-05 Thread Robert Dodier
it is certainly 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, send an email to sage-support

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

2014-06-05 Thread Robert Dodier
as 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

[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

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 . ch4r...@hotmail.com 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

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

2014-05-29 Thread Robert Dodier
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: bug in comparison of function field elements

2014-04-30 Thread Robert Bradshaw
On Tue, Apr 29, 2014 at 10:57 PM, Robert Bradshaw rober...@math.washington.edu wrote: On Tue, Apr 29, 2014 at 9:07 AM, Volker Braun vbraun.n...@gmail.com wrote: On Tuesday, April 29, 2014 3:58:14 PM UTC+1, Simon King wrote: Yes there is! The hook is the hash function. CPython implementation

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

2014-04-29 Thread Robert Bradshaw
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 issues. - Robert -- You received this message because you are subscribed

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

2014-04-23 Thread Robert Dodier
)*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 receiving emails from it, send an email to sage-support+unsubscr...@googlegroups.com

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

2014-04-16 Thread Robert Dodier
, 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...@googlegroups.com. To post to this group, send

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 gtzana...@gmail.com wrote: On Fri, Mar 14, 2014 at 4:49 PM, Robert Bradshaw rober...@math.washington.edu wrote: Note that intL[i][introws[i]] + j %w == 0: would probably be just (or nearly) as fast as ((int(tupleL[i])[int(rows[i])])+j

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

2014-03-14 Thread Robert Bradshaw
Note that intL[i][introws[i]] + j %w == 0: would probably be just (or nearly) as fast as ((int(tupleL[i])[int(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,

[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 Groups

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

2014-01-28 Thread Robert Jacobson
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, January 28, 2014 8:29:16 PM UTC-5

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 jsku...@gmail.com 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

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
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. On 01/03/2014 02:48 AM, Robert

[sage-support] Re: simplification of symbolic expressions

2013-11-24 Thread Robert Dodier
; (%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 luck, have fun, hope this helps. Robert Dodier -- You received

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,

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

2013-08-21 Thread Robert Bradshaw
On Mon, Aug 19, 2013 at 6:30 PM, Dima Pasechnik dimp...@gmail.com wrote: On 2013-08-19, Vincent Knight knigh...@cf.ac.uk 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

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

[sage-support] Re: Integration of formal functions

2013-07-27 Thread Robert Dodier
), 'integrate(...)). '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 from this group

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

2013-07-22 Thread Robert Dodier
On 2013-07-18, Ed Scheinerman edward.scheiner...@gmail.com 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

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

2013-07-19 Thread Robert Dodier
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.com. To post to this group, send email to sage

[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

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

2013-04-23 Thread Robert Dodier
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.com. To post

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 juan...@gmail.com wrote: I have the

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

2013-03-17 Thread Robert Bradshaw
The syntax R.A,d = 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.A,d=QQ[] sage: R Multivariate Polynomial Ring in A, d over Rational Field On Sun, Mar 17, 2013 at 12:51 AM, Rolandb

Re: [sage-support] sage vs matlab

2013-03-05 Thread Robert Bradshaw
might not have much choice but if you're doing research than you can use whatever tools allow you to work/collaborate 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

[sage-support] Re: Bug in solve?

2013-02-24 Thread Robert Dodier
(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 encountered a Lisp error: BINDING-STACK

[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

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 lalwani.ak...@gmail.com 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

[sage-support] Re: Limit and variable

2013-02-01 Thread Robert Dodier
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 this message because you are subscribed to the Google Groups sage-support group

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

2013-01-23 Thread Robert Bradshaw
, 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 lfahlb...@gmail.com wrote: Hi - I would appreciate if someone could point me

[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

[sage-support] solving Diophantine equations in Sage

2012-12-08 Thread Robert Dodier
as 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 unsubscribe from

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

2012-12-08 Thread Robert Dodier
. Are there implementations 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@googlegroups.com

[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 subscribed to the Google

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.a = 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] algebraic number to integer conversion?

2012-10-17 Thread Robert Bradshaw
On Wed, Oct 17, 2012 at 9:37 PM, Eric Kangas eric.c.kan...@gmail.com 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 =

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 irv...@naml.us wrote: On Wed, Sep 26, 2012 at 10:42 PM, Robert Bradshaw rober...@math.washington.edu wrote: On Wed, Sep 26, 2012 at 8:54 PM, Geoffrey Irving irv...@naml.us wrote: On Wed, Sep 26, 2012 at 6:03 PM, Robert Bradshaw rober

Re: [sage-support] Re: Notebook server very slow. What are typical bottlenecks? I have ver5.2 and 1Gb RAM

2012-09-26 Thread Robert Bradshaw
I'd go for several 100M of RAM each, likely 0.5G to be comfortable, plus some memory for the OS and server itself. Throw 4G at it and it should behave much better, 8G and you should be good to go. Something like https://github.com/jasongrout/sage-forker would likely greatly reduce this

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

2012-09-26 Thread Robert Bradshaw
On Wed, Sep 26, 2012 at 4:28 PM, Geoffrey Irving irv...@naml.us wrote: Hello, I recently used sage to write a code generation script for exact geometric predicates: https://github.com/otherlab/simplicity Since it's a python script that imports sage, the simplicity script is GPL. Not

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

2012-09-26 Thread Robert Bradshaw
On Wed, Sep 26, 2012 at 8:54 PM, Geoffrey Irving irv...@naml.us wrote: On Wed, Sep 26, 2012 at 6:03 PM, Robert Bradshaw rober...@math.washington.edu wrote: On Wed, Sep 26, 2012 at 4:28 PM, Geoffrey Irving irv...@naml.us wrote: Hello, I recently used sage to write a code generation script

Re: [sage-support] randint - Maybe one bug !

2012-09-19 Thread Robert Bradshaw
backwards incompatible (considering the number of Python packages, including internally in our own libraryies, that understand and expect Python ints). - Robert On Wed, Sep 19, 2012 at 9:42 AM, Christophe BAL projet...@gmail.com wrote: What I think very confusing is that 1/4 is the Sage division

[sage-support] Multivariate polynomials: libsingular vs. polydict

2012-09-05 Thread Robert Samal
! Robert -- 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 unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com. Visit this group at http

[sage-support] Re: Linear programming in sage

2012-08-30 Thread Robert Samal
parameters are documented in GLPKBackend.solver_parameter which I didn't find. (I suppose one has to import it first.) Also I didn't find it in the sage tutorials describing MixedIntegerLinearProgram(). Anyway, thanks again! Robert -- You received this message because you are subscribed

  1   2   3   4   5   6   7   8   9   10   >