In [ ]:

    

import sys
sys.path.append('../')

%matplotlib notebook


    

In [ ]:

    

from pycrest.mesh import *
from pycrest.refinement import *
from pycrest.visualize import *
from math import pi


    

In [ ]:

    

# A simple triangulation of the unit square, using two triangles
vertices = [
    (0.0, 0.0),
    (1.0, 0.0),
    (1.0, 1.0),
    (0.0, 1.0)
]
elements = [
    (3, 0, 1),
    (1, 2, 3)
]
initial_mesh = Mesh2d(vertices, elements)


    

In [ ]:

    

# Configure parameters for refinement
tolerance = 0.05
mesh = bisect_to_tolerance(initial_mesh, tolerance)


    

In [ ]:

    

plot_triangulation(mesh)


    

In [ ]:

    

mesh.num_vertices


    

In [ ]:

    

mesh.num_elements


    

In [ ]:

    

(coarse, fine) = threshold(initial_mesh, tolerance, [0, 2], [3 * pi/4, pi / 2])


    

In [ ]:

    

plot_triangulation(coarse)


    

In [ ]:

    

plot_triangulation(fine)


    

In [ ]:

    

fine.num_elements


    

In [ ]:

    

coarse.num_elements


    

In [ ]: