Nick

@ksoni-1258 : I don't think I ever removed the imaginary part of the complex number to do complex-valued quadrature.

I think instead of iterating through the matrix you need to make sure the overloads are mapping appropriately. I think an `if constexpr` branch on the matrix type might help...

L1 generally indicates an L1 norm. The L1 norm allows one to compute the condition number of integration and hence give a correct estimate of the error that can be...

This shows that you have a misuse of the Eigen library. You might have bit off a bit more than you can chew with this particular task; could there be...

Yeah I'm fine with the `if constexpr` on the type.

I'm confused; do you want documentation that conventions are different or do you think this is a bug?

We could use `std::random` with a bunch of casting, though we wouldn't be able to query abscissas that are not representable as `float`, `double`, and `long double`. Is it super...

Ok, if this is gonna make a mess, I'll go ahead and try to clean it up. This hasn't been released yet, so we can make backwards-incompatible changes. Let me...

@jzmaddock : Is there any way we could do a module while preserving what we have for C++11 users?

Looks like a great start. I'll go through your other items on your checklist later.