Coordinate calculation of the Ricci and Weyl tensors from a basis

> restart:

> grtw();

`GRTensorII Version 1.79 (R6)`

`2 February 2001`

`Developed by Peter Musgrave, Denis Pollney and Kay...

`Copyright 1994-2001 by the authors.`

`Latest version available from: http://grtensor.phy...

`e:/Grtii(6)/Metrics`

We load the bases.

> qload(mix1);

`Default spacetime` = mix1

`For the mix1 spacetime:`

Coordinates

x(up)

`x `^a = vector([Theta, Phi, Psi, T])

`Basis inner product`

eta(bup,bup)

eta^`(a)`*``^`(b)` = matrix([[1, 0, 0, 0], [0, 1, 0...

`Basis (covariant components)`

w1(dn)

omega1[a] = vector([-exp(a(T))*sqrt(1-Psi^2)/(sqrt(...

w2(dn)

omega2[a] = vector([-exp(b(T))*Psi/(sqrt(1-Theta^2)...

w3(dn)

omega3[a] = vector([0, Theta*exp(c(T)), exp(c(T))/(...

w4(dn)

omega4[a] = vector([0, 0, 0, exp(a(T)+b(T)+c(T))])

`Mixmaster metric (e.g. MTW Box 30.1, Theta = cos(t...

From the bases we calculate the metric components.

> grcalc(g(dn,dn,pdn),g(up,up),detg);

`CPU Time ` = .50e-1

> gralter(_,6,7);

Component simplification of a GRTensorII object:

Applying routine expand to object g(dn,dn,pdn)

Applying routine expand to object g(up,up)

Applying routine expand to object detg

Applying routine factor to object g(dn,dn,pdn)

Applying routine factor to object g(up,up)

Applying routine factor to object detg

`CPU Time ` = .140

Though not needed for the following calculations, it is sometimes convenient to see the line element.

> grcalc(ds);

`CPU Time ` = .10e-1

> grdisplay(_);

`For the mix1 spacetime:`

`Line element`

` ds`^2 = (-exp(a(T))^2+exp(a(T))^2*Psi^2-exp(b(T))...
` ds`^2 = (-exp(a(T))^2+exp(a(T))^2*Psi^2-exp(b(T))...
` ds`^2 = (-exp(a(T))^2+exp(a(T))^2*Psi^2-exp(b(T))...

We now proceed to calculate the coordinate components of the Ricci and Weyl tensors.

> grcalc(Chr(dn,dn,up));

`CPU Time ` = .141

> gralter(_,6,7);

Component simplification of a GRTensorII object:

Applying routine expand to object Chr(dn,dn,up)

Applying routine factor to object Chr(dn,dn,up)

`CPU Time ` = .190

> grcalcalter(R(dn,dn),6);

Simplification will be applied during calculation.

Applying routine expand to object R(dn,dn)

`CPU Time ` = .40e-1

> gralter(_,6,7);

Component simplification of a GRTensorII object:

Applying routine expand to object R(dn,dn)

Applying routine factor to object R(dn,dn)

`CPU Time ` = .110

> grmap(_,autoAlias,`x`);

Applying routine autoAlias to R(dn,dn)

Warning, alias or macro b[T] defined in terms of b

Warning, alias or macro a[T] defined in terms of a

Warning, alias or macro c[T] defined in terms of c

> grdisplay(_);

`For the mix1 spacetime:`

`Covariant Ricci`

R(dn,dn)

`R `[Theta*Theta] = 1/2*(exp(a)^2*exp(c)^4+exp(b)^2...
`R `[Theta*Theta] = 1/2*(exp(a)^2*exp(c)^4+exp(b)^2...

`R `[Theta*Phi] = -1/2*Psi*(Psi-1)^2*(Psi+1)^2*(exp...
`R `[Theta*Phi] = -1/2*Psi*(Psi-1)^2*(Psi+1)^2*(exp...

`R `[Phi*Phi] = -1/2*(-exp(b)^2*Theta^2*exp(c)^4+2*...
`R `[Phi*Phi] = -1/2*(-exp(b)^2*Theta^2*exp(c)^4+2*...
`R `[Phi*Phi] = -1/2*(-exp(b)^2*Theta^2*exp(c)^4+2*...
`R `[Phi*Phi] = -1/2*(-exp(b)^2*Theta^2*exp(c)^4+2*...
`R `[Phi*Phi] = -1/2*(-exp(b)^2*Theta^2*exp(c)^4+2*...

`R `[Phi*Psi] = 1/2*Theta*(Psi-1)*(Psi+1)*(exp(b)^4...

`R `[Psi*Psi] = 1/2*(exp(b)^4-exp(c)^4+exp(a)^4-2*c...

`R `[T*T] = -a[T,T]-b[T,T]-c[T,T]+2*a[T]*b[T]+2*c[T...

> grcalcalter(Ricciscalar,6);

Simplification will be applied during calculation.

Applying routine expand to object tRicciscalar

`CPU Time ` = 0.

> gralter(_,6,7);

Component simplification of a GRTensorII object:

Applying routine expand to object Ricciscalar

Applying routine factor to object Ricciscalar

`CPU Time ` = .30e-1

> grcalcalter(C(dn,dn,dn,dn),6);

Simplification will be applied during calculation.

Applying routine expand to object R(dn,dn,dn,dn)

Applying routine expand to object tRicciscalar

Applying routine expand to object C(dn,dn,dn,dn)

`CPU Time ` = .220

> gralter(_,6,7);

Component simplification of a GRTensorII object:

Applying routine expand to object C(dn,dn,dn,dn)

Applying routine factor to object C(dn,dn,dn,dn)

`CPU Time ` = 1.692

> grmap(_,autoAlias,`x`);

Applying routine autoAlias to C(dn,dn,dn,dn)

> grdisplay(_);

`For the mix1 spacetime:`

`Covariant Weyl`

`C `[Theta*Phi*Theta*Phi] = 1/6*(-4*exp(a)^2*a[T]*e...
`C `[Theta*Phi*Theta*Phi] = 1/6*(-4*exp(a)^2*a[T]*e...
`C `[Theta*Phi*Theta*Phi] = 1/6*(-4*exp(a)^2*a[T]*e...
`C `[Theta*Phi*Theta*Phi] = 1/6*(-4*exp(a)^2*a[T]*e...
`C `[Theta*Phi*Theta*Phi] = 1/6*(-4*exp(a)^2*a[T]*e...
`C `[Theta*Phi*Theta*Phi] = 1/6*(-4*exp(a)^2*a[T]*e...
`C `[Theta*Phi*Theta*Phi] = 1/6*(-4*exp(a)^2*a[T]*e...
`C `[Theta*Phi*Theta*Phi] = 1/6*(-4*exp(a)^2*a[T]*e...
`C `[Theta*Phi*Theta*Phi] = 1/6*(-4*exp(a)^2*a[T]*e...
`C `[Theta*Phi*Theta*Phi] = 1/6*(-4*exp(a)^2*a[T]*e...
`C `[Theta*Phi*Theta*Phi] = 1/6*(-4*exp(a)^2*a[T]*e...

`C `[Theta*Phi*Theta*Psi] = -1/6*(Psi-1)*(Psi+1)*(-...
`C `[Theta*Phi*Theta*Psi] = -1/6*(Psi-1)*(Psi+1)*(-...
`C `[Theta*Phi*Theta*Psi] = -1/6*(Psi-1)*(Psi+1)*(-...
`C `[Theta*Phi*Theta*Psi] = -1/6*(Psi-1)*(Psi+1)*(-...
`C `[Theta*Phi*Theta*Psi] = -1/6*(Psi-1)*(Psi+1)*(-...
`C `[Theta*Phi*Theta*Psi] = -1/6*(Psi-1)*(Psi+1)*(-...

`C `[Theta*Phi*Theta*T] = 1/2*Psi*Theta*(Psi-1)*(Ps...
`C `[Theta*Phi*Theta*T] = 1/2*Psi*Theta*(Psi-1)*(Ps...

`C `[Theta*Phi*Phi*Psi] = -1/6*(-3*exp(a)^2*exp(b)^...
`C `[Theta*Phi*Phi*Psi] = -1/6*(-3*exp(a)^2*exp(b)^...
`C `[Theta*Phi*Phi*Psi] = -1/6*(-3*exp(a)^2*exp(b)^...

`C `[Theta*Phi*Phi*T] = -1/2*Theta*(-3*exp(b)^2*Psi...
`C `[Theta*Phi*Phi*T] = -1/2*Theta*(-3*exp(b)^2*Psi...

`C `[Theta*Phi*Psi*T] = -1/2*(-exp(b)^2*a[T]-b[T]*e...

`C `[Theta*Psi*Theta*Psi] = -1/6*(-exp(a)^2*exp(b)^...
`C `[Theta*Psi*Theta*Psi] = -1/6*(-exp(a)^2*exp(b)^...
`C `[Theta*Psi*Theta*Psi] = -1/6*(-exp(a)^2*exp(b)^...
`C `[Theta*Psi*Theta*Psi] = -1/6*(-exp(a)^2*exp(b)^...
`C `[Theta*Psi*Theta*Psi] = -1/6*(-exp(a)^2*exp(b)^...

`C `[Theta*Psi*Theta*T] = -1/2*Psi*(-3*exp(b)^2*b[T...

`C `[Theta*Psi*Phi*Psi] = 1/6*(Psi-1)*(Psi+1)*(-3*e...
`C `[Theta*Psi*Phi*Psi] = 1/6*(Psi-1)*(Psi+1)*(-3*e...
`C `[Theta*Psi*Phi*Psi] = 1/6*(Psi-1)*(Psi+1)*(-3*e...

`C `[Theta*Psi*Phi*T] = -1/2*(Psi^2*b[T]*exp(c)^2+2...
`C `[Theta*Psi*Phi*T] = -1/2*(Psi^2*b[T]*exp(c)^2+2...

`C `[Theta*T*Theta*T] = 1/6*(-2*exp(a)^2*exp(b)^2*e...
`C `[Theta*T*Theta*T] = 1/6*(-2*exp(a)^2*exp(b)^2*e...
`C `[Theta*T*Theta*T] = 1/6*(-2*exp(a)^2*exp(b)^2*e...
`C `[Theta*T*Theta*T] = 1/6*(-2*exp(a)^2*exp(b)^2*e...
`C `[Theta*T*Theta*T] = 1/6*(-2*exp(a)^2*exp(b)^2*e...

`C `[Theta*T*Phi*Psi] = -1/2*(exp(c)^2*c[T]+3*exp(a...
`C `[Theta*T*Phi*Psi] = -1/2*(exp(c)^2*c[T]+3*exp(a...

`C `[Theta*T*Phi*T] = -1/6*Psi*(Psi-1)^2*(Psi+1)^2*...
`C `[Theta*T*Phi*T] = -1/6*Psi*(Psi-1)^2*(Psi+1)^2*...
`C `[Theta*T*Phi*T] = -1/6*Psi*(Psi-1)^2*(Psi+1)^2*...

`C `[Phi*Psi*Phi*Psi] = 1/6*(Theta-1)*(Theta+1)*(ex...
`C `[Phi*Psi*Phi*Psi] = 1/6*(Theta-1)*(Theta+1)*(ex...
`C `[Phi*Psi*Phi*Psi] = 1/6*(Theta-1)*(Theta+1)*(ex...
`C `[Phi*Psi*Phi*Psi] = 1/6*(Theta-1)*(Theta+1)*(ex...
`C `[Phi*Psi*Phi*Psi] = 1/6*(Theta-1)*(Theta+1)*(ex...

`C `[Phi*Psi*Phi*T] = 1/2*Psi*(Theta-1)*(Theta+1)*(...

`C `[Phi*T*Phi*T] = -1/6*Theta^2*exp(c)^4*exp(a)^2-...
`C `[Phi*T*Phi*T] = -1/6*Theta^2*exp(c)^4*exp(a)^2-...
`C `[Phi*T*Phi*T] = -1/6*Theta^2*exp(c)^4*exp(a)^2-...
`C `[Phi*T*Phi*T] = -1/6*Theta^2*exp(c)^4*exp(a)^2-...
`C `[Phi*T*Phi*T] = -1/6*Theta^2*exp(c)^4*exp(a)^2-...
`C `[Phi*T*Phi*T] = -1/6*Theta^2*exp(c)^4*exp(a)^2-...
`C `[Phi*T*Phi*T] = -1/6*Theta^2*exp(c)^4*exp(a)^2-...
`C `[Phi*T*Phi*T] = -1/6*Theta^2*exp(c)^4*exp(a)^2-...
`C `[Phi*T*Phi*T] = -1/6*Theta^2*exp(c)^4*exp(a)^2-...
`C `[Phi*T*Phi*T] = -1/6*Theta^2*exp(c)^4*exp(a)^2-...
`C `[Phi*T*Phi*T] = -1/6*Theta^2*exp(c)^4*exp(a)^2-...

`C `[Phi*T*Psi*T] = 1/6*exp(c)^2*(Psi-1)*(Psi+1)*(e...
`C `[Phi*T*Psi*T] = 1/6*exp(c)^2*(Psi-1)*(Psi+1)*(e...
`C `[Phi*T*Psi*T] = 1/6*exp(c)^2*(Psi-1)*(Psi+1)*(e...

`C `[Psi*T*Psi*T] = 1/6*exp(c)^2*(exp(b)^4+exp(b)^2...
`C `[Psi*T*Psi*T] = 1/6*exp(c)^2*(exp(b)^4+exp(b)^2...
`C `[Psi*T*Psi*T] = 1/6*exp(c)^2*(exp(b)^4+exp(b)^2...

>