Function screen 

pub fn screen<T: Scalar>(
    ncells: usize,
    cell_type: ReferenceCellType,
) -> SingleElementGrid<T, CiarletElement<T, IdentityMap, T>>

Create a grid of a square screen

Create a grid of the square [0,1]^2. The input ncells is the number of cells along each side of the square.