I have a simple SeqAIJ example where MatTransposeMatMult is 30x slower than MatTranspose followed by MatMatMult. This is embarrassing when people try to compare performance. Are there any cases where these custom implementations provide value rather than just confusing people or turning them away from PETSc?