def METHOD_NAME(self, comm, dtype):
    """
    Setup mesh and tacs assembler for problem we will be testing.
    """
    # Create the stiffness object
    props = constitutive.MaterialProperties(rho=2570.0, E=70e9, nu=0.3, ys=350e6)
    stiff = constitutive.PlaneStressConstitutive(props, t=0.1, tNum=0)
    # Set up the basis function
    model = elements.LinearElasticity2D(stiff)
    basis = elements.LinearQuadBasis()
    elem = elements.Element2D(model, basis)
    # Allocate the TACSCreator object
    vars_per_node = model.getVarsPerNode()
    creator = TACS.Creator(comm, vars_per_node)
    if comm.rank == 0:
        num_elems = nx * ny
        num_nodes = (nx + 1) * (ny + 1)
        # discretize plate
        x = np.linspace(0, Lx, nx + 1, dtype)
        y = np.linspace(0, Ly, ny + 1, dtype)
        xyz = np.zeros([nx + 1, ny + 1, 3], dtype)
        xyz[:, :, 0], xyz[:, :, 1] = np.meshgrid(x, y, indexing="ij")
        node_ids = np.arange(num_nodes).reshape(nx + 1, ny + 1)
        # Set connectivity for each element
        conn = []
        for i in range(nx):
            for j in range(ny):
                conn.append(
                    [
                        node_ids[i, j],
                        node_ids[i + 1, j],
                        node_ids[i, j + 1],
                        node_ids[i + 1, j + 1],
                    ]
                )
        conn = np.array(conn, dtype=np.intc).flatten()
        ptr = np.arange(0, 4 * num_elems + 1, 4, dtype=np.intc)
        comp_ids = np.zeros(num_elems, dtype=np.intc)
        creator.setGlobalConnectivity(num_nodes, ptr, conn, comp_ids)
        # Set up the boundary conditions (fixed at bottom left corner, roller on top right)
        bcnodes = np.array([node_ids[0, 0], node_ids[0, -1]], dtype=np.intc)
        bcvars = np.array([0, 1, 0], dtype=np.intc)
        # Set the boundary condition pointers
        bcptr = np.array([0, 2, 3], dtype=np.intc)
        creator.setBoundaryConditions(bcnodes, bcptr, bcvars)
        # Set the node locations
        creator.setNodes(xyz.flatten())
    # Set the elements for each (only one) component
    element_list = [elem]
    creator.setElements(element_list)
    # Create the tacs assembler object
    assembler = creator.createTACS()
    return assembler