Hello Everyone,

I want to implement a model in FEniCS which uses a so called “local newton iteration”, where the equation for the internal (tensor-valued) variable \boldsymbol{a} is to be solved at the Gaus point level. From the example of plasticity from Jeremy Bleyer I understand that I should use `Quadrature`

-Elements for the representation at the Gauss-Points.

The equation has the form \boldsymbol{a} = \exp{\boldsymbol{b}} \cdot \boldsymbol{c} where \boldsymbol{b} = \boldsymbol{b}(\boldsymbol{a}) is a non-linear function of \boldsymbol{a} and \boldsymbol{c} = \boldsymbol{c}(\boldsymbol{a} ,\boldsymbol{u}) a nonlinear function of \boldsymbol{a} and the displacement \boldsymbol{u}.

What I do not know is how to solve equations in FEniCS at the Gauss points without defining a weak form. Is this possible?

Thanks in advance!