Abstract. In this paper we present a computationally efficient algorithm utilizing a fully or seminonlocal graph Laplacian for solving a wide range of learning problems in binary data classification and image processing. In their recent work [Multiscale Model. Simul., 10 (2012) 1. Introduction. This work develops a fast algorithm for a recent variational method in a graph setting. The method is inspired by diffuse interface models that have been used in a variety of problems, such as those in fluid dynamics and materials science. We consider data represented as nodes in a weighted graph, and each edge is assigned a numerical value describing the similarity between the nodes. In spectral graph theory, this approach is successfully used to perform various learning tasks in imaging and data clustering. The standard techniques of the theory are thoroughly described in [13,45], and the graph Laplacian, which is discussed in more detail in section 2.2, is introduced as one of the fundamental concepts. In imaging, spectral methods are often used in image segmentation applications, as shown in [54,33,14].We are particularly interested in nonlocal total variation methods, as they are a link between spectral graph theory and diffuse interface models and thus can be used as a motivation for our algorithm. These methods are used in numerous image processing applications. They were initially developed as methods for image denoising [9,29] but were successfully applied to many other image processing problems such as inpainting and reconstruction in [30,63,49], image deblurring in [40], and manifold processing in [18].Bertozzi and Flenner introduce a graph-based model based on the Ginzburg-Landau func-