search

trilinos / Trilinos

Primary repository for the Trilinos Project
https://trilinos.org/
Other
1.22k stars 568 forks source link

Tpetra: Distributed dense matrix type with Layout right #13565

Open Adrian-Diaz opened 2 weeks ago

Adrian-Diaz commented 2 weeks ago

@trilinos/tpetra

Greetings Everyone,

Is there a Tpetra type (or other package with a similar distributed type) that can serve as a dense matrix using Layout Right? In theory the multivector type could work but it seems to me thats hard coded to always use Layout Left. I also saw CRS Matrix which allows Layout Right with a 1D storage view, but that class seems to carry the overhead of the graph matrix for a dense application. Thanks in advance for any additional info on the matter.

Adrian Diaz

cgcgcg commented 1 week ago

@trilinos/tpetra

csiefer2 commented 1 week ago

You are correct. MultiVector is hardcoded to LayoutLeft. AFAIK, we don't have a LayoutRight distributed object.

Adrian-Diaz commented 1 week ago

Hmmm, would anything break if Multivector or a renamed copy of it was compiled with layout right for that argument? At least anything major that comes to mind, naturally there's lots of details and it may not even compile.

csiefer2 commented 4 days ago

Vector extraction would break in a "LayoutRightMultiVector"... that'll break some of the Kyrlov solvers, for instance.

That said, if there's a strong driver for a layout right multivector, we could definitely consider adding one.