Ask Your Question
0

manifolds : ricci_scalar() has no attribute 'at' [closed]

asked 2020-05-20 15:52:31 +0100

LPsFR gravatar image

gP3.ricci_scalar().at(p) does not produce value of ricci_scalar at point p. (gP3 is my metric). whilst ricci().at(p) is working. Error message : AttributeError: 'DiffScalarFieldAlgebra_with_category.element_class' object has no attribute 'at'

edit retag flag offensive reopen merge delete

Closed for the following reason the question is answered, right answer was accepted by LPsFR
close date 2020-05-21 11:41:47.627017

1 Answer

Sort by ยป oldest newest most voted
0

answered 2020-05-20 21:38:02 +0100

eric_g gravatar image

updated 2020-05-21 10:48:25 +0100

To get the value of a scalar field at a point, simply use the call method, i.e. the parenthesis operator. In your case:

gP3.ricci_scalar()(p)

Actually, at() is reserved to tensor fields of valence $>0$, since for them the call method has a different meaning. For instance, if g is the metric tensor and u and v are two vector fields, the call method of g is used to denote the bilinear form action of g on the pair (u,v) as g(u,v). The following identity, which involves the various call methods and at(), holds:

g.at(p)(u.at(p), v.at(p)) == g(u, v)(p)

Here is a full example:

sage: E.<x,y> = EuclideanSpace()
sage: g = E.metric()
sage: g.display()
g = dx*dx + dy*dy
sage: u = E.vector_field(-y, x)
sage: v = E.vector_field(x+y, x-y)
sage: p = E((2, 3)); p
Point on the Euclidean plane E^2
sage: bool( g.at(p)(u.at(p), v.at(p)) == g(u, v)(p) )
True
edit flag offensive delete link more

Comments

As simple as that. Many thanks.

LPsFR gravatar imageLPsFR ( 2020-05-21 11:11:39 +0100 )edit

Question Tools

Stats

Asked: 2020-05-20 15:52:31 +0100

Seen: 287 times

Last updated: May 21 '20