40 Reputation

One Badge

2 years, 25 days

MaplePrimes Activity

These are questions asked by prldb

Suppose I specify a metric in say the differential geometry package or Physics package that has an arbitary function f(x,y) of some of the coordinates as a component. Can I specify Maple to treat the function as real valued when Maple is asked to compute curvature etc for the metric? In some examples I have performed Maple sometimes returns expressions in f(x,y) for components that are in fact zero once one treats f(x,y) as real valued. It would be preferable to have Maple actually return '0' in such cases.


I tried the 'assume' command but Maple complained it was being applied to a protected name and I tried using

use RealDomain in simplify (5) end use

after Maple's computations, where (5) lablled Maple's output, but Maple didn't perform any simplification.

Given a metric, to compute quantities in the NP formalism one needs to specify a null tetrad. In the various examples in the help pages, sometimes the tetrad is specified simply as a list of 4 vectors, e.g., NT := [...] and sometimes evalDG is applied as in NT := evalDG({...]). Using the first format, Maple accepted NT as argument in NPSpinCoefficients but  NPCurvatureScalars(SpinCoefficients,NT) complained that the second argument wasn't a list of four vectors. When I used the second format, both commands returned the expected results. Why the difference?

If I call for the metric (27.27) in Stephani et al. in the Physics package, I expected the null tetrad employed to compute say the Weyl scalars would be the one given in conjunction with (27.27) (equation (27.22) in Stephani et al.). But, Weyl[scalars] returned long expressions for each of the five scalars, whereas with respect to the null tetrad in Stephani et al one expects the first two (or last two) scalars to be zero since the null vector k in their tetrad is the multiple pnd of Weyl. e_[nullvectors] and Setup(tetradmetric=null) followed by e_[ ] seem to output the same null tetrad, which does not appear to be that of Stephani et al. I assume this null tetrad is the null tetrad associated to the orthonomral tetrad that e_[ ]  returns if one hasn't used Setup(tetradmetric=null). How does Maple select this default orthonormal tetrad? What is the best way to set the null tetrad of Stephani et al as the null tetrad to compute Weyl[scalars]?

Here is a simpler example. Calling [27, 37, 1]. Stephani et al give the null tetrad in terms of the spacetime coordinates along with the metric in their equaiton (27.37). After using Setup(tetradmetic=null), e_[ ] returns a tetrad that might, I suppose, be that of Stephani et al. in disguise. Specifically, the tetrad vector defining the null congruence should have only a component with respect to the coordinate r, yet the Maple output gives expressions for the components with respct to the u and r coordinates for both the k and l elements of the null tetrad (while the expressions for the complex element m is exactly what one would expect). As in the previous example, all Weyl scalars are given by nontrivial expressions, even though two of them should be zero (since k is a multiple pnd). So is it the case that the experssions in the Maple output where one expects zero are in fact zero in disguise? The experssions are complicated enough that it is not obvious.  I have uploaded the Maple document for these calculations.

The Help page Physics/tensors-a complete guide states that spacetime metrics from Kramer et al. are referenced by chapter, section, and equation number, e.g., g_[[12, 16, 1]]. But there is no section 16 in Ch 12 and equations within each chapter are numbered sequentially without reference to section. By playing around it seems that in fact the first number is chapter, the second number is equation number, and the third number refers to subcases of the metric, when they are specified in the text. Is that correct?

Also, the output I get from say g_[[27, 27, 1]], or any other attempt  made, is just the metric, without any specification of the coordinates etc, which the Help pages susggest should be part of the output.

I have been computing curvature etc of a four-dimensional metric that invloves three arbitary functions using the DifferentialGeometry package. I am now intereted in whether Maple can compute the Killing vectors. I have used the code of Example 3 from the Maple help page on KillingVectors in the  Tensor subpackage of the package DifferentialGeoemtry. If I set the arbitrary functions to be constants using the' auxiliaryequations' part of the code, then Maple outputs Killing vectors, so all good so far. If I leave the three arbitary functions to be arbitrary, however, specifying only that they be nonzero (as in Example 3) Maple had not produced any output after five hours, though when I moved the cursor over the tool bars at the top of the screen it turned to an hour glass, indicating that Maple was busy. Does anyone have experience with whether I am asking too much of Maple or whether I just need to give it more time to produce answers? If the latter, how long? The Example 3 indicates Maple will consider several possibilties including the completely arbitrary option, though that example only involves one arbitary function, so having multiple arbitary functions obviously increases the complexity.

Page 1 of 1