[sage-support] Question related to Sagemath permission on Github

2024-01-21 Thread Juan Grados
I have a question about how you deal with deletion policies in Sagemath. Specifically, I plan to create a repository with many people involved. This repo MUST have several admins (who can potentially delete the repo). To protect the repo from arbitrary deletion, I need to know if Github allows addi

[sage-support] Question about solving polynomial system over RR in sage

2021-08-22 Thread Juan Grados
I would like to solve a polynomial system over RR in sage. See the toy example below >>> var_names = ['xor_0_7_x0', 'xor_0_7_x12', 'xor_0_7_y0'] >>> nvars = 3 >>> R=PolynomialRing(RealField(53), nvars, var_names, order="degneglex") >>> a=R.gens()[0]+R.gens()[1]+R.gens()[2] >>>

Re: [sage-support] About Polyhedron

2021-02-16 Thread Juan Grados
w.lncc.br/) http://juaninf.blogspot.com - El lun, 8 feb 2021 a las 15:35, Juan Grados () escribió: > Thanks you to reply. The paper does not show the equations. > - > D.Sc. Juan del Carmen Grados Vásquez > Labor

Re: [sage-support] About Polyhedron

2021-02-08 Thread Juan Grados
inequality (-1, -1, -1, 3, 0, 3, -2) x + 2 >= 0 >>> > An inequality (0, 1, 0, 0, 0, 0, 0) x + 0 >= 0 >>> > An inequality (1, 0, 0, 0, 1, -1, -1) x + 1 >= 0 >>> > An inequality (0, 0, 0, 0, 0, 0, 1) x + 0 >= 0 >>> > An inequality (1,

Re: [sage-support] About Polyhedron

2021-02-07 Thread Juan Grados
Nacional de Computação Científica Tel: +55 21 97633 3228 (http://www.lncc.br/) http://juaninf.blogspot.com - El dom, 7 feb 2021 a las 22:51, Juan Grados () escribió: > Yes, but according to that paper it will be 65, and not 37.

Re: [sage-support] About Polyhedron

2021-02-07 Thread Juan Grados
> An inequality (1, 0, 0, 0, 0, 0, 0) x + 0 >= 0 > > An inequality (0, 0, 0, 0, 1, 0, 0) x + 0 >= 0 > > An inequality (0, 0, 0, 0, 0, 1, 0) x + 0 >= 0 > > An inequality (0, -1, 0, 1, 0, 1, -1) x + 1 >= 0 > > An inequality (-1, 0, 0, 1, 0, 1, -1) x + 1 >=

[sage-support] About Polyhedron

2021-02-07 Thread Juan Grados
Dear members, I am trying to reproduce page 9 of https://eprint.iacr.org/2016/407.pdf but until now is not possible to find the 65 inequalities that paper says. I am thinking that maybe this is because the version of SAGE I am using (this is 9). Do you think that there is any chance to obtain 65 i

[sage-support] Re: [sage-algebra] Improving the next code

2018-02-25 Thread Juan Grados
way to sum particular_soln + homogeneous_soln in a fast way? 2018-02-25 19:23 GMT-03:00 David Joyner : > On Sun, Feb 25, 2018 at 4:11 PM, Juan Grados wrote: > > How can I improve the time for the next code?. Basically, I want to > solve a > > large undetermined binary line

[sage-support] Improving the next code

2018-02-25 Thread Juan Grados
How can I improve the time for the next code?. Basically, I want to solve a large undetermined binary linear system and then I need to calculate its hamming weight. A = random_matrix(GF(2), 10, 12, density=0.55) b = random_vector(GF(2), 10) particular_soln = A.solve_right(b, check=False) A_right_k

[sage-support] Re: Fast solving of undetermined system over finite fields

2018-02-23 Thread Juan Grados
sorry, dimensions 1 x 1 El 23 de febrero de 2018, 5:14, Juan Grados escribió: > How can I solving a undetermined system with dimension 1 x 1000 over > finite fields (GF(2)) in SAGE in fast way. > > -- >

[sage-support] Fast solving of undetermined system over finite fields

2018-02-23 Thread Juan Grados
How can I solving a undetermined system with dimension 1 x 1000 over finite fields (GF(2)) in SAGE in fast way. -- - MSc. Juan del Carmen Grados Vásquez Laboratório Nacional de Computação Científica Tel: +55 21 97633 3228 (ht

Re: [sage-support] Re: Patterson Algorithm

2017-06-11 Thread Juan Grados
Panos, I had write patterson algorithm here http://juaninf.blogspot.com.br/2013/04/function-make-div-with-id-mycell-sage.html . 2017-03-07 9:47 GMT-03:00 Johan S. H. Rosenkilde : > Hi Panos, > > > I finally implement the decoder using lattice basis reduction (using LLL) > > I presume you mean F[

Re: [sage-support] Re: Matrix of operations after gaussian elimination

2017-06-10 Thread Juan Grados
More specifically, I need transform a partity check matrix into a systematic form, but I need the transform matrix 2017-06-10 7:06 GMT-03:00 Juan Grados : > I need an implementaton to perfom a Gaussian elimination on a matrix with > entries in GF(2). This implementation should re

Re: [sage-support] Re: Matrix of operations after gaussian elimination

2017-06-10 Thread Juan Grados
I need an implementaton to perfom a Gaussian elimination on a matrix with entries in GF(2). This implementation should return a transformation matrix and the result of the gassual elimination. 2017-06-10 6:51 GMT-03:00 Simon King : > Hi Juan, > > On 2017-06-10, Juan Grados wrote: >

Re: [sage-support] Re: Matrix of operations after gaussian elimination

2017-06-09 Thread Juan Grados
Thank by I get ValueError: too many values to unpack I think because my matrix has entries in GF(2). 2017-06-09 4:01 GMT-03:00 Kwankyu : > Is this what you want? > > sage: m = matrix(2,[1,2,3,4]) > sage: H,U = m.echelon_form(transformation=True) > sage: H > [1 0] > [0 2] > sage: U > [-2 1] > [

[sage-support] Matrix of operations after gaussian elimination

2017-06-08 Thread Juan Grados
I need the Matrix of operations after gaussian elimination in sage. How Can I get that? -- - MSc. Juan del Carmen Grados Vásquez Laboratório Nacional de Computação Científica Tel: +55 21 97633 3228 (http://www.lncc.br/) http://ju

[sage-support] Generating all solutions of a system over finite field in SAGE

2016-08-04 Thread Juan Grados
I need a operator, or function, to found all solutions of a linear system over the finite field F_2. For example in the next code I have generated a random linear system of the form Ax=b, where A=[[1,1,1], [0,1,1]] and b = [0,1], but when I used the operator \ I have only one solution [1,1,0], but

[sage-support] SAT Random Solutions

2016-07-25 Thread Juan Grados
Hello I am trying to get different random solution for every loop in my code sage from sage.sat.boolean_polynomials import solve as solve_sat # optional - cryptominisat sr = mq.SR(1,1,1,4,gf2=True,polybori=True) F,s = sr.polynomial_system() iter = 0 while(iter < 5): s

[sage-support] Degree of Regularity

2016-05-13 Thread Juan Grados
Dears, Is there any fucntion to calculate the Degree of Regularity of Quadratic Semi-Regular Determined Polynomial Systems in SAGE. (Looking for that, I found only for overdetermined systems the method degree_of_semi_regularity()) thanks by your attention -- ---

[sage-support] Question about Boolean Formulas

2015-12-09 Thread Juan Grados
Dears members, I have two question about the package sage.logic.booleval of SAGE. 1.-In the documentation of that package I have found a similar formula of this "t = ['&', ['|', 'a', 'b'], ['|', 'a', 'c']]". My question is it is posible to write three literals in each clause? 2.-In that documentat

[sage-support] Question about Casting over Finite Field

2015-08-09 Thread Juan Grados
I have declared the finite Field "K" and the equation "e". How I can cast the equation "e" to finite field "K".? K=GF(2) e = 0==r1405 + 97557948918122409/154603142862411371*r1406 - 26113332388244438/463809428587234113*r1408 + 180996540766820643/154603142862411371*r1409 thanks -- ---

[sage-support]

2015-05-06 Thread Juan Grados
I want to plot an histogram in SAGE from a list, my_data, but I want to plot all values of the axis X, in the below example this must be X=[0,1,2,3] but I get 2,4,6,8,10, Could you help me please? my_data = [0,1,2,10] BC=bar_chart(my_data, width=1, color="lightblue", figsize=4) sho

[sage-support] How convert from PolynomialSequence to DIMAC format

2015-01-28 Thread Juan Grados
Dears members, I have a object PolynomialSequence in ANF. PF = PolynomialSequence(total_polynomials, cipherSimon.RR) I need output this PF object in DIMAC format file. Could you help me please? -- - MSc. Juan del Carmen Grado

Re: [sage-support] Re: Question about subs in PolynomialSequence

2014-08-31 Thread Juan Grados
thanks! 2014-08-31 6:26 GMT-03:00 Dima Pasechnik : > On 2014-08-31, Juan Grados wrote: > > Dears members, > > > > I'm trying to use subs method in PolynomialSequence, but when I doing > that > > I get the same PolynomialSequence. For example in code below

[sage-support] Question about subs in PolynomialSequence

2014-08-30 Thread Juan Grados
Dears members, I'm trying to use subs method in PolynomialSequence, but when I doing that I get the same PolynomialSequence. For example in code below my polynomial sequence PS has only one polynomial (P) ... I'm trying to use subs but I get the same polynomial. Where is wrong? reset() from sage.

[sage-support] monomials of polynomial_zz_pex.Polynomial_ZZ_pEX

2014-08-29 Thread Juan Grados
Dears members, How I will be able to extract the monomials of univariate polynomial ring "P" of type polynomial_zz_pex.Polynomial_ZZ_pEX? K = GF(2^128,'t') PR. = PolynomialRing(K,"X") P = PR.random_element() -- - MSc. Juan del Ca

[sage-support] Question About Finite Field propertie in polynomial ring

2014-08-28 Thread Juan Grados
Dears members, Let be the field q = 2 K. = GF(q^n) and the Polynomial Ring PR = PolynomialRing(K,"X") Let be a random monomial of PR for example P = t*X^(q^a). Is there any method in sage to reduce X degree of polynomial P, such that equivalent polynomial is t*X^(q^b) where b = mod(P.degree(X)-1,n

Re: [sage-support] Question about solve_sat

2014-08-27 Thread Juan Grados
thanks 2014-08-27 14:05 GMT-03:00 Martin Albrecht : > It means that the system is believed to be unsolvable. > > On Wednesday 27 Aug 2014 13:54:47 Juan Grados wrote: > > Dear members, > > > > I'm trying to solve the attach formulas but I get FALSE, anyone know w

[sage-support] Question about solve_sat

2014-08-27 Thread Juan Grados
Dear members, I'm trying to solve the attach formulas but I get FALSE, anyone know what means that, or Why I get FALSE? sage: solve_sat(PolynomialSequence(F,RR),n=infinity) False thanks -- - MSc. Juan del Carmen Grados Vásquez

[sage-support] Re: GF(2n) to GF(2) ANF

2014-08-22 Thread Juan Grados
I'm see that substitution is very slow in SAGE for many many. Is there any way to get best perfomance to substitution? thanks 2014-08-22 7:42 GMT-03:00 Juan Grados : > I'm wrote a method to get the symbolic coefficients of polynomial p \in > GF(2^pp)[X0,...,X_{nvars-1}] transla

[sage-support] GF(2n) to GF(2) ANF

2014-08-22 Thread Juan Grados
I'm wrote a method to get the symbolic coefficients of polynomial p \in GF(2^pp)[X0,...,X_{nvars-1}] translate to GF(2). The method was wrotte using companion matrix [1] of irreducible polynomial GF(2^pp). For example for pp = 2. The polynomial p = X0*X1 = (x_00+x_01*t)*(x_10+x_11*t) = (x_00*x_10+x

[sage-support] Collect Coefficients MPolynomial_libsingular

2014-08-21 Thread Juan Grados
I'm want a method to get the symbolic coefficients of polynomial p \in GF(2^pp)[X0,...,X_{nvars-1}] translate to GF(2)[t]. For example for pp = 2. The polynomial p = X0*X1 = (x_00+x_01*t)*(x_10+x_11*t) = (x_00*x_10+x_01*x_11) + (x_01*x_11+x_01*x_10 + x_00*x_11)*t then the method get [(x_00*x_10+x_0

Re: [sage-support] Re: Several Substitution

2014-08-18 Thread Juan Grados
Dear Emmanuel, Thank, ... one last question ... How I will be able to extract the coefficients of t^0,t^1,...,t^(p-1) 2014-08-18 5:49 GMT-03:00 Emmanuel Charpentier < emanuel.charpent...@gmail.com>: > If I follow your cpde correctly, you are working on polynoms in X0,..,X3. > What you attempt to

[sage-support] Several Substitution

2014-08-17 Thread Juan Grados
Dears members I need doing substitution of values vars_gf2subs in vars_GF variables on p polynomial. I'm trying p.subs(vars_GF[i]=vars_gf2subs[i]) bu tI get error. How I will be able to doing that? p = 3 k. = GF(2^p) nvars = 4 vars_GF = [] vars_gf2 = [] vars_gf2subs = [] l = 0 for i in range(nvars

Re: [sage-support] Question about SAT solver

2014-08-13 Thread Juan Grados
nd the logic in solve() doesn't handle this case. Please open > a > ticket and CC me. > > On Tuesday 12 Aug 2014 09:01:33 Juan Grados wrote: > > Dears members, > > > > I'm trying to understand why the solver not found solution in this code > > >

[sage-support] Question about SAT solver

2014-08-12 Thread Juan Grados
Dears members, I'm trying to understand why the solver not found solution in this code R. = BooleanPolynomialRing() S = PolynomialSequence([x*y+z,x+y]) sol = S.solve(); sol [] For me the solution is x=1;y=1 and z=1, or I'm wrong? thanks -- -

[sage-support] Cryptominisat Error

2014-08-12 Thread Juan Grados
Dear members, I'm installed the cryptominisat solver using: ./sage -i cryptominisat-2.9.6 in Sage Version 6.2, Release Date: 2014-05-06 on Ubuntu 14.04 LTS but when I trying to execute the command from sage.sat.solvers.cryptominisat import CryptoMiniSat I get the error: Traceback (most recent c

[sage-support] Question about ortogonal subspaces

2014-07-15 Thread Juan Grados
Dears members, I trying to solve the next problem. Let be the chain of subspaces J1 \subset J2 \subset J3 \subset J4 over the finite field GF(3), where dim(J1) = 2, dim(J2) = 4, dim(J3) = 6 and dim(J4)=8. I want extract the basis vector of the subspace J4-J3, J3-J2 and J2-J1. For J4-J3 I get using

[sage-support] Question about Mpi4py

2014-06-04 Thread Juan Grados
As I can use the "Reduce" function to add 4 matrices copies (diff_table), which were declared using the "zeros" function?. Here my code in SAGE from mpi4py import MPI from numpy import* comm = MPI.COMM_WORLD rank = comm.rank rol = lambda val, r_bits, max_bits: \ (val <

Re: [sage-support] Re: pretty print question

2014-04-01 Thread Juan Grados
thanks but I get $$x= \left(\begin{array}{rr} 1 & 2 \\ 3 & 4 \end{array}\right) $$\newcommand{\Bold}[1]{\mathbf{#1}}\hbox{x=\left(\begin{array}{rr} 1 & 2 \\ 3 & 4 \end{array}\right)} the latex code and not pretty print 2014-04-01 17:59 GMT-03:00 Dominique Laurain : > x = var('x') > > A=ma

[sage-support] pretty print question

2014-04-01 Thread Juan Grados
Dears members, I want make a pretty_print in the follow code but I get false, and I want get the pretty expression of $$x=\left(\begin{array}{rr} 1 & 2 \\ 3 & 4 \end{array}\right)$$ var('x') >A=matrix([[1,2],[3,4]]) >pretty_print(x==A) How I will be able to fix that? -- -

[sage-support] Re: Factorial Carry Value in Python

2013-11-13 Thread Juan Grados
0<= u_i <= l-i 2013/11/13 Juan Grados > Let be s between 1 and l!-1 an integer value then s can expressed uniquely > than: > > s = u1*(l-1)! + u2*(l-2)!+ ... ul*0 > > Is there any function to find the values u1, u2, .

[sage-support] Factorial Carry Value in Python

2013-11-13 Thread Juan Grados
Let be s between 1 and l!-1 an integer value then s can expressed uniquely than: s = u1*(l-1)! + u2*(l-2)!+ ... ul*0 Is there any function to find the values u1, u2, ..., ul in SAGE or python? -- - MSc. Juan del Carmen Grados V

Re: [sage-support] Linear System Equations Polynomials

2013-10-24 Thread Juan Grados
1,p) 2013/10/23 Ivan Andrus > Do you mean something like: > > R. = PR.quotient(PR.ideal(p)) > # change your formulas to Xbar instead of X > A \ b > # ==> (a^3 + a, a^2, (a^3 + a^2)*Xbar^2 + (a + 1)*Xbar + a^3 + a) > > -Ivan > > On Oct 23, 2013, at 1:14 PM, Juan Grado

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

2013-10-23 Thread Juan Grados
+ 2*x + 4) > sage: b > (4*x^2 + x + 4, x^2 + 2*x + 4, 4*x^2 + 2*x + 3, x^2 + 2*x + 4) > > > On Wed, Oct 23, 2013 at 11:36 AM, Juan Grados wrote: > > thanks, and if the system coefficients are polynomials over finite > field? Is > > there PolynomialModRing(p(x))?

Re: [sage-support] Linear System Equations Polynomials

2013-10-23 Thread Juan Grados
ed matrix(PolynomialModRing(p),A)\vector(PolynomialModRing(p),b) but PolynomialModRing not exist ... 2013/10/23 John Cremona > On 23 October 2013 19:50, Juan Grados wrote: > > Is there in sage, any instruction to solve a linear system equations > > module p(x) (polynomial over finite

[sage-support] Linear System Equations Polynomials

2013-10-23 Thread Juan Grados
Is there in sage, any instruction to solve a linear system equations module p(x) (polynomial over finite field), where the system coefficients are polynomials over finite field?. I know that for integers exists, example (thanks Purkayastha) sage: I6 = IntegerModRing(6) sage: M = random_matrix(I6,

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

2013-10-23 Thread Juan Grados
0, 2, 1) > > > On 10/23/2013 09:11 PM, Juan Grados wrote: > >> Exists in python, any instruction to solve a linear system equations >> module n (integer). In mathematica for example If exists >> <http://mathematica.**stackexchange.com/questions/** >>

[sage-support]

2013-10-23 Thread Juan Grados
Exists in python, any instruction to solve a linear system equations module n (integer). In mathematica for example If exists . thanks -- -

Re: [sage-support] Re: Embed Cells

2013-04-19 Thread Juan Grados
Dears members, Now I get other problem When I trying m = 4; K_.<a> = GF(2); F.<a> = GF(2^m,modulus=a^4+a+1);#creado el campo donde vivirá el código I get ValueError: Modulus parameter not understood. 2013/4/19 Juan Grados <juan...@gmail.com> > Dear Jason, > >

Re: [sage-support] Re: Embed Cells

2013-04-19 Thread Juan Grados
Dear Jason, I get, the problem is that code $(function () { ... must be before head tag. Thanks by your attention 2013/4/19 Juan Grados > I am trying put this in my entry blog the complete code in the blogger > editor is below. The website is > http://juaninf.blogspot.com.b

Re: [sage-support] Re: Embed Cells

2013-04-19 Thread Juan Grados
" Hola() 2013/4/19 Jason Grout > On 4/19/13 12:07 PM, Juan Grados wrote: > >> Thanks by your reply, I trying, but yet not work. >> >> $(function () { >> // Make the div with id 'mycell' a Sage cell >> sagecell.makeSagecell({**inputL

Re: [sage-support] Re: Embed Cells

2013-04-19 Thread Juan Grados
x27;}); }); def Hola(): print "Hola" Hola() 2013/4/19 Jason Grout > On 4/18/13 11:33 PM, Juan Grados wrote: > >> Dears members, >> >> I like use the SAGE to embed any code (example below). But I want linked >> two cells using this >>

[sage-support] Embed Cells

2013-04-18 Thread Juan Grados
Dears members, I like use the SAGE to embed any code (example below). But I want linked two cells using this instruction but this not work. How I will be able to fix this?. For example I trying call from "compute" cel

Re: [sage-support] right_kernel()

2013-04-18 Thread Juan Grados
matrix in other form? H: is a parity check matrix. T: Toepliz V: Vandermonde D: Diagonal 2013/4/18 William Stein > On Thu, Apr 18, 2013 at 1:50 PM, Juan Grados wrote: > > Dears members, > > > > I executed right_kernel() for any matrix H. Why always I get a basis > ma

[sage-support] right_kernel()

2013-04-18 Thread Juan Grados
Dears members, I executed right_kernel() for any matrix H. Why always I get a basis matrix in this form [I|A] (I:identity)? -- - MSc. Juan del Carmen Grados Vásquez Laboratório Nacional de Computação Científica Tel: +55 24 2233-

Re: [sage-support] Solve Binary Linear System

2013-04-18 Thread Juan Grados
ge: %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 all matrix entries and all entrie of > vector b >

[sage-support] Solve Binary Linear System

2013-04-17 Thread Juan Grados
I have the equation Ax=b where all matrix entries and all entrie of vector b are in GF(2). How I will be able to solve this linear system equation over GF(2) in SAGE software? -- - MSc. Juan del Carmen Grados Vásquez Laboratório

Re: [sage-support] Fraction Fields coefficients

2012-07-01 Thread Juan Grados
Hi Robert, thanks by your response, but i get the same error AttributeError: 'FractionFieldElement_1poly_field' object has no attribute 'list' 2012/7/1 Robert Bradshaw > Perhaps your looking for something like > > sage: R. = QQ[] > sage: f = 2 * x^2 / (x^3 + 5); f > 2*x^2/(x^3 + 5) > sage: f.n

[sage-support] Question

2012-06-30 Thread Juan Grados
Dears members, How I will can get the coefficients list in Fraction Field. When I write p.lis() i obtained this message: 'FractionFieldElement_1poly_field' object has no attribute 'list' -- - MSc. Juan del Carmen Grados Vásquez

Re: [sage-support] Question about polynomials roots

2012-03-02 Thread Juan Grados
My question is for the support algorithm, for example: Chien, Berlekamp-Trace ... etc 2012/3/2 Martin Albrecht > Have you read the documentation? It's massive as I just found out: > > sage: F = GF(3) > sage: P. = F[] > sage: p = P.random_element() > sage: p.roots? > > Or check the source: > > sa

Re: [sage-support] Re: Question about plot

2012-02-24 Thread Juan Grados
I fix this, thanks by your time 2012/2/24 Juan Grados > Yes in console terminal, when i make this > > p1 = plot(spline(pattersonCpb), > (__builtin__.max(pattersonCpb)[0],__builtin__.min(pattersonCpb)[0]))+points(pattersonCpb) > show(p1) > the plot work fine, but when i make

Re: [sage-support] Re: Question about plot

2012-02-24 Thread Juan Grados
n 2/24/12 2:36 PM, Juan Grados wrote: > >> Jason Grout idea is good for me, but I need empiler several "plots", i >> after use show(plot1+plot2+...). How do this?. >> > > Exactly like you said. Make each plot: > > plot1=plot(spline(pts), (0,5)) > ... &g

Re: [sage-support] Question about plot

2012-02-24 Thread Juan Grados
Jason Grout idea is good for me, but I need empiler several "plots", i after use show(plot1+plot2+...). How do this?. 2012/2/24 Dan Drake > On Fri, 24 Feb 2012 at 01:39AM -0200, Juan Grados wrote: > > I need smooth line, (interpolation the points), line(vec) plot line > wi

Re: [sage-support] Question about plot

2012-02-23 Thread Juan Grados
I need smooth line, (interpolation the points), line(vec) plot line without smoth 2012/2/24 D. S. McNeil > > I am a two-tuple vector, "vet=[(1,2),(3,4),(5,6),..]", i want plot > > this data with a line aproximation, (interpolation this points), exist > > any parameter in list_plot function for t

Re: [sage-support] Question about perfomance

2012-02-19 Thread Juan Grados
William thanks for your answer, this information where I can read? ... I'm doing a good job and I would reference it 2012/2/19 William Stein > On Sun, Feb 19, 2012 at 11:24 AM, Juan Grados wrote: > > Dear members. > > > > is there a way to know how much is lost, in co

[sage-support] Question about perfomance

2012-02-19 Thread Juan Grados
Dear members. is there a way to know how much is lost, in computing time, using python instead of C/ansi to program an algorithm? -- - Juan del Carmen Grados Vásquez Laboratório Nacional de Computação Científica Tel: +55 24 2233-

Re: [sage-support] Re: Question about rdtsc()

2012-02-19 Thread Juan Grados
thanks 2012/2/19 Oleksandr Kazymyrov > Dear Juaninf, > > You can call "rdtsc" in SAGE (Accessing C Functions in Separate > Files) > or function "%time". > > Best regards, > Oleksandr > > -- > To post to this group, send email to sage-suppor

Re: [sage-support] Question about cycles per second in SAGE

2012-02-14 Thread Juan Grados
sorry, the measure is cycles/byte. sorry again, how obtain this measure in SAGE? 2012/2/14 David Kirkby > On 14 February 2012 03:40, juaninf wrote: > >> Hi everybody. >> >> Exists any command in SAGE to get the spent cpu-cycles per second in >> SAGE? > > > What do you mean by "spent cycles per

Re: [sage-support] Re: Question about download SAGE

2012-02-14 Thread Juan Grados
thanks! 2012/2/14 emil > > > On Feb 14, 6:57 pm, Juan Grados wrote: > > I dont permissions for this in my server > > Well, if there is no better solution (like asking sysadmin, or > somebody uploads other package ) you could build lzma (or successor xz-

Re: [sage-support] Re: Question about download SAGE

2012-02-14 Thread Juan Grados
I dont permissions for this in my server 2012/2/14 emil > > > On Feb 14, 5:21 pm, juaninf wrote: > > Hi everybody > > > > Please I want download the SAGE math for fedora-32bits, but i need the > > file with extension *.tar not tar.lzmaa, because my tar program no > > have support fot this

Re: [sage-support] Re: Question about minimal polynomial to build GF

2011-12-11 Thread Juan Grados
... thanks, by your attention, I am reading a documentation to Implement Galois Field in other place, ... and declare in C this static unsigned prim_poly[MAX_EXT_DEG + 1] = {◦1 , ◦3 , /∗ not used ∗/ ◦7 , ◦13 , ◦23 , ◦45 , ◦103 , ◦203 , ◦435 , ◦1041 , ◦2011 , ◦4005 , ◦10123 , ◦20033 , ◦42103 , ◦10

Re: [sage-support] Re: Question about minimal polynomial to build GF

2011-12-11 Thread Juan Grados
but when I put this polynomial x^8+x^7+x^4+x^3+x+1 I get ValueError: finite field modulus must be irreducible but it is not 2011/12/11 achrzesz > > > On Dec 12, 12:10 am, juaninf wrote: > > Hi everybody > > > > I want choose different minimal polynomial to build a Galois Field > > 2^m, how? >

Re: [sage-support] Questión About factor() function

2011-11-04 Thread Juan Grados
for all, factors ... 2011/11/5 juaninf > Dears. > > In this code > > F. = GF(2^(4)) > PR = PolynomialRing(F,'X') > M = vector(PR,8) > #get a g1 > g1 = (X + 1) * (X + a) * (X + a^3) * (X + a^5) * (X + a^5 + a^4 + a^3 > + a) * (X + a^6 + a^5 + a^4 + a^2) * (X + a^7 + a^2 + a + 1) * (X + > a^7 + a^

Re: [sage-support] Re: Question about MPI in SAGE

2011-10-17 Thread Juan Grados
Dear Maarten sorry 2011/10/16 Maarten Derickx > Dear Juan Grados, > > I don't think people have forgotten it, but merely not answering since you > are asking in the wrong place. This is the sage mailinglist not the mpi4py one > and the problems you are experiencing are p

Re: [sage-support] Question about sagenb.org

2011-10-12 Thread Juan Grados
help me please 2011/10/12 Juan Grados > yes I trying but ... > > > http://groups.google.com/group/sage-support/browse_thread/thread/b581f46aad988472# > > > 2011/10/12 Robert Bradshaw > >> No. Also note that sagenb.org is on a single computer shared among >&g

Re: [sage-support] Question about sagenb.org

2011-10-12 Thread Juan Grados
yes I trying but ... http://groups.google.com/group/sage-support/browse_thread/thread/b581f46aad988472# 2011/10/12 Robert Bradshaw > No. Also note that sagenb.org is on a single computer shared among > tens of thousands of users, so if you want to do something > computationally expensive you're

Re: [sage-support] Re: Question about MPI in SAGE

2011-10-10 Thread Juan Grados
please dont forget 2011/10/5 Juan Grados > I maked this thread > http://groups.google.com/group/mpi4py/browse_thread/thread/2e08c06f271c3069 > now I get this error > > > ImportError: > /home/juaninf/Escritorio/juaninf/Programas/sage/local/lib/python2.6/site-packages/mpi4p

Re: [sage-support] Question about legend place in plot

2011-10-10 Thread Juan Grados
thanks 2011/10/10 D. S. McNeil > > I want that legend is in left side, ... how? > > Try something like: > > p = plot(sin,legend_label="sin") > p.set_legend_options(loc='upper left') > p.show() > > After making some plot p, you can type help(p.set_legend_options) to > see more information about t

Re: [sage-support] Re: Question about Convert Polynomial

2011-10-09 Thread Juan Grados
thanks, fix my problem 2011/10/9 Maarten Derickx > In sage the types are not that important (it is just an implementation > detail), but the parent of an object is since that says what the > mathematical meaning of an object is. Try: > > q=p.parent()(q) > > -- > To post to this group, send email

Re: [sage-support] Re: Question about MPI in SAGE

2011-10-05 Thread Juan Grados
I maked this thread http://groups.google.com/group/mpi4py/browse_thread/thread/2e08c06f271c3069 now I get this error ImportError: /home/juaninf/Escritorio/juaninf/Programas/sage/local/lib/python2.6/site-packages/mpi4py/MPI.so: undefined symbol: _PyByteArray_empty_string 2011/10/3 Juan Grados

Re: [sage-support] Re: Question about Plot

2011-10-04 Thread Juan Grados
thanks, Jason, John, I've already done with your help 2011/10/4 Jason Grout > On 10/4/11 5:28 PM, juaninf wrote: > >> Hi everyone >> >> I am have 4 vectors, and I want plot 2 graphics in the same place, >> >> vector1 vs vector2 color = red >> vector3 vs vector4 color = yelow >> > > You might fin

Re: [sage-support] Re: Question about MPI in SAGE

2011-10-03 Thread Juan Grados
... 2011/10/2 Juan Grados > My openmpi > > openmpi 1.4.1-3 > > > 2011/10/2 Volker Braun > >> On Sunday, October 2, 2011 4:25:54 PM UTC+2, juaninf wrote: >>> >>> My MPI version is 0.6, >> >> >> I guess thats not OpenMPI, then? >

Re: [sage-support] Re: Question about MPI in SAGE

2011-10-02 Thread Juan Grados
My openmpi openmpi 1.4.1-3 2011/10/2 Volker Braun > On Sunday, October 2, 2011 4:25:54 PM UTC+2, juaninf wrote: >> >> My MPI version is 0.6, > > > I guess thats not OpenMPI, then? > > >> -- > To post to this group, send email to sage-support@googlegroups.com > To unsubscribe from this group,

Re: [sage-support] Re: Question about MPI in SAGE

2011-10-02 Thread Juan Grados
My MPI version is 0.6, Ubuntu version Ubuntu 11.04 Sage version 4.7.1 2011/10/2 Volker Braun > I'd guess lt_dlexit is a symbol that libtool defines that has something to > do with getting out of a dynamic library. Seems like the mpi4py package > misses some linker flags. You should post some in

[sage-support] Question about MPI in SAGE

2011-10-01 Thread Juan Grados
Hi everyone, I want use MPI, but when I put this line from mpi4py import MPI I get this error Traceback (most recent call last): File "", line 1, in File "_sage_input_5.py", line 10, in exec compile(u'open("___code___.py","w").write("# -*- coding: utf-8 -*-\\n" + _support_.preparse_w

Re: [sage-support] Question about xgcd method

2011-10-01 Thread Juan Grados
fixed 2011/9/30 juaninf > Hi everybody, > > I want implement a modified extend Euclidean Algorithm, (egcd > function), but this give wrong results, below my egcd, please help me > to fix ... > def egcd(p1,p2): >if p2 == PR(0): >return (p1,1,0) >else: >(q1, r1) = (p1).quo_

Re: [sage-support] Re: Question about congruence

2011-09-30 Thread Juan Grados
I found this method in polynomial_element.pyx file thanks 2011/9/30 Juan Grados > Hi David, > > Yes already, but I want build other new, because I need make modification > in EEA (add other stop parade), but I first that my "EEA Algorithm" work > fine ... > &g

Re: [sage-support] Re: Question about congruence

2011-09-30 Thread Juan Grados
Hi David, Yes already, but I want build other new, because I need make modification in EEA (add other stop parade), but I first that my "EEA Algorithm" work fine ... thanks 2011/9/30 David Joyner > On Fri, Sep 30, 2011 at 5:11 AM, juaninf wrote: > > Hi everybody, > > > > I want implement a

Re: [sage-support] Question about Patterson Algorithm Implementation

2011-09-28 Thread Juan Grados
I have already sent, but I dont answer ... because I expect please only if anyelse can help me iff a time ... 2011/9/28 David Joyner > On Wed, Sep 28, 2011 at 6:12 PM, Juan Grados wrote: > > Hi David, > > > > Yes I understand, but now I think that have a logic problem in

Re: [sage-support] Question about Patterson Algorithm Implementation

2011-09-28 Thread Juan Grados
Hi David, Yes I understand, but now I think that have a logic problem in algorithm, but I don't know where ... i "copying lines" from [Ict2011], ... 2011/9/28 David Joyner > On Wed, Sep 28, 2011 at 5:58 PM, Juan Grados wrote: > > help please! > > > They did s

Re: [sage-support] Question about Patterson Algorithm Implementation

2011-09-28 Thread Juan Grados
help please! 2011/9/28 Juan Grados > in the end line > > print sigma.roots(), > > always give empty vector, here sigma.roots() should nonzero vector > > 2011/9/28 Juan Grados > >> Hi thanks for your answers, >> >> I used _inverter_, _mul_, _add_ etc

Re: [sage-support] Question about Patterson Algorithm Implementation

2011-09-28 Thread Juan Grados
in the end line print sigma.roots(), always give empty vector, here sigma.roots() should nonzero vector 2011/9/28 Juan Grados > Hi thanks for your answers, > > I used _inverter_, _mul_, _add_ etc, because apparently > the implementation work fine but only "apparently&quo

Re: [sage-support] Question about Patterson Algorithm Implementation

2011-09-28 Thread Juan Grados
Hi thanks for your answers, I used _inverter_, _mul_, _add_ etc, because apparently the implementation work fine but only "apparently", i think that the essencial problem is with _invert_ method, but now I used inverse_mod , but I dont where are the error, I implemented Berlekamp Algorithm too, fr

[sage-support] Questión about primitive element

2011-09-28 Thread Juan Grados
Hi How I will get a primite element ... F = GF(2) PRF. = PolynomialRing(F); print PRF Phi = PRF.quotient(z^4+z+1); Phi.primitive_element() . ? -- - Juan del Carmen Grados Vásquez Laboratório Nacional de Computação