deal with dirichlet boundary condition
deal with dirichlet boundary condition and code cleaning
update README.md