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