blocked index sets

Barry Smith bsmith at mcs.anl.gov
Wed Aug 27 13:08:25 CDT 2008


On Aug 27, 2008, at 12:13 PM, Matthew Knepley wrote:

> On Wed, Aug 27, 2008 at 12:06 PM, Lisandro Dalcin  
> <dalcinl at gmail.com> wrote:
>> So, Do all us agree my proposed fix should be pushed? I'll wait for
>> Matt comments/complaints...
>
> I complain that IS is a fake parallel object.

   In what way is IS a fake parallel object? Maybe we can resolve
your concern about it being a fake parallel object.

    Barry


> However, if
> GetSize/GetLocalSize already
> do this, then yes we should change the ISBlock version as well.
>
>   Matt
>
>> On Wed, Aug 27, 2008 at 1:13 PM, Barry Smith <bsmith at mcs.anl.gov>  
>> wrote:
>>>
>>> On Aug 27, 2008, at 10:23 AM, Matthew Knepley wrote:
>>>
>>>> There is no concept of global for IS. They are purely serial. AO  
>>>> is the
>>>> only
>>>> global construct with indices.
>>>
>>>  This is kind of true, and maybe used to be completely true. But  
>>> IS does
>>> have a communicator and that communicator can be MPI_COMM_WORLD or
>>> any parallel communicator.  In other words the IS is evolving to  
>>> be an
>>> object
>>> that can be parallel in the same sense as vecs or mats
>>>
>>>  There are already ISGetSize() and ISGetLocalSize() so it sure  
>>> makes sense
>>> to have the same paradgm for the ISGetBlockSize().
>>>
>>>
>>>  Barry
>>>
>>> Originally IS had no parallel concept, then we added the ISGetSize/ 
>>> LocalSize
>>> but forgot to do it for the ISBlock...
>>>
>>>
>>>>
>>>>
>>>> Matt
>>>>
>>>> On Wed, Aug 27, 2008 at 10:09 AM, Lisandro Dalcin <dalcinl at gmail.com 
>>>> >
>>>> wrote:
>>>>>
>>>>> I believe we have to review the interface of ISBlock. Currently,
>>>>> ISBlockGetSize() return the number of LOCAL block indices. This  
>>>>> is not
>>>>> consistent with other naming conventions for getting local and  
>>>>> glocal
>>>>> sizes. I propose to change this to the following
>>>>>
>>>>> 1) change: ISBlockGetSize() returns the number global blocks
>>>>> 2) addition:  ISBlockGetLocalSize() return the number of local  
>>>>> blocks
>>>>>
>>>>> Comments?
>>>>>
>>>>>
>>>>> --
>>>>> 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
>>>>>
>>>>>
>>>>
>>>>
>>>>
>>>> --
>>>> 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
>>
>>
>
>
>
> -- 
> 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
>




More information about the petsc-dev mailing list