Many thanks Daniel for your help. On Tue, Nov 29, 2016 at 2:25 PM, Daniel Arndt < [email protected]> wrote:
> Anup, > > Looking at the resulting mesh after using grid_y_transform, the mesh is > still rectangular. > Therefore, it is not surprising that make_periodicity_constraints works > without offset, but not with: > > Tensor<1, dim> offset; > offset[dim-1] = width*tan_theta*scale; > GridTools::collect_periodic_faces(dof_handler_eta, 1, 11, 0, > periodicity_vector1, offset); > > However, the above code works if you turn grid_y_transform into > > Point<dim> pt_out = pt_in; > pt_out[2] += x*tan_theta; > return pt_out; > > Best, > Daniel > > -- > 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. > -- 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.
