added bd->Moperator and made all bd's share J and ksp so that LU factorization sparsity structure only done once. However, code runs but no longer converges quickly. Not sure why