Chaste Commit::baa90ac2819b962188b7562f2326be23c47859a7
|
#include <ElectrodesStimulusFactory.hpp>
Public Member Functions | |
ElectrodesStimulusFactory (std::vector< std::pair< AbstractChasteRegion< DIM > *, AbstractChasteRegion< DIM > * > > &rElectrodePairs, std::vector< double > &rStimulusMagnitudes, std::vector< double > &rDurations, std::vector< double > &rPeriods, std::vector< double > &rStarts, std::vector< double > &rEnds) | |
boost::shared_ptr< AbstractStimulusFunction > | CreateStimulusForNode (Node< DIM > *pNode) |
~ElectrodesStimulusFactory () | |
void | SetCompatibleExtracellularStimulus () |
void | GroundSecondElectrode (bool grounded) |
Public Member Functions inherited from AbstractStimulusFactory< DIM > | |
virtual boost::shared_ptr< AbstractStimulusFunction > | CreateStimulusForNode (Node< ELEMENT_DIM > *pNode) |
virtual unsigned | GetNumberOfCells () |
AbstractStimulusFactory () | |
virtual | ~AbstractStimulusFactory () |
void | SetMesh (AbstractTetrahedralMesh< ELEMENT_DIM, ELEMENT_DIM > *pMesh) |
AbstractTetrahedralMesh< ELEMENT_DIM, ELEMENT_DIM > * | GetMesh () |
std::vector< AbstractChasteRegion< ELEMENT_DIM > * > | GetRegionsToBeGrounded () |
Protected Types | |
typedef LinearBasisFunction< DIM > | BasisFunction |
Protected Member Functions | |
void | CheckForElectrodesIntersection () |
double | ComputeElectrodeTotalFlux (AbstractChasteRegion< DIM > *pRegion, double stimulusMagnitude) |
Protected Attributes | |
std::vector< std::pair< AbstractChasteRegion< DIM > *, AbstractChasteRegion< DIM > * > > & | mrElectrodePairs |
std::vector< double > & | mrMagnitudes |
std::vector< double > & | mrDurations |
std::vector< double > & | mrPeriods |
std::vector< double > & | mrStarts |
std::vector< double > & | mrEnds |
std::vector< double > | mMagnitudesElectrode1 |
std::vector< double > | mMagnitudesElectrode2 |
bool | mGroundSecondElectrode |
Protected Attributes inherited from AbstractStimulusFactory< DIM > | |
AbstractTetrahedralMesh< ELEMENT_DIM, ELEMENT_DIM > * | mpMesh |
std::vector< AbstractChasteRegion< ELEMENT_DIM > * > | mGroundedRegions |
Friends | |
class | TestStimulusFactory |
This class implements the specification of tow electrodes with a RegularStimulus applied to them. It makes sure the compatibility conditions is verified by scaling the magnitude of the second electrode according to the flux flowing thorugh each of them.
User needs to pass in a vector of pairs of electrodes and, for each pair, the stimulation parameters. This class will then implement the CreateStimulusForNode accordingly by creating a RegularStimulus object for each node. Each node belonging to any "first electrode" will have an extracellular volume stimulus equals to the value that is passed in, while each node belonging to any "second elecrode" will have an extracellular volume stimulus corrected by the flux (proportional to electrode volume) in such a way that compatibility conditions of the extended bidomain equations are followed.
It is possible to ground all second electrodes by calling GroundSecondElectrode. In this case, the member variable mGroundSecondElectrode will have size > 0 and other classes can behave accordingly (e.g., problem class will set Dirichlet boundary conditions accordingly).
Definition at line 67 of file ElectrodesStimulusFactory.hpp.
|
protected |
typedef for basis function
Definition at line 121 of file ElectrodesStimulusFactory.hpp.
ElectrodesStimulusFactory< DIM >::ElectrodesStimulusFactory | ( | std::vector< std::pair< AbstractChasteRegion< DIM > *, AbstractChasteRegion< DIM > * > > & | rElectrodePairs, |
std::vector< double > & | rStimulusMagnitudes, | ||
std::vector< double > & | rDurations, | ||
std::vector< double > & | rPeriods, | ||
std::vector< double > & | rStarts, | ||
std::vector< double > & | rEnds | ||
) |
Constructor. Electrodes and stimulation parameters need to be passed in.
rElectrodePairs | the pairs of electrodes |
rStimulusMagnitudes | the magnitudes of the stimuli (microA / cm^3). First electrode will have magnitude value and second electrode will have -magnitude (before being corrected to ensure equal flux). |
rDurations | the duration of each stimulus (ms) |
rPeriods | the period of each stimulus (ms) |
rStarts | the start time of each stimulus (ms). |
rEnds | the end of each stimulation (ms) |
Definition at line 45 of file ElectrodesStimulusFactory.cpp.
References EXCEPTION, ElectrodesStimulusFactory< DIM >::mMagnitudesElectrode1, ElectrodesStimulusFactory< DIM >::mMagnitudesElectrode2, and ElectrodesStimulusFactory< DIM >::mrMagnitudes.
ElectrodesStimulusFactory< DIM >::~ElectrodesStimulusFactory | ( | ) |
Destructor
Definition at line 73 of file ElectrodesStimulusFactory.cpp.
|
protected |
Helper method to check if any of the electrodes intersects with the other. it needs to check whether any element that conatins any node of electrode 1 does not conatin also nodes of eledtrode 2
Definition at line 78 of file ElectrodesStimulusFactory.cpp.
References PetscTools::Barrier(), and EXCEPTION.
|
protected |
Helper method to compute electrode flux. This method is the main functionality of this class.
It essentially computes the contribution of each electrode to the RHS vector of the linear system. The loops over elements and gauss points mimicks the ones in AssembleOnElement. Interpolation is done for the value of the magnitude of the electrode.
NOTE that this method assumes the use of a 2nd order Gauss quadrature rule (the default value of Chaste FE assemblers).
pRegion | the electrode |
stimulusMagnitude | the stimulus magnitude |
Definition at line 186 of file ElectrodesStimulusFactory.cpp.
References AbstractChasteRegion< SPACE_DIM >::DoesContain(), AbstractElement< ELEMENT_DIM, SPACE_DIM >::GetIndex(), AbstractElement< ELEMENT_DIM, SPACE_DIM >::GetNumNodes(), GaussianQuadratureRule< ELEMENT_DIM >::GetNumQuadPoints(), GaussianQuadratureRule< ELEMENT_DIM >::GetWeight(), and GaussianQuadratureRule< ELEMENT_DIM >::rGetQuadPoint().
boost::shared_ptr< AbstractStimulusFunction > ElectrodesStimulusFactory< DIM >::CreateStimulusForNode | ( | Node< DIM > * | pNode | ) |
pNode | the node for which to create the stimulus |
Definition at line 162 of file ElectrodesStimulusFactory.cpp.
References Node< SPACE_DIM >::GetPoint().
void ElectrodesStimulusFactory< DIM >::GroundSecondElectrode | ( | bool | grounded | ) |
Allow the user to ground the second electrode
grounded | : true if second electrode is grounded |
Definition at line 117 of file ElectrodesStimulusFactory.cpp.
|
virtual |
This method checks whether the users wanted a grounded electrode. If so, it fills in mGroundedRegions accordingly.
If not, it calls ComputeElectrodeTotalFlux and scales the magnitudes of the electrodes stimuli accordingly
Reimplemented from AbstractStimulusFactory< DIM >.
Definition at line 123 of file ElectrodesStimulusFactory.cpp.
References PetscTools::ReplicateException().
|
friend |
Definition at line 71 of file ElectrodesStimulusFactory.hpp.
|
protected |
Whether the second electrode is grounded
Definition at line 111 of file ElectrodesStimulusFactory.hpp.
|
protected |
Used to store temporarily magnitudes of electrode 1 for correction
Definition at line 102 of file ElectrodesStimulusFactory.hpp.
Referenced by ElectrodesStimulusFactory< DIM >::ElectrodesStimulusFactory().
|
protected |
Used to store temporarily magnitudes of electrode 2 for correction
Definition at line 107 of file ElectrodesStimulusFactory.hpp.
Referenced by ElectrodesStimulusFactory< DIM >::ElectrodesStimulusFactory().
|
protected |
Vector of stimuli durations (ms). Must be the same size as mrElectrodePairs.
Definition at line 85 of file ElectrodesStimulusFactory.hpp.
|
protected |
Vector of pairs, each pair representing a pair of electrodes.
Definition at line 76 of file ElectrodesStimulusFactory.hpp.
|
protected |
Vector of stimuli end times (ms). Must be the same size as mrElectrodePairs.
Definition at line 97 of file ElectrodesStimulusFactory.hpp.
|
protected |
Vector of stimuli magnitudes (microA / cm^3). Must be the same size as mrElectrodePairs.
Definition at line 81 of file ElectrodesStimulusFactory.hpp.
Referenced by ElectrodesStimulusFactory< DIM >::ElectrodesStimulusFactory().
|
protected |
Vector of stimuli periods (ms). Must be the same size as mrElectrodePairs.
Definition at line 89 of file ElectrodesStimulusFactory.hpp.
|
protected |
Vector of stimuli start times (ms). Must be the same size as mrElectrodePairs.
Definition at line 93 of file ElectrodesStimulusFactory.hpp.