Solver in the latest stage. The scheme works good but the fixed point iterations seems to have problem. That is the result diverges only in the fixed point iterations. It should be checked.