really cache nodes and weights up to nGQ=60, i.e. up to gauss.quad(n=61)