65 c_vector<double, DIM+1>& rPhi,
66 c_matrix<double, DIM, DIM+1>& rGradPhi,
68 c_vector<double, 1>& rU,
69 c_matrix<double, 1, DIM>& rGradU,
72 c_matrix<double, DIM, DIM> pde_diffusion_term = this->mpEllipticPde->ComputeDiffusionTerm(rX);
75 if (mLinearInUCoeffInSourceTerm != 0)
77 return prod( trans(rGradPhi), c_matrix<double, DIM, DIM+1>(prod(pde_diffusion_term, rGradPhi)) )
78 - mLinearInUCoeffInSourceTerm * outer_prod(rPhi,rPhi);
82 return prod( trans(rGradPhi), c_matrix<double, DIM, DIM+1>(prod(pde_diffusion_term, rGradPhi)) );
virtual c_matrix< double, 1 *(DIM+1), 1 *(DIM+1)> ComputeMatrixTerm(c_vector< double, DIM+1 > &rPhi, c_matrix< double, DIM, DIM+1 > &rGradPhi, ChastePoint< DIM > &rX, c_vector< double, 1 > &rU, c_matrix< double, 1, DIM > &rGradU, Element< DIM, DIM > *pElement)