B0 = funB(r1.ravel())