46 c_vector<double, ELEMENT_DIM+1> &rPhi,
47 c_matrix<double, SPACE_DIM, ELEMENT_DIM+1> &rGradPhi,
49 c_vector<double,2> &rU,
50 c_matrix<double, 2, SPACE_DIM> &rGradU ,
53 c_vector<double,1> empty_u;
54 c_matrix<double,1,SPACE_DIM> empty_grad_u;
56 c_matrix<double,ELEMENT_DIM+1,ELEMENT_DIM+1> normal_monodomain_mat
57 = mMonodomainAssembler.ComputeMatrixTerm(rPhi,rGradPhi,rX,empty_u,empty_grad_u,pElement);
59 c_matrix<
double,2*(ELEMENT_DIM+1),2*(ELEMENT_DIM+1)> ret
60 = zero_matrix<double>(2*(ELEMENT_DIM+1),2*(ELEMENT_DIM+1));
63 matrix_slice<c_matrix<double, 2*ELEMENT_DIM+2, 2*ELEMENT_DIM+2> >
64 slice00(ret, slice(0, 2, ELEMENT_DIM+1), slice(0, 2, ELEMENT_DIM+1));
65 slice00 = normal_monodomain_mat;
MonodomainPurkinjeVolumeAssembler(AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM > *pMesh, MonodomainTissue< ELEMENT_DIM, SPACE_DIM > *pTissue)
c_matrix< double, 2 *(ELEMENT_DIM+1), 2 *(ELEMENT_DIM+1)> ComputeMatrixTerm(c_vector< double, ELEMENT_DIM+1 > &rPhi, c_matrix< double, SPACE_DIM, ELEMENT_DIM+1 > &rGradPhi, ChastePoint< SPACE_DIM > &rX, c_vector< double, 2 > &rU, c_matrix< double, 2, SPACE_DIM > &rGradU, Element< ELEMENT_DIM, SPACE_DIM > *pElement)