follow dev40, finish bistabcg on all coarse grid in ./test/test-dev41.ipynb, however, the cg version does not converge, and a preliminary judgment is that r_matvec_p destroys the positive definiteness of the original matvec, moreover, the tolerance degree of the result error on the coarse grid needs further discussion...... orz