In my opinion, the method in here
|
void GaussQuadrature::GetPoints( |
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 ?
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 ?