@@ -0,0 +1,877 @@
+#!/usr/bin/env python%0A# -*- coding: utf-8 -*-%0A%0Afrom math import pi%0Aimport pygmsh%0A%0Afrom helpers import compute_volume%0A%0A%0Adef test(lcar=0.05):%0A    geom = pygmsh.built_in.Geometry()%0A%0A    # Draw a square%0A    poly = geom.add_polygon(%5B%0A        %5B+0.5, +0.0, 0.0%5D,%0A        %5B+0.0, +0.5, 0.0%5D,%0A        %5B-0.5, +0.0, 0.0%5D,%0A        %5B+0.0, -0.5, 0.0%5D,%0A        %5D,%0A        lcar=lcar,%0A        )%0A%0A    axis = %5B0, 0, 1.0%5D%0A%0A    geom.extrude(%0A        poly,%0A        translation_axis=axis,%0A        rotation_axis=axis,%0A        point_on_axis=%5B0.0, 0.0, 0.0%5D,%0A        angle=0.5*pi,%0A        num_layers=5,%0A        recombine=True%0A        )%0A%0A    ref = 3.98156496566%0A    points, cells, _, _, _ = pygmsh.generate_mesh(geom)%0A    assert abs(compute_volume(points, cells) - ref) %3C 1.0e-2 * ref%0A    return points, cells%0A%0A%0Aif __name__ == '__main__':%0A    import meshio%0A    meshio.write('rotated_layers.vtu', *test())%0A

add test for the rotated layers extrusion