Re: [ublas] bug on compressed_matrix ?

Subject: Re: [ublas] bug on compressed_matrix ?
From: Hoang Giang Bui (hgbk2008_at_[hidden])
Date: 2013-01-16 14:06:44

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
> ublas_at_[hidden]
> http://lists.boost.org/mailman/listinfo.cgi/ublas
> Sent to: hgbk2008_at_[hidden]

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