artorg-unibe-ch / spline_mesher

Meshing tool for homogenised finite elements based on HR-pQCT images
https://artorg-unibe-ch.github.io/spline_mesher/
MIT License
0 stars 0 forks source link

move to upper scope when testing is done #48

Closed github-actions[bot] closed 1 year ago

github-actions[bot] commented 1 year ago

https://github.com/simoneponcioni/spline-mesher/blob/8a72294712a3793c7b88d0048a5fbe1cbd087ec6/src/spline_mesher/quad-refinement.py#L818


            connected_points.append(tag_unique)
        points_in_surf.append(connected_points)

    line_tags_intersurf = trab_refinement.create_intersurface_connection(points_in_surf)

    for _iter in range(1, len(points_in_surf)):
        (
            lines_lower_dict,
            lines_upper_dict,
            lines_intersurf_dict,
        ) = trab_refinement.create_line_dict(
            lines_lower_surf=line_tags[_iter - 1],
            lines_upper_surf=line_tags[_iter],
            lines_intersurf=line_tags_intersurf[_iter - 1],
        )

        # fmt: off
        start_time = time.time()
        print("Finding closed curve loops...")
        curve_loops_tags = fcc.find_closed_curve_loops(lines_lower_dict, lines_upper_dict, lines_intersurf_dict)
        end_time = time.time()
        exec_time = end_time - start_time
        print(f"Time to find closed curve loops: {exec_time:.2f} seconds")
        # fmt: on
        gmsh.model.occ.synchronize()
        intersurface_surfaces = []
        for cl in curve_loops_tags:
            curve_loop = trab_refinement.add_curve_loop(cl)
            surf = trab_refinement.gmsh_add_surface([curve_loop])
            intersurface_surfaces.append(surf)

    # TODO: move to upper scope when testing is done
    gmsh.model.occ.synchronize()
    for intersurf_tags in line_tags_intersurf:
        for line in intersurf_tags:
            gmsh.model.mesh.setTransfiniteCurve(line, 2, "Progression", 1.0)
    for surf in intersurface_surfaces:
        gmsh.model.mesh.setTransfiniteSurface(surf)

    gmsh.option.setNumber("Mesh.RecombineAll", 1)
    gmsh.option.setNumber("Mesh.RecombinationAlgorithm", 1)
    gmsh.option.setNumber("Mesh.Recombine3DLevel", 2)
    gmsh.option.setNumber("Mesh.ElementOrder", 2)

    # * 10. Create 2D mesh
    gmsh.model.occ.synchronize()
    gmsh.model.mesh.generate(2)
github-actions[bot] commented 1 year ago

Closed in 74e1d5a33575593185f19f7bfa837e7ade1f39fd