c/c++ sparse matrix library

Hello folks :)
I am working on a certain piece of code and i've been searching for a sparse matrix library for c/c++ for quite a while.
I have been freelance searching in the web and for older threads on this board, but the people usually searching for such a lib don't seem to know what they exactly want and don't really seem to have specific. On the other hand, I know perfectly well what I am looking for:

1) I need a sparse CRS matrix format for complex matrices.

2) I want to have an easy access like A->val(i,j) or simply A(i,j) or likewise to accesses that element.

3) Most importantly, I want to be able to set any entry (i,j) . This is were most libs i have tested so far fail.

I have tried Blitz++ (impressive library), but blitz only provides support for dense data structures. I don't really see how anybody can craft such a powerful lib and not include the most important HPC application of all ( large linear systems )

I have tried sparselib++ v1.7 ... i have rewritten the entire comprow_double.h for support complex values and then i realized that sparselib++ actually does not support the option to set the entry (j,k) if it wasn't previously a entry in your sparse structure. I want some library in which i can set a previously zero entry and some functionality should shift all the rowPtrs and colInds and values and resize all arrays... sparselib does not do this.

very important note: i do not really need any other functionality - if the lib has a method to transpose a dense matrix its fine, if it contains matrix*vector or some simple gauss-seidel relaxation or krylov method its cool, but i don't really need it.
I am aware that this is such a small piece of code, i could easily write it myself - however ... i don't really like to do things that have surely be done like 1242341 times before over again :)

it would be very cool, if someone would know such librarys and could suggest something :)

regards, sascha

Feb 11 '10 #1

Subscribe Post Reply

2701

Similar topics

Sparse matrices

by: George Sakkis | last post by:

Is there any sparse matrix package compatible with Numeric/Numarray ? Ideally, the implementation of a matrix (dense/sparse) should be transparent to the application. However the APIs of the only...

Python

sparse matrix class implementation

by: mariaczi | last post by:

Hi, I code class to storage sparse matrix row compressed and i have a problem with implements method to setVal and addVal. I will replace later this methods overloaded operator(). Please, can You...

C / C++

Matrix library where's the speed ?

by: Frank-O | last post by:

Hi , Recently I have been commited to the task of "translating" some complex statistical algorithms from Matlab to C++. The goal is to be three times as fast as matlab ( the latest) . I've...

C / C++

Finding Nonzero Elements in a Sparse Matrix

by: deLenn | last post by:

Hi, Does scipy have an equivalent to Matlab's 'find' function, to list the indices of all nonzero elements in a sparse matrix? Cheers.

Python

large and sparse matrices

by: mediratta | last post by:

Hi, I want to allocate memory for a large matrix, whose size will be around 2.5 million x 17000. Three fourth of its rows will have all zeroes, but it is not known which will be those rows. If I...

C / C++

Ax=b Jacobi 100000x100000 sparse matrix

by: hvmclrhu | last post by:

Hi I have a big problem. When we compile serial.c with gcc, I get this error program is generating the sparse matrix Segmentation fault I think ı have to use malloc() but I don't know how to...

C / C++

Sparse Matrix implemented as a doubly-linked list

by: adam.kleinbaum | last post by:

Hi there, I'm a novice C programmer working with a series of large (30,000 x 30,000) sparse matrices on a Linux system using the GCC compiler. To represent and store these matrices, I'd like to...

C / C++

C++ source code to generate sparse matrices.

by: Poseidon13 | last post by:

Hi everyone,

C / C++

implemeting sparse and non-sparse matrixes

by: ishakteyran | last post by:

hello to all.. i have a realy tough assignment which requires me to add, substract, multiply, and get inverse of non-sparse and sparse matrixes.. in a more clear way it wants me to to the...

Java

Sparse Matrix

by: Mark | last post by:

I'm trying to figure out the best way to go about doing this. I have a "map" for a game, which I'd like to store in a matrix. Some cells will be empty (NULL), and some will hold objects. I need a...

C / C++

How to turn on java script in a villaon keypad mobile phone

by: Charles Arthur | last post by:

How do i turn on java script on a villaon, callus and itel keypad mobile phone

Java

Merging data from multiple Excel files

by: ryjfgjl | last post by:

In our work, we often receive Excel tables with data in the same format. If we want to analyze these data, it can be difficult to analyze them because the data is spread across multiple Excel files...

Data Management

Navigating the Data Structures and Algorithms (DSA)

by: BarryA | last post by:

What are the essential steps and strategies outlined in the Data Structures and Algorithms (DSA) roadmap for aspiring data scientists? How can individuals effectively utilize this roadmap to progress...

Algorithms / Advanced Math

Is that possible of reading the .csv file in column wise and the column have different lengths ?

by: Sonnysonu | last post by:

This is the data of csv file 1 2 3 1 2 3 1 2 3 1 2 3 2 3 2 3 3 the lengths should be different i have to store the data by column-wise with in the specific length. suppose the i have to...

C / C++

How to build RAID in BIOS?

by: Hystou | last post by:

There are some requirements for setting up RAID: 1. The motherboard and BIOS support RAID configuration. 2. The motherboard has 2 or more available SATA protocol SSD/HDD slots (including MSATA, M.2...

Computer Hardware

What is ONU?

by: marktang | last post by:

ONU (Optical Network Unit) is one of the key components for providing high-speed Internet services. Its primary function is to act as an endpoint device located at the user's premises. However,...

General

Problem With Comparison Operator <=> in G++

by: Oralloy | last post by:

Hello folks, I am unable to find appropriate documentation on the type promotion of bit-fields when using the generalised comparison operator "<=>". The problem is that using the GNU compilers,...

C / C++

Maximizing Business Potential: The Nexus of Website Design and Digital Marketing

by: jinu1996 | last post by:

In today's digital age, having a compelling online presence is paramount for businesses aiming to thrive in a competitive landscape. At the heart of this digital strategy lies an intricately woven...

Online Marketing

Discussion: How does Zigbee compare with other wireless protocols in smart home applications?

by: tracyyun | last post by:

Dear forum friends, With the development of smart home technology, a variety of wireless communication protocols have appeared on the market, such as Zigbee, Z-Wave, Wi-Fi, Bluetooth, etc. Each...

General