void ConstructGridVC(Triang2D & G,
Geom & destGeom,
G2 const& srcG,
Geom2 const& srcGeom,
VertexMap & VertexCorr, // maps src::handles to
CellMap & CellCorr) // complex2d::handles
{
// construct combinatorial grid
ConstructGrid0(G,srcG,VertexCorr,CellCorr);
// copy geometry
destGeom.set_grid(G);
typedef typename grid_types<G2>::VertexIterator src_vertex_it;
typedef typename grid_types<Triang2D>::Vertex Vertex;
for(src_vertex_it src_v = srcG.FirstVertex(); ! src_v.IsDone(); ++src_v) {
Vertex V(G,VertexCorr(src_v.handle()));
assign_point(destGeom.coord(V) , srcGeom.coord(*src_v));
}
}
sums()
{
assign_point(sum_m, T());
assign_point(sum_ms, T());
}