debugged. You may turn on the intralayer interaction by uncommenting the compute_coulomb_form line in function run