CellwiseSourcePde.cpp
00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036 #include "CellwiseSourcePde.hpp"
00037 #include "ApoptoticCellProperty.hpp"
00038 #include "Exception.hpp"
00039
00040 template<unsigned DIM>
00041 CellwiseSourcePde<DIM>::CellwiseSourcePde(AbstractCellPopulation<DIM,DIM>& rCellPopulation, double coefficient)
00042 : mrCellPopulation(rCellPopulation),
00043 mCoefficient(coefficient)
00044 {
00045 }
00046
00047 template<unsigned DIM>
00048 const AbstractCellPopulation<DIM,DIM>& CellwiseSourcePde<DIM>::rGetCellPopulation() const
00049 {
00050 return mrCellPopulation;
00051 }
00052
00053 template<unsigned DIM>
00054 double CellwiseSourcePde<DIM>::GetCoefficient() const
00055 {
00056 return mCoefficient;
00057 }
00058
00059 template<unsigned DIM>
00060 double CellwiseSourcePde<DIM>::ComputeConstantInUSourceTerm(const ChastePoint<DIM>& rX, Element<DIM,DIM>* pElement)
00061 {
00062 return 0.0;
00063 }
00064
00065 template<unsigned DIM>
00066 double CellwiseSourcePde<DIM>::ComputeLinearInUCoeffInSourceTerm(const ChastePoint<DIM>& rX, Element<DIM,DIM>* pElement)
00067 {
00068 NEVER_REACHED;
00069 return 0.0;
00070 }
00071
00072 template<unsigned DIM>
00073 double CellwiseSourcePde<DIM>::ComputeLinearInUCoeffInSourceTermAtNode(const Node<DIM>& rNode)
00074 {
00075 double coefficient = 0.0;
00076
00077 if (mrCellPopulation.IsCellAttachedToLocationIndex(rNode.GetIndex()))
00078 {
00079 CellPtr p_cell = mrCellPopulation.GetCellUsingLocationIndex(rNode.GetIndex());
00080
00081 bool cell_is_apoptotic = p_cell->HasCellProperty<ApoptoticCellProperty>();
00082
00083 if (!cell_is_apoptotic)
00084 {
00085 coefficient = mCoefficient;
00086 }
00087 }
00088
00089 return coefficient;
00090 }
00091
00092 template<unsigned DIM>
00093 c_matrix<double,DIM,DIM> CellwiseSourcePde<DIM>::ComputeDiffusionTerm(const ChastePoint<DIM>& rX)
00094 {
00095 return identity_matrix<double>(DIM);
00096 }
00097
00099
00101
00102 template class CellwiseSourcePde<1>;
00103 template class CellwiseSourcePde<2>;
00104 template class CellwiseSourcePde<3>;
00105
00106
00107 #include "SerializationExportWrapperForCpp.hpp"
00108 EXPORT_TEMPLATE_CLASS_SAME_DIMS(CellwiseSourcePde)