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.
Use larger L1 cache and PTX annotations by default in CUDA code.
Add an option to disable the L1 cache.
Add support for the propagate on-read scheme.