Dear all,
by reading related documents, now I understand that the terms containing log(s) come from R^{lub}.
The only left question is about the variables {r2bcorr_para_self, r2bcorr_para_mix, r2bcorr_perp_self, r2bcorr_perp_self} in function "sd_compute_resistance_matrix_sparse()".
The first two variables corresponds to \alpha_{s} and \alpha_{m} in equations (3.17), (3.21) in David's thesis.
But how about the latter two ?
Would anyone like to do me a favour, and to explain a little bit where they come from ?
I get confused because they are different from my expectation,
\beta_{s} or \beta_{m} in equations (3.19) and (3.23).