dielec_mat.f90

 (Walter&Cohen obtain 11.3 with 15 bands)

      kmesh         bands      cutoff      epsilon
                                Ry    (G=G'=0,q=[0.01 0 0])

  4  4  4 1 1 1      15          2         19.41619
  4  4  4 0 0 0      15          2         12.12227
  5  5  5 1 1 1      15          2         11.65866 
  5  5  5 0 0 0      15          2         15.81732
  6  6  6 1 1 1      15          2         13.89199
  6  6  6 0 0 0      15          2         11.51096
 10 10 10 1 1 1      15          2         11.69469
 10 10 10 0 0 0      15          2         11.44243
 20 20 20 1 1 1      15          2         11.44239
 20 20 20 0 0 0      15          2         11.44160 

 10 10 10 0 0 0      15          2         11.44243
 10 10 10 0 0 0      30          2         11.44743
 10 10 10 0 0 0      50          2         11.44773 

 10 10 10 0 0 0      15          2         11.44243
 10 10 10 0 0 0      15          5         11.44249 
 10 10 10 0 0 0      15         20         11.44246




 coulomb.f90 [ 1/epsi^{-1}(G=0,G'=0,q) ]

      kmesh        cutoff      epsilon
                    Ry    (G=G'=0,q=[0.01 0 0])

  4  4  4 1 1 1      2          17.6616 
  4  4  4 0 0 0      2          10.8284
  5  5  5 1 1 1      2          10.4232
  5  5  5 0 0 0      2          14.3534 
  6  6  6 1 1 1      2          12.5597
  6  6  6 0 0 0      2          10.3008
 10 10 10 1 1 1      2          10.4756                    
 10 10 10 0 0 0      2          10.2375      
 20 20 20 1 1 1      2          10.2323   
 20 20 20 0 0 0      2          10.2323       

 10 10 10 0 0 0      2          10.2375 
 10 10 10 0 0 0      5          10.2218
 10 10 10 0 0 0     20          10.2291