ufl team mailing list archive

["Garth N. Wells" <gnw20@xxxxxxxxx>]

Is there any particular reason why 'rot' is not implemented,

    def rot(self, o, a):
        raise NotImplementedError # FIXME


and 'curl' is,

    def curl(self, o, a):
        # o = curl a = "cross(nabla, a)"
        def c(i, j):
            return a[j].dx(i) - a[i].dx(j)
        return as_vector((c(1,2), c(2,0), c(0,1))) # FIXME: Verify this


Should there be a check in curl that the vector function is in R^3?

Garth