Sorry I didn't state my problem more accurate. The entries in the input matrixes and vectors​ are correct -- they have right values and correct nonzero structure. The KSP was setup so I can solve the lienar systems for an unsteady state convective problem. I use KSPBCGS + PCILU and so far I only try sequential on a very small problem. I wonder if blocked ilu is actually used in MATBAIJ2MATAIJ case.
Thanks,
