MTL: (no subject)

[jaime.silva]

Hello,
i think that is correct is a bug. I debug the code and i think if you change
this in the mtl.h file it corrects the bug. The method is called:

template <class Matrix, class VecX, class VecZ>
inline void
mult_symm__(const Matrix& A, const VecX& x, VecZ& z, row_tag)

and the diff file:
line 1590
<     if (A.is_upper()){
<       tmp += *j * x[j.row()];
---
>     if (A.is_upper()) {
>       tmp += *j * x[j.column()];
line 1599
<       tmp += *j * x[j.row()];
---
>       tmp += *j * x[j.column()];
line 1604
<       tmp += *j * x[j.row()];
---
>       tmp += *j * x[j.column()];
 
Regards,

Jaime Silva

*****************************************************************
IPC - Institute for Polymer And Composites

Department of Polymer Engineering,

Universidade do Minho,Campus de Azurém
4800 - 058 Guimarães,Portugal



Phone [phone]
Fax: [phone]

URL: http://ipc.uminho.pt/
_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/