Barry Smith wrote:

Good point, I have removed it.

I put it in because I wanted an easy way to test that PETSc double arrays are always 8 byte aligned (and the unaligned struct values were giving me lots of false alarms).

A related issue is PetscMallocN which (in optimized mode) gives unaligned arrays even if malloc always returns aligned memory. Consider PetscMalloc2(3,PetscInt,&ai,3,PetscScalar,&a); I think it is preferable for PetscMallocN to align all the pointers. Jed