In Sage, the following code produces error messages ending with "TypeError: cannot evaluate symbolic expression numerically":
from mpmath import *
var('a')
print sqrt(a)
If I remove the line "from mpmath import *", it prints "sqrt(a)" as expected.
How can I get mpmath to work with sage variables?
To complement @vdelecroix answer, the main issue is that, when you write:
from mpmath import *
you import hundreds of functions into the namespace, overwriting the ones that are imported by Sage, this includes sqrt.
If you want to use both Sage's symbolics and mpmath numerics, you have to be more parsimonious and only import what you need. For example, if the functions you want to import from mpmath are acos and atan, you can write:
from mpmath import acos, atan
Also, perhaps a more long-term solution, if there is something that works better with mpmath than with the Sage's default (or that seems absent from Sage), you should advertise it so that we use mpmath for such a feature in Sage.
As a complement: If you ant to use both mpmath's sqrt and Sage's one, you have three options: Either you use from mpmath import sqrt as mpsqrt and call mpsqrt, or you use import mpmath and call mpmath.sqrt, or you use import mpmath as mp and call mp.sqrt. (The as ... is the way to define an alias, so you can use anything in place of mpsqrt or mp in my examples.) In all cases, sqrt still refers to Sage's function.
mpmath is a library intended to compute functions to high precision. If you try to use it on a symbolic variable it will fail. The input must be numeric (and Sage numbers work well)
sage: import mpmath
sage: mpmath.sqrt(5.3)
mpf('2.3021728866442674')
sage: mpmath.sqrt(pi)
mpf('1.7724538509055159')
I'm attempting to do a mix of symbolic and high-precision-numeric computations in a loop. It appears (from my example above) that merely importing mpmath precludes doing any symbolic computations involving sqrt anywhere in that session.
Specifically, it appears that the "ordinary" sqrt works with symbolic variables, but that the mpmath sqrt does not. Is there a workaround?
Yes: do not use from mpmath import * but say explicitely what you need to import.
Asked: 2017-12-03 00:10:43 +0100
Seen: 680 times
Last updated: Dec 04 '17
Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.
