_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>
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
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>
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
next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110623/22a428a4/attachment.htm>
hment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110623/8b30f9d6/attachment-0001.htm>
"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
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
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>
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
ttachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110623/edb210f3/attachment-0001.htm>
:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110623/453b6563/attachment.htm>
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>
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>
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>
L attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110623/ed891675/attachment.htm>
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
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>
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
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
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
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
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
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
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>
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>
>
--
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>
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>
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>
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>
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>
-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>
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
t.
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110623/79fde707/attachment.pgp>
34 matches
Mail list logo