More comments
Add ewald2() function from ABINIT
Test direct sum, polish
Add test_ewald2()