def METHOD_NAME(
    self, assembler, force_vec, dv_pert_vec, ans_pert_vec, xpts_pert_vec
):
    """
    Setup user-defined vectors for analysis and fd/cs sensitivity verification
    """
    local_num_nodes = assembler.getNumOwnedNodes()
    vars_per_node = assembler.getVarsPerNode()
    # The nodes have been distributed across processors now
    # Let's find which nodes this processor owns
    xpts0 = assembler.createNodeVec()
    assembler.getNodes(xpts0)
    xpts0_array = xpts0.getArray()
    # Split node vector into numpy arrays for easier parsing of vectors
    local_xyz = xpts0_array.reshape(local_num_nodes, 3)
    local_x, local_y, local_z = local_xyz[:, 0], local_xyz[:, 1], local_xyz[:, 2]
    # No need to populate force vector, since we're using inertial forces
    # Create temporary dv vec for doing fd/cs
    dv_pert_array = dv_pert_vec.getArray()
    dv_pert_array[:] = 1.0
    # Create temporary state variable vec for doing fd/cs
    ans_pert_array = ans_pert_vec.getArray()
    # Define perturbation array that uniformly moves all nodes on right edge of plate to the right
    ans_pert_array = ans_pert_array.reshape(local_num_nodes, vars_per_node)
    ans_pert_array[local_y == 0.0, 0] = 1.0
    # Define perturbation array that uniformly moves all nodes on right edge of plate to the right
    xpts_pert_array = xpts_pert_vec.getArray()
    xpts_pert_array = xpts_pert_array.reshape(local_num_nodes, 3)
    # Define perturbation array that uniformly moves all nodes on right edge of plate to the right
    xpts_pert_array[local_x == Lx, 0] = 1.0
    return