Teko: Enable dof-split for block matrix

[maxfirmbach]

Question

@trilinos/teko @MalachiTimothyPhillips @cgcgcg

Hello,

currently it's possible in Teko to split an operator based on DOFs e.g. fluid matrix goes in and with

  <Parameter name="Strided Blocking" type="string" value="2 1"/>

I'm able to generate a block operator splitting velocity and pressure DOFs.

For certain examples in our code I would need to apply this procedure on a block in a block matrix e.g. I have 3 blocks resembling structure, thermo and a scalar transport fields respectively and I want to split the scalar transport problem into [1 1].

E.g. [3 1 2] -> [3 1 [1 1]]

Is this something that could be easily implemented?

Thanks!

[jhux2]

I have 3 blocks resembling structure, thermo and a scalar transport fields respectively and I want to split the scalar transport problem into [1 1].

Is the 3-block matrix stored as a single point matrix, 3 point matrices, or something else?

[maxfirmbach]

It's a Thyra::PhysicallyBlockedLinearOp with A = [structure, thermo, scatra]. So A(0,0) would give you the structure stiffness matrix and so forth. The A(i,j) blocks are just simple Crs_Matrices.

[maxfirmbach]

I would like to be able to do a dof-split e.g. on A(2,2).

[jhux2]

@cgcgcg @MalachiTimothyPhillips Can you advise Max?

[MalachiTimothyPhillips]

@maxfirmbach Is the idea that you'd pair the block-recursive structure with the hierarchical block Gauss-Seidel? I want to make sure I understand the use-case.

[maxfirmbach]

@MalachiTimothyPhillips I think that would be something we could do, yes.