//! Idea of implementation: define an ordering for the nodes (with respect to the axis of anisotropy)
//! and then use a merge sort algorithm to renumerate the nodes along the lines which are parallel to this axis.
//! to do so, we introduce a 'NodeList' class.
//! everything is handeled by the PermutationManager
//! everything is handled by the PermutationManager
// type of nodes:
...
...
@@ -151,7 +151,7 @@ public:
// the information about the axis of anisotropy is in NodeList, since this list contains a methode "less_equal", which describes the ordering of the nodes.
// the information about the axis of anisotropy is in NodeList, since this list contains a method "less_equal", which describes the ordering of the nodes.
// this ordering depends on the axis of anisotropy
template<classGridViewImp>
classPermutationManager{
...
...
@@ -430,9 +430,9 @@ public:
for(inti=0;i<node_list_.size();++i)
{
std::cout<<"Original node["<<get_inverse_permuted_index(i)<<"] = ("<<node_list_[i].get_global_coordinate()<<") recieved permuted index "<<i<<std::endl;
std::cout<<"Original node["<<get_inverse_permuted_index(i)<<"] = ("<<node_list_[i].get_global_coordinate()<<") received permuted index "<<i<<std::endl;
// alternative:
//std::cout << "Original node[" << get_inverse_permuted_index( i ) << "] = (" << original_node_list_[get_inverse_permuted_index( i )].get_global_coordinate() << ") recieved permuted index " << i << std::endl;
//std::cout << "Original node[" << get_inverse_permuted_index( i ) << "] = (" << original_node_list_[get_inverse_permuted_index( i )].get_global_coordinate() << ") received permuted index " << i << std::endl;