42 c_vector<double, DIM+1> &rPhi,
43 c_matrix<double, DIM, DIM+1> &rGradPhi,
45 c_vector<double,3> &rU,
46 c_matrix<double,3,DIM> &rGradU ,
50 c_matrix<
double,3*(DIM+1),3*(DIM+1)> ret = zero_matrix<double>(3*(DIM+1),3*(DIM+1));
54 c_matrix<double, DIM+1, DIM+1> basis_outer_prod = outer_prod(rPhi, rPhi);
57 matrix_slice<c_matrix<double, 3*DIM+3, 3*DIM+3> >
58 slice00(ret, slice (0, 3, DIM+1), slice (0, 3, DIM+1));
59 slice00 = basis_outer_prod;
62 matrix_slice<c_matrix<double, 3*DIM+3, 3*DIM+3> >
63 slice11(ret, slice (1, 3, DIM+1), slice (1, 3, DIM+1));
64 slice11 = basis_outer_prod;
67 matrix_slice<c_matrix<double, 3*DIM+3, 3*DIM+3> >
68 slice22(ret, slice (2, 3, DIM+1), slice (2, 3, DIM+1));
69 slice22 = basis_outer_prod;
virtual c_matrix< double, 3 *(DIM+1), 3 *(DIM+1)> ComputeMatrixTerm(c_vector< double, DIM+1 > &rPhi, c_matrix< double, DIM, DIM+1 > &rGradPhi, ChastePoint< DIM > &rX, c_vector< double, 3 > &rU, c_matrix< double, 3, DIM > &rGradU, Element< DIM, DIM > *pElement)