[petsc-users] Transpose of Block Matrix with aij type

[Karl Rupp]

Hi,

why would you expect that the transpose of a 3x8 matrix is not a 8x3-matrix?

Best regards,
Karli


On 04/23/2013 03:27 PM, Joon hee Choi wrote:
> Hello,
>
> I tried to get transpose of block matrix(just with aij type), but the result was not a block matrix. For example,
>
> A =
> 1 2 3 4 | 4 3 2 1
> 2 3 4 5 | 5 4 3 2
> 3 4 5 6 | 6 5 4 3
>
>
> AT(expected) =
> 1 2 3 4
> 2 3 4 5
> 3 4 5 6
> -------
> 4 3 2 1
> 5 4 3 2
> 6 5 4 3
>
>
> AT(result) =
> 1 2 3
> 2 3 4
> 3 4 5
> 4 5 6
> 4 5 6
> 3 4 5
> 2 3 4
> 1 2 3
>
> If someone knows about this problem, please let me know it.
>
> Thank you
>