23 for(
const auto& region :
poly.mesh.regions())
25 for(const auto& face : region.second->faces())
27 if(face.second->periodicFacePair.second)
29 if( int(region.second->regionID)==face.second->periodicFacePair.second->regionIDs.first
30 ||int(region.second->regionID)==face.second->periodicFacePair.second->regionIDs.second)
32 const auto& grain(poly.grains.at(region.second->regionID));
33 std::cout<<
"MeshFace "<<face.first<<
" checking if periodicShift is a LatticeVector "<<std::flush;
34 const auto lv(grain->latticeVector(face.second->periodicFacePair.first));
35 std::cout<<
"LatticeVector"<<
"="<<lv.base().transpose()<<std::endl;