def n_sip_terms(sip_order): assert(sip_order >= 2) n = 0 for order in range(2, sip_order+1): n += 2 * (order + 1) return n # for parallelization of matching in 'intrabrickshift'...