gauss quadrature #92
Description
In my opinion, the method in here
tempestremap/src/GaussQuadrature.cpp
Line 23 in 60b99ca
can be optimized.
Those positions and weights in the reference element are fixed, the arrays dG and dW should not be allocated for every triangle. Some static arrays can be initialized with actual weights and gauss point locations, why alocate and set for every triangle ?