Processing math: 100%
Ask Your Question
0

How to assume the sum of some variables is equal to a constant?

asked 7 years ago

ensaba gravatar image

Suppose I have this expression:

a(t1+t2+t31)b

This can be simplified to 0 if we assume t1+t2+t3=1.

How can I accomplish this in Sage? I've tried:

var('t1', 't2', 't3', 'a', 'b')
expr1 = ((t1 + t2 + t3 - 1)*a)/b
expr1.full_simplify()

assume(t1 + t2 + t3 == 1)
expr1.full_simplify()

I expected the second call tofull_simplify() to return 0 but it returned the same result as the first call, which is at1+at2+at3ab

Preview: (hide)

Comments

Note that you can work ( in the sense of algebraic geometry) as follows:

sage: R.<t1,t2,t3,a,b> = PolynomialRing(QQ)
sage: Q = R.quotient( t1+t2+t3-1 ).fraction_field()
sage: expr1 = ((t1 + t2 + t3 - 1)*a)/b
sage: expr1
(t1*a + t2*a + t3*a - a)/b
sage: expr1.parent()
Fraction Field of Multivariate Polynomial Ring in t1, t2, t3, a, b over Rational Field
sage: Q( expr1 )
0

Working in the quotient ring modulo the "assumed relations" will always give the reduced expressions without contorsions.

dan_fulea gravatar imagedan_fulea ( 7 years ago )

Thanks for the suggestion. Do you have any idea why assume() didn't work?

ensaba gravatar imageensaba ( 7 years ago )
1

Asking for ?assume we see for instance in

   EXAMPLES:

   Assumptions are typically used to ensure certain relations are
   evaluated as true that are not true in general.

that the command is not designed to work with simplifications hand in hand, but rather with special, unstructured conditions that make an integral, a square root, etc. be defined, or to have evaluations of special functions like sin and cos. To give an unrelated example where we would also "expect more", let us consider...

sage: var('k'); assume( k, 'integer' );
sage: cos( k*pi ) == (-1)^k
cos(pi*k) == (-1)^k
sage: bool( cos( k*pi ) == (-1)^k )
True

We have to insist to get the whole information.

In your case just insist to get the information: bool( expr1 == 0 ) .

dan_fulea gravatar imagedan_fulea ( 7 years ago )

1 Answer

Sort by » oldest newest most voted
1

answered 7 years ago

ndomes gravatar image
var('t1', 't2', 't3', 'a', 'b')
expr1 = ((t1 + t2 + t3 - 1)*a)/b
equ = t1 + t2 + t3 == 1
expr1.subs(solve(equ,t1))
Preview: (hide)
link

Comments

Any idea why assume() didn't work?

ensaba gravatar imageensaba ( 7 years ago )

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

Question Tools

1 follower

Stats

Asked: 7 years ago

Seen: 449 times

Last updated: Aug 05 '17