else
{
/* @whattodo */
- /* (*os) << ' ' << step; */
- /* for (unsigned int i=0; i<vectors.size(); ++i) */
- /* for (unsigned int j=0; j<vectors(i)->size(); ++j) */
- /* (*os) << ' ' << (*vectors(i))(j); */
- /* (*os) << std::endl; */
+ (*os) << ' ' << step;
+ for (unsigned int i=0; i<vectors.size(); ++i)
+ for (unsigned int j=0; j<vectors(i)->size(); ++j)
+ (*os) << ' ';// << (*vectors(i))(j);
+ (*os) << std::endl;
Assert ((false), ExcMessage ("This function is corrupt: @whattodo"));
}
return *this;
* Assuming that the stored vector is a cell vector, extract the given
* element from it.
*/
- virtual
+ virtual
double
get_cell_data_value (const unsigned int cell_number) const = 0;
for (unsigned int dataset=0; dataset<this->cell_data.size(); ++dataset)
{
+ // @whattodo - this next line does not cause a problem (yet), but probably will one day.
const double value
= this->cell_data[dataset]->get_cell_data_value (cell_and_index->second);
for (unsigned int q=0; q<n_q_points; ++q)
get_vector_element (const VectorType &vector,
const unsigned int cell_number)
{
- // @whattodo Note, there should be a way to get the value type from a VectorType
+ // @whattodo Note, there should be a way to get the value type from a VectorType!
// return vector[cell_number];
Assert ((false), ExcMessage ("This function is corrupt: @whattodo"));
}
= std::distance (this->triangulation->begin_active(),
typename Triangulation<DH::dimension,DH::space_dimension>::active_cell_iterator(cell_and_face->first));
+ // @whattodo - this next line does not cause a problem (yet), but probably will one day.
const double value
= this->cell_data[dataset]->get_cell_data_value (cell_number);
for (unsigned int q=0; q<n_q_points; ++q)
const unsigned int cell_number
= std::distance (this->triangulation->begin_active(),
typename Triangulation<DH::dimension,DH::space_dimension>::active_cell_iterator(*cell));
+
+ // @whattodo - this next line does not cause a problem (yet), but probably will one day.
const double value
= this->cell_data[dataset]->get_cell_data_value (cell_number);
switch (DH::dimension)