Dear all, i was trying to use Dealii to solve a 3-dimensional elastic-plastic question on a cubic by applying periodic boundary condition like in step-45 (v. 8.5.0). std::vector<GridTools::PeriodicFacePair<typename DoFHandler<dim>:: cell_iterator> > periodicity_vector; GridTools::collect_periodic_faces(dof_handler, 1, 2, 2, periodicity_vector);//1,2,2: the first two numbers refer to the cell faces at the boundary, the third number to the direction DoFTools::make_periodicity_constraints<DoFHandler<dim> >(periodicity_vector, constraint_matrix);
But i got different results as i have created the mesh using GridGenerator::subdivided_hyper_rectangle (triangulation, rep, p1, p2);//rep is the repetition method 1 to get the mesh: p1 (0, 0, 0) p2(1, 1, 1) rep=(8, 8, 8) no global refinement method 2 to get the mesh: p1 (0, 0, 0) p2(1, 1, 1) rep=(1, 1, 1) global refinement = 3 I just wonder why there come two different results when these two meshes are actually the same. Thank you guys for the help! Chuyi -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. For more options, visit https://groups.google.com/d/optout.
