Assumptions¶
The GenericDeclaration
class provides assumptions about a symbol or
function in verbal form. Such assumptions can be made using the assume()
function in this module, which also can take any relation of symbolic
expressions as argument. Use forget()
to clear all assumptions.
Creating a variable with a specific domain is equivalent with making an
assumption about it.
There is only rudimentary support for consistency and satisfiability checking in Sage. Assumptions are used both in Maxima and Pynac to support or refine some computations. In the following we show how to make and query assumptions. Please see the respective modules for more practical examples.
In addition to the global assumptions()
database, assuming()
creates reusable, stackable context managers allowing for temporary
updates of the database for evaluation of a (block of) statements.
EXAMPLES:
The default domain of a symbolic variable is the complex plane:
sage: var('x')
x
sage: x.is_real()
False
sage: assume(x,'real')
sage: x.is_real()
True
sage: forget()
sage: x.is_real()
False
Here is the list of acceptable features:
sage: ", ".join(map(str, maxima("features")._sage_()))
'integer, noninteger, even, odd, rational, irrational, real, imaginary,
complex, analytic, increasing, decreasing, oddfun, evenfun, posfun,
constant, commutative, lassociative, rassociative, symmetric,
antisymmetric, integervalued'
Set positive domain using a relation:
sage: assume(x>0)
sage: x.is_positive()
True
sage: x.is_real()
True
sage: assumptions()
[x > 0]
Assumptions also affect operations that do not use Maxima:
sage: forget()
sage: assume(x, 'even')
sage: assume(x, 'real')
sage: (-1)^x
1
sage: (-gamma(pi))^x
gamma(pi)^x
sage: binomial(2*x, x).is_integer()
True
Assumptions are added and in some cases checked for consistency:
sage: assume(x>0)
sage: assume(x<0)
Traceback (most recent call last):
...
ValueError: Assumption is inconsistent
sage: forget()
- class sage.symbolic.assumptions.GenericDeclaration(var, assumption)¶
Bases:
sage.structure.unique_representation.UniqueRepresentation
This class represents generic assumptions, such as a variable being an integer or a function being increasing. It passes such information to Maxima’s declare (wrapped in a context so it is able to forget) and to Pynac.
INPUT:
var
– the variable about which assumptions are being madeassumption
– a string containing a Maxima feature, either user defined or in the list given bymaxima('features')
EXAMPLES:
sage: from sage.symbolic.assumptions import GenericDeclaration sage: decl = GenericDeclaration(x, 'integer') sage: decl.assume() sage: sin(x*pi) 0 sage: decl.forget() sage: sin(x*pi) sin(pi*x) sage: sin(x*pi).simplify() sin(pi*x)
Here is the list of acceptable features:
sage: ", ".join(map(str, maxima("features")._sage_())) 'integer, noninteger, even, odd, rational, irrational, real, imaginary, complex, analytic, increasing, decreasing, oddfun, evenfun, posfun, constant, commutative, lassociative, rassociative, symmetric, antisymmetric, integervalued'
Test unique representation behavior:
sage: GenericDeclaration(x, 'integer') is GenericDeclaration(SR.var("x"), 'integer') True
- assume()¶
Make this assumption.
- contradicts(soln)¶
Return
True
if this assumption is violated by the given variable assignment(s).INPUT:
soln
– Either a dictionary with variables as keys or a symbolic relation with a variable on the left hand side.
EXAMPLES:
sage: from sage.symbolic.assumptions import GenericDeclaration sage: GenericDeclaration(x, 'integer').contradicts(x==4) False sage: GenericDeclaration(x, 'integer').contradicts(x==4.0) False sage: GenericDeclaration(x, 'integer').contradicts(x==4.5) True sage: GenericDeclaration(x, 'integer').contradicts(x==sqrt(17)) True sage: GenericDeclaration(x, 'noninteger').contradicts(x==sqrt(17)) False sage: GenericDeclaration(x, 'noninteger').contradicts(x==17) True sage: GenericDeclaration(x, 'even').contradicts(x==3) True sage: GenericDeclaration(x, 'complex').contradicts(x==3) False sage: GenericDeclaration(x, 'imaginary').contradicts(x==3) True sage: GenericDeclaration(x, 'imaginary').contradicts(x==I) False sage: var('y,z') (y, z) sage: GenericDeclaration(x, 'imaginary').contradicts(x==y+z) False sage: GenericDeclaration(x, 'rational').contradicts(y==pi) False sage: GenericDeclaration(x, 'rational').contradicts(x==pi) True sage: GenericDeclaration(x, 'irrational').contradicts(x!=pi) False sage: GenericDeclaration(x, 'rational').contradicts({x: pi, y: pi}) True sage: GenericDeclaration(x, 'rational').contradicts({z: pi, y: pi}) False
- forget()¶
Forget this assumption.
- has(arg)¶
Check if this assumption contains the argument
arg
.EXAMPLES:
sage: from sage.symbolic.assumptions import GenericDeclaration as GDecl sage: var('y') y sage: d = GDecl(x, 'integer') sage: d.has(x) True sage: d.has(y) False
- sage.symbolic.assumptions.assume(*args)¶
Make the given assumptions.
INPUT:
*args
– a variable-length sequence of assumptions, each consisting of:any number of symbolic inequalities, like
0 < x, x < 1
a subsequence of variable names, followed by some property that should be assumed for those variables; for example,
x, y, z, 'integer'
would assume that each ofx
,y
, andz
are integer variables, andx, 'odd'
would assume thatx
is odd (as opposed to even).
The two types can be combined, but a symbolic inequality cannot appear in the middle of a list of variables.
OUTPUT:
If everything goes as planned, there is no output.
If you assume something that is not one of the two forms above, then an
AttributeError
is raised as we try to call itsassume
method.If you make inconsistent assumptions (for example, that
x
is both even and odd), then aValueError
is raised.Warning
Do not use Python’s chained comparison notation in assumptions. Python literally translates the expression
0 < x < 1
to(0 < x) and (x < 1)
, but the value ofbool(0 < x)
isFalse
whenx
is a symbolic variable. Therefore, by the definition of Python’s logical “and” operator, the entire expression is equal to0 < x
.EXAMPLES:
Assumptions are typically used to ensure certain relations are evaluated as true that are not true in general.
Here, we verify that for \(x>0\), \(\sqrt{x^2}=x\):
sage: assume(x > 0) sage: bool(sqrt(x^2) == x) True
This will be assumed in the current Sage session until forgotten:
sage: bool(sqrt(x^2) == x) True sage: forget() sage: bool(sqrt(x^2) == x) False
Another major use case is in taking certain integrals and limits where the answers may depend on some sign condition:
sage: var('x, n') (x, n) sage: assume(n+1>0) sage: integral(x^n,x) x^(n + 1)/(n + 1) sage: forget()
sage: var('q, a, k') (q, a, k) sage: assume(q > 1) sage: sum(a*q^k, k, 0, oo) Traceback (most recent call last): ... ValueError: Sum is divergent. sage: forget() sage: assume(abs(q) < 1) sage: sum(a*q^k, k, 0, oo) -a/(q - 1) sage: forget()
An integer constraint:
sage: n,P,r,r2 = SR.var('n, P, r, r2') sage: assume(n, 'integer') sage: c = P*e^(r*n) sage: d = P*(1+r2)^n sage: solve(c==d,r2) [r2 == e^r - 1] sage: forget()
Simplifying certain well-known identities works as well:
sage: n = SR.var('n') sage: assume(n, 'integer') sage: sin(n*pi) 0 sage: forget() sage: sin(n*pi).simplify() sin(pi*n)
Instead of using chained comparison notation, each relationship should be passed as a separate assumption:
sage: x = SR.var('x') sage: assume(0 < x, x < 1) # instead of assume(0 < x < 1) sage: assumptions() [0 < x, x < 1] sage: forget()
If you make inconsistent or meaningless assumptions, Sage will let you know:
sage: assume(x<0) sage: assume(x>0) Traceback (most recent call last): ... ValueError: Assumption is inconsistent sage: assume(x<1) Traceback (most recent call last): ... ValueError: Assumption is redundant sage: assumptions() [x < 0] sage: forget() sage: assume(x,'even') sage: assume(x,'odd') Traceback (most recent call last): ... ValueError: Assumption is inconsistent sage: forget()
You can also use assumptions to evaluate simple truth values:
sage: x, y, z = var('x, y, z') sage: assume(x>=y,y>=z,z>=x) sage: bool(x==z) True sage: bool(z<x) False sage: bool(z>y) False sage: bool(y==z) True sage: forget() sage: assume(x>=1,x<=1) sage: bool(x==1) True sage: bool(x>1) False sage: forget()
- class sage.symbolic.assumptions.assuming(*args, **kwds)¶
Bases:
object
Temporarily modify assumptions.
Create a context manager in which temporary assumptions are added (or substituted) to the current assumptions set.
The set of possible assumptions and declarations is the same as for
assume()
.This can be useful in interactive mode to discover the assumptions necessary to a given integration, or the exact solution to a system of equations.
It can also be used to explore the branches of a
cases()
expression.As with
assume()
, it is an error to add an assumption either redundant or inconsistent with the current assumption set (unlessreplace=True
is used). See examples.INPUT:
*args
– assumptions (same format as forassume()
).replace
– a boolean (defaultFalse
).Specifies whether the new assumptions are added to (default) or replace (if
replace=True
) the current assumption set.
OUTPUT:
A context manager useable in a
with
statement (see examples).EXAMPLES:
Basic functionality : inside a
with assuming:()
block, Sage uses the updated assumptions database. After exit, the original database is restored.sage: var("x") x sage: forget(assumptions()) sage: solve(x^2 == 4,x) [x == -2, x == 2] sage: with assuming(x > 0): ....: solve(x^2 == 4,x) ....: [x == 2] sage: assumptions() []
The local assumptions can be stacked. We can use this functionality to discover incrementally the assumptions necessary to a given calculation (and by the way, to check that Sage’s default integrator (Maxima’s, that is), sometimes nitpicks for naught).
sage: var("y,k,theta") (y, k, theta) sage: dgamma(y,k,theta)=y^(k-1)*e^(-y/theta)/(theta^k*gamma(k)) sage: integrate(dgamma(y,k,theta),y,0,oo) Traceback (most recent call last): ... ValueError: Computation failed since Maxima requested additional constraints; using the 'assume' command before evaluation *may* help (example of legal syntax is 'assume(theta>0)', see `assume?` for more details) Is theta positive or negative? sage: a1=assuming(theta>0) sage: with a1:integrate(dgamma(y,k,theta),y,0,oo) Traceback (most recent call last): ... ValueError: Computation failed since Maxima requested additional constraints; using the 'assume' command before evaluation *may* help (example of legal syntax is 'assume(k>0)', see `assume?` for more details) Is k positive, negative or zero? sage: a2=assuming(k>0) sage: with a1,a2:integrate(dgamma(y,k,theta),y,0,oo) Traceback (most recent call last): ... ValueError: Computation failed since Maxima requested additional constraints; using the 'assume' command before evaluation *may* help (example of legal syntax is 'assume(k>0)', see `assume?` for more details) Is k an integer? sage: a3=assuming(k,"noninteger") sage: with a1,a2,a3:integrate(dgamma(y,k,theta),y,0,oo) 1 sage: a4=assuming(k,"integer") sage: with a1,a2,a4:integrate(dgamma(y,k,theta),y,0,oo) 1
As mentioned above, it is an error to try to introduce redundant or inconsistent assumptions.
sage: assume(x > 0) sage: with assuming(x > -1): "I won't see this" Traceback (most recent call last): ... ValueError: Assumption is redundant sage: with assuming(x < -1): "I won't see this" Traceback (most recent call last): ... ValueError: Assumption is inconsistent
- sage.symbolic.assumptions.assumptions(*args)¶
List all current symbolic assumptions.
INPUT:
args
– list of variables which can be empty.
OUTPUT:
list of assumptions on variables. If args is empty it returns all assumptions
EXAMPLES:
sage: var('x, y, z, w') (x, y, z, w) sage: forget() sage: assume(x^2+y^2 > 0) sage: assumptions() [x^2 + y^2 > 0] sage: forget(x^2+y^2 > 0) sage: assumptions() [] sage: assume(x > y) sage: assume(z > w) sage: sorted(assumptions(), key=lambda x: str(x)) [x > y, z > w] sage: forget() sage: assumptions() []
It is also possible to query for assumptions on a variable independently:
sage: x, y, z = var('x y z') sage: assume(x, 'integer') sage: assume(y > 0) sage: assume(y**2 + z**2 == 1) sage: assume(x < 0) sage: assumptions() [x is integer, y > 0, y^2 + z^2 == 1, x < 0] sage: assumptions(x) [x is integer, x < 0] sage: assumptions(x, y) [x is integer, x < 0, y > 0, y^2 + z^2 == 1] sage: assumptions(z) [y^2 + z^2 == 1]
- sage.symbolic.assumptions.forget(*args)¶
Forget the given assumption, or call with no arguments to forget all assumptions.
Here an assumption is some sort of symbolic constraint.
INPUT:
*args
– assumptions (default: forget all assumptions)
EXAMPLES:
We define and forget multiple assumptions:
sage: forget() sage: var('x,y,z') (x, y, z) sage: assume(x>0, y>0, z == 1, y>0) sage: sorted(assumptions(), key=lambda x:str(x)) [x > 0, y > 0, z == 1] sage: forget(x>0, z==1) sage: assumptions() [y > 0] sage: assume(y, 'even', z, 'complex') sage: assumptions() [y > 0, y is even, z is complex] sage: cos(y*pi).simplify() 1 sage: forget(y,'even') sage: cos(y*pi).simplify() cos(pi*y) sage: assumptions() [y > 0, z is complex] sage: forget() sage: assumptions() []
- sage.symbolic.assumptions.preprocess_assumptions(args)¶
Turn a list of the form
(var1, var2, ..., 'property')
into a sequence of declarations(var1 is property), (var2 is property), ...
EXAMPLES:
sage: from sage.symbolic.assumptions import preprocess_assumptions sage: preprocess_assumptions([x, 'integer', x > 4]) [x is integer, x > 4] sage: var('x, y') (x, y) sage: preprocess_assumptions([x, y, 'integer', x > 4, y, 'even']) [x is integer, y is integer, x > 4, y is even]