I found a bug in the implementation of LCBOP I submitted some time ago.
Als,o there are some lines are used for checking of correctness of
some functions which should be commented out in the final version. The
diff is below.
< double Nij = N[i]-f_c_ij;
pair_lcbop.cpp (41.9 KB)