Zeroth order solution in c_test is incorrect
The zeroth order solution in c_test should be the final solution if the elastic constants are uniform throughout the system. However, this is not true. This is also the case in the original Fortran code so it isn't a porting issue.