sage: import numpy
sage: numpy.roots?
The values in the rank-1 array p are coefficients of a
polynomial.
If the length of p is n+1 then the polynomial is
p[0] * x**n + p[1] * x**(n-1) + ... + p[n-1]*x + p[n]
sage: a=numpy.array([1,0,1],dtype=float)
sage: numpy.ro
On 8/5/07, Robert Miller <[EMAIL PROTECTED]> wrote:
>
> The following handles factoring over RR:
>
> elif is_RealField(R):
> n = pari.set_real_precision(int(3.5*R.prec()) + 1)
> G = list(self._pari_('x').factor())
>
> Would this work for RDF too? Is there a faster w
The following handles factoring over RR:
elif is_RealField(R):
n = pari.set_real_precision(int(3.5*R.prec()) + 1)
G = list(self._pari_('x').factor())
Would this work for RDF too? Is there a faster way?
On Aug 4, 12:06 pm, "William Stein" <[EMAIL PROTECTED]> wrote
On 8/4/07, Robert Miller <[EMAIL PROTECTED]> wrote:
>
> Here is the full traceback.
>
> sage: M = random_matrix(RDF, 4, 4)
> sage: M.eigenspaces()
> ---
>Traceback (most recent call
> last)
> /Volumes/HOME/robert/sage-2.7.
Here is the full traceback.
sage: M = random_matrix(RDF, 4, 4)
sage: M.eigenspaces()
---
Traceback (most recent call
last)
/Volumes/HOME/robert/sage-2.7.2/ in ()
/Volumes/HOME/robert/sage-2.7.2/matrix2.pyx in
matrix2.Matri
On 8/3/07, Joshua Kantor <[EMAIL PROTECTED]> wrote:
> The eigen routine uses numpy. If you do
> p,e=m.eigen()
>
> then p is an array of eigenvalues and e is a matrix whose columns are
> eigenvectors.
> The documentation works fine in my local install not sure why it fails
> on sage.math
This is p
The eigen routine uses numpy. If you do
p,e=m.eigen()
then p is an array of eigenvalues and e is a matrix whose columns are
eigenvectors.
The documentation works fine in my local install not sure why it fails
on sage.math
Josh
On Aug 3, 2:43 pm, Robert Bradshaw <[EMAIL PROTECTED]>
wrote:
>
Linear algebra over RR is generic and, in many cases, really bad so
far. Use RDF or even scipy unless you really need arbitrary
precision. (If the RDF interface is inconsistent, fix it and send a
patch :-)
- Robert
On Aug 1, 2007, at 6:28 AM, Robert Miller wrote:
>
> Hi everyone,
>
> I ha
