efficient matrix block assignment
Umut Tabak
u.tabak at tudelft.nl
Thu Aug 6 16:37:35 CDT 2009
> Indeed. As barry pointed out, you can put a rectangular m x n *set* of
> *blocks*...
>
> For example:
>
> array([[0, 0, 1, 1, 2, 2],
> [0, 0, 1, 1, 2, 2],
> [3, 3, 4, 4, 5, 5],
> [3, 3, 4, 4, 5, 5]])
>
> is a 2x3 set of blocks where each (square) block has bs=2, i.e. each
> block is 2x2 square.
>
Was trying to read Barry's email again because I missed that
*set* of
*blocks*...
and to sketch on paper some matrices, right on time arrived your email
;) , thanks for the example.
Best,
Umut
More information about the petsc-users
mailing list