void printColoring (const Epetra_MapColoring & ColorMap, Epetra_CrsGraph * A, bool verbose) {
int NumColors = ColorMap.NumColors();
int * ListOfColors = ColorMap.ListOfColors();
EpetraExt::CrsGraph_MapColoringIndex MapColoringIndexTransform( ColorMap );
vector<Epetra_IntVector> & ColIndices = MapColoringIndexTransform( *A );
if( verbose )
{
cout << endl;
cout << "***************************************\n";
cout << "Column Indexing by Color:\n";
cout << "***************************************\n";
cout << endl;
for( int i = 0; i < NumColors; ++i )
{
cout << "COLOR: " << ListOfColors[i] << endl;
cout << ColIndices[i];
}
cout << endl;
}
return;
}
