5. 1. 5 Application

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