finsberg / cardiac_benchmark Goto Github PK

Checking the line of code for pressure I see that its being used ufl.cofac, but the presssure boundary condition is given as:

Could it be that a transpose T is missing at the end? I would think something like:

def _pressure_term(self, F, w):
        return dolfin.inner(self.parameters["p"] * ufl.det(F) * ufl.inv(ufl.grad(F)).T * self.N, w) * self.ds(
            self.endo,
        )