On Thu, 20 May 2010 16:17:20 -0500, Yujie <recrusader at gmail.com> wrote: > At least, there are "0" and "MatSetSizes" between MatConjugate and > MatRealPart in seqdense.c; > there are three "0" between MatConjugate and MatRealPart (why not be > two "0") in mpidense.c?
The indexing is wrong, run in a debugger and see which function gets called. It won't be the one you expect (but the skew should tell you where you are off). Or just use petsc-3.1 which has this feature. Jed
