First implementation of forces due to deriv of gsi, compiles, no difference in output
since additional matrices were set to zero