Boost logo

Ublas :

Subject: Re: [ublas] simple question on CSR matrix
From: Nasos Iliopoulos (nasos_i_at_[hidden])
Date: 2010-08-12 10:01:35


you can use the find_element(i,j) member function, that returns a pointer to the element when it exists and 0 otherwise.

Also for filling sparse matrices, I would suggest taking a look at Gunter's wiki:

and playing a bit with the various sparse matrix types.

I find that for FEA, it is usually more efficient to use a generalized_vector_of_vector to assemble the matrix and then copy that into a compressed_matrix, to send it to a solver. Of course this process heavily depends on problem size.


> From: rrossi_at_[hidden]
> To: ublas_at_[hidden]
> Date: Thu, 12 Aug 2010 14:21:58 +0200
> Subject: [ublas] simple question on CSR matrix
> Dear list,
> i have a CSR matrix, of which i preallocated the zeros doing
> push_back
> is there a way to ensure that if i access
> A(i,j) and i,j is not in the graph an error is thrown? i mean... i would
> prefer an error (even segfault) than ublas silently attempting to resize
> the matrix
> any hint is very welcome
> thank you
> Riccardo
> --
> ________________________________________________________________
> Riccardo Rossi, Ph.D, Civil Engineer
> member of the Kratos Group:
> Centro Internacional de Métodos Numéricos en Ingeniería (CIMNE)
> Universidad Politécnica de Cataluña (UPC)
> Edificio C-1, campus Norte UPC
> Gran Capitan, s/n
> 08034 Barcelona, España
> Tel. (+34) 93 401 73 99
> ________________________________________________________________
> Los datos de carácter personal contenidos en el mensaje, se registrarán
> en un fichero para facilitar la gestión de las comunicaciones de CIMNE.
> Se pueden ejercitar los derechos de acceso, rectificación, cancelación y
> oposición por escrito, dirigiéndose a nuestras oficinas de CIMNE, Gran
> Capitán s/n, Edificio C1 - Campus Norte UPC, 08034 Barcelona, España.
> Les dades de caràcter personal contingudes en aquest missatge es
> registraran en un fitxer per facilitar la gestió de les comunicacions
> del CIMNE. Es poden exercir els drets d'accés, rectificació,
> cancel·lació i oposició, per escrit a les nostres oficines del CIMNE,
> Gran Capità s/n, Edifici C1, Campus Nord UPC, 08034 Barcelona, Espanya.
> All personal data contained in this mail will be processed
> confidentially and stored in a file property of CIMNE in order to manage
> corporate communications. You may exercise the right of access,
> rectification, deletion and objection by letter sent to CIMNE, Gran
> Capitán, Edificio C1 - Campus Norte UPC, 08034 Barcelona, Spain.
> _______________________________________________
> ublas mailing list
> ublas_at_[hidden]
> Sent to: nasos_i_at_[hidden]