Maple Validate Your CM Invariants with Maple: Open Source Code Available for Testing

  • Thread starter Thread starter bcrowell
  • Start date Start date
  • Tags Tags
    Maple
AI Thread Summary
The discussion centers on the Carminati-McLenaghan (CM) invariants, a set of 16 polynomial curvature invariants implemented in a proprietary Maple package. An open-source version of these invariants has been developed for the Maxima system, with tests conducted to verify their accuracy in various spacetimes. However, there is a lack of available reference values for the CM invariants in finite cases, making it difficult to validate the results. The author seeks assistance from anyone with access to Maple to compare outputs and confirm the correctness of the calculations. Collaboration is encouraged, especially with users familiar with Maple, to enhance the validation process.
bcrowell
Staff Emeritus
Science Advisor
Insights Author
Messages
6,723
Reaction score
431
There is a set of 16 polynomial curvature invariants called the Carminati-McLenaghan invariants, described here: https://en.wikipedia.org/wiki/Carminati–McLenaghan_invariants . They've been implemented (I think by Carminati and McLenaghan themselves) in a free Maple package described here: http://grtensor.phy.queensu.ca/Griihelp/cmscalar.help . Maple itself, however, is proprietary. I've implemented the CM invariants as open-source code https://github.com/bcrowell/cm_invariants that works in the open-source computer algebra system Maxima. I've written up a bunch of tests, e.g., calculating the invariants in spacetimes where I know that they should vanish, or spacetimes where I know that some of them should diverge at a curvature singularity. However, I haven't found any tabulations online of what the CM invariants are *supposed* to be in cases where they're finite. For example, there is an invariant called ##W_1##, and for the Schwarzschild spacetime I get ##W_1=6m^2/r^6##, but although this seems reasonable, I don't have any way to check whether it's right (e.g., whether the numerical coefficient should really be 6).

Would anyone who has a copy of Maple be willing to run the Maple implementation of the CM invariants and tell me some results that I could use to check whether my code is calculating correct output? The spacetimes that I have used so far for tests are in this test suite: https://github.com/bcrowell/cm_invariants/tree/master/tests .

Any help would be much appreciated!
 
Physics news on Phys.org
Unfortunately, I don't own Maple. If I can help with some Mathematica package, let me know.
 
In another thread, I saw that @Ray Vickson uses Maple. He may be willing to help you.
 
  • Like
Likes bcrowell
Back
Top