Great. Normally I have no use for Gauss rules beyond order 4 (at most cubic finite elements in the library). But it is good to have access to arbitrary order.
Thanks! P On Wednesday, December 17, 2014 1:35:34 PM UTC-8, Steven G. Johnson wrote: > > In particular, it looks like you just need: > > param_coords, weights = Base.QuadGK.gauss(Float64, order) > > in your GaussRule function. >