20 std::ofstream file(filename.c_str());
23 for (
size_t i=0; i<_names.size(); i++)
24 file << _names[i] <<
"\t";
31 typedef typename GridView::template Codim<0>::Iterator CellIterator;
32 CellIterator it = _gv.template begin<0>();
33 CellIterator end = _gv.template end<0>();
34 for (; it != end; ++it)
36 int i = _is.index(*it);
38 assert (i == counter++);
40 assert(it->geometry().corners() == 2);
41 const FieldVector<ctype,dimworld>& left = it->geometry().corner(0);
42 const FieldVector<ctype,dimworld>& right = it->geometry().corner(1);
43 assert(left[0] < right[0]);
45 writeRow(file, left, _data[2*i]);
46 writeRow(file, right, _data[2*i+1]);
51 typedef typename GridView::template Codim<dimworld>::Iterator VertexIterator;
52 VertexIterator it = _gv.template begin<dimworld>();
53 VertexIterator end = _gv.template end<dimworld>();
54 for (; it != end; ++it) {
57 writeRow(file, it->geometry().corner(0), _data[_is.index(*it)]);
68 const FieldVector<ctype,dimworld>& position,
69 const std::vector<float> & data)
const
71 assert (data.size() == _names.size());
73 file << position <<
"\t";
75 for (
size_t j=0; j<data.size(); j++)
76 file << data[j] <<
"\t";
85 template<
class Gr
idView>
86 template<
class DataContainer>
90 assert((t == cellData && _is.size(0) == data.size())
92 _names.push_back(name);
100 int shift = (t==vertexData ? 1 : 0);
101 for (
size_t i=0; i<_is.size(0); i++)
103 _data[c++].push_back(data[i]);
104 _data[c++].push_back(data[i+shift]);
110 for (
size_t i=0; i<_is.size(dimworld); i++)
111 _data[i].push_back(data[i]);