Re: bug on compressed_matrix ?

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]