[petsc-dev] DAGetMatrix

Barry Smith bsmith at mcs.anl.gov
Tue Sep 28 21:18:34 CDT 2010


On Sep 28, 2010, at 5:11 PM, Lisandro Dalcin wrote:

> On 28 September 2010 17:39, Barry Smith <bsmith at mcs.anl.gov> wrote:
>> 
>> 
>> 
>> On Sep 27, 2010, at 2:52 PM, Lisandro Dalcin wrote:
>> 
>>> DAGetMatrix(da,&mat,"shell") or DAGetMatrix(da,&mat,"python") or any
>>> other matrix type the DA does not know generate an error. Could this
>>> be changed to returning a matrix with the proper local&global sizes
>>> and type set ?
>> 
>>  Done but not tested, let me know if you have any problems. Note that the local sizes are set but for unknown matrix times I am not sure if they will always automatically compute the global sizes. Please report problems.
>> 
> 
> I'm wondering why you do not set the global sizes ...

   No particular reason I guess. It never did since the individual Mat's always figured it out.

   I have changed it so DAGetMat() now sets the global sizes also. After you hg pull.

   Barry

> 
> 
> -- 
> Lisandro Dalcin
> ---------------
> CIMEC (INTEC/CONICET-UNL)
> Predio CONICET-Santa Fe
> Colectora RN 168 Km 472, Paraje El Pozo
> Tel: +54-342-4511594 (ext 1011)
> Tel/Fax: +54-342-4511169




More information about the petsc-dev mailing list