Ask Your Question
3

Using matrix elements as arguments

asked 2010-11-26 10:44:53 +0200

DisneySage gravatar image

updated 2010-11-28 13:44:48 +0200

niles gravatar image

I have a rather easy question, or so it would seem. I have looked for an answer but was unable to find one anywhere so I'm asking it here.

I am making a very simple iterative algorithm for which the input as well as the output at the end of every iteration is a vector (or matrix for that matter). What I want to do is use these elements as arguments for several functions during each of the iteration. So for example

x=var('x')
y=var('y')
z=matrix(2,1,[ [1],[1] ]
f=x^2+y^3
H=f.hessian()

Then what I would like to do is say

H(z[0],z[1])

or

H(z)

But no matter what I try I can't seem to get it to work. Ideas?

edit retag flag offensive close merge delete

3 Answers

Sort by ยป oldest newest most voted
4

answered 2010-11-28 13:42:41 +0200

niles gravatar image

There is a subtle difference between "symbolic expressions" and "callable symbolic expressions", which are also termed "functions". Writing

sage: f=x^2+y^3

makes f a symbolic expression, which displays as

sage: f
x^2 + y^3

On the other hand, writing

sage: g(x,y)=x^2+y^3

makes g a callable symbolic expression, which displays as

sage: g
(x, y) |--> x^2 + y^3

With f, you need to use f.substitute to substitute values, but with g, since you have already informed sage of the variable order, you can use it like a function

sage: g(1,3)
28

The way you define your finction determines what kind of expression the corresponding Hessian is; note the difference in syntax below:

First method:

sage: x=var('x')
sage: y=var('y')
sage: z=matrix(2,1,[ [1],[1] ])
sage: f=x^2+y^3
sage: H=f.hessian()
sage: H.substitute(x=z[0,0],y=z[1,0])
[2 0]
[0 6]

sage: f
x^2 + y^3
sage: H
[  2   0]
[  0 6*y]

Second method, making g a callable symbolic expression:

sage: x=var('x')
sage: y=var('y')
sage: z=matrix(2,1,[ [1],[1] ])
sage: g(x,y)=x^2+y^3
sage: H=g.hessian()
sage: H(z[0,0],z[1,0])
[2 0]
[0 6]

sage: g
(x, y) |--> x^2 + y^3
sage: H
[  (x, y) |--> 2   (x, y) |--> 0]
[  (x, y) |--> 0 (x, y) |--> 6*y]
edit flag offensive delete link more
2

answered 2010-11-26 12:06:51 +0200

updated 2010-11-26 12:08:04 +0200

This works for me:

sage: H.substitute(x=z[0,0],y=z[1,0])

(z is a matrix, so it requires two indices to specify an element.)

edit flag offensive delete link more
0

answered 2015-01-06 00:08:51 +0200

wonder gravatar image

I guess

H = f.hessian().function(x,y)
H(z[0,0],z[1,0])

would also work

edit flag offensive delete link more

Your Answer

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

Add Answer

Question Tools

Stats

Asked: 2010-11-26 10:44:53 +0200

Seen: 2,197 times

Last updated: Jan 06 '15