Hi Paul, thanks, I'll take care of the patch. You'll get comments tomorrow.
Best regards, Karli On 03/09/2013 02:30 PM, Paul Mullowney wrote: > Hi, > > I've gotten GPU based icc preconditioners working with the following set > of patches. One uses ICC by loading an upper triangular matrix into the > cusparse class, and then setting the symmetry option, i.e. > > ierr = > MatCreateSeqAIJCUSPARSE(comm,n,m,PETSC_NULL,num_entries_per_row,&A);CHKERRQ(ierr); > > ierr = MatSetOption(A,MAT_SYMMETRIC,PETSC_TRUE); > > I've organized the implementations into 2 files: > 1) icc-cholesky-cusparse-organization.patch : This patch has the basic > infrastructure for calling the symbolic and numeric factorizations (both > ICC and Cholesky). Currently these routines make scoping calls to the > implementations in aij/seq/aijfact.c > > 2) icc-cholesky-implementation.patch : This patch contains the matrix > construction on the GPU. This is done in the method > MatSeqAIJCUSPARSEBuildICCTriMatrices. > > This patch has significant changes so comments or suggestions are > definitely welcome. > Thanks, > -Paul
