Compute gamma1 instead of loading it from a common block.

Fixes #447.

This will break tests since PyClaw examples need to be updated too.