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.
