Lisandro, Thanks, I have incorporated these and a few other blocksize fixes. Barry On Apr 10, 2012, at 4:36 AM, Lisandro Dalcin wrote:
Barry, could you review the attached patch? Basically, many Mat preallocation routines should not set the block size to 1, and instead honor what the user could have set calling MatSetBlockSize(). Please note that in this patch I'm not addressing all instances of this issue in the codebase, I've just worked on some of the matrix types (seq/mpidense, seq/mpiaij, shell) that I know for sure the fix seems to make sense.
-- Lisandro Dalcin --------------- CIMEC (INTEC/CONICET-UNL) Predio CONICET-Santa Fe Colectora RN 168 Km 472, Paraje El Pozo 3000 Santa Fe, Argentina Tel: +54-342-4511594 (ext 1011) Tel/Fax: +54-342-4511169 <matblocksize.diff>