Re: bug on compressed_matrix ?
Hoang Giang Bui <[email protected]>
| Newsgroups | gmane.comp.lib.boost.ublas |
|---|---|
| Message-ID | <[email protected]> |
On 01/16/13 20:00, Thomas Klimpel wrote: > Hoang Giang Bui wrote: >> Yes, I need the correct csr format as standard since we rely on ublas >> for computing element stiffness matrix and global stiffness matrix stem >> from structural problem. Furthermore, I need to take the transpose of >> ublas matrix to feed in a subroutine to compute the restriction operator >> in our multilevel solver. This subroutine uses the compressed data for >> speed purpose. This behaviour causes segmentation false error in my code >> since it read a negative number of entries in the last row. > Can you try calling the "complete_index1_data()" member function and see whether it fixes your problem? (You still have to take care that "index2_data()" and "value_data()" return references to arrays whose actual size is "meaningless".) > > Regards, > Thomas > _______________________________________________ > ublas mailing list > [email protected] > http://lists.boost.org/mailman/listinfo.cgi/ublas > Sent to: [email protected] Yes, it does the trick. And I can handle the index2_data and value_data without problem if index1_data is correct. Thank you very much. But this confusion could be avoided if it is mentioned somewhere in the documentation. Btw, is there any performance penalty if complete_index1_data is called? BR Giang Bui _______________________________________________ ublas mailing list [email protected] http://lists.boost.org/mailman/listinfo.cgi/ublas Sent to: [email protected]