efficient matrix block assignment

Lisandro Dalcin dalcinl at gmail.com
Thu Aug 6 16:27:33 CDT 2009


>> Barry Smith wrote:
>>> MatSetValuesBlock() you can put in a rectangular set of blocks.
>>>

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.




>
>>
>> Thanks for the clarifications,
>> Umut
>
>
>
> --
> What most experimenters take for granted before they begin their experiments
> is infinitely more interesting than any results to which their experiments
> lead.
> -- Norbert Wiener
>



-- 
Lisandro Dalcín
---------------
Centro Internacional de Métodos Computacionales en Ingeniería (CIMEC)
Instituto de Desarrollo Tecnológico para la Industria Química (INTEC)
Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET)
PTLC - Güemes 3450, (3000) Santa Fe, Argentina
Tel/Fax: +54-(0)342-451.1594


More information about the petsc-users mailing list