faster? implicit integration, reduced max nb of iter of cg to 20, implemented solver without matrix assembly