Subprogram multidip_levin::nested_cgreen_correct_levin (Z, a, b, c, N, sa, sb, m, l, k, v): Generalize for arbitrary dimension!
Subprogram multidip_romberg::nested_cgreen_correct_romberg (Z, a, b, c, N, sa, sb, m, l, k, v): Generalize for arbitrary dimension!
Subprogram multidip_romberg::nested_cgreen_eval (Z, c, N, sa, sb, m, l, k, rs, asy): Generalize to arbitrary dimension!
Subprogram multidip_routines::multiint_chain (moldat, r0, Ei, esc, omega, ie, c, N, state, ichanf, sb, k, l, m): We do not need to calculate the integrals for all final-initial channel pairs as done now, because when the initial channels have the same quantum number 'm', the integrals will differ only by the coefficient 'ap'. This would reduce the run time for calculations with large partial wave expansion.
Subprogram multidip_routines::setup_initial_state (states, moldat, irr, lubnd, Ei): Take into account modified bound state energy calculated by BOUND.