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.