From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 120826 invoked by alias); 19 Jan 2016 16:43:10 -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 120814 invoked by uid 89); 19 Jan 2016 16:43:10 -0000 Authentication-Results: sourceware.org; auth=none X-Virus-Found: No X-Spam-SWARE-Status: No, score=2.2 required=5.0 tests=AWL,BAYES_20,FREEMAIL_FROM,RCVD_IN_DNSWL_LOW,SPF_PASS autolearn=ham version=3.3.2 spammy=Storage, Row, ccs, crs X-HELO: mail-ig0-f178.google.com Received: from mail-ig0-f178.google.com (HELO mail-ig0-f178.google.com) (209.85.213.178) by sourceware.org (qpsmtpd/0.93/v0.84-503-g423c35a) with (AES128-GCM-SHA256 encrypted) ESMTPS; Tue, 19 Jan 2016 16:43:08 +0000 Received: by mail-ig0-f178.google.com with SMTP id t15so81611907igr.0 for ; Tue, 19 Jan 2016 08:43:08 -0800 (PST) X-Google-DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=1e100.net; s=20130820; h=x-gm-message-state:mime-version:in-reply-to:references:date :message-id:subject:from:to:content-type; bh=jD6CGGspWEep2D+aDxjqXlZbclbOESDkFupVK4EXJWY=; b=E61iua/Y6uwmVSasWZ43ogNjf5+iY/PmxMRIxif+Y8qx3r+bPqO3xQfhrfEzns2lqB Y/GzXBw+14Ml3kJRNRetqTn06Nhig7DE6tFG3b/QUyjGpvYdZPJXEWK0IimmidtcxWip jB0SgQ63uYnbotUpG/JtOyGaWMpsL7Qmi32CRirh+kt8mReYLPyv0HY1IZXSXk9ANBXG 1GC44vQiUYVjaQgi+Vit87S+Ok3SMiZYL6IR4Ow4PLiDIIjwqslwz6KrsTGn4aK2q78u VHixWkZfzOJinp0aWOqUWrOaFfPYJnoAfDSJXLcRmqc2MMBupDtBjq8JENQiAo849RZE VJwA== X-Gm-Message-State: AG10YOT+RiaecQjh8YVq6S7gKvVuxMDxVs0wAaFI/EPt++lsdgDuUpEmkucBPvwKN20M4QoUWEmrLHfNIsBsNQ== MIME-Version: 1.0 X-Received: by 10.50.79.196 with SMTP id l4mr17025890igx.11.1453221786462; Tue, 19 Jan 2016 08:43:06 -0800 (PST) Received: by 10.79.118.156 with HTTP; Tue, 19 Jan 2016 08:43:06 -0800 (PST) In-Reply-To: References: Date: Tue, 19 Jan 2016 16:43:00 -0000 Message-ID: Subject: Re: Sparse matrix extension From: Alexis Tantet To: gsl-discuss@sourceware.org Content-Type: text/plain; charset=UTF-8 X-SW-Source: 2016-q1/txt/msg00000.txt.bz2 Dear GSLers, As a scientist rather than a developer, I have developed an extension of the sparse matrix module (CRS, I/O, manipulation, see below), which I have tested. These modifications conserve the structure of the original module and be useful for a large number of sparse matrices users. I'm not familiar with the contributing process here. My repository can be found there: https://github.com/atantet/gslsp Unfortunately, I did not know of the gsl.git repository and I forked it froml: https://github.com/drjerry/gslsp , which seems to be a bit older than gsl.git. How can I push/merge to gsl.git ? Should it be as an update or another extension? Is it necessary to adapt to the newest version of the code ? Best regards, Alexis Tantet CHANGES.md: Extension of the sparse matrix module of GSL =================================== Introduction ------------ Usages of sparse matrices are numerous in scientific computing. When numerical linear algebra problems become large, sparse matrices become necessary to avoid memory overload and unnecessary computations, at the cost of element access and matrix construction. As a result, most large scale linear solvers or eigen solvers perform on sparse matrices. Fortunately, a very useful sparse matrix module has recently been introduced to GSL. However, important features are still lacking, such has Compressed Row Storage (CRS) matrices, input/output functions and other matrix properties and manipulation functions. This new version attempts to address this, conserving the original structure of the module and conventions. Major changes ------------- * Add CRS format and update functions manipulating compressed matrices : - additional flag GSL_SPMATRIX_CRS and macro GSLSP_ISMATRIX ( gsl_spmatrix.h ) - additional members innerSize and outerSize used to iterate matrix elements ( gsl_spmatrix.h ) - rename some variables for coherence ( gsl_spmatrix.h , *.c ) - update all functions on compressed matrices ( *.c ) * Allow to sum duplicate elements when compressing ( spcompress.c ) : - modify gsl_spmatrix_compress - add gsl_spmatrix_sum_duplicate * CCS <-> CRS and fast transpose inplace in spswap.c : - add gsl_spmatrix_switch_major - add gsl_spmatrix_transpose * Add printing and scanning functions in spio.c : - add gsl_spmatrix_fprintf - add gsl_spmatrix_fscanf * Add manipulation functions in spmanip.c (particularly useful for Markov chain transition matrices) : - add gsl_spmatrix_get_rowsum : get vector of sum over row elements - add gsl_spmatrix_get_colsum : get vector of sum over column elements - add gsl_spmatrix_div_rows : divide all elements of each row by a vector element - add gsl_spmatrix_div_cols : divide all elements of each column by a vector element * Add test functions in atprop.c : - add gsl_spmatrix_gt_elements : greater than test for each matrix element - add gsl_spmatrix_ge_elements : greater or equal than test for each matrix element - add gsl_spmatrix_lt_elements : lower than test for each matrix element - add gsl_spmatrix_le_elements : lower or equal than test for each matrix element - add gsl_spmatrix_any : test if any non-zero element in matrix Other minor changes have been made, such as error tests. test.c has also been updated to test new features.