evaluate sum_ij (qL_i (G_ij(q_G)) qR_j) ,
where qL=(dX^1/dalpha,dX^2/dalpha ,dzeta/dalpha) and qR=(dX^1/dbeta,dX^2/dbeta ,dzeta/dbeta) and
dzeta_dalpha then known to be either 0 of ds and dtheta and 1 for dzeta
evaluate sum_k sum_ij (qL_i d/dq^k(G_ij(q_G)) qR_j) q_vec^k, k=1,2,3
where qL=(dX^1/dalpha,dX^2/dalpha ,dzeta/dalpha) and qR=(dX^1/dbeta,dX^2/dbeta ,dzeta/dbeta) and
dzeta_dalpha then known to be either 0 of ds and dtheta and 1 for dzeta
Evaluate value and n derivatives at x of all basis functions with support in local cell
derivs[i,j] = (d/dx)^i B_j(x) for 0 <= i <= n and jmin <= j <= jmin+degree