41 c_vector<double, DIM+1> &rPhi,
42 c_matrix<double, DIM, DIM+1> &rGradPhi,
44 c_vector<double,2> &rU,
45 c_matrix<double,2,DIM> &rGradU ,
48 c_matrix<
double,2*(DIM+1),2*(DIM+1)> ret = zero_matrix<double>(2*(DIM+1), 2*(DIM+1));
52 c_matrix<double, DIM+1, DIM+1> basis_outer_prod = outer_prod(rPhi, rPhi);
55 matrix_slice<c_matrix<double, 2*DIM+2, 2*DIM+2> >
56 slice00(ret, slice (0, 2, DIM+1), slice (0, 2, DIM+1));
57 slice00 = basis_outer_prod;
c_matrix< double, 2 *(DIM+1), 2 *(DIM+1)> ComputeMatrixTerm(c_vector< double, DIM+1 > &rPhi, c_matrix< double, DIM, DIM+1 > &rGradPhi, ChastePoint< DIM > &rX, c_vector< double, 2 > &rU, c_matrix< double, 2, DIM > &rGradU, Element< DIM, DIM > *pElement)