c_vector< double, 1 *(ELEM_DIM+1)> ComputeVectorTerm(c_vector< double, ELEM_DIM+1 > &rPhi, c_matrix< double, SPACE_DIM, ELEM_DIM+1 > &rGradPhi, ChastePoint< SPACE_DIM > &rX, c_vector< double, 1 > &rU, c_matrix< double, 1, SPACE_DIM > &rGradU, Element< ELEM_DIM, SPACE_DIM > *pElement)