[petsc-users] Fwd: superlu_dist with same_nonzero_pattern

2011-06-23 Thread Hong Zhang
_PASTIX are >>> undeclared either. Am I missing some header files? I used the same >>> procedure as I complied regular petsc, which had no such problems. ?(I >>> tried to attach the configure log in a previous mail, but it is >>> bounced back.) >>> >>> Can you give me some suggestions on using superlu_dist with >>> same_nonzero_pattern? Thank you. >>> >>> Best, >>> Xiangdong >>> >> > -- next part -- A non-text attachment was scrubbed... Name: superlu_dist_test2.tar.gz Type: application/x-gzip Size: 3750 bytes Desc: not available URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110623/e4b5c57c/attachment.bin>

[petsc-users] Ghost communication

2011-06-23 Thread Milan Mitrovic
I thought at some point that I maybe have the order of the arguments wrong... The first few lines in particles_dcop_apply_petsc are: ! get ghosted versions of vectors CALL VecGhostGetLocalForm(from,xl,info) CALL VecSet(to,0.0_MK,info) ! get arrays from pets

[petsc-users] Ghost communication

2011-06-23 Thread Jed Brown
d,info) >END SUBROUTINE dcop_mult > Where do you call VecGhostGetLocalForm()? -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110623/b7cfed59/attachment.htm>

[petsc-users] Ghost communication

2011-06-23 Thread Milan Mitrovic
SUBROUTINE dcop_mult(A,x,y,info) ! args Mat A Vec x Vec y INTEGER, INTENT( OUT) :: info ! vars INTEGER :: id info = 0 ! update ghosts CALL VecGhostUpdateBegin(x,INSERT_VALUES,SCATTER_FORWARD

[petsc-users] Size of off-diagonal matrix

2011-06-23 Thread Jed Brown
next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110623/22a428a4/attachment.htm>

[petsc-users] Ghost communication

2011-06-23 Thread Jed Brown
hment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110623/8b30f9d6/attachment-0001.htm>

[petsc-users] Size of off-diagonal matrix

2011-06-23 Thread Jed Brown
"sparse vectors" (i.e. an additional layer of indirection) or the whole dense global vector would need to be gathered to every process which would be horribly non-scalable. -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipe

[petsc-users] Ghost communication

2011-06-23 Thread Milan Mitrovic
Well, I found where the problem is. I wanted to us KSP to solve a system for which I created a matrix-free matrix, and I provide the MatMult method. I created the vectors the same way I showed before, but when I ran KSPSolve I found that the ghost communication did not happen automatically during t

[petsc-users] Size of off-diagonal matrix

2011-06-23 Thread Jed Brown
defining the sequential matrix operating on those values. Unless you already happen to use the same format, you're better off using MatSetValues(). -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110623/441df6e3/attachment.htm>

[petsc-users] Fwd: superlu_dist with same_nonzero_pattern

2011-06-23 Thread Xiangdong Liang
Thanks, Hong. I tried the runtime option -mat_superlu_dist_equil NO. However, the problem is still there. I've upload my short codes here: http://math.mit.edu/~xdliang/superlu_dist_test.zip In this code, I first generate my sparse matrix M with dimension N-by-N, then solve Mx=J for a few times. E

[petsc-users] Size of off-diagonal matrix

2011-06-23 Thread Jed Brown
ttachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110623/edb210f3/attachment-0001.htm>

[petsc-users] PETSs Mesh Examples

2011-06-23 Thread Jed Brown
: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110623/453b6563/attachment.htm>

[petsc-users] SNES_DIVERGED_LS_FAILURE

2011-06-23 Thread Jed Brown
m for the Bratu problem on page 2 here: http://www.nada.kth.se/kurser/kth/2D1290/comp03_3.pdf -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110623/cee69d03/attachment.htm>

[petsc-users] SNES_DIVERGED_LS_FAILURE

2011-06-23 Thread Juha Jäykkä
achment was scrubbed... Name: not available Type: application/pgp-signature Size: 836 bytes Desc: This is a digitally signed message part. URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110623/b47f0a67/attachment.pgp>

[petsc-users] SNES_DIVERGED_LS_FAILURE

2011-06-23 Thread Juha Jäykkä
available Type: application/pgp-signature Size: 836 bytes Desc: This is a digitally signed message part. URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110623/417a20ff/attachment.pgp>

[petsc-users] PCCOMPOSITE and ApplyTranspose

2011-06-23 Thread Jed Brown
L attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110623/ed891675/attachment.htm>

[petsc-users] Size of off-diagonal matrix

2011-06-23 Thread John Fettig
On Thu, Jun 23, 2011 at 11:17 AM, Jed Brown wrote: > On Thu, Jun 23, 2011 at 17:09, John Fettig wrote: >> >> Is there any way to access the map from global numbering to the >> compressed numbering? > > 1. Certainly not before you create the matrix. > 2. You have this information in some form if y

[petsc-users] SNES_DIVERGED_LS_FAILURE

2011-06-23 Thread Juha Jäykkä
ed... Name: not available Type: application/pgp-signature Size: 836 bytes Desc: This is a digitally signed message part. URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110623/f439add8/attachment.pgp>

[petsc-users] SNES_DIVERGED_LS_FAILURE

2011-06-23 Thread Klaij, Christiaan
you can have valgrind suppress the false errors from openmpi, see http://www.open-mpi.de/faq/?category=debugging#valgrind_clean > Date: Thu, 23 Jun 2011 11:49:44 +0100 > From: Juha J?ykk? > Subject: Re: [petsc-users] SNES_DIVERGED_LS_FAILURE > To: PETSc users list > Message-ID: <201106231149.505

[petsc-users] Size of off-diagonal matrix

2011-06-23 Thread John Fettig
On Thu, Jun 23, 2011 at 11:06 AM, John Fettig wrote: > On Thu, Jun 23, 2011 at 11:03 AM, Jed Brown wrote: >> On Thu, Jun 23, 2011 at 16:59, John Fettig wrote: >>> >>> I understand everything you say here except what "compressing" means. >>> It is compressed by virtue of being in CSR format, but

[petsc-users] Size of off-diagonal matrix

2011-06-23 Thread John Fettig
On Thu, Jun 23, 2011 at 11:03 AM, Jed Brown wrote: > On Thu, Jun 23, 2011 at 16:59, John Fettig wrote: >> >> I understand everything you say here except what "compressing" means. >> It is compressed by virtue of being in CSR format, but then you do >> some additional mapping with the column indic

[petsc-users] Size of off-diagonal matrix

2011-06-23 Thread John Fettig
On Thu, Jun 23, 2011 at 10:50 AM, Jed Brown wrote: > On Thu, Jun 23, 2011 at 16:47, John Fettig wrote: >> >> I'm trying to use MatCreateMPIAIJWithSplitArrays, but it isn't clear >> to me how the off-diagonal portion is supposed to be constructed. > > Okay, it's defined by compressing all the sort

[petsc-users] Size of off-diagonal matrix

2011-06-23 Thread John Fettig
On Thu, Jun 23, 2011 at 10:25 AM, Jed Brown wrote: > On Thu, Jun 23, 2011 at 14:53, John Fettig wrote: >> ?Is there >> anyplace in the documentation that would help me understand this? >> Alternatively, where would I look in mpiaij.c? > > The MPIAIJ format works with blocks that are in SeqAIJ for

[petsc-users] Size of off-diagonal matrix

2011-06-23 Thread John Fettig
On Wed, Jun 22, 2011 at 10:30 PM, Barry Smith wrote: > > On Jun 22, 2011, at 9:02 AM, John Fettig wrote: > >> On Wed, Jun 22, 2011 at 9:56 AM, Matthew Knepley >> wrote: >>> It is compressed, so it only has columns that have a nonzero. >> >> Thanks for answering so quickly. ?I had just assumed th

[petsc-users] PETSs Mesh Examples

2011-06-23 Thread Matthew Knepley
armareddy84 at gmail.com >>> Phone: +1-512-350-9082 >>> United States of America. >>> >>> >> >> >> -- >> 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 >> > > > > -- > - > Dharmendar Reddy Palle > Graduate Student > Microelectronics Research center, > University of Texas at Austin, > 10100 Burnet Road, Bldg. 160 > MER 2.608F, TX 78758-4445 > e-mail: dharmareddy84 at gmail.com > Phone: +1-512-350-9082 > United States of America. > > -- 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-users/attachments/20110623/9d557f38/attachment.htm>

[petsc-users] PETSs Mesh Examples

2011-06-23 Thread Dharmendar Reddy
ER 2.608F, TX 78758-4445 >> e-mail: dharmareddy84 at gmail.com >> Phone: +1-512-350-9082 >> United States of America. >> >> > > > -- > 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 > -- - Dharmendar Reddy Palle Graduate Student Microelectronics Research center, University of Texas at Austin, 10100 Burnet Road, Bldg. 160 MER 2.608F, TX 78758-4445 e-mail: dharmareddy84 at gmail.com Phone: +1-512-350-9082 United States of America. -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110623/7c4951d6/attachment.htm>

[petsc-users] SNES_DIVERGED_LS_FAILURE

2011-06-23 Thread Matthew Knepley
> -- 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-users/attachments/20110623/c7120d58/attachment-0001.htm>

[petsc-users] PETSs Mesh Examples

2011-06-23 Thread Matthew Knepley
croelectronics Research center, > University of Texas at Austin, > 10100 Burnet Road, Bldg. 160 > MER 2.608F, TX 78758-4445 > e-mail: dharmareddy84 at gmail.com > Phone: +1-512-350-9082 > United States of America. > > -- 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-users/attachments/20110623/2bb2f805/attachment.htm>

[petsc-users] PETSs Mesh Examples

2011-06-23 Thread Dharmendar Reddy
harmendar Reddy Palle Graduate Student Microelectronics Research center, University of Texas at Austin, 10100 Burnet Road, Bldg. 160 MER 2.608F, TX 78758-4445 e-mail: dharmareddy84 at gmail.com Phone: +1-512-350-9082 United States of America. -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110623/997ea4f0/attachment.htm>

[petsc-users] PETSs Mesh Examples

2011-06-23 Thread Matthew Knepley
gin 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-users/attachments/20110623/f32733fa/attachment.htm>

[petsc-users] SNES_DIVERGED_LS_FAILURE

2011-06-23 Thread Matthew Knepley
e 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-users/attachments/20110623/07ff6c1e/attachment-0001.htm>

[petsc-users] PETSs Mesh Examples

2011-06-23 Thread Dharmendar Reddy
-mail: dharmareddy84 at gmail.com Phone: +1-512-350-9082 United States of America. -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110623/fa662b9b/attachment.htm>

[petsc-users] PCCOMPOSITE and ApplyTranspose

2011-06-23 Thread Andrew T Barker
I've made a PCCOMPOSITE that contains two PCSHELLS, both of which have a transpose set with PCShellSetApplyTranspose(). However the composite PC does not seem to support the transpose: [0]PETSC ERROR: No support for this operation for this object type! [0]PETSC ERROR: PC does not have apply t

[petsc-users] SNES_DIVERGED_LS_FAILURE

2011-06-23 Thread Juha Jäykkä
t. URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110623/79fde707/attachment.pgp>