Don't compile in clang

Some more testing is needed, but looks promising so far

Test that all derivatives match with those of finite differentiation to numerical precision (they do!)