Provide diag - Githubissues

Algebraic-Programming / ALP

Home of ALP/GraphBLAS and ALP/Pregel, featuring shared- and distributed-memory auto-parallelisation of linear algebraic and vertex-centric programs. Soon with more to come!

Apache License 2.0

24 stars 4 forks source link

Provide diag #238

Open anyzelman opened 10 months ago

anyzelman commented 10 months ago

These are ideally provided as an algorithm (though may performance-wise require to be a primitive instead):

RC diag( Matrix< T > &out, const Matrix< T > &in );
RC diag( Matrix< T > &out, const Vector< T > &in );
RC diag( Vector< T > &out, const Matrix< T > &in );

anyzelman commented 5 months ago

Crossref issue #233 -- the matrix factory presently defines a diag function, but found that an iterator-only approach in conjunction with buildMatrixUnique cannot effect parallel I/O. One solution for that is to instead rely on the functionality proposed under this issue, such that input values for the diagonal can be ingested into a Vector, which is then in turn promoted to a diagonal matrix.