For problems with zero-conductivity regions, project the vector onto the compatible subspace:, where
is the discrete gradient restricted to the interior nodes of the regions with zero conductivity.
For problems with zero-conductivity regions, project the vector onto the compatible subspace:, where
is the discrete gradient restricted to the interior nodes of the regions with zero conductivity. This ensures that x is orthogonal to the gradients in the range of
.
This function is typically called after the solution iteration is complete, in order to facilitate the visualization of the computed field. Without it the values in the zero-conductivity regions contain kernel components.