Precomputed weights and nodes? about integratorxx HOT 1 CLOSED

I'm not a fan of keeping around precomputed grids which can take an arbitrary npts input (e.g. GC, GL, etc). For some (e.g. Gill's MultiExp) it can't be (practically) avoided, so I'm fine tabulating those. The reason we tabulate the angular quadratures is that they're hard to compute. Computing these quadratures is trivial, so this wouldn't be a necessary performance optimization and would increase compile times, etc.

Your point regarding precision is well taken - I would opt to simply have a path that computes the nodes/weights in multiple / arbitrary precision programatically. There are C/C++ libraries that can do this (GMP comes to mind, but its GPL...). All the functions are templated, so it should be relatively straightforward to do.

We can definitely track that if there's an immediate use case (maybe for the unit tests as you suggest?)

from integratorxx.