The application of the discretization statements on the respective elements such as vertices can be written in the following manner:

R_FVM = sum<VE>(_v)[O(_1, _v)*A/d*sum<EV>(_e)[O(_e, _1) * q]]; vertex_iterator vit = C.vertex_begin(); while (vit.valid()) { R = R_FVM(*vit); std::cout << *vit << R << std::endl; }

Using an iterator, all vertices of a given cell complex can be traversed. The traversal loop, terminates because the iterator becomes invalid. The residual expression `R` is calculated for each vertex. Afterwards, the vertex and the respective residuum are listed.

Michael 2008-01-16