LDG with multiple variable with different boundary conditions
When using the LDG operator with multiple variables, only the first variable is solved for correctly.
This is due to the fact that m_traceNormals is initialized as follows
pFields[0]->GetTrace()->GetNormals(m_traceNormals);i.e., it is based on the first variable. However the ordering of the DOFs in the trace is not unique and it depends on the boundary conditions.
Attached is a minimal example to reproduce the issue.