removed mallocs from inner loop of Mat_CheckInode_FactorLU, also added a bunch of const to array variables in that file
move Mat_CheckInode_FactorLU to symbolic factorization so it doesn't get called in EVERY numeric factorization

set matsolve function pointer for SeqAIJ in the numeric factorization, not in the checkinode so matsolve cannot be called on matrix that has not been numerically factored
optimize MatSetValues in fdcoloring using an array of spidx - see MatTransposeColoringCreate_SeqAIJ()
-------------------
new impl of MatFDColoringxxx_SeqAIJ works now :-)
mv MatGetColumnIJ_SeqAIJ_Color() from file matmatmult.c to aij.c
optimize MatFDColoringApply_SeqAIJ using local const array pointers
mv checking !mat->assembled from MatFDColoringCreate_xxx() to interface MatFDColoringCreate
mv redundant valications in MatFDColoringApply_xxx()
cleanup and optimize MatFDColoringApply_SeqAIJ()
cleanup and optimize
add timer for profiling segments of MatFDColoringApply_SeqAIJ
experimenting different ways of storing row, col and den2sp -- rowcolden2sp3 gives the best performance
cleanup MatFDColoringCrete and MatFDColoringApply for SeqAIJ
cleanup
fix memory leak
mv MatFDColoringApply_SeqAIJ() from aij.c to fdaij.c; rm MatFDColoringApply_SeqAIJ_old; MatFDColoringCreate_SeqAIJ_den2sp() supports SeqBAIJ matrix.
add MatFDColoringApply_SeqBAIJ
bugfix of MatFDColoringApply_BAIJ
implement version of STOREVALS for MatFDColoringApply_SeqBAIJ, but it does not improve the version of !STOREVALS
rm impl of STOREVALS