The cmetric()
function of the ctensor
package does not compute the inverse metric when the metric has non-diagonal components. This can be remedied by explicitly computing the inverse matrix before any calls to scurvature()
or ricci()
, which seems to indicate that the computation of the inverse metric inside cmetric()
is buggy.
The attached file reproduces the error with a minimal two-dimensional non-diagonal metric and exhibits the above mentioned workaround.
wxbug_report()
output:
wxMaxima version: 14.12.1 Maxima version: 5.35.1.2 Maxima build date: 2015-01-15 21:27:01 Host type: i686-pc-mingw32 System type Win32 6.0.6001 (Service Pack 1) X86 Lisp implementation type: SBCL Lisp implementation version: 1.2.7
You are using a version of Maxima in which the matrix inversion code was broken (unrelated to ctensor). Please upgrade to a later version.
The problem has gone with Maxima version 5.36.1.
Thank you very much!