ML symbols in petscksp.so

symbols, PETSC_ARCH=ompi is only missing ML symbols. Jed -- next part -- A non-text attachment was scrubbed... Name: signature.asc Type: application/pgp-signature Size: 261 bytes Desc: OpenPGP digital signature URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20091029/b1b2558a/attachment.pgp>

changes for PetscMap

Ok, I have pushed all these changes, including PetscMap --> PetscLayout Barry On Oct 27, 2009, at 4:56 PM, Dima Karpeyev wrote: > What about a different name? I vote for PetscLayout, because that's > what it is. > Dmitry. > > On Tue, Oct 27, 2009 at 4:44 PM, Matthew Knepley > wrot

remove 'x11' test from 'make test'

Thanks for looking that up (I am wondering why it's not default) -- I did that and hopefully X won't need restarting. Boyana -- Boyana Norris, Computer Scientist, Argonne National Laboratory norris at mcs.anl.gov, +1.630.252.7908, http://www.mcs.anl.gov/~norris/ On Oct 28, 2009, at 10:09 PM, B

ML symbols in petscksp.so

crubbed... Name: signature.asc Type: application/pgp-signature Size: 261 bytes Desc: OpenPGP digital signature URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20091029/a899b7e8/attachment.pgp>

most efficient way to use local CSR matrix?

Mark, Barry, Thanks for the help. For now I'm sticking with the approach of copying the csr matrix and using the csr data structures to do the preallocations. I'll eventually get around to writing the code for assembling directly into the petsc matrix. I have two more questions. 1) On 1 p

most efficient way to use local CSR matrix?

,irow,rowptr[i+1]-rowptr[i],&colind[rowptr[i]],&a[rowptr[i]],INSERT_VALUES); >>>>>> } >>>>>> >>>>>> where rowptr and colind are the CSR indexing arrays in the local >>>>>> numbering. >>>>>> >>>>>> I would like to eliminate the duplicate memory (and the copy above). >>>>>> Is there a recommended way to let PETSc reference array 'a' directly or a >>>>>> way to get rid of 'a' and translate our CSR assembly routines to use low >>>>>> level PETSc data structures? >>>>>> >>>>>> So far I've added MatMPIAIJSetPreallocationCSR to try to speed up the >>>>>> first solve, but the documentation is clear that 'a' is not a pointer to >>>>>> the >>>>>> AIJ storage so I still need the copy. I tried using >>>>>> MatCreateMPIAIJWithSplitArrays with the 'o' arrays set to zero but I got >>>>>> indexing errors, and the documentation really seems to imply that the >>>>>> local >>>>>> matrix is really not a standard CSR anyway. If that's true then it seems >>>>>> like the options I have are to write a shell matrix for parallel CSR >>>>>> storage >>>>>> or write some new assembly routines that use MatSetValuesLocal. I'm more >>>>>> concerned about the memory than the overhead of MatSetValuesLocal, but it >>>>>> would certainly be easier on me to use the CSR-based assembly routines we >>>>>> already have. >>>>>> >>>>>> Thanks, >>>>>> Chris >>>>>> >>>>>> >>>>>> >>>>>> >>>>> >>>> >>> >>> >> > -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20091029/7a7f7f4e/attachment.html>

missing CHKERRQ in MatNullSpaceRemove()

While trying to update petsc4py to latest changes in petsc-dev repo, I've noticed that at the very end of MatNullSpaceRemove() we have: if (sp->remove){ ierr = (*sp->remove)(sp,vec,sp->rmctx); } PetscFunctionReturn(0); } Any good reason for not CHKERRQ(ierr) after calling the userd-defin

missing CHKERRQ in MatNullSpaceRemove()

results to which their experiments lead. -- Norbert Wiener -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20091029/22b1856b/attachment.html>

most efficient way to use local CSR matrix?

rors, and the documentation really seems to imply > that the local matrix is really not a standard CSR anyway. If that's > true then it seems like the options I have are to write a shell > matrix for parallel CSR storage or write some new assembly routines > that use MatSetValuesLocal. I'm more concerned about the memory than > the overhead of MatSetValuesLocal, but it would certainly be easier > on me to use the CSR-based assembly routines we already have. > > Thanks, > Chris > > > > > > > > > > > > > -- > What most experimenters take for granted before they begin their > experiments is infinitely more interesting than any results to which > their experiments lead. > -- Norbert Wiener -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20091029/f88ff12f/attachment.html>

ML symbols in petscksp.so

On Thu, Oct 29, 2009 at 4:45 PM, Jed Brown wrote: > > I realize that the real problem here was OpenCV's libml and the fact > that linkers don't resolve symbols by starting with the most recent -L > path [*], but we should at least remember that putting -L{PETSC_LIB_DIR} > at the beginning of the l

most efficient way to use local CSR matrix?

Hmm, this should not happen. The matrix should be identical in both cases (note the initial residual is the same in both cases so they may be identical). Here's one more thing you can try. In the same routine create TWO matrices; one each way, then use MatAXPY() to take the differen

ML symbols in petscksp.so

On Oct 29, 2009, at 3:52 PM, Lisandro Dalcin wrote: > On Thu, Oct 29, 2009 at 4:45 PM, Jed Brown wrote: >> >> I realize that the real problem here was OpenCV's libml and the fact >> that linkers don't resolve symbols by starting with the most recent >> -L >> path [*], but we should at least re

ML symbols in petscksp.so

On Thu, 29 Oct 2009, Barry Smith wrote: > > On Oct 29, 2009, at 3:52 PM, Lisandro Dalcin wrote: > > > On Thu, Oct 29, 2009 at 4:45 PM, Jed Brown wrote: > > > > > > I realize that the real problem here was OpenCV's libml and the fact > > > that linkers don't resolve symbols by starting with the

ML symbols in petscksp.so

ntually gets heinous, and everyone always goes back to all headers. I am for keeping each package in its own libraries. Matt > Satish > -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20091029/5d58c8da/attachment.html>

ML symbols in petscksp.so

On Thu, Oct 29, 2009 at 7:01 PM, Satish Balay wrote: > On Thu, 29 Oct 2009, Barry Smith wrote: > >> >> On Oct 29, 2009, at 3:52 PM, Lisandro Dalcin wrote: >> >> > On Thu, Oct 29, 2009 at 4:45 PM, Jed Brown wrote: >> > > >> > > I realize that the real problem here was OpenCV's libml and the fact >

ML symbols in petscksp.so

On Thu, 29 Oct 2009, Lisandro Dalcin wrote: > >> > IMHO, I think that what Jed suggested in previous mail about using > >> > -Wl,-whole-archive ${PETSC_EXTERNAL_LIB} -Wl,-no-whole-archive when > >> > --with-shared is in action could be a VERY good idea... Then PETSc > >> > link lines will not need

ML symbols in petscksp.so

(a) makes the linking API uniform, -L/path/to/petsc/lib -lpetsc* (b) causes a problem if they update that lib, don't update PETSc, and then expect their code to be calling the updated version Jed -- next part -- A non-text attachment was scrubbed... Name: signatu

most efficient way to use local CSR matrix?

I think I added the correct test code and got a difference of 0, but spooles is still not returning the correct result, at least not until the 2nd iteration. I'm attaching some of the code and output. Chris code for mat creates --- MatCreateMPIAIJ(Py_PETSC_COM

most efficient way to use local CSR matrix?

On Oct 29, 2009, at 5:02 PM, Chris Kees wrote: > I think I added the correct test code and got a difference of 0, but > spooles is still not returning the correct result, at least not > until the 2nd iteration. I'm attaching some of the code and output. We have not been actively updating sp

most efficient way to use local CSR matrix?

Have you run this with the debug version of PETSc? It does additional argument testing that might catch a funny value passed in. Barry On Oct 29, 2009, at 5:02 PM, Chris Kees wrote: > I think I added the correct test code and got a difference of 0, but > spooles is still not returni

most efficient way to use local CSR matrix?

On Oct 29, 2009, at 8:43 PM, Hong Zhang wrote: > > On Oct 29, 2009, at 5:02 PM, Chris Kees wrote: > >> I think I added the correct test code and got a difference of 0, >> but spooles is still not returning the correct result, at least not >> until the 2nd iteration. I'm attaching some of the

most efficient way to use local CSR matrix?

I'm just using a simple configuration: ./config/configure.py --with-clanguage=C --with-cc='/usr/bin/mpicc - arch x86_64' --with-cxx='/usr/bin/mpicxx -arch x86_64' --without- fortran --with-mpi-compilers --without-shared --without-dynamic -- download-parmeti s=ifneeded --download-spooles=ifneede