Added the possibility to impose the dirichlet boundary condition explicitly through the api (as in default.py->problem). This change avoids random results obtained with the vii methodology and could be reverted in the future if a better fix is found.