A new and rigorous SPN theory – Part IV: Numerical qualification of GSP3(0) and the generalized transverse integration nodal method