minor updates
 * changed outer and inner loop (j->i) vs (i->j) for better memory locality
 * disabled suptitle and printmat functions
 * moved constant coefficient computation out of a double loop