update: fulled_lda1b update with a variable to control which orbital density to fix
update: fix a small bug into search_mu call in lda1b