Jt

where J is a sparse matrix, A is diagonal and Jt is the transpost of J.

In Ublas documentation I can't found any reference to a

diagonal_matrix type even

if it seems to me that such a type exists. If there is support for

this type how should it be used?

Does also anyone know what's the state of SuperLU bindings?

Thanks in advance for any help you will provide.

Regards,

Luca