Re: [sage-support] Re: persistent homology?

2021-05-30 Thread Guillaume Rousseau
I have not yet benchmarked it but am working on it, however it is definitely slower in comparison to Ripser / gudhi and co. (although if I understand correctly Ripser only computes homology for Vietoris-Rips complexes?), since it is a very direct python implementation of Zomorodian's algo with

Re: [sage-support] Re: persistent homology?

2021-05-26 Thread Guillaume Rousseau
Hello all, I know I am a bit late to the party, but I have written a Sage module for computing persistent homology, and I just submitted a ticket ! -- You received this message because you are subscribed to the Google Groups "sage-support" group. To unsubscribe from this group and stop receivin

[sage-support] Bug report : function Matrix.set_block() may not work correctly with extension fields

2016-05-23 Thread Guillaume Holler
Hello, I don't know if this is the right place to submit bugs. I have consistently reproduced the following bug: using set_block when matrix space is on a given NumberField or an extension field (at least when there are multiple generators) works as expected the first time (in a sage sessio

Re: [sage-support] symbolic bitwise

2015-04-14 Thread Guillaume CONNAN
-BEGIN PGP SIGNED MESSAGE- Hash: SHA256 Le 10/04/2015 09:42, Jori Mantysalo a écrit : > On Wed, 8 Apr 2015, Guillaume CONNAN wrote: > >>>> p1 = lambda z: ((3*z^3 | 2*z^2) ^^ 5*z) + 2 >>> >>> What is the supposed result of >>> >>> 3

Re: [sage-support] symbolic bitwise

2015-04-08 Thread Guillaume CONNAN
-BEGIN PGP SIGNED MESSAGE- Hash: SHA256 Le 08/04/2015 09:26, Jori Mantysalo a écrit : > On Sat, 4 Apr 2015, Guillaume CONNAN wrote: > >> p1 = lambda z: ((3*z^3 | 2*z^2) ^^ 5*z) + 2 > > What is the supposed result of > > 3*z^3 | 2*z^2 p1(2) is ((24 | 8) ^^ 10) +

Re: [sage-support] symbolic bitwise

2015-04-04 Thread Guillaume CONNAN
-BEGIN PGP SIGNED MESSAGE- Hash: SHA256 Le 03/04/2015 19:19, Jori Mantysalo a écrit : > On Fri, 3 Apr 2015, Guillaume CONNAN wrote: > >> it's to do something like ( P1 * P2 | a ) & P3 + 2 > > Can you give a concrete example of "P1 * P2 | a"? somethi

Re: [sage-support] symbolic bitwise

2015-04-03 Thread Guillaume CONNAN
-BEGIN PGP SIGNED MESSAGE- Hash: SHA256 Le 03/04/2015 11:23, Jori Mantysalo a écrit : > On Thu, 2 Apr 2015, Guillaume CONNAN wrote: > >> is there a library to do some symbolic calculations over integer >> but with bitwise operations ? I need to work on symbolic >

[sage-support] symbolic bitwise

2015-04-02 Thread Guillaume CONNAN
-BEGIN PGP SIGNED MESSAGE- Hash: SHA256 Hello, is there a library to do some symbolic calculations over integer but with bitwise operations ? I need to work on symbolic polynomials mixing arithmetic + and * and bitwise >> << & | ^ . Thanks, - -- Guillaume Conna

[sage-support] Re: bug desolve ?

2010-05-08 Thread Guillaume
OK, but what about the assume(k>0) which is not assumed ? sage: x=var('x') sage: f=function('f',x) sage: k=var('k') sage: assume(k>0) sage: desolve(diff(f(x),x,2)/f(x)==k,[f,x]) TypeError: Computation failed since Maxima requested additional constraints (try the command 'assume(k>0)' before i

[sage-support] bug desolve ?

2010-05-07 Thread Guillaume
Hello, I want to solve f''/f=k with k in R sage: x=var('x') sage: f=function('f',x) sage: k=var('k') sage: desolve(diff(f(x),x,2)/f(x)==k^2+1,[f,x]) k1*e^(I*sqrt(-k^2 - 1)*x) + k2*e^(-I*sqrt(-k^2 - 1)*x) I put k^2+1 since sage would'nt solve even with assume(k>0) -- To post to this group, se

[sage-support] edo

2010-05-07 Thread Guillaume
mmand 'assume(k>0)' before integral or limit evaluation, for example): Is k positive, negative, or zero? sage: -- Guillaume -- 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: convert trigonometric/hyperbolic functions to exponentials

2010-03-12 Thread Guillaume
and weirder... sage: a*sinh(log(a)).substitute_function(sinh,sh).simplify_full() 1/2*a*log(a) - 1/2 sage: t*sinh(log(t)).substitute_function(sinh,sh).simplify_full() 1/2*t^2 - 1/2 this 1/2*a*log(a) - 1/2 is a bit unexpected... well, I'm discovering sage but these results look strange. che