Use omega instead of dividing by tau0 in the relaxation code.
Move some functionality into the sym_codegen module.

Change KernelCodePrinter to minimize the number of division operations.
Add an explanation for why we modify _print_Mul.