bool SaveGridHierarchyTransformed(MultiGrid& mg, const char* filename,
number offset)
{
PROFILE_FUNC_GROUP("grid");
// cast away constness
SubsetHandler& sh = mg.get_hierarchy_handler();
APosition aPos;
// uses auto-attach
Grid::AttachmentAccessor<Vertex, APosition> aaPos(mg, aPos, true);
// copy the existing position to aPos. We take care of dimension differences.
// Note: if the method was implemented for domains, this could be implemented
// in a nicer way.
if(mg.has_vertex_attachment(aPosition))
ConvertMathVectorAttachmentValues<Vertex>(mg, aPosition, aPos);
else if(mg.has_vertex_attachment(aPosition2))
ConvertMathVectorAttachmentValues<Vertex>(mg, aPosition2, aPos);
else if(mg.has_vertex_attachment(aPosition1))
ConvertMathVectorAttachmentValues<Vertex>(mg, aPosition1, aPos);
// iterate through all vertices and apply an offset depending on their level.
for(size_t lvl = 0; lvl < mg.num_levels(); ++lvl){
for(VertexIterator iter = mg.begin<Vertex>(lvl);
iter != mg.end<Vertex>(lvl); ++iter)
{
aaPos[*iter].z() += (number)lvl * offset;
}
}
// finally save the grid
bool writeSuccess = SaveGridToFile(mg, sh, filename, aPos);
// clean up
mg.detach_from_vertices(aPos);
return writeSuccess;
}
void TestSubdivision(const char* fileIn, const char* fileOut, int numRefs)
{
PROFILE_FUNC_GROUP("grid");
//todo: Callbacks have to make sure that their attachment is accessible in the grid.
// even if they were initialized before the attachment was attached to the grid.
MultiGrid mg;
SubsetHandler sh(mg);
SubdivisionLoopProjector<APosition> refCallback(mg, aPosition, aPosition);
GlobalMultiGridRefiner ref(mg, &refCallback);
if(LoadGridFromFile(mg, sh, fileIn)){
for(int lvl = 0; lvl < numRefs; ++lvl){
ref.refine();
}
ProjectToLimitPLoop(mg, aPosition, aPosition);
SaveGridToFile(mg, mg.get_hierarchy_handler(), fileOut);
}
else{
UG_LOG("Load failed. aborting...\n");
}
}
bool SaveGridHierarchy(MultiGrid& mg, const char* filename)
{
PROFILE_FUNC_GROUP("grid");
return SaveGridToFile(mg, mg.get_hierarchy_handler(), filename);
}