#include <BidomainNeumannSurfaceTermAssembler.hpp>
Inherits AbstractFeSurfaceIntegralAssembler< ELEMENT_DIM, SPACE_DIM, 2 >.

Public Member Functions | |
| BidomainNeumannSurfaceTermAssembler (AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM > *pMesh, BoundaryConditionsContainer< ELEMENT_DIM, SPACE_DIM, 2 > *pBoundaryConditions, unsigned numQuadPoints=2) | |
Protected Member Functions | |
| virtual c_vector< double, 2 *ELEMENT_DIM > | ComputeVectorSurfaceTerm (const BoundaryElement< ELEMENT_DIM-1, SPACE_DIM > &rSurfaceElement, c_vector< double, ELEMENT_DIM > &rPhi, ChastePoint< SPACE_DIM > &rX) |
Assembler which sets up the surface integral integrals for the bidomain equations, assuming that the boundary conditions are written: div(sigma_i grad phi_i) . n = g1 and div(sigma_e grad phi_e) dot n = g2.
These are not 'natural' boundary conditions for the para-elliptic bidomain equations (natural BCs for the second
Hence we don't use the NaturalNeumannSurfaceTermAssembler and have a special class here. It means that any BCs specified for bidomain and put in a BoundaryConditionsContainer should be for div(sigma_i grad phi_i) . n and div(sigma_e grad phi_e) . n.
Definition at line 49 of file BidomainNeumannSurfaceTermAssembler.hpp.
| BidomainNeumannSurfaceTermAssembler< ELEMENT_DIM, SPACE_DIM >::BidomainNeumannSurfaceTermAssembler | ( | AbstractTetrahedralMesh< ELEMENT_DIM, SPACE_DIM > * | pMesh, | |
| BoundaryConditionsContainer< ELEMENT_DIM, SPACE_DIM, 2 > * | pBoundaryConditions, | |||
| unsigned | numQuadPoints = 2 | |||
| ) | [inline] |
Constructor
| pMesh | The mesh | |
| pBoundaryConditions | The boundary conditions container | |
| numQuadPoints | Number of quad points (per dimension) to use |
Definition at line 76 of file BidomainNeumannSurfaceTermAssembler.hpp.
| c_vector< double, 2 *ELEMENT_DIM > BidomainNeumannSurfaceTermAssembler< ELEMENT_DIM, SPACE_DIM >::ComputeVectorSurfaceTerm | ( | const BoundaryElement< ELEMENT_DIM-1, SPACE_DIM > & | rSurfaceElement, | |
| c_vector< double, ELEMENT_DIM > & | rPhi, | |||
| ChastePoint< SPACE_DIM > & | rX | |||
| ) | [inline, protected, virtual] |
This method returns the vector to be added to full vector for a given Gauss point in BoundaryElement, ie, essentially the INTEGRAND in the boundary integral part of the definition of the vector. The arguments are the bases, x and current solution computed at the Gauss point.
| rSurfaceElement | the element which is being considered. | |
| rPhi | The basis functions, rPhi(i) = phi_i, i=1..numBases | |
| rX | The point in space |
Reimplemented from AbstractFeSurfaceIntegralAssembler< ELEMENT_DIM, SPACE_DIM, 2 >.
Definition at line 87 of file BidomainNeumannSurfaceTermAssembler.hpp.
References AbstractFeSurfaceIntegralAssembler< ELEMENT_DIM, SPACE_DIM, 2 >::mpBoundaryConditions.
1.6.3