I have changed some of the data structures for the various AIJ matrix formats. On you next pull you will need to recompile all the matrix code. cd src/mat make ACTION=libfast tree Let me know if the behavior of the matrices changes due to me update to the code. Barry I'm starting to get the various AIJ matrices to shared more common data structures and code.