53 std::string pvj_file_name;
54 bool file_specified =
true;
63 file_specified =
false;
67 if (!file_specified || !junction_file.
Exists() )
70 WARNING(
"No Purkinje-Ventricular junction (.pvj) file found. Junctions must be specified manually.");
76 if (!junction_stream.good())
82 while (junction_stream.good())
84 std::string this_line;
85 getline(junction_stream, this_line);
87 if (this_line==
"" || this_line==
"\r")
89 if (junction_stream.eof())
98 std::stringstream line(this_line);
105 if (mpMixedDimensionMesh->rGetNodePermutation().size() != 0)
107 unsigned mapped_node_id = mpMixedDimensionMesh->rGetNodePermutation()[node_id];
109 mJunctionMap[mapped_node_id] = resistance;
113 mJunctionMap[node_id] = resistance;
127 assert(mpMixedDimensionMesh);
129 CableRangeAtNode cable_range = mpMixedDimensionMesh->GetCablesAtNode(pNode);
130 double total_cross_sectional_area = 0.0;
132 iter != cable_range.second;
137 total_cross_sectional_area += M_PI*cable_radius*cable_radius;
144 p_pvj_purkinje_stim->SetAppliedToPurkinjeCellModel();
145 p_pvj_ventricular_stim->SetVentricularCellModel(pCardiacCell);
146 p_pvj_ventricular_stim->SetPurkinjeCellModel(pPurkinjeCell);
147 p_pvj_purkinje_stim->SetVentricularCellModel(pCardiacCell);
148 p_pvj_purkinje_stim->SetPurkinjeCellModel(pPurkinjeCell);
151 boost::shared_ptr<MultiStimulus> p_multi_stim_ventricular(
new MultiStimulus);
152 p_multi_stim_ventricular->AddStimulus(p_pvj_ventricular_stim);
156 boost::shared_ptr<MultiStimulus> p_multi_stim_purkinje(
new MultiStimulus);
157 p_multi_stim_purkinje->AddStimulus(p_pvj_purkinje_stim);