# How to set a metric tensor inverse?

I have a complicated metric tensor. In order to compute the Laplacian in this metric, I need the metric tensor inverse. But this exceeds my memory capacity after ~36 hours of computing.

I would like instead to set the inverse with what I believe to be the correct entries. Is that possible? Reviewing the documentation and the source code, I see no indication of such an operation.

I wonder if g.inverse().comp(XI.frame())[:] = Ginv[:] would work, where g denotes the metric tensor, and Ginv the matrix of expressions which should be correct.

BTW, I have the metric and its inverse in one frame. It is after mapping into a second frame [XI] that the problem arises.

edit retag close merge delete

Sort by » oldest newest most voted

g._inverse.set_comp(XI.frame())[:] = Ginv[:] should work (notice that it starts with g._inverse, not g.inverse(), and that set_comp must be used instead of comp).

Probably we should implement a method set_inverse() to allow easily for this.

more