On 27.05.19 20:40, Tracy Hall wrote:
> If anyone is curious to see an example
> [...]
> sage: for ii in range(1):
> : alarm(10)
> : try:
> : if prod(P.gens())^2 in anideal:
> : cancel_alarm()
> [...]
> : except AlarmInterrupt:
> :
On 27.05.19 18:08, Jeroen Demeyer wrote:
> On 2019-05-27 16:11, Tracy Hall wrote:
>> in particular it seems to have undefined
>> behavior in subsequent calls after it has been interrupted by alarm().
>
> This is simply a fact of life and not really considered to be a bug. You
> should not rely on
On 20.05.19 05:29, jianrong wrote:
> I am trying to write a function in sage to permute a list of elements.
An alternative might me something along the following:
|
Permutation("(1,2,3)(4,5)").action(['a', 'b', 'c', 'd', 'e'])
|
--
You received this message because you are subscribed to the
On 04.05.19 06:04, Andrew wrote:
> The `Permutation` function is more general. For example, the folllowing
> all work: [...]
Thank you. But sorry, this does not answer my question. Maybe I should
be more precise:
What is the idea behind
|
sage: Permutation([5,4,3,2,1]).parent()
Standard
sage: Permutation([5,4,3,2,1]).parent()
Standard permutations
sage: Permutation('(1,5)(2,4)(3)').parent()
Standard permutations of 5
Why is the first not restricted to permutations of length 5?
--
You received this message because you are subscribed to the Google Groups
"sage-support" group.
On 27.02.19 16:12, Daniel Krenn wrote:
> On 27.02.19 15:35, Dima Pasechnik wrote:
>> So you get your normal vectors in the subspace parallel to the affine hull
>> of P.
>
> Thank you, looks easy :) (I am now just using orthogonal=True in my
> case, as I do not want to
On 27.02.19 15:35, Dima Pasechnik wrote:
>> Yes, this is the interesting case. The problem then is going back from
>> the projection. I guess that orthogonality is ususally destroyed here...
> One can ensure it is orthonormal:
>
> sage: P = polytopes.simplex(2)
> sage:
>
On 27.02.19 14:34, Dima Pasechnik wrote:
> On Wed, Feb 27, 2019 at 1:02 PM Daniel Krenn wrote:
>> Is there an easy way in SageMath to compute the in- or outward surface
>> normal vector of these faces of P? (in contrast to doing it all from
>> scratch). If not, are there me
Say we have
sage: P = polytopes.simplex(2)
sage: P.faces(1)
(<0,1>, <0,2>, <1,2>)
Is there an easy way in SageMath to compute the in- or outward surface
normal vector of these faces of P? (in contrast to doing it all from
scratch). If not, are there methods that might help, so that not
The H-representation consists of equations and inequalities and the
equations seem to be in some canonical form. Is there a method that
returns the non-free variables (or indices), i.e. that are the variables
completely determined by the equations meaning once a value for the
other variables is
On 2019-02-08 10:07, Daniel Krenn wrote:
> Let I be an ideal. Then I might want to compute something involving
> Groebner basis, e.g. computing I.variety().
> Now suppose one wants to select a particular algorithm for the
> computation of the Groebner basis. Then (due to caching) I u
On 2018-12-05 11:07, Dima Pasechnik wrote:
For integrating a polynomial over a polyhedron LattE is used but if the
dimension is not full, then it is not implemented, see
sage: x, y = polygens(QQ, 'x, y')
sage: P = Polyhedron(vertices=[[0,0],[1,1]])
On 2019-02-08 10:33, Simon King wrote:
> On 2019-02-08, Daniel Krenn wrote:
>> Let I be an ideal. Then I might want to compute something involving
>> Groebner basis, e.g. computing I.variety().
>> Now suppose one wants to select a particular algorithm for the
>> comput
Let I be an ideal. Then I might want to compute something involving
Groebner basis, e.g. computing I.variety().
Now suppose one wants to select a particular algorithm for the
computation of the Groebner basis. Then (due to caching) I use something
along the lines of
GB =
On 2018-12-04 15:36, Dima Pasechnik wrote:
> On Tue, Dec 4, 2018 at 1:47 PM Daniel Krenn wrote:
>> [I might not be that familiar with integration over polyhedra, but
>> shouldn't that basically be a somehow "nice" transformation where some
>> Jacobi-determinant
On 2018-11-13 06:20, Michael Beeson wrote:
> def demo():
> var('N,x')
> test = ((N*(3*I*sqrt(3) + 9) + N*(3*I*sqrt(3) - 3)))*x
> print("test = ")
> print(test)
> print("test.full_simplify() = ")
> print(test.full_simplify())
>
> Here is the output
>
> sage: demo()
>
> test =
>
>
On 2018-10-31 12:21, deSitter wrote:
> My text terminal uses a black background. The default text is dark blue
> on jet black. This is not a winning start to 8.4. How do I change the
> text color. Seriously, this is the sort of little thing that is
> absolutely exasperating, to get a perfect build
Is there an algebraic structure for the group of roots of unity in
SageMath?
To clearify: I think of something whose elements are all roots of unity
or all n-th root of unity for a given n and the group operation is
multiplication.
Best, Daniel
--
You received this message because you are
On 2017-10-22 17:27, Jianrong Li wrote:
> Let $r=x_{1}^{4} + 2 \, x_{1}^{3} x_{2} + 4 \, x_{1}^{2} x_{2}^{2} + 2
> \, x_{1}
> x_{2}^{3} + x_{2}^{4} + 2 \, x_{1}^{3} x_{3} + 2 \, x_{2}^{3} x_{3} + 4
> \, x_{1}^{2} x_{3}^{2} + 4 \, x_{2}^{2} x_{3}^{2} + 2 \, x_{1} x_{3}^{3}
> + 2 \, x_{2} x_{3}^{3}
On 2017-08-12 12:10, Daniel Krenn wrote:
> On 2017-08-12 11:14, Volker Braun wrote:
>> Caches from @cached_method are pickled by default, but the class can
>> opt-out of this. Sometimes this is necessary to make pickling work (not
>> every Cython object is pickle-able).
&
On 2017-08-12 11:14, Volker Braun wrote:
> Caches from @cached_method are pickled by default, but the class can
> opt-out of this. Sometimes this is necessary to make pickling work (not
> every Cython object is pickle-able).
Hmmmshouldn't the below work then?
sage: class A(SageObject):
:
I have an ideal and computed its Groebner basis (6 hours). Now I want to
store this result, which is not the problem itself, but when loading
this data, it seems not possible to reconstruct the ideal (meaning
create it with the original generators and "bootstrap" the groebner
basis found earlier).
On 2017-07-20 16:34, Vincent Delecroix wrote:
> sage: algdep(RealField(20)(1.4142), 2)
> x^2 - 2
>
> Note that you have to specify the maximal degree. Be careful about the
> input
>
> sage: algdep(1.4142, 2)
> 5000*x - 7071
>
> This is using the PARI/GP algdep command (whose algorithm is a LLL
Dear all,
I vaguely remember a discussion on this list some years ago about the
following: Given a floating point number, what is a good guess for the
simplest possible algebraic number that approximates this?
E.g. The input 1.4142 gives back sqrt(2) or x^2-2 and says that this is
approximates
On 2017-04-28 22:11, William Stein wrote:
> In SMC in a **Sage worksheet**, do
>
> show(g)
Thanks. (I've tried g.show() in the SMC, and then one gets only the
usual static representation).
Best
Daniel
--
You received this message because you are subscribed to the Google Groups
Is there a way to click on a vertex of a SageMath-graph (in some
notebook, Jupyter or SMC-notebook) and drag it around, thus manually
positioning the vertices?
(I think there was a graph_editor in the old Sage Notebook which could
do this, but this does not seem to work in one of the new
Shouldn't the following two both work?
sage: '%.2f' % (pi.n(),)
'3.14'
sage: '{:.2f}'.format(pi.n())
---
ValueErrorTraceback (most recent call last)
in ()
> 1 '{:.2f}'.format(pi.n())
On 2017-01-07 11:09, Vincent Delecroix wrote:
> Le 05/01/2017 à 11:33, Daniel Krenn a écrit :
>> It seems like the embedding is simply ignored in some sense:
>> sage: QQbar(QuadraticField(-1, 'I', embedding=-CC.gen()).gen())
>> I
>
> To my mind, this one should a
On 2017-01-05 11:29, Daniel Krenn wrote:
> On 2017-01-05 10:55, Daniel Krenn wrote:
>> as there is no *canonical* coercion as no embedding of the number field
>> is specified.
>>
>> How can I specify this embedding such that it is used e.g. for the
>> sym
On 2017-01-05 10:55, Daniel Krenn wrote:
> as there is no *canonical* coercion as no embedding of the number field
> is specified.
>
> How can I specify this embedding such that it is used e.g. for the
> symbolic I?
This looks weird: I is defined in sage.libs.pynac.
I want to take an number field element embedded in the symbolic ring
like the imaginary I and add it to an algebraic number:
sage: I.pyobject() + QQbar(sqrt(2))
This results in
TypeError: unsupported operand parent(s) for '+': 'Number Field in I
with defining polynomial x^2 + 1' and
On 2016-12-13 08:06, Ralf Stephan wrote:
> On Monday, December 12, 2016 at 10:55:32 AM UTC+1, Daniel Krenn wrote:
> and the source in src/sage/symbolic/integration/integral.py
> you must walk the expression tree and apply subs to the first operand of
> all instances of
>
Let f be a sum of symbolic expressions. One of the summands is
sage: a = ((z^3 - 10*z^2 + 17*z - 8)/(z^4 + z^3 + z^2 + z + 1)).integrate(z)
sage: a
integrate((z^3 - 10*z^2 + 17*z - 8)/(z^4 + z^3 + z^2 + z + 1), z)
which cannot be integrated. However, doing
sage: a.subs({z^4 + z^3 + z^2 + z + 1:
On 2016-12-07 10:00, Dima Pasechnik wrote:
> I don't see how using git worktree would reduce recompilation, for the
> merging workflow you recommend touches a lot of files, potentially.
In your SageMath root, assuming you are at branch develop:
git worktree add ../merge
cd ../merge
git trac
On 2016-12-07 08:58, Dima Pasechnik wrote:
> For the purposes of reviewing you do not care about history being
> clean in your local branches, and what I describe does not lead to
> much recompilation.
For reviewing (without any new commits) this works.
> Whereas in particular with old tickets,
On 2016-12-06 22:02, Justin C. Walker wrote:
> I am starting in a new, empty directory, and 'git' seems to want a repository
> specified.
>
> I have a "global" .gitconfig file set up.
>
> A couple of questions:
>
> Should I check out the 'develop' branch first, and then incorporate (how?)
>
Pressing TAB after a dot on some object gives (in the recent 7.4.beta1)
this:
sage: M = Matrix([1])
sage: M.
M.act_on_polynomial M.anticommutator
M.add_multiple_of_column M.antitranspose
M.add_multiple_of_row M.apply_map >
The tex-lines
\documentclass{tikzposter}
\usepackage{sagetex}
\begin{document}
Bla.
\end{document}
fail with
! Undefined control sequence.
\@enddocumenthook ... \jobname .tex'}
\fi \ST@wsf
{_st_.endofdoc()}\@ifundef...
l.5 \end{document}
Any ideas?
Best
On 2016-04-01 21:31, Emmanuel Charpentier wrote:
> * revert to by original patch against 7.1rc0, update it and push this
> to Trac for review ? Or
> * push by jumbo patch ?
If the merge of 7.2.beta2 into your 7.1.rc0+whatever works without a
conflict, then there is no need to do this
On 2016-02-17 17:25, Dima Pasechnik wrote:
> it's not related to RIF, it seems. Compilation ends with
>
> ImportError: cannot import name ZZ
>
> and even the following does not work:
>
> import sage.all
> from sage.all import *
> from sage.rings.integer_ring import ZZ
> print ZZ(1000)
>
>
On 2016-02-17 10:16, Dima Pasechnik wrote:
>
>
> On Wednesday, February 17, 2016 at 7:26:18 AM UTC, Daniel Krenn wrote:
>
> On 2016-02-17 08:05, Jeroen Demeyer wrote:
> > On 2016-02-17 07:33, Daniel Krenn wrote:
> >> Calling this in my working d
On 2016-02-17 08:05, Jeroen Demeyer wrote:
> On 2016-02-17 07:33, Daniel Krenn wrote:
>> Calling this in my working directory with
>>sage problem.spyx
>> fails with
>>Traceback (most recent call last):
>>...
>>ImportError: cannot import name
My problem boils down to:
problem.spyx
from sage.rings.real_interval_field import RealIntervalField # fails
RealIntervalField(100)(4.2) # not needed in minimal non-working example
Calling this in my working directory with
sage problem.spyx
fails with
Traceback
On 2015-11-09 08:06, P Purkayastha wrote:
> [...] so that everything gets coerced to Sage
> integers, irrespective of the input:
>
> Q = Integer(q)
FYI, this is a conversion; use
Q = ZZ.coerce(q)
instead.
Daniel
--
You received this message because you are subscribed to the Google Groups
Am 2015-04-14 um 21:34 schrieb Volker Braun:
Inline plots work in Sage now, you just need a sufficently new version.
In 6.6
show(x+1)
still shows the result externally, not inline.
Daniel
--
You received this message because you are subscribed to the Google Groups
sage-support group.
To
Can someone explain me the difference between SR.symbol and SR.var ?
(BTW: SR.symbol does not have a description)
Daniel
--
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
Am 2015-04-10 um 23:12 schrieb Michael Orlitzky:
On 04/10/2015 05:09 PM, Daniel Krenn wrote:
Can someone explain me the difference between SR.symbol and SR.var ?
(BTW: SR.symbol does not have a description)
SR.symbol is faster, but only lets you declare one variable at a time.
SR.var
Am 2015-04-11 um 00:03 schrieb Volker Braun:
Common confusion. See also: http://trac.sagemath.org/ticket/17958
Thanks.
Daniel
--
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,
Am 2015-02-22 um 19:31 schrieb Volker Braun:
I think thats fixed in pyopenssl-0.14. I opened a ticket
at http://trac.sagemath.org/ticket/17831
Easiest solution is probably pip install pyopenssl
Seems to work now.
Thanks.
Daniel
On Saturday, February 21, 2015 at 12:47:31 PM UTC+1, John
Am 2014-10-27 um 19:08 schrieb mjs:
With the recent release of Firefox 33 (now shipped to Fedora, coming
soon to an OS near you), I can no longer connect to my Sage 6.3 server.
The error I get is:
An error occurred during a connection to sage.math.clemson.edu:34567.
The key does not
Am 2014-09-09 um 05:03 schrieb Miguel Yorro:
Hi! I'm using Craig Finch Sage Beginner's Guide to learn to use Sage.
However, when I ran the following code (both in the notebook and
terminal mode)
|
sage:var('x')
x
sage:sinc(x)=sin(x)/x
sage:plot(sinc,(x,-10,10))
Try one of the following:
Am 2014-08-29 um 21:25 schrieb Daniel Krenn:
I want to solve polynomial equations and in order to do so, I do
something like:
sage: R.x,y = PolynomialRing(QQ, order='lex')
sage: I = R.ideal([x*y-1, x^2-y^2])
sage: I.groebner_basis()
[x - y^3, y^4 - 1]
Meanwhile, I found, which seems to do
I want to solve polynomial equations and in order to do so, I do
something like:
sage: R.x,y = PolynomialRing(QQ, order='lex')
sage: I = R.ideal([x*y-1, x^2-y^2])
sage: I.groebner_basis()
[x - y^3, y^4 - 1]
Now I have to take the equation with only one variable, find the
solutions for it (over
Am 2014-08-13 um 08:18 schrieb Dima Pasechnik:
On 2014-08-12, Nils Bruin nbr...@sfu.ca wrote:
--=_Part_5256_754630738.1407883007654
Content-Type: text/plain; charset=UTF-8
On Tuesday, August 12, 2014 7:29:01 AM UTC-7, Daniel Krenn wrote:
It doesn't look like the results above help
Am 2014-08-12 um 17:42 schrieb Nils Bruin:
On Tuesday, August 12, 2014 7:29:01 AM UTC-7, Daniel Krenn wrote:
(type 'sage.symbolic.expression.Expression', 28257)
This could be a large number (although expressions are recursive data
structures, so one complicated expression can cause
Am 2014-08-11 um 16:19 schrieb Nils Bruin:
On Monday, August 11, 2014 4:46:45 AM UTC-7, Daniel Krenn wrote:
As expected. If you want to get some idea of what is taking memory at
the sage side you could do something like
import gc
from collections import Counter
gc.collect()
pre={id(c
Am 2014-08-10 um 15:30 schrieb Dima Pasechnik:
On 2014-08-08, Nils Bruin nbr...@sfu.ca wrote:
--=_Part_203_2052776100.1407511132909
Content-Type: text/plain; charset=UTF-8
Would at least restarting ecl release the memory?
(it's not clear how the latter
can be done; would
Am 2014-08-11 um 09:24 schrieb Nils Bruin:
On Sunday, August 10, 2014 11:16:19 PM UTC-7, Daniel Krenn wrote:
In my case, this (maxima_calculus.reset()) gives the following:
[...] I wasn't able to
reproduce the segfault. Would you have a concise code fragment that can
reliably
Am 2014-08-10 um 17:38 schrieb Nils Bruin:
On Sunday, August 10, 2014 6:31:03 AM UTC-7, Dima Pasechnik wrote:
(it's not clear how the latter can be done; would
maxima_calculus.reset() do the job?)
Tried again (since I had troubles before; see other posting), but memory
is not released.
Am 2014-08-08 um 17:18 schrieb Nils Bruin:
On Friday, August 8, 2014 3:02:03 AM UTC-7, Dima Pasechnik wrote:
so doing
sage: import sage.libs.ecl
sage: sage.libs.ecl.ecl_eval((ext:set-limit 'ext:heap-size 0))
might solve your problem.
Many thanks, this solved my original problem; I can
I'm tying to do a length computation, which unfortunately terminates by
TypeError: ECL says: Memory limit reached. Please jump to an outer
pointer, quit program and enlarge the
memory limits before executing the program again.
How can do this enlarging?
FWI: What I have are a lot of (not too
How do I solve the symbolic equation x^2 == 16*sin(x) symbolically? I
don't even get 0 as solution.
Daniel
--
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
On Thursday, 5 December 2013 09:22:09 UTC+1, Daniel Krenn wrote:
How do I solve the symbolic equation x^2 == 16*sin(x) symbolically? I
don't even get 0 as solution.
Daniel
--
You received this message because you are subscribed to the Google
Groups sage-support group.
To unsubscribe
Am 2013-10-19 04:41, schrieb meaning:
I have a problem when calling Mathematica 9 in Sage 5.12.
I installed Mathematica 9, Sage 5.12 and claim their direction in
.bash_profile file.
The direction of Mathematica is: /usr/local/bin/mathematica
The direction of Sage is:
Am 2013-08-22 01:12, schrieb 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.
Ok, thanks for the clearification.
Your check might be the
bottleneck,
I need an array of Elements of RealIntervalField and I want to iterate
(a lot of times) through it. How can I do that in a fast way in Sage's
cython (i.e. in the notebook or in a .spyx-file)?
An (minimal) example what I basically want is given below (PS).
I tried the following things (and a lot
Am 2013-08-01 01:09, schrieb Oscar Lazo:
Hello dear sage users!
I would like to view a 3d plot from a specific viewpoint using tachyon.
For instance
save: var('s')
sage: spiral=parametric_plot3d((cos(s),sin(s),s/8),
(s,-16*pi,16*pi),aspect_ratio=1,plot_points=200)
sage: show(spiral,
The worksheet [1] contains the command 'pi.n?'. Nothing else is in the
worksheet. But Sage hangs up evaluating that line. Creating a new worksheet
and inserting 'pi.n?' there does not lead to that behavior. So what is
wrong with [1]?
Daniel
PS: I didn't create a trac-tricket for it, since it
Is there something to solve recursions (e.g. linear recursions, but also
others) in Sage? Or, formulated in another way: Is there something in Sage
like RSolve in Mathematica?
--
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to
Am Donnerstag, 7. Juni 2012 22:54:15 UTC+2 schrieb luisfe:
What is the propeer way to rotate an object without suffering zooming
effect? This is related to the more general question. Is there a way to
deal with the camera using tachyon viewer? For instance, it seems that for
implicit3d
This is now #12967 http://trac.sagemath.org/sage_trac/ticket/12967.
Am Freitag, 11. Mai 2012 10:58:53 UTC+2 schrieb Robert Samal:
Hi!
By some random experiments I discovered the following weirdness:
sage: bool(piInfinity)
False
sage: bool(piInfinity)
True
So far it seems that pi
This is now #12968 http://trac.sagemath.org/sage_trac/ticket/12968. (I
couldn't find a ticket for it, but the problem was documented in the source
code (see other answers to that posting). Does someone know that there is a
ticket for it already?)
Am Sonntag, 13. Mai 2012 18:50:11 UTC+2 schrieb
Am Freitag, 18. Mai 2012 21:16:08 UTC+2 schrieb kcrisman:
On Friday, May 18, 2012 2:49:49 PM UTC-4, arshpreet singh wrote:
If you just append the right directory to your PATH in a .bashrc or
.profile
that should work. I have
export PATH=$PATH:'/Applications/MathApps/'
(which
Am 2012-04-12 17:03, schrieb Jan Groenewald:
Hi
On 12 April 2012 16:36, P Purkayastha ppu...@gmail.com
mailto:ppu...@gmail.com wrote:
Did you plot it and see the graph? Maybe this is not the graph you
wanted? The following command will plot the graph for you (it looks
The label of an edge of an graph is positioned (by default) at the midpoint
of an edge (the center of the label is equal to that midpoint). That
means, that the text lays directly over the edge, so it cannot be read that
good.
Is there a way to position the label in another way? E.g. position
Am 2012-03-21 22:30, schrieb Starx:
Check out sage.sets.family.Family, it might be what you're looking for.
Many thanks, that is really what I'm looking for.
Daniel
--
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to
Is there something like an infinite dictionary in Sage/Python?
More precisely, is there something where
- i can put in values like in a dictionary,
- but maybe also a function which tells me how to map a key to a value,
- and maybe also something that maps a key to a (finite) set of keys and
Am 2012-03-21 22:15, schrieb Simon King:
It seems unclear to me what you realy want to do.
On the one hand, you say that you want to provide a function f to the
infinite dictionary D, such that D[key]==f(key) for all keys - but
then, what would you gain compared with directly calling f?
How do I change the default precision used? E.g. I want to enter {{{a
= 1.2}}} and want that a is an element of RealField(100) without
explicitly telling to use that field each time.
What does not work is:
sage: RR = RealField(100)
sage: a = 1.2
sage: a.parent()
Real Field with 53 bits of
79 matches
Mail list logo