[petsc-dev] Fwd: A few remaining issues with Mat block sizes

[Satish Balay]

pushed now.

satish

On Tue, 10 Apr 2012, Barry Smith wrote:

> 
>    Satish,
> 
>      Please add this patch.
> 
>     Thanks
> 
>    Barry
> 
> 
> Begin forwarded message:
> 
> > From: Lisandro Dalcin <dalcinl at gmail.com>
> > Date: April 10, 2012 12:43:39 PM CDT
> > To: For users of the development version of PETSc <petsc-dev at mcs.anl.gov>
> > Subject: Re: [petsc-dev] A few remaining issues with Mat block sizes
> > Reply-To: For users of the development version of PETSc <petsc-dev at mcs.anl.gov>
> > 
> > On 10 April 2012 17:58, Barry Smith <bsmith at mcs.anl.gov> wrote:
> >> 
> >>   Lisandro,
> >> 
> >>    Thanks, I have incorporated these and a few other blocksize fixes.
> >> 
> > 
> > Many thanks, Barry! However, I still think you are missing this one (I
> > found the issue after running petsc4py testsuite). Pasted below and
> > also attached. Please review and push.
> > 
> > 
> > diff -r 2d1c2b03ef3a src/mat/impls/aij/seq/aij.c
> > --- a/src/mat/impls/aij/seq/aij.c	Tue Apr 10 09:57:31 2012 -0500
> > +++ b/src/mat/impls/aij/seq/aij.c	Tue Apr 10 20:36:48 2012 +0300
> > @@ -3477,9 +3477,12 @@
> >   PetscErrorCode ierr;
> > 
> >   PetscFunctionBegin;
> > +  if (Ii[0]) SETERRQ1(PETSC_COMM_SELF,PETSC_ERR_ARG_OUTOFRANGE,
> > "Ii[0] must be 0 it is %D", Ii[0]);
> > +
> > +  ierr = PetscLayoutSetUp(B->rmap);CHKERRQ(ierr);
> > +  ierr = PetscLayoutSetUp(B->cmap);CHKERRQ(ierr);
> > +
> >   ierr = MatGetSize(B, &m, &n);CHKERRQ(ierr);
> > -
> > -  if (Ii[0]) SETERRQ1(PETSC_COMM_SELF,PETSC_ERR_ARG_OUTOFRANGE,
> > "Ii[0] must be 0 it is %D", Ii[0]);
> >   ierr = PetscMalloc((m+1) * sizeof(PetscInt), &nnz);CHKERRQ(ierr);
> >   for(i = 0; i < m; i++) {
> >     nz     = Ii[i+1]- Ii[i];
> > 
> > 
> > -- 
> > Lisandro Dalcin
> > ---------------
> > CIMEC (INTEC/CONICET-UNL)
> > Predio CONICET-Santa Fe
> > Colectora RN 168 Km 472, Paraje El Pozo
> > 3000 Santa Fe, Argentina
> > Tel: +54-342-4511594 (ext 1011)
> > Tel/Fax: +54-342-4511169
> [see attached file: matblocksize-seqaij-csr.diff]
>