的非线性微分方程，我可以用以下形式表示： ，类似于Fenics 教程中给出的“非线性泊松”示例，除了非线性现在是一个矩阵。$u:\mathbb{R}^2\rightarrow\mathbb{R}^3$$\nabla (g(u)_{ij} \nabla(u)^j) = 0$$g$

但是我想添加一个额外的表单约束：

$u_x \cdot u_y = 0$。

我怀疑有一种方法可以使用拉格朗日乘数来做到这一点，但我不确定如何表达这一点。