code: def spa_tmplt_engine(htilde,  kmin,  phase_order,
                      delta_f,  piM,  pfaN,
                      pfa2,  pfa3,  pfa4,  pfa5,  pfl5,
                      pfa6,  pfl6,  pfa7, amp_factor):
      
      taylorf2_kernel(htilde.data,  kmin,  phase_order,
                      delta_f,  piM,  pfaN,
                      pfa2,  pfa3,  pfa4,  pfa5,  pfl5,
                      pfa6,  pfl6,  pfa7, amp_factor)
docs: Calculate the spa tmplt phase