(no subject)

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 +=3D *j * x[j.row()];
---
>     if (A.is_upper()) {
>       tmp +=3D *j * x[j.column()];
line 1599
<       tmp +=3D *j * x[j.row()];
---
>       tmp +=3D *j * x[j.column()];
line 1604 =

<       tmp +=3D *j * x[j.row()];
---
>       tmp +=3D *j * x[j.column()];
 =

Regards,

Jaime Silva

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

Department of Polymer Engineering,

Universidade do Minho,Campus de Azur=E9m
4800 - 058 Guimar=E3es,Portugal



Phone (+351) 253 510 337/338
Fax: (+351) 253510339

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