(I hope this isn't a RTFM. I saw something in the manual about hermitian 
matrices related to Lapack, so it only seemed to be dense.)

Is there provision for hermitian *sparse* matrices in Julia, in the sense 
that if A is "declared" hermitian, only its lower (or upper) triangle is 
stored and A*x does the right thing? That's always been missing from 
Matlab, and would be an important feature to have. In smooth optimization, 
most operators end up being symmetric (e.g., Hessians, or KKT matrices).

Thanks.

Reply via email to