46 "Circular stimulus.");
64 const TiXmlElement* pXml)
69 m_nq = pField->GetTotPoints();
78 const TiXmlElement *pXmlparameter;
80 pXmlparameter = pXml->FirstChildElement(
"p_x1");
81 m_px1 = atof(pXmlparameter->GetText());
83 pXmlparameter = pXml->FirstChildElement(
"p_y1");
84 m_py1 = atof(pXmlparameter->GetText());
86 pXmlparameter = pXml->FirstChildElement(
"p_z1");
87 m_pz1 = atof(pXmlparameter->GetText());
89 pXmlparameter = pXml->FirstChildElement(
"p_r1");
90 m_pr1 = atof(pXmlparameter->GetText());
92 pXmlparameter = pXml->FirstChildElement(
"p_is");
93 m_pis = atof(pXmlparameter->GetText());
95 pXmlparameter = pXml->FirstChildElement(
"p_strength");
115 int dim =
m_field->GetCoordim(0);
118 int nq =
m_field->GetNpoints();
132 for(
int j=0; j<nq; j++)
134 outarray[0][j] += v_amp
141 for(
int j=0; j<nq; j++)
143 outarray[0][j] += v_amp
152 for(
int j=0; j<nq; j++)
154 outarray[0][j] += v_amp
void Initialise()
Initialise the protocol storage and set initial conditions.
virtual void v_Update(Array< OneD, Array< OneD, NekDouble > > &outarray, const NekDouble time)
NekDouble m_chiCapMembrane
std::vector< std::pair< std::string, std::string > > SummaryList
ProtocolSharedPtr m_Protocol
Stimulus protocol to apply.
LibUtilities::SessionReaderSharedPtr m_session
Session.
boost::shared_ptr< SessionReader > SessionReaderSharedPtr
StimulusCirc(const LibUtilities::SessionReaderSharedPtr &pSession, const MultiRegions::ExpListSharedPtr &pField, const TiXmlElement *pXml)
int m_nq
Number of physical points.
boost::shared_ptr< ExpList > ExpListSharedPtr
Shared pointer to an ExpList object.
virtual void v_GenerateSummary(SolverUtils::SummaryList &s)
static std::string className
Name of class.
MultiRegions::ExpListSharedPtr m_field
Transmembrane potential field from PDE system.
StimulusFactory & GetStimulusFactory()
static StimulusSharedPtr create(const LibUtilities::SessionReaderSharedPtr &pSession, const MultiRegions::ExpListSharedPtr &pField, const TiXmlElement *pXml)
Creates an instance of this class.
tKey RegisterCreatorFunction(tKey idKey, CreatorFunction classCreator, tDescription pDesc="")
Register a class with the factory.