Seg fault when using defintiion of quadrature points when 'c' direction is specified larger than 'a' and b' directions
I have noted some segmentation faults when I specify the quadrature points on prisms of the form:
<E COMPOSITE=“C[0]” BASISTYPE=“Modified_A,Modified_A,Modified_B” NUMMODES="3,3,3" POINTSTYPE="GaussLobattoLegendre,GaussLobattoLegendre,GaussRadauMAlpha1Beta0" NUMPOINTS="5,5,6" FIELDS="u,v,w" />
instead of
<E COMPOSITE=“C[0]" BASISTYPE="Modified_A,Modified_A,Modified_B" NUMMODES="3,3,3" POINTSTYPE="GaussLobattoLegendre,GaussLobattoLegendre,GaussRadauMAlpha1Beta0" NUMPOINTS="6,6,5" FIELDS="u,v,w" />
From the core dump it would appear to be arising the evaluation of the derivative factors within collections. I am not sure you should be able to specify the 'c' drection to be larger than the 'a' and 'b' directions but if not we should have an ASSERT in to stop it.