[petsc-dev] changes coming to PETSc master branch

[Barry Smith]

  Updated in barry/feature-vec-lock and next

  Thanks

  Barry

> On Jan 31, 2015, at 10:53 AM, Satish Balay <balay at mcs.anl.gov> wrote:
> 
> Barry,
> 
> ftp://ftp.mcs.anl.gov/pub/petsc/nightlylogs/archive/2015/01/30/filtered-make_next_arch-pardiso_cg.log
> 
> Perhaps the following fix?
> 
> Satish
> 
> ------------
> diff --git a/src/mat/impls/aij/seq/mkl_pardiso/mkl_pardiso.c b/src/mat/impls/aij/seq/mkl_pardiso/mkl_pardiso.c
> index f76267e..52e8a4f 100644
> --- a/src/mat/impls/aij/seq/mkl_pardiso/mkl_pardiso.c
> +++ b/src/mat/impls/aij/seq/mkl_pardiso/mkl_pardiso.c
> @@ -182,8 +182,8 @@ PetscErrorCode MatSolve_MKL_PARDISO(Mat A,Vec b,Vec x)
> 
> 
>   mat_mkl_pardiso->nrhs = 1;
> -  ierr = VecGetArrayRead(x,&xarray);CHKERRQ(ierr);
> -  ierr = VecGetArray(b,&barray);CHKERRQ(ierr);
> +  ierr = VecGetArray(x,&xarray);CHKERRQ(ierr);
> +  ierr = VecGetArrayRead(b,&barray);CHKERRQ(ierr);
> 
>   /* solve phase */
>   /*-------------*/
> @@ -207,8 +207,8 @@ PetscErrorCode MatSolve_MKL_PARDISO(Mat A,Vec b,Vec x)
> 
> 
>   if (mat_mkl_pardiso->err < 0) SETERRQ1(PETSC_COMM_SELF,PETSC_ERR_LIB,"Error reported by MKL_PARDISO: err=%d. Please check manual\n",mat_mkl_pardiso->err);
> -  ierr = VecRestoreArrayRead(x,&xarray);CHKERRQ(ierr);
> -  ierr = VecRestoreArray(b,&barray);CHKERRQ(ierr);
> +  ierr = VecRestoreArray(x,&xarray);CHKERRQ(ierr);
> +  ierr = VecRestoreArrayRead(b,&barray);CHKERRQ(ierr);
>   mat_mkl_pardiso->CleanUp = PETSC_TRUE;
>   PetscFunctionReturn(0);
> }