Ab initio calculation of orthogonal parameters for complex configurations using B-splines