[sage-support] Re: performance of applying a map to an integer vector

2013-11-20 Thread Felix Breuer
n using machine precision arithmetic (be it via Cython or NumPy), is there a way to tell whether at some point throughout a long computation there were numerical overflows - without paying a huge speed penalty for making this check? Thanks, Felix On Wednesday, November 20, 2013 8:00:41 PM UTC+1, Nil

[sage-support] Re: performance of applying a map to an integer vector

2013-11-20 Thread Felix Breuer
y big integers are very few and far between, but they do exist? Thanks, Felix On Wednesday, November 20, 2013 7:27:11 PM UTC+1, Felix Breuer wrote: > > Hi Simon, hi Maarten, > > thank you for your answers! > > @Simon: > > Yes, I figured that calling vector(...) would imp

[sage-support] Re: performance of applying a map to an integer vector

2013-11-20 Thread Felix Breuer
this fits to Cython's types? Cheers, Felix On Wednesday, November 20, 2013 7:10:06 PM UTC+1, Maarten Derickx wrote: > > It seems like your code is mostly doing easy integer operations in tight > for loops (although in your case the for loop is hidden in v.apply_map). If > you car

[sage-support] performance of applying a map to an integer vector

2013-11-20 Thread Felix Breuer
apply_map). (Do something like return vector((foo(vi,d) for vi in v)) instead of apply_map does not make much a difference.) Now, my questions are: 1) Why is this so slow? Am I missing something here? 2) Is there anything I can do to improve performance? Thank you very much, Felix -- You

[sage-support] Re: Polyhedron construction fails on specific list of vertices

2013-04-23 Thread Felix Wellen
With "backend='ppl', base_ring=QQ" it works! Thanks again ;) Am Dienstag, 23. April 2013 14:10:26 UTC+2 schrieb Volker Braun: > > This is now http://trac.sagemath.org/14479 (needs review) > > > -- You received this message because you are subscribed to the Google Groups "sage-support" group. T

[sage-support] Polyhedron construction fails on specific list of vertices

2013-04-23 Thread Felix Wellen
Hi, while trying to render the intersection of a 4-cube and a 3-cube (I guess that really doesn't matter), I stopped on trouble with the following: sage: point_list = [[0.132, -1.028, 0.028],[0.5, 0.5, -1.5],[-0.5, 1.5, -0.5],[0.5, 0.5, 0.5],[1.5, -0.5, -0.5],[-0.332, -0.332, -0.668],[-1.332,

[sage-support] Re: computing the cup-product in the cohomology of a simplicial complex

2012-01-04 Thread Felix Breuer
I'd be happy to do some more work on it. I would need some pointers, though, on what still needs to be done. I can write up some documentation next week. Cheers, Felix PS: I am at the Joint Meetings right now. Are you at JMM? Or are there any Sage-related events at JMM? -- To post to

[sage-support] Re: computing the cup-product in the cohomology of a simplicial complex

2011-12-05 Thread Felix Breuer
re of first :) Cheers, Felix -- 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 For more options, visit this group at http://groups.google.com/group/sage-support URL: http://www.sagemath.org

[sage-support] Re: computing the cup-product in the cohomology of a simplicial complex

2011-12-04 Thread Felix Breuer
help! Felix -- 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 For more options, visit this group at http://groups.google.com/group/sage-support URL: http://www.sagemath.org

[sage-support] Re: computing the cup-product in the cohomology of a simplicial complex

2011-12-04 Thread Felix Breuer
Hi John! Thanks for your detailed answer! I will try to figure this out and get back to you! Felix -- 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 For more options, visit this group

[sage-support] computing the cup-product in the cohomology of a simplicial complex

2011-12-03 Thread Felix Breuer
d to have? Sorry, if these questions are somewhat fuzzy, but I am new to sage. Thank you for any help you can give! Best, Felix -- 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 F

[sage-support] Re: OpenGL, WebGL and Sage

2011-05-03 Thread Felix Lawrence
On May 3, 1:00 pm, Jason Grout wrote: > At one point, some people were experimenting also with using GPU to > speed up linear algebra, but I don't think any code made it into Sage. There are trac tickets for this at #10009 and #10010 for OpenCL and CUDA respectively. They use Andreas Kloeckner's

[sage-support] Re: Scipy acts differently for *i and j

2010-11-05 Thread Felix Lawrence
ignoring the I. Not as convenient as it could be, but at least it doesn't give you an incorrect answer! Cheers, Felix On Nov 4, 11:57 pm, kcrisman wrote: > On Nov 3, 10:41 pm, Felix Lawrence wrote: > > > Here's why there's a discrepancy: > > sage: import

[sage-support] Re: Scipy acts differently for *i and j

2010-11-03 Thread Felix Lawrence
nverters, and you send it this one: lambda s: complex(s.replace('i','j')) Does anyone know if this is fixed in newer versions of numpy/scipy? Cheers, Felix On Nov 3, 2:43 am, kcrisman wrote: > Dear support, > > We get some inconsistencies with using I and j in sc

[sage-support] Re: building fortran extensions on Mac OS X

2010-06-01 Thread Felix Lawrence
else to review the patch before it gets included in a new version of sage. Cheers, Felix On Jun 1, 11:45 pm, Chris Kees wrote: > Hi, > > I'm trying to build a fortran extension module with f2py using and spkg I've > built for sage 4.4.2 on Mac OS X 10.6. When I just run &

Re: [sage-support] Re: [sage-edu] GR package

2010-02-05 Thread Felix Lenders
Hi, you could also check out reduce[1] with the redten[2] package, which worked very well and fast(!) for me, not comparable to the Maple tensor package. Regards, Felix [1] http://www.reduce-algebra.com/ [2] http://www.utsc.utoronto.ca/~harper/redten.html -- To post to this group, send email

Re: [sage-support] bug in numerical approximation?

2010-01-24 Thread Felix Lenders
>0) sage: g = integrate(t^(x-1)*exp(-t),(t,0,infinity)) sage: g(3).n() --- TypeError Traceback (most recent call last) /home/felix/ in () /home/felix/Temp/sage-4.3/local/lib/python2.6/site-packages/sage/symbolic/expression.so in sa

[sage-support] Re: f2py on mac

2010-01-20 Thread Felix Lawrence
It turned out that there was another bug in there. At least on my machine, f2py was still broken on 4.3.1.rc1. I tracked down the problem (mac gcc doesn't support -shared) and there's a patch ready for review at ticket #8010. Cheers, Felix On Jan 15, 10:42 am, gsw wrote: > >

[sage-support] Re: How to list (and remove) loaded objects in a sage session

2010-01-17 Thread Felix Lawrence
more sensible not to display these omnipresent identifiers? It's hard to find my own variables in the mess! (I'm on 10.6, core 2 duo, running 4.3, if this is a bug) Cheers, Felix -- 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 For more options, visit this group at http://groups.google.com/group/sage-support URL: http://www.sagemath.org

[sage-support] f2py on mac

2010-01-13 Thread Felix Lawrence
t 1)), it gives me a much longer error message, which I have posted at http://www.physics.usyd.edu.au/~felix/sage/sage_fortran_error.txt Is there a way to get sage to prefer /usr/bin/gcc to /usr/local/bin/ gcc that is easier than temporarily editing my path globally? Does sage's fortran normally work on mac

[sage-support] sage-mode

2010-01-08 Thread Felix Lenders
ugh space for the next command as well. What can I do about this? Thank you very much for your help! Felix -- 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 For more options, visit this

[sage-support] what is wrong?

2009-12-10 Thread felix
Hello, When I try: var('t') r=vector((t^2,5*t,t^2-16)) v=derivative(r(t),t) norm_v=v.norm() print r print v pos=r(0) print pos vel=v(0) i receive an error for the last command vel(0)...i suppose i'm doing something wrong see error below Thank you for help... FG Traceback (most recent ca

[sage-support] lattice polytopes and sage

2009-12-10 Thread Felix Breuer
Sage. However barvinok has not been packaged yet and I could not find out whether other such software, e.g. LattE, has already been included in Sage. Thanks, Felix -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+uns

[sage-support] symbolic linear algebra: beginner questions

2009-11-24 Thread Felix Breuer
],[c2])') M = maxima('matrix(a-c,b-c)') M = maxima('matrix(transpose(a-c),transpose(b-c))') In the last two lines I get errors about the arguments being "invalid rows". Any suggestions? Thanks, Felix -- To post to this group, send email to sage-support@goo

[sage-support] Re: Weave broken in Sage 4.1 Mac intel 64 bit

2009-10-22 Thread felix
Hi, just came back to check for any progress in this thread. Tried your patch and it seems to work fine for me. Thanks a lot, Felix --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email

[sage-support] Re: Sage 4.1.1 VMWare image

2009-10-11 Thread Felix Breuer
different host. Correct? Thanks, Felix On 11 Okt., 21:20, William Stein wrote: > On Sun, Oct 11, 2009 at 12:03 PM, Felix Breuer wrote: > > > Hello! > > > I am trying to run the VMWare image of Sage 4.1.1. I downloaded and > > unzipped the image and ran it with the vmw

[sage-support] Sage 4.1.1 VMWare image

2009-10-11 Thread Felix Breuer
ce is that running notebook does not give any error message, it just returns to the prompt immediately. Any suggestions? Thanks, Felix --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send emai

[sage-support] Re: Weave broken in Sage 4.1 Mac intel 64 bit

2009-09-10 Thread felix
Hi, I just got a trac account and create the ticket at http://trac.sagemath.org/sage_trac/ticket/6914 . By the Way: I verified, that the problem is also present with the 32 bit mac intel binary. - Felix --~--~-~--~~~---~--~~ To post to this group, send email to

[sage-support] Re: weave corrupt

2009-09-09 Thread felix
Actually, the code from your last post works perfectly under sage 4.0.1 (OSX intel). [39mcreating /var/folders/hF/hFYXAHyRF680EWYF08TbAE+++TI/-Tmp-/felix/ python25_interme\ diate/compiler_56e514baa7dba73b5ded18f9a64c0373 [0m Time: CPU 0.29 s, Wall: 6.09 s 179997 Time: CPU 0.01 s, Wall: 0.01

[sage-support] Re: weave corrupt

2009-09-07 Thread felix
Hi, What about: import scipy.weave from scipy.weave import converters ? - Felix On 4 Sep., 21:56, Patrick Hammer wrote: > Hi! > > import weave > from weave import converters > > gives: > > Traceback (most recent call last): >   File "", line 1, in >

[sage-support] Re: Weave broken in Sage 4.1 Mac intel 64 bit

2009-09-07 Thread felix
g to sage 4.0.1 when I intend to use weave (I only use it to test code fragments before I build proper extensions) which is currently good enough for me to get my work done. - Felix On 25 Aug., 09:09, Robert Bradshaw wrote: > I wonder if this broke in the upgrade to 2.6. The errors below mean

[sage-support] Weave broken in Sage 4.1 Mac intel 64 bit

2009-08-20 Thread felix
r Mac OS X 10.5.8 on my MacBook Pro Unibody rev. A. It used to work fine before. Weave doesn't even try to compile anything anymore it just produces a KeyError: >>> from scipy import weave >>> weave.inline("""printf("hello");""")

[sage-support] solve()

2008-03-24 Thread felix
Hello, I have a system of 5 equations: x=2*a+1 x=3*b+2 x=4*c+3 x=5*d+4 x=6*e+5 a,b,c,d,e,x are integers variables. I would solve the system with solve(). 1)How to declare that solve have to assume x,a,b,c,d,e as integers? The answer is k*60+59 k integer. Can Sage with the solve function r

[sage-support] restrict the solutions of solve() to ZZ

2008-03-24 Thread felix
Hello, A very simple question: I need to resolve a system of equations with the solve(..) function. How to declare that : 1) the variables belong to N (positives integers ) 2) filter the answers of the resolved system to N. Thank you. Felix