16098dc10SBarry Smith 2b45d2f2cSJed Brown #include <petsc-private/matimpl.h> /*I "petscmat.h" I*/ 3b45d2f2cSJed Brown #include <petsc-private/vecimpl.h> 46098dc10SBarry Smith 54a2ae208SSatish Balay #undef __FUNCT__ 64a2ae208SSatish Balay #define __FUNCT__ "MatConvert_Shell" 719fd82e9SBarry Smith PetscErrorCode MatConvert_Shell(Mat oldmat, MatType newtype,MatReuse reuse,Mat *newmat) 8dfbe8321SBarry Smith { 9676c34cdSKris Buschelman Mat mat; 106098dc10SBarry Smith Vec in,out; 11dfbe8321SBarry Smith PetscErrorCode ierr; 12b24ad042SBarry Smith PetscInt i,M,m,*rows,start,end; 136098dc10SBarry Smith MPI_Comm comm; 1487828ca2SBarry Smith PetscScalar *array,zero = 0.0,one = 1.0; 156098dc10SBarry Smith 163a40ed3dSBarry Smith PetscFunctionBegin; 17*ce94432eSBarry Smith ierr = PetscObjectGetComm((PetscObject)oldmat,&comm);CHKERRQ(ierr); 186098dc10SBarry Smith 196098dc10SBarry Smith ierr = MatGetOwnershipRange(oldmat,&start,&end);CHKERRQ(ierr); 206098dc10SBarry Smith ierr = VecCreateMPI(comm,end-start,PETSC_DECIDE,&in);CHKERRQ(ierr); 216098dc10SBarry Smith ierr = VecDuplicate(in,&out);CHKERRQ(ierr); 226098dc10SBarry Smith ierr = VecGetSize(in,&M);CHKERRQ(ierr); 236098dc10SBarry Smith ierr = VecGetLocalSize(in,&m);CHKERRQ(ierr); 247c307921SBarry Smith ierr = PetscMalloc((m+1)*sizeof(PetscInt),&rows);CHKERRQ(ierr); 252205254eSKarl Rupp for (i=0; i<m; i++) rows[i] = start + i; 266098dc10SBarry Smith 27f69a0ea3SMatthew Knepley ierr = MatCreate(comm,&mat);CHKERRQ(ierr); 28f69a0ea3SMatthew Knepley ierr = MatSetSizes(mat,m,M,M,M);CHKERRQ(ierr); 294a6a16e5SBarry Smith ierr = MatSetType(mat,newtype);CHKERRQ(ierr); 30d0f46423SBarry Smith ierr = MatSetBlockSize(mat,oldmat->rmap->bs);CHKERRQ(ierr); 314a6a16e5SBarry Smith ierr = MatSetUp(mat);CHKERRQ(ierr); 326098dc10SBarry Smith 336098dc10SBarry Smith for (i=0; i<M; i++) { 342dcb1b2aSMatthew Knepley ierr = VecSet(in,zero);CHKERRQ(ierr); 356098dc10SBarry Smith ierr = VecSetValues(in,1,&i,&one,INSERT_VALUES);CHKERRQ(ierr); 366098dc10SBarry Smith ierr = VecAssemblyBegin(in);CHKERRQ(ierr); 376098dc10SBarry Smith ierr = VecAssemblyEnd(in);CHKERRQ(ierr); 386098dc10SBarry Smith 396098dc10SBarry Smith ierr = MatMult(oldmat,in,out);CHKERRQ(ierr); 406098dc10SBarry Smith 416098dc10SBarry Smith ierr = VecGetArray(out,&array);CHKERRQ(ierr); 42676c34cdSKris Buschelman ierr = MatSetValues(mat,m,rows,1,&i,array,INSERT_VALUES);CHKERRQ(ierr); 436098dc10SBarry Smith ierr = VecRestoreArray(out,&array);CHKERRQ(ierr); 446098dc10SBarry Smith 456098dc10SBarry Smith } 46606d414cSSatish Balay ierr = PetscFree(rows);CHKERRQ(ierr); 476bf464f9SBarry Smith ierr = VecDestroy(&in);CHKERRQ(ierr); 486bf464f9SBarry Smith ierr = VecDestroy(&out);CHKERRQ(ierr); 49676c34cdSKris Buschelman ierr = MatAssemblyBegin(mat,MAT_FINAL_ASSEMBLY);CHKERRQ(ierr); 50676c34cdSKris Buschelman ierr = MatAssemblyEnd(mat,MAT_FINAL_ASSEMBLY);CHKERRQ(ierr); 51676c34cdSKris Buschelman 52ceb03754SKris Buschelman if (reuse == MAT_REUSE_MATRIX) { 538ab5b326SKris Buschelman ierr = MatHeaderReplace(oldmat,mat);CHKERRQ(ierr); 54c3d102feSKris Buschelman } else { 55676c34cdSKris Buschelman *newmat = mat; 56c3d102feSKris Buschelman } 573a40ed3dSBarry Smith PetscFunctionReturn(0); 586098dc10SBarry Smith } 596098dc10SBarry Smith 606098dc10SBarry Smith 616098dc10SBarry Smith 62