interim
test_gradient_fwd gives J values 1/2 of what they should be