From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 14458 invoked by alias); 30 Jan 2014 23:40:34 -0000 Mailing-List: contact gsl-discuss-help@sourceware.org; run by ezmlm Precedence: bulk List-Id: List-Subscribe: List-Archive: List-Post: List-Help: , Sender: gsl-discuss-owner@sourceware.org Received: (qmail 14448 invoked by uid 89); 30 Jan 2014 23:40:34 -0000 Authentication-Results: sourceware.org; auth=none X-Virus-Found: No X-Spam-SWARE-Status: No, score=-1.9 required=5.0 tests=AWL,BAYES_00,RP_MATCHES_RCVD,SPF_PASS autolearn=ham version=3.3.2 X-HELO: ipmx5.colorado.edu Received: from ipmx5.colorado.edu (HELO ipmx5.colorado.edu) (128.138.128.235) by sourceware.org (qpsmtpd/0.93/v0.84-503-g423c35a) with ESMTP; Thu, 30 Jan 2014 23:40:33 +0000 From: Patrick Alken Received: from bonanza.ngdc.noaa.gov ([140.172.179.41]) by smtp.colorado.edu with ESMTP/TLS/DHE-RSA-AES256-SHA; 30 Jan 2014 16:40:31 -0700 Message-ID: <52EAE2EF.1020205@colorado.edu> Date: Thu, 30 Jan 2014 23:40:00 -0000 User-Agent: Mozilla/5.0 (X11; Linux x86_64; rv:24.0) Gecko/20100101 Thunderbird/24.2.0 MIME-Version: 1.0 To: "gsl-discuss@sourceware.org" Subject: Sparse matrix linear solver added to GSL Content-Type: text/plain; charset=ISO-8859-1; format=flowed Content-Transfer-Encoding: 7bit X-SW-Source: 2014-q1/txt/msg00024.txt.bz2 Hello all, The sparse matrix routines are in pretty good shape now and they are all located in a branch called 'sparse' on the gsl git. Furthermore I've added a sparse linear solver (for A x = b) based on the Generalized Residual Minimum Method (GMRES) iterative solver. GMRES is considered a very good/robust solver for a wide class of matrices, though it does have convergence issues on some matrices. In these cases it helps to precondition the system but I have not currently implemented any preconditioners. Maybe one day I'll get around to that. Everything is documented, with an example of solving the 1D Poisson equation sparse system, and I've added lots of automated tests so things are working well. Anyone who uses sparse matrices in their work I'd be interested in any feedback. I plan to merge the sparse branch into the master branch in the next week or so. Patrick