Use dirichlet constraints instead of setting values to zero by hand